Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  573 1726 (N32E48 L=284)   470     44        1 ( 549/1701- 599/1747) -23.2 E52-E45 N30-N33
      2  498 1640 (N25E49 L=283)   418     41        0 ( 471/1614- 519/1662) -15.9 E52-E45 N23-N26
      3  283 1410 (N08E58 L=273)  2684    319        9 ( 233/1275- 329/1523) -45.0 E68-E51 N00-N15
      4 1105 1410 (N18E01 L=331)   498    210        0 (1003/1396-1208/1427)   9.5 E07-W05 N16-N19
      5 1504 1281 (N17W26 L=357)  2130    304       10 (1385/1178-1592/1375) -11.3 W17-W32 N10-N21
      6 1116 1306 (N11W01 L=333)   205     93        0 (1087/1263-1156/1333)  59.0 W00-W03 N08-N13
      7  193 1301 (N01E65 L=267)   662     83        1 ( 165/1276- 232/1322)   8.1 E69-E59 N00-N02
      8 1377 1246 (N12W18 L=349)   998    235        4 (1305/1165-1445/1300) -21.8 W12-W22 N07-N15
      9 1211 1253 (N10W08 L=339)   298     85        2 (1157/1225-1238/1280) -11.3 W04-W09 N08-N11
     10  310 1120 (S09E49 L=283)   269     71        0 ( 270/1105- 345/1132)   0.0 E52-E45 S09-S08
     11 1465 1094 (N05W25 L=357)   248    109        0 (1410/1081-1519/1106)  23.2 W21-W28 N03-N06
     12 1706  904 (S02W46 L= 17)   565    123        1 (1642/ 883-1771/ 923)  16.7 W40-W50 S03-N00
     13  428  829 (S25E37 L=294)   996    286        1 ( 333/ 755- 560/ 871)  22.8 E46-E27 S30-S21
     14 1223  656 (S25W19 L=350)  3578    951        7 ( 931/ 578-1641/ 820)   8.1 W00-W49 S29-S17
     15  449  658 (S36E37 L=295)   264     31        2 ( 426/ 643- 476/ 678)   0.0 E39-E34 S37-S34
     16 1781  603 (S16W64 L= 36)  1330    183        2 (1713/ 555-1860/ 648)  36.9 W56-W73 S20-S11
     17  487  535 (S44E35 L=296)  2128    547        4 ( 252/ 486- 745/ 563)   2.0 E69-E11 S46-S40
     18  929  463 (S42W04 L=335)   310    113        0 ( 879/ 446- 989/ 485)   0.0 W00-W08 S43-S41
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y