Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1125 1833 (N49E22 L=153)   571     91        0 (1075/1810-1166/1853)   0.0 E26-E17 N47-N50
      2  615 1703 (N24E48 L=126)   711    107        0 ( 590/1647- 639/1762) -56.3 E52-E45 N19-N28
      3 1315 1366 (N20W08 L=183)   278     92        1 (1285/1330-1350/1405) -26.6 W05-W11 N18-N21
      4  903 1341 (N08E15 L=160)   265     27        1 ( 878/1323- 924/1361)  56.3 E16-E13 N06-N09
      5 1246 1217 (N09W08 L=182)   472    155        0 (1225/1151-1268/1297) -53.1 W04-W10 N05-N13
      6 1741 1016 (N12W45 L=219)   402     76        0 (1721/ 983-1756/1050) -56.3 W43-W46 N09-N14
      7  179  937 (S30E64 L=111)  1298    168        1 ( 114/ 906- 266/ 966)  11.8 E76-E52 S31-S26
      8  520  773 (S34E27 L=147)   430    151        0 ( 471/ 731- 561/ 839)  56.3 E30-E24 S37-S28
      9  246  802 (S37E56 L=119)   589    112        0 ( 198/ 780- 300/ 829)  18.4 E63-E48 S39-S34
     10  443  769 (S35E34 L=141)   269     51        0 ( 424/ 737- 463/ 800) -63.4 E35-E31 S37-S33
     11 1279  700 (S18W24 L=198)   398    103        4 (1226/ 679-1329/ 720)  26.6 W20-W26 S20-S16
     12  777  401 (S52W02 L=176)   536    109        0 ( 732/ 370- 823/ 442) -12.5 E02-W07 S53-S49
  

    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