Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  549  682 (S24E33 L=217)   734    154        2 ( 515/ 630- 601/ 744) -76.0 E35-E29 S27-S19
      2 1755  689 (S22W55 L=305)   493     71        0 (1725/ 664-1776/ 725)  56.3 W52-W57 S23-S20
      3 1709  566 (S31W56 L=305)   545     51        1 (1680/ 542-1740/ 582)   0.0 W52-W59 S32-S29
      4  420  388 (S43E60 L=189)   460     64        0 ( 391/ 363- 447/ 413)  11.3 E68-E53 S44-S41
      5  692  313 (S51E33 L=216)   516     67        1 ( 668/ 285- 718/ 349)  40.6 E37-E29 S53-S47
      6  802  314 (S51E22 L=228)  1025    178        1 ( 729/ 289- 880/ 343) -12.1 E28-E13 S53-S48
      7  610  310 (S51E43 L=206)   476     53        1 ( 587/ 290- 640/ 335)   7.1 E48-E39 S52-S48
  

    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