Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1021 1539 (N31W07 L=158)   536    159        0 ( 960/1487-1080/1592) -48.8 W03-W10 N27-N35
      2 1867 1312 (N07W68 L=219)   979     75        2 (1842/1282-1902/1350) -33.7 W64-W74 N04-N09
      3  384 1245 (N21E41 L=110)  1483    240        3 ( 291/1186- 485/1294)   3.4 E49-E32 N17-N23
      4 1874  909 (S17W62 L=213)  4427    348       15 (1791/ 762-1931/1046) -83.3 W52-W72 S26-S09
      5 1586  883 (S15W34 L=185)  1519    370        1 (1542/ 713-1617/1045)  81.5 W31-W36 S25-S05
      6 1474  769 (S21W26 L=177)   249     92        0 (1446/ 730-1496/ 812)  51.3 W23-W27 S23-S18
      7 1708  596 (S36W49 L=200)   958    159        0 (1644/ 542-1769/ 638)  22.6 W42-W55 S38-S33
      8  589  580 (S22E35 L=116)   319     87        0 ( 571/ 537- 619/ 614)  38.7 E37-E32 S24-S19
      9  305  571 (S18E61 L= 90)   467     59        3 ( 288/ 544- 326/ 602) -90.0 E63-E58 S20-S16
     10 1518  487 (S41W31 L=182)   364     87        0 (1468/ 472-1559/ 504)  -6.3 W26-W35 S41-S39
  

    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