Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  655 1742 (N52E25 L= 69)   869     62        1 ( 598/1703- 712/1775)  24.0 E30-E20 N48-N54
      2  275 1462 (N34E54 L= 40)   699    117        3 ( 253/1425- 298/1496) -51.3 E57-E51 N31-N36
      3  342 1436 (N31E47 L= 48)   489    108        1 ( 311/1399- 373/1456) -33.7 E49-E43 N28-N32
      4  817 1334 (N19E10 L= 84)   252    109        0 ( 766/1289- 846/1361) -38.7 E13-E08 N16-N21
      5 1484 1281 (N10W31 L=126)   472    105        1 (1428/1235-1532/1328) -39.8 W28-W34 N07-N13
      6  589  812 (S10E29 L= 65)   246     43        1 ( 566/ 788- 614/ 836) -33.7 E30-E27 S11-S08
      7  475  467 (S30E47 L= 47)  2599    546        2 ( 375/ 365- 558/ 598)  28.4 E62-E37 S36-S22
      8 1695  516 (S40W49 L=144)  1306    264        1 (1593/ 458-1782/ 579)  21.8 W40-W56 S43-S35
  

    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