Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1463 1740 (N51W25 L= 75)   388     25        1 (1443/1717-1489/1759)   9.5 W22-W28 N48-N53
      2  337 1603 (N22E66 L=345)  2087    188        3 ( 267/1552- 410/1661)  36.0 E74-E59 N18-N26
      3  627 1491 (N16E33 L= 17)   365    113        1 ( 613/1445- 646/1543) -63.4 E34-E30 N13-N19
      4 1878 1131 (N16W61 L=111)  3695    306        8 (1795/1043-1945/1225) -17.5 W50-W72 N10-N21
      5  592  863 (S22E25 L= 25)   502    149        1 ( 534/ 842- 661/ 893)  24.0 E29-E20 S23-S18
      6 1417  365 (S39W47 L= 97)   455     50        1 (1390/ 342-1438/ 400) -76.0 W44-W50 S40-S36
      7  471  370 (S58E48 L=  2)   943     98        1 ( 423/ 351- 527/ 395)   0.0 E56-E38 S59-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