Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1164 1817 (N57E14 L= 36)  1796    256        2 (1057/1765-1270/1848)  39.5 E22-E04 N48-N64
      2 1779 1131 (N21W51 L=101)  1867    265        4 (1679/1064-1850/1195)  -3.0 W40-W59 N17-N24
      3 1675 1108 (N17W41 L= 91)   984    203        4 (1612/1018-1726/1167) -32.5 W35-W46 N12-N20
      4  244 1123 (S11E57 L=353)   895    173        1 ( 198/1057- 301/1164)  48.8 E61-E51 S15-S07
      5 1770 1033 (N14W51 L=101)   785    134        5 (1745/ 969-1807/1085)  90.0 W48-W55 N10-N17
      6 1229  610 (S22W22 L= 72)   338    120        0 (1178/ 580-1284/ 648)  50.2 W19-W24 S24-S18
      7 1708  498 (S18W67 L=117)   840    133        0 (1673/ 443-1738/ 552) -66.8 W64-W68 S21-S14
      8 1492  437 (S27W49 L= 99)  1160    157        4 (1426/ 410-1572/ 462)  33.7 W43-W55 S29-S23
      9 1302  373 (S36W39 L= 89)   432     76        1 (1265/ 359-1342/ 396)  45.0 W35-W41 S37-S33
     10 1165  233 (S50W40 L= 90)   456     46        1 (1141/ 205-1186/ 256) -24.0 W36-W45 S52-S48
     11 1324  192 (S49W63 L=113)   785     83        2 (1289/ 166-1356/ 210)  -2.7 W55-W76 S49-S46
  

    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