Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  668 1744 (N27E47 L=114)   567     91        0 ( 646/1700- 695/1785) -60.3 E49-E44 N23-N30
      2 1019 1701 (N33E21 L=140)   646     97        1 ( 988/1665-1060/1748) -35.5 E24-E17 N30-N36
      3  390 1517 (N07E55 L=106)  1681    291        3 ( 275/1475- 515/1554)  22.8 E65-E46 N02-N11
      4  801 1375 (N07E22 L=140)   353    114        0 ( 765/1338- 836/1426)  80.5 E22-E20 N04-N10
      5  360 1338 (S04E50 L=111)   823    215        1 ( 290/1255- 414/1405) -26.6 E58-E43 S08-N00
      6 1631 1347 (N29W29 L=191)   434     94        0 (1598/1296-1657/1384) -38.7 W26-W32 N26-N30
      7  183 1205 (S15E65 L= 97)  1613    320        0 (  84/1084- 260/1316)  38.3 E81-E58 S22-S07
      8 1409 1204 (N13W17 L=178)   203    109        1 (1379/1174-1448/1222)   0.0 W15-W19 N11-N14
      9  360 1174 (S14E45 L=116)   597     65        2 ( 335/1132- 391/1214) -68.2 E48-E42 S16-S11
     10  456 1145 (S14E37 L=124)   493    125        1 ( 425/1104- 518/1174)  45.0 E39-E33 S16-S11
     11 1367  750 (S13W27 L=188)  1024    290        2 (1241/ 718-1504/ 792)  32.7 W20-W34 S16-S07
     12 1688  416 (S19W71 L=232)   519     70        0 (1659/ 387-1719/ 439)  68.2 W68-W73 S21-S16
  

    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