Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1301 1462 (N22W26 L=141)   277     62        1 (1273/1444-1334/1480) -56.3 W24-W28 N20-N23
      2  427 1354 (N32E36 L= 79)   274     96        1 ( 386/1326- 471/1372)   7.1 E39-E31 N30-N32
      3 1890 1020 (S15W63 L=178)   897    169        1 (1832/ 957-1921/1087)  23.2 W55-W69 S17-S11
      4  183  827 (N04E65 L= 50)   548     64        1 ( 160/ 799- 208/ 857)  20.6 E69-E60 N02-N05
      5  857  759 (S11E15 L=100)   680    222        0 ( 805/ 660- 917/ 840)  36.0 E20-E09 S15-S06
      6  331  699 (S06E53 L= 62)   398     51        1 ( 306/ 677- 359/ 727)  26.6 E55-E50 S06-S03
      7 1742  642 (S36W48 L=163)  1338    260        0 (1667/ 576-1847/ 722) -24.3 W38-W69 S43-S28
      8  350  567 (S14E57 L= 58)  1220    151        2 ( 309/ 519- 397/ 640)  23.2 E64-E50 S16-S09
  

    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