Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1219 1820 (N51E11 L= 83)   488     24        0 (1192/1802-1247/1837)  36.9 E13-E09 N48-N53
      2  791 1748 (N32E37 L= 57)   527     77        0 ( 763/1709- 823/1785) -29.1 E41-E32 N29-N34
      3  422 1651 (N18E61 L= 33)  1252    110        2 ( 381/1616- 480/1691)  18.4 E68-E54 N15-N21
      4  764 1600 (N21E31 L= 63)   368     41        0 ( 739/1575- 792/1625)  76.0 E32-E29 N19-N23
      5  186  897 (S31E63 L= 31)   446     65        0 ( 170/ 864- 197/ 932) -90.0 E65-E60 S32-S28
      6 1548  375 (S29W61 L=155)   886    139        0 (1499/ 318-1602/ 422) -77.5 W58-W63 S33-S24
      7  971  234 (S59W30 L=124)  1337    208        2 ( 903/ 156-1035/ 313) -14.5 W17-W48 S62-S54
  

    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