Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  640 1395 (N27E21 L=303)   316     89        1 ( 605/1370- 685/1411) -26.6 E23-E17 N25-N28
      2  735 1262 (N17E15 L=309)  1757    433        6 ( 603/1139- 833/1394)  57.3 E24-E08 N11-N25
      3  220 1234 (N22E55 L=269)   613    200        2 ( 157/1204- 265/1259)   4.8 E62-E50 N21-N23
      4 1348  836 (S15W18 L=341)   675    218        0 (1280/ 746-1404/ 951)  47.7 W12-W22 S20-S09
      5  912  780 (S13E10 L=314)   186     84        1 ( 881/ 757- 966/ 797) -11.3 E11-E06 S14-S12
      6  228  677 (S11E64 L=260)  3710    565       11 ( 162/ 578- 284/ 785)  38.7 E76-E56 S16-S04
      7 1734  669 (S31W47 L= 11)  1547    230        2 (1655/ 612-1784/ 760)  60.9 W40-W52 S34-S24
      8  700  615 (S21E26 L=298)   993    238        2 ( 631/ 524- 749/ 719)  42.5 E33-E21 S25-S14
      9 1523  592 (S33W30 L=354)   308     32        0 (1501/ 568-1544/ 617)  45.0 W27-W31 S35-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