casa6
Public
Commits
| 978 978 | beamcs = ia.coordsys() |
| 979 979 | beam_array = ia.getchunk(axes=[beamcs.findcoordinate("spectral")['pixel'][0],beamcs.findcoordinate("stokes")['pixel'][0]],dropdeg=True) |
| 980 980 | nn = beam_array.shape |
| 981 981 | xextent = nn[0]*cell_asec*0.5 |
| 982 982 | xextent = [xextent,-xextent] |
| 983 983 | yextent = nn[1]*cell_asec*0.5 |
| 984 984 | yextent = [-yextent,yextent] |
| 985 985 | flipped_array = beam_array.transpose() |
| 986 986 | ttrans_array = flipped_array.tolist() |
| 987 987 | ttrans_array.reverse() |
| 988 - | pl.imshow(ttrans_array,interpolation='bilinear',cmap=pl.cm.jet,extent=xextent+yextent,origin="bottom") |
| 988 + | pl.imshow(ttrans_array,interpolation='bilinear',cmap=pl.cm.jet,extent=xextent+yextent,origin="lower") |
| 989 989 | psfim.replace(project+"/","") |
| 990 990 | pl.title(psfim,fontsize="x-small") |
| 991 991 | try: |
| 992 992 | b = qa.convert(beam['major'],'arcsec')['value'] |
| 993 993 | pl.xlim([-3*b,3*b]) |
| 994 994 | pl.ylim([-3*b,3*b]) |
| 995 995 | ax = pl.gca() |
| 996 996 | pl.text(0.05,0.95,"bmaj=%7.1e\nbmin=%7.1e" % (beam['major']['value'], |
| 997 997 | beam['minor']['value']), |
| 998 998 | transform = ax.transAxes,bbox=dict(facecolor='white', alpha=0.7),size="x-small",verticalalignment="top") |