Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1109 1738 (N39E17 L= 62)   649     66        1 (1063/1713-1152/1758)  -8.1 E20-E13 N36-N40
      2 1332 1537 (N31W05 L= 83)  2424    717        5 (1236/1240-1415/1716) -53.1 E06-W16 N14-N42
      3 1521 1632 (N46W18 L= 96)   468     65        0 (1484/1607-1554/1665)  -8.1 W13-W21 N43-N47
      4  913 1526 (N19E20 L= 59)   843    220        1 ( 854/1437- 986/1612) -31.6 E26-E13 N14-N23
      5  263 1472 (N03E67 L= 11)  1509    213        6 ( 230/1403- 303/1542) -63.4 E72-E62 S01-N07
      6 1782 1208 (N25W44 L=123)  1557    236        1 (1697/1124-1855/1296) -16.4 W35-W52 N22-N27
      7 1559 1029 (N07W31 L=109)   399    146        0 (1512/ 978-1614/1057)   0.0 W27-W35 N05-N08
      8  256  925 (S30E50 L= 28)   673    160        1 ( 207/ 898- 310/ 944)   5.7 E55-E45 S31-S27
      9  739  447 (S50E02 L= 76)   683    115        2 ( 686/ 411- 805/ 489)  -5.2 E07-W04 S51-S46
     10  861  211 (S64W26 L=104)   433     41        1 ( 836/ 194- 886/ 228)  -7.1 W21-W29 S66-S62
  

    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