Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  329 1084 (S14E47 L=291)  1813    429        7 ( 196/1025- 406/1175)  26.6 E59-E40 S17-S07
      2 1609 1086 (N08W35 L= 13)   361    111        0 (1593/1026-1623/1147) -82.9 W33-W35 N04-N12
      3  122  920 (S25E66 L=272)  1747    249        0 (  93/ 835- 151/ 986)  71.6 E69-E63 S29-S20
      4 1773  931 (N04W51 L= 29)   667    132        3 (1749/ 876-1795/ 973)  90.0 W49-W53 N00-N06
      5  622  881 (S22E22 L=316)   404    129        1 ( 565/ 842- 672/ 930)  40.6 E26-E19 S24-S18
      6 1301  725 (S19W23 L=  1)   483    156        0 (1223/ 701-1363/ 770)  36.9 W17-W26 S20-S14
      7 1624  716 (S12W45 L= 23)   578    154        0 (1564/ 684-1677/ 736)   0.0 W40-W50 S13-S10
      8  425  618 (S41E38 L=300)   700    290        1 ( 349/ 564- 499/ 694)  34.5 E46-E30 S46-S35
      9 1272  617 (S25W24 L=  2)   384    103        0 (1217/ 593-1333/ 647)  39.8 W21-W27 S27-S22
     10 1692  593 (S17W55 L= 33)   554    145        2 (1655/ 551-1734/ 635) -29.7 W50-W58 S19-S14
     11 1494  468 (S29W45 L= 23)  3979    699        9 (1249/ 423-1724/ 528)  17.5 W28-W66 S35-S20
  

    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