Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1330 1672 (N35W45 L=158)  2710    458        7 (1134/1603-1524/1704) -14.8 W27-W61 N28-N39
      2 1121 1516 (N32W21 L=135)   678    199        3 (1084/1433-1167/1602)  65.6 W18-W25 N26-N37
      3  455 1072 (N23E35 L= 79)   267     72        0 ( 424/1050- 490/1099) -38.7 E37-E32 N21-N25
      4 1863 1019 (S18W58 L=172)   401     94        0 (1847/ 981-1875/1063)  53.1 W55-W60 S20-S15
      5 1684  877 (S20W37 L=151)   293     69        0 (1667/ 841-1700/ 907)  36.9 W35-W39 S21-S18
      6 1391  855 (S12W17 L=130)   403    136        2 (1333/ 825-1463/ 886) -35.5 W13-W20 S15-S09
      7  229  818 (N12E61 L= 52)   493     84        2 ( 197/ 780- 259/ 854) -59.0 E64-E57 N09-N14
      8 1724  730 (S31W40 L=153)  2249    371        7 (1628/ 589-1847/ 852)  21.4 W29-W54 S37-S25
      9 1614  800 (S22W31 L=144)   321     79        0 (1572/ 787-1653/ 815) -26.6 W27-W33 S23-S20
     10 1271  665 (S20W05 L=119)   421    132        1 (1244/ 601-1295/ 725)  54.5 W01-W07 S23-S16
     11  581  423 (S18E47 L= 67)   743    221        2 ( 546/ 355- 628/ 491)  56.3 E51-E42 S22-S13
     12  638  324 (S26E49 L= 65)   357     53        0 ( 614/ 301- 661/ 351)  23.2 E52-E45 S27-S24
  

    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