Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  906 1756 (N36E27 L= 90)   786    173        3 ( 837/1702- 986/1789)  -8.1 E33-E19 N33-N37
      2  304 1571 (N14E67 L= 50)   683    120        2 ( 273/1525- 335/1617)  81.9 E69-E64 N10-N17
      3 1083 1585 (N28E09 L=108)   270     67        0 (1063/1555-1109/1617)  90.0 E10-E07 N25-N30
      4 1426 1568 (N35W14 L=132)   731    146        2 (1350/1531-1497/1592)  14.0 W08-W20 N32-N36
      5  532 1241 (S04E34 L= 84)   350    107        0 ( 510/1193- 547/1288) -68.2 E35-E32 S06-N00
      6 1701 1159 (N16W40 L=157)   675    144        0 (1644/1110-1747/1203) -11.3 W34-W44 N14-N17
      7  217 1152 (S13E59 L= 59)   578     98        1 ( 172/1120- 258/1175) -10.3 E64-E53 S14-S11
      8  795 1060 (S09E14 L=103)   660    291        0 ( 687/ 994- 891/1114)  28.6 E19-E08 S14-S05
      9 1919 1050 (N17W64 L=182)   706     91        2 (1901/1000-1938/1092)  80.5 W62-W67 N13-N19
     10  226  995 (S22E56 L= 61)   881    114        0 ( 178/ 955- 276/1043)  24.4 E61-E50 S24-S18
     11 1755  927 (N05W48 L=166)  1026    228        0 (1733/ 822-1773/1033) -65.6 W44-W51 S01-N10
     12  936  762 (S24W00 L=118)   563    212        0 ( 839/ 722-1022/ 816)  -4.8 E06-W06 S25-S22
     13 1759  717 (S06W56 L=173)   423     50        1 (1739/ 684-1775/ 748)  90.0 W53-W58 S08-S04
     14 1710  645 (S12W55 L=172)  1101    177        1 (1656/ 575-1742/ 704) -24.4 W47-W60 S15-S07
     15 1786  632 (S10W63 L=181)  1060    139        3 (1742/ 585-1822/ 694)  81.9 W60-W66 S13-S06
     16 1656  570 (S17W54 L=171)   349     36        1 (1627/ 548-1679/ 592)  14.0 W51-W56 S18-S16
     17 1679  460 (S22W64 L=181)  1361    148        3 (1605/ 432-1742/ 499)  32.5 W57-W69 S24-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