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()
