Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1286 1776 (N41W52 L=169)   408     79        0 (1271/1734-1301/1812)  20.0 W46-W57 N38-N43
      2  191 1270 (N41E65 L= 53)  1007    120        2 ( 173/1203- 203/1324) -69.4 E70-E62 N36-N44
      3  362 1043 (N23E43 L= 75)  1210    323        3 ( 249/ 928- 450/1123) -40.9 E54-E37 N15-N29
      4  153 1041 (N27E68 L= 50)  1913    133        5 ( 119/ 983- 190/1097)  40.6 E74-E61 N23-N29
      5  191  923 (N19E63 L= 55)  2448    292        5 ( 140/ 793- 236/1042)  65.0 E72-E56 N11-N26
      6 1753  925 (S20W44 L=162)   836    180        1 (1712/ 855-1797/1010)  35.0 W39-W49 S23-S16
      7  256  945 (N20E55 L= 63)   594    135        2 ( 239/ 870- 275/ 998)  60.3 E56-E52 N15-N22
      8 1502  901 (S13W24 L=142)  1069    306        0 (1381/ 841-1618/ 998)   3.0 W15-W34 S15-S11
      9 1799  703 (S37W48 L=166)   437     44        1 (1772/ 682-1825/ 724)   0.0 W44-W51 S38-S35
     10  862  391 (S26E30 L= 88)  3388    893        2 ( 688/ 252-1145/ 624)  11.8 E47-E03 S34-S18
     11 1209  509 (S28E04 L=114)   623    237        0 (1140/ 439-1243/ 619)  39.3 E10-W01 S31-S22
     12 1669  567 (S41W33 L=151)   727    132        0 (1615/ 538-1727/ 607)   0.0 W27-W39 S42-S39
  

    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