Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  807 1689 (N30E35 L=230)  1276    217        2 ( 724/1618- 880/1741) -10.6 E42-E26 N26-N33
      2  634 1649 (N23E44 L=221)  1488    231        4 ( 548/1560- 706/1724)  77.0 E47-E40 N15-N29
      3 1303 1423 (N26W06 L=271)  2230    478        3 (1188/1207-1364/1596) -48.0 E05-W14 N14-N35
      4  879 1395 (N13E18 L=247)   207     65        1 ( 852/1367- 907/1423) -26.6 E20-E16 N12-N14
      5 1191 1215 (N11W05 L=270)   296    157        0 (1112/1167-1240/1262) -14.0 W00-W08 N08-N12
      6 1620 1047 (N13W35 L=300)   387     97        0 (1606/1001-1632/1105) -59.0 W33-W36 N10-N15
      7  281 1057 (S19E51 L=214)   662    130        0 ( 233/1015- 331/1086)  33.7 E55-E45 S22-S16
      8  833  998 (S10E10 L=255)   307    151        0 ( 809/ 944- 891/1034)  63.4 E11-E07 S13-S07
      9  647  924 (S19E20 L=245)   477    220        0 ( 613/ 862- 700/ 981) -39.8 E22-E15 S21-S14
     10 1777  834 (N06W54 L=319)  1046    220        1 (1714/ 754-1816/ 935)  84.8 W50-W57 N00-N11
     11  281  918 (S27E49 L=216)   462     91        1 ( 241/ 904- 334/ 933)   5.7 E53-E43 S29-S25
     12  608  456 (S50E14 L=251)  1683    278        2 ( 507/ 388- 720/ 547) -10.3 E25-E02 S53-S46
  

    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