Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1223 1633 (N43W24 L= 95)   569    152        1 (1174/1590-1276/1660)  15.3 W18-W29 N40-N44
      2  468 1335 (N28E38 L= 33)  2427    762        9 ( 254/1244- 737/1415) -15.6 E59-E16 N21-N33
      3  206  955 (N05E62 L=  9)   510    120        3 ( 173/ 936- 243/ 979)   6.3 E66-E57 N03-N06
      4  303  894 (N01E52 L= 19)   461     87        2 ( 270/ 870- 330/ 926)  26.6 E55-E48 N00-N03
      5 1650  800 (S15W41 L=112)   444    117        0 (1603/ 769-1696/ 826)  14.0 W36-W44 S16-S13
      6 1312  287 (S50W18 L= 89)   366     57        1 (1290/ 267-1333/ 310) -26.6 W15-W21 S52-S47
      7  871  214 (S52E26 L= 45)  1037    211        1 ( 772/ 193- 970/ 248) -12.5 E34-E15 S55-S47
  

    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