Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  423 1590 (N25E52 L=153)   446     70        1 ( 388/1566- 459/1612)  31.0 E55-E49 N23-N27
      2  764 1422 (N16E20 L=185)   947    274        2 ( 717/1324- 853/1528) -52.4 E24-E13 N10-N23
      3 1831 1234 (N18W57 L=262)   369     34        1 (1807/1219-1851/1255)  18.4 W53-W59 N16-N19
      4 1694  916 (S04W45 L=250)   394     70        0 (1674/ 879-1716/ 952) -63.4 W43-W46 S05-S01
      5  699  740 (S25E18 L=187)   526    231        0 ( 610/ 697- 802/ 773)  17.1 E24-E11 S28-S23
      6  334  603 (S37E53 L=152)   719    156        0 ( 286/ 542- 384/ 664) -53.1 E56-E48 S40-S32
      7  877  285 (S57E03 L=202)  2561    442        6 ( 699/ 256-1084/ 317)   4.2 E24-W18 S60-S52
  

    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