gains=numpy.conj(calxSB) * numpy.transpose(calxSB) polAcc=numpy.multiply(polAcc,gains) earth_radius=(xyz[0,0]**2+xyz[0,1]**2+xyz[0,2]**2)**.5 for fid,fn in enumerate(acc_files): for sid,sb in enumerate(sbs): obs.date=ephemTime[fid,sid] src=ephem.FixedBody() src._ra=obs.sidereal_time() src._dec=obs.lat src.compute(obs) uvw=lss.xyz2uvw(xyz, src, obs, freqs[sid]) #uvw plot #pylab.plot(uvw[:,:,0],uvw[:,:,1],'.') stokesAcc=numpy.array([xxAcc[:,:,fid,sid],xyAcc[:,:,fid,sid],yxAcc[:,:,fid,sid],yyAcc[:,:,fid,sid]]) xxIm=lss.fftImage(stokesAcc[0],uvw,px,res,mask=False,conv='fast') xyIm=lss.fftImage(stokesAcc[1],uvw,px,res,mask=False,conv='fast') yxIm=lss.fftImage(stokesAcc[2],uvw,px,res,mask=False,conv='fast') yyIm=lss.fftImage(stokesAcc[3],uvw,px,res,mask=False,conv='fast') iIm=(xxIm+yyIm).real qIm=(xxIm-yyIm).real uIm=(xyIm+yxIm).real vIm=(xyIm-yxIm).imag pylab.subplot(2,2,1) pylab.imshow(iIm) pylab.xlabel('Pixels (E-W)') pylab.ylabel('Pixels (N-S)') pylab.title('I') pylab.colorbar() pylab.subplot(2,2,2)
##instpolMat=numpy.zeros((px[0],px[0],2,2),dtype=complex) ##for xCoord in range(px[0]): ## for yCoord in range(px[1]): ## bMat=numpy.matrix([[xxIm[xCoord,yCoord],xyIm[xCoord,yCoord]],[yxIm[xCoord,yCoord],yyIm[xCoord,yCoord]]]) ## U,s,V=numpy.linalg.svd(bMat) ## J=numpy.dot(U,numpy.dot(numpy.eye(2)*numpy.sqrt(s),V)) ## instpolMat[xCoord,yCoord]=J ###remove system noise factor ##noiseFactor=numpy.sum((instpolMat[:,:,0,0][imMask]+instpolMat[:,:,1,1][imMask])/2.)/len(instpolMat[:,:,0,0][imMask]) ##instpolMat-=noiseFactor*numpy.eye(2) ##instpolMat/=numpy.max(instpolMat[:,:,0,0]) #normalization pylab.figure() #pylab.imshow((xxIm+yyIm).real,cmap=pylab.get_cmap('gray')) ##Qx = pylab.quiver(instpolMat[:,:,0,0],instpolMat[:,:,0,1],color='b',pivot='middle') ##Qy = pylab.quiver(instpolMat[:,:,1,0],instpolMat[:,:,1,1],color='r',pivot='middle') ##iIm=(xxIm+yyIm).real ##iIm[imMask]=0. ##pylab.imshow(iIm,cmap=pylab.get_cmap('gray')) #pylab.colorbar() xxIm=lss.fftImage(xxAcc[:,:,fid,sid],uvw,px,res) pylab.imshow(xxIm.real) pylab.colorbar() if opts.savefig: pylab.savefig('dft_s%i.png'%sbs[0]) else: pylab.show()