Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  607 1679 (N42E32 L= 79)   405     87        1 ( 592/1634- 623/1724) -66.8 E33-E29 N38-N45
      2 1664 1629 (N35W57 L=167)   836     77        2 (1622/1598-1696/1652) -15.9 W51-W61 N32-N36
      3 1225 1542 (N29W15 L=126)  1003    304        6 (1092/1485-1353/1596) -19.4 W06-W23 N24-N32
      4  313 1360 (N21E50 L= 61)   444     68        1 ( 300/1325- 332/1403) -68.2 E51-E47 N18-N23
      5 1500 1360 (N16W32 L=142)   383    119        1 (1453/1332-1539/1394)  23.2 W28-W35 N14-N18
      6 1114 1300 (N13W06 L=117)   200    110        0 (1107/1240-1124/1353) -90.0 W05-W06 N09-N16
      7 1915  984 (S06W67 L=177)   962    158        2 (1865/ 924-1944/1051) -30.3 W60-W72 S09-S02
      8  509  958 (S05E32 L= 78)   498    169        1 ( 440/ 909- 550/1000) -45.0 E36-E29 S07-S01
      9 1779  932 (S10W51 L=162)  1222    208        7 (1695/ 884-1874/ 994)  16.4 W43-W61 S12-S05
     10  606  907 (S08E26 L= 84)   427     89        2 ( 583/ 864- 640/ 950) -78.7 E27-E23 S10-S05
     11  214  813 (S11E60 L= 51)  3071    392        6 ( 149/ 679- 286/ 916)  60.3 E67-E50 S19-S05
     12 1822  833 (S15W57 L=167)  2806    263        8 (1738/ 740-1889/ 910)  48.4 W48-W66 S21-S10
     13  712  860 (S11E20 L= 91)   685    197        3 ( 627/ 833- 816/ 883)   0.0 E25-E13 S12-S10
     14  394  719 (S18E44 L= 66)  1804    431        1 ( 296/ 611- 463/ 865)  82.4 E57-E38 S24-S09
     15 1703  815 (S17W46 L=156)   666     99        2 (1671/ 772-1732/ 864)  71.6 W43-W48 S19-S13
     16 1233  669 (S25W12 L=123)   495    180        0 (1182/ 601-1282/ 748)  60.9 W09-W15 S28-S19
     17  719  383 (S42E28 L= 82)  1660    328        3 ( 581/ 282- 828/ 483)  22.1 E48-E16 S49-S35
     18 1122  210 (S60W06 L=117)   932     98        3 (1082/ 172-1174/ 238)  21.8 W01-W12 S64-S57
  

    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