Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1091 1770 (N52E11 L=162)   636    131        0 (1031/1745-1150/1791)  33.7 E15-E05 N48-N54
      2  745 1564 (N29E29 L=144)  2055    386        4 ( 586/1466- 906/1642) -13.1 E44-E14 N25-N32
      3 1621 1560 (N44W43 L=215)   818     98        0 (1569/1535-1668/1585)  12.1 W35-W49 N41-N47
      4  924 1297 (N15E10 L=162)   380    142        2 ( 875/1232- 964/1344) -40.6 E13-E06 N11-N17
      5  306 1284 (N05E54 L=118)   347     45        1 ( 287/1255- 331/1309) -36.9 E56-E51 N03-N06
      6 1730 1170 (N19W47 L=219)   283     48        0 (1706/1145-1753/1188) -21.8 W44-W49 N17-N20
      7 1897 1038 (N13W67 L=240)  4082    208        3 (1860/ 939-1939/1124) -61.4 W61-W76 N07-N18
      8 1797  937 (N06W55 L=228)   706    139        1 (1779/ 868-1813/1003) -82.9 W53-W57 N01-N09
      9  940  844 (S12E02 L=170)   355    124        0 ( 889/ 812- 988/ 889)  50.2 E05-W00 S14-S08
     10 1017  665 (S22W05 L=177)  1669    523        1 ( 789/ 621-1182/ 782)  34.3 E09-W13 S27-S12
     11 1143  215 (S54W33 L=205)   508     57        1 (1114/ 195-1165/ 240)  63.4 W29-W36 S55-S51
  

    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