Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  798 1803 (N35E46 L=257)  2432    384        4 ( 662/1704- 879/1865)  25.0 E56-E33 N28-N40
      2  486 1656 (N17E58 L=245)   672    123        3 ( 462/1594- 512/1707) -53.1 E61-E53 N12-N20
      3  579 1565 (N13E44 L=258)   440    225        0 ( 540/1483- 631/1660) -72.9 E46-E42 N07-N20
      4 1294 1319 (N18W07 L=310)   340    143        0 (1256/1262-1335/1384) -35.5 W03-W10 N16-N21
      5 1718  926 (N10W46 L=348)   571    141        1 (1674/ 868-1770/ 987)  66.0 W42-W48 N05-N14
      6 1820  859 (N10W58 L=  0)   797    164        1 (1794/ 791-1850/ 943)  90.0 W56-W59 N05-N15
      7  161  870 (S36E68 L=235)   949    109        0 ( 123/ 819- 200/ 920)  21.3 E80-E60 S39-S32
      8  316  872 (S33E45 L=257)   323     44        1 ( 288/ 853- 340/ 889)  -9.5 E48-E42 S33-S31
      9  541  727 (S36E23 L=280)   305     64        0 ( 505/ 715- 579/ 742)  18.4 E25-E19 S37-S34
     10 1805  664 (S00W66 L=  8)  1460    158        2 (1770/ 624-1854/ 735)  -6.3 W58-W76 S02-N02
     11  988  532 (S34W14 L=317)   507    132        1 ( 930/ 516-1057/ 561)  26.6 W10-W18 S36-S32
     12 1511  282 (S30W71 L= 13)  1335     95        3 (1471/ 256-1568/ 314)  18.4 W64-W80 S32-S27
     13 1349  237 (S39W62 L=  5)   555     66        0 (1324/ 219-1384/ 268)   0.0 W55-W67 S40-S37
  

    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