Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  491 1684 (N40E45 L=143)   781    146        2 ( 446/1658- 536/1705)  -6.3 E50-E40 N37-N41
      2  803 1688 (N40E18 L=170)   304     76        0 ( 769/1666- 830/1705) -33.7 E21-E15 N37-N41
      3  233 1209 (N08E56 L=132)   522    113        1 ( 218/1173- 263/1254) -51.3 E58-E52 N05-N10
      4  175 1127 (N03E62 L=126)  1992    295        6 ( 135/1029- 213/1214) -73.3 E68-E57 S02-N08
      5 1881 1038 (N00W63 L=251)   801    134        1 (1847/ 974-1916/1087) -32.0 W58-W67 S03-N03
      6 1013  856 (S14E01 L=188)   349    116        1 ( 978/ 790-1051/ 914)  66.8 E02-W01 S17-S10
      7  522  560 (S32E36 L=152)   414     93        2 ( 480/ 544- 565/ 581)   0.0 E39-E32 S33-S31
  

    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