Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1611 1365 (N33W29 L=124)   477     72        0 (1592/1327-1629/1406) -68.2 W27-W31 N30-N35
      2 1796 1113 (N23W50 L=145)  2340    361        4 (1681/ 986-1899/1193)  42.0 W41-W65 N12-N30
      3 1391  873 (S03W25 L=120)   225     51        0 (1369/ 852-1424/ 888)  56.3 W23-W26 S04-S01
      4  376  759 (S37E39 L= 56)  2456    269        7 ( 302/ 672- 455/ 866)  46.8 E48-E32 S44-S28
      5 1746  659 (S04W60 L=155)   585     67        1 (1734/ 622-1760/ 702) -38.7 W56-W62 S05-S01
      6 1363  611 (S18W33 L=128)   651    202        1 (1291/ 574-1442/ 680)  63.4 W28-W34 S21-S11
      7 1137  614 (S24W19 L=114)   366     55        2 (1105/ 589-1170/ 635)  11.3 W16-W21 S25-S22
      8 1443  570 (S18W40 L=135)   959    174        1 (1377/ 534-1531/ 604)  48.8 W35-W44 S21-S13
      9 1058  573 (S29W16 L=111)   298     65        1 (1035/ 551-1092/ 600)  56.3 W13-W17 S30-S27
     10  670  506 (S46E09 L= 86)   270     71        0 ( 642/ 476- 688/ 539)  71.6 E10-E07 S48-S42
     11  951  498 (S37W12 L=107)   299     44        2 ( 930/ 475- 978/ 524)  63.4 W10-W13 S39-S35
  

    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