Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  715 1631 (N30E28 L=127)  1192    387        1 ( 666/1553- 831/1768) -32.3 E37-E18 N24-N40
      2 1714 1398 (N25W46 L=201)  1189    183        2 (1645/1337-1769/1458) -26.6 W40-W50 N22-N28
      3  915 1410 (N17E10 L=145)   209     71        1 ( 879/1397- 955/1428)   0.0 E12-E07 N16-N18
      4  307 1365 (N11E53 L=102)   559    120        1 ( 268/1312- 341/1418)  81.9 E55-E50 N07-N14
      5 1222 1353 (N16W09 L=164)   239     50        1 (1188/1333-1244/1376)   0.0 W07-W10 N15-N17
      6 1320 1321 (N15W15 L=170)   376     78        2 (1289/1285-1348/1364)  78.7 W13-W17 N13-N18
      7 1112 1214 (N07W04 L=159)  1285    382        2 ( 968/1060-1189/1334) -45.0 E06-W08 S01-N12
      8 1742 1191 (N13W46 L=201)   266    109        0 (1700/1171-1774/1217)   0.0 W42-W49 N11-N14
      9  552 1132 (S02E30 L=125)   534    146        0 ( 493/1088- 619/1162) -14.0 E33-E25 S04-N00
     10 1015 1045 (S04E01 L=154)   212     65        0 ( 982/1026-1057/1063)   0.0 E02-W01 S04-S03
     11  855  973 (S10E10 L=145)   257    104        0 ( 811/ 943- 901/1005)  38.7 E12-E07 S11-S07
     12  143  858 (S20E68 L= 87)  1785    150        2 ( 113/ 792- 171/ 942) -90.0 E74-E64 S23-S14
     13 1407  886 (S10W25 L=180)   352     49        1 (1377/ 853-1431/ 914) -56.3 W22-W26 S11-S07
     14  279  858 (S20E51 L=104)   307     53        1 ( 259/ 833- 294/ 891)  63.4 E53-E49 S21-S17
     15 1096  785 (S19W07 L=162)  1426    710        0 ( 843/ 719-1401/ 855)   8.4 E09-W25 S23-S13
     16  385  656 (S32E45 L=110)   765    174        0 ( 308/ 622- 454/ 700)  14.0 E53-E37 S34-S28
     17 1408  464 (S35W36 L=191)   816    127        1 (1341/ 444-1482/ 488)   0.0 W29-W42 S36-S34
     18 1412  364 (S42W42 L=197)   748    140        1 (1371/ 335-1483/ 412)  -4.1 W36-W50 S44-S38
  

    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