Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  894 1763 (N36E29 L= 89)   706    127        1 ( 834/1727- 955/1790)   0.0 E33-E22 N34-N38
      2 1079 1590 (N28E09 L=108)   280     68        1 (1057/1559-1112/1623)  78.7 E10-E07 N25-N30
      3  295 1565 (N13E67 L= 50)  1080    132        4 ( 253/1518- 331/1615)  74.1 E72-E64 N09-N16
      4 1414 1573 (N35W13 L=131)   697    153        0 (1342/1538-1480/1595)  16.7 W08-W18 N32-N37
      5 1136 1348 (N14E00 L=117)   204    105        0 (1105/1301-1173/1388) -45.0 E02-W02 N11-N15
      6  524 1244 (S04E34 L= 83)   369    127        0 ( 501/1185- 539/1291) -81.9 E36-E32 S07-N00
      7 1696 1161 (N16W39 L=157)   706    136        0 (1639/1106-1740/1206) -26.6 W34-W43 N13-N17
      8  421 1157 (S10E40 L= 78)   225     84        0 ( 411/1115- 429/1202) -78.7 E40-E39 S12-S07
      9  212 1153 (S13E59 L= 59)   556     89        1 ( 171/1122- 250/1181) -11.3 E64-E54 S14-S11
     10  793 1061 (S09E14 L=104)   814    321        0 ( 674/ 987- 899/1117)  24.8 E20-E07 S15-S05
     11 1915 1048 (N16W64 L=181)   602     79        1 (1900/1000-1932/1091)  90.0 W62-W66 N13-N19
     12  221  995 (S22E56 L= 61)   857    107        1 ( 172/ 956- 270/1039)  24.4 E62-E51 S24-S19
     13 1749  922 (N04W48 L=166)  1066    219        1 (1723/ 817-1766/1027) -74.7 W44-W50 S02-N09
     14  931  765 (S24E00 L=118)   559    232        1 ( 832/ 725-1025/ 818)   0.0 E07-W06 S25-S22
     15 1722  671 (S10W54 L=172)  1529    211        5 (1649/ 579-1775/ 751) -73.3 W47-W59 S14-S04
     16 1780  632 (S10W63 L=180)   945    114        2 (1742/ 584-1817/ 689)  81.9 W59-W65 S13-S06
     17 1655  574 (S17W54 L=172)   367     59        1 (1626/ 548-1674/ 609) -21.8 W51-W56 S18-S14
     18 1672  464 (S22W63 L=181)  1311    187        2 (1575/ 436-1748/ 502)  31.6 W55-W68 S26-S17
  

    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