Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  890 1897 (N52E39 L=292)   934     99        1 ( 850/1873- 938/1919)  -8.7 E46-E32 N49-N53
      2  650 1798 (N38E50 L=281)   838     89        4 ( 626/1754- 675/1849) -42.0 E56-E44 N33-N42
      3 1216 1452 (N22W04 L=335)   300    112        1 (1175/1415-1260/1481) -23.2 W00-W07 N20-N23
      4 1854 1088 (N15W57 L= 28)   657     78        2 (1820/1060-1878/1132)  51.3 W53-W60 N12-N17
      5  417  942 (S20E38 L=293)   606    180        1 ( 361/ 882- 467/1002)  48.8 E43-E34 S24-S16
      6 1701  651 (S14W54 L= 26)   353     51        1 (1673/ 637-1734/ 668)  31.0 W51-W57 S15-S12
      7  332  583 (S45E52 L=279)   632     46        2 ( 311/ 553- 353/ 617)  38.7 E56-E48 S46-S42
      8 1695  561 (S19W59 L= 30)  1083    183        3 (1620/ 511-1764/ 607)  32.5 W51-W64 S22-S14
      9 1116  542 (S33W16 L=347)   460    200        1 (1003/ 526-1194/ 568)  22.6 W09-W21 S36-S30
     10 1629  479 (S25W58 L= 29)   700     85        1 (1600/ 445-1678/ 501)  51.3 W54-W62 S27-S22
     11 1037  246 (S57W26 L=358)   424     62        0 (1007/ 229-1069/ 266)  18.4 W23-W30 S58-S55
  

    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