Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  711 1841 (N62E24 L=358)   585     40        1 ( 682/1813- 735/1864)   0.0 E29-E20 N58-N64
      2  813 1736 (N49E09 L= 13)   383     64        1 ( 796/1701- 834/1767) -80.5 E10-E07 N45-N51
      3 1350 1522 (N26W27 L= 49)   998    227        1 (1238/1473-1439/1546) -24.8 W19-W32 N22-N28
      4 1046 1432 (N23W05 L= 27)   198     75        1 (1011/1418-1083/1450) -26.6 W03-W07 N22-N24
      5 1694 1398 (N15W51 L= 73)   975    166        2 (1637/1367-1745/1435) -26.6 W45-W56 N13-N17
      6  669 1384 (N24E20 L=  2)   935    260        3 ( 568/1323- 794/1427) -26.6 E27-E11 N18-N27
      7  346 1410 (N29E47 L=335)   390     69        0 ( 309/1393- 387/1425)   7.1 E50-E42 N28-N30
      8 1315 1272 (N11W20 L= 42)   211     77        1 (1286/1238-1340/1296) -63.4 W18-W21 N08-N12
      9 1260 1174 (N06W15 L= 37)   551    215        0 (1193/1132-1360/1233) -33.7 W11-W21 N02-N09
     10  372 1177 (N14E42 L=340)   277    106        0 ( 355/1127- 395/1222)  68.2 E43-E39 N11-N16
     11  996 1117 (N04E01 L= 21)   266     91        1 ( 982/1072-1011/1169)  80.5 E01-W00 N01-N07
     12 1319 1092 (N00W18 L= 40)   188     76        0 (1281/1076-1350/1106) -14.0 W16-W20 N00-N01
     13 1833  940 (S13W57 L= 79)   533     68        1 (1796/ 920-1865/ 965)  18.4 W53-W61 S14-S11
     14 1729  825 (S19W47 L= 69)   509     81        4 (1696/ 792-1764/ 866) -59.0 W44-W50 S21-S16
     15 1651  737 (S24W41 L= 63)   912    219        4 (1629/ 634-1683/ 842) -85.6 W38-W44 S30-S17
     16  848  617 (S24E16 L=  6)   681    133        2 ( 808/ 556- 880/ 686)  63.4 E18-E12 S27-S19
     17  745  504 (S30E25 L=357)   405     77        1 ( 707/ 472- 786/ 530)  15.9 E28-E21 S32-S29
  

    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