Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  632 1744 (N53E27 L=162)   429     65        0 ( 609/1717- 659/1765)  45.0 E29-E24 N50-N54
      2  472 1557 (N39E37 L=153)   389     80        0 ( 438/1539- 509/1570)  -9.5 E39-E33 N37-N40
      3 1628 1504 (N22W48 L=238)   896    238        1 (1523/1483-1705/1523) -11.3 W39-W54 N20-N23
      4  237 1328 (N26E55 L=135)   462    105        1 ( 217/1302- 274/1352) -18.4 E57-E50 N24-N27
      5  335  789 (S08E49 L=141)   292     79        0 ( 302/ 774- 367/ 804) -11.3 E51-E45 S08-S06
      6  698  592 (S24E26 L=163)   304     63        1 ( 666/ 570- 735/ 614)  11.3 E28-E23 S24-S22
      7  788  336 (S43E29 L=161)   344     70        0 ( 767/ 308- 812/ 367) -90.0 E30-E27 S45-S40
  

    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