Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1165 1878 (N57W59 L=200)   468     75        0 (1140/1840-1189/1911)   7.4 W48-W71 N55-N58
      2 1090 1417 (N28W14 L=156)   390     74        2 (1062/1379-1117/1454)  36.9 W12-W17 N26-N30
      3 1507 1388 (N16W40 L=182)   357     44        0 (1475/1364-1538/1412) -45.0 W38-W41 N14-N18
      4 1747  874 (S20W45 L=186)   371     60        2 (1715/ 854-1782/ 892) -23.2 W41-W48 S21-S18
      5 1441  743 (S19W19 L=161)   260     53        0 (1421/ 719-1465/ 772)  26.6 W17-W21 S20-S17
      6 1705  643 (S34W40 L=182)   408     31        1 (1687/ 618-1723/ 671) -76.0 W38-W42 S36-S32
      7  673  454 (S20E37 L=105)  1787    375        2 ( 570/ 318- 764/ 588)  29.1 E52-E25 S28-S12
      8  431  480 (S15E55 L= 86)   780    172        0 ( 403/ 408- 452/ 561)  35.5 E63-E49 S19-S09
      9  922  382 (S30E24 L=118)   286     61        0 ( 889/ 364- 947/ 402)  11.3 E26-E21 S31-S29
     10  545  292 (S30E61 L= 81)  1131    113        1 ( 488/ 258- 579/ 326) -51.3 E68-E53 S33-S27
  

    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