Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1612 1272 (N06W43 L=151)   387     91        1 (1592/1225-1632/1306)  45.0 W39-W45 N03-N07
      2  230 1207 (N34E54 L= 54)   378     51        0 ( 214/1188- 252/1230) -45.0 E55-E51 N32-N35
      3  247  990 (N21E54 L= 54)   480    142        1 ( 216/ 959- 280/1024) -45.0 E57-E51 N18-N22
      4 1800  744 (S32W47 L=156)  2321    380        3 (1696/ 621-1906/ 875)  21.4 W36-W60 S37-S27
      5 1616  664 (S30W30 L=138)   682    135        1 (1578/ 606-1650/ 737)  54.5 W26-W33 S33-S26
      6  304  580 (S05E63 L= 45)   655    103        2 ( 267/ 542- 332/ 621) -71.6 E65-E59 S07-S01
      7  646  350 (S26E45 L= 64)   760    222        2 ( 602/ 271- 700/ 427)  72.9 E47-E40 S32-S19
      8  717  222 (S37E49 L= 59)   519    147        2 ( 685/ 202- 750/ 260)  21.8 E53-E43 S38-S34
  

    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