def psmapcurrent(cdata,chan='ch2',cmode='T',nbins=360):
"""
function to combine already read science files and pointing
files, and create pseudomap. chan is 'ch2', cmode is 'T'
"""
combined = combine_cofe_h5_pointing(cdata['dd'],cdata['pp'])
psmap=cu.phasebin(nbins,combined['az'],combined['sci_data'][chan][cmode])
pcolormesh(psmap.T)
return psmap
def psmapcurrent(cdata,chan='ch2',cmode='T',nbins=360):
"""
function to combine already read science files and pointing
files, and create pseudomap. chan is 'ch2', cmode is 'T'
"""
combined = combine_cofe_h5_pointing(cdata['dd'],cdata['pp'])
psmap=cu.phasebin(nbins,combined['az'],combined['sci_data'][chan][cmode])
pcolormesh(psmap.T)
return psmap
#saves combined file as well as one with phasebin subtracted from each channel to remove
# an azimuth template (first order subtraction of sidelobe effects))
import sys
sys.path.append('/smbmount/labuse/software_github_repos/cofe-python-analysis-tools/utils_meinhold/')
sys.path.append('/smbmount/labuse/software_github_repos/cofe-python-analysis-tools/utils_zonca/')
sys.path.append('/usr/local/lib/python2.7/site-packages/')
import cPickle
from glob import glob
import matplotlib
matplotlib.use('agg')
import cofe_util as cu
import realtime as rt
import numpy as np
yrmodays=['20130926','20130927','20130928','20130929','20130930','20131001','20131002','20131003']
for yrmoday in yrmodays:
#yrmoday=sys.argv[1]
cdata=rt.getdatanow(yrmoday)
chans=[]
for i in range(6):
chans.append('ch'+str(i))
for chan in chans:
for mode in ['T','Q','U']:
psmap=cu.phasebin(360,cdata['az'],cdata['sci_data'][chan][mode],degrees=True)
template=np.mean(psmap,axis=1)
alltemplate=np.interp(cdata['az'],np.arange(360),template)
cdata['sci_data'][chan][mode]=cdata['sci_data'][chan][mode]-alltemplate
cPickle.dump(cdata,open('combined_data_'+yrmoday+'_sub.pkl','w'),protocol=-1)
import pyfits
import cPickle
import cofe_util as util
get_ipython().magic(u"cd /cofe/flight_data/Level1/1.2")
get_ipython().magic(u"autocall 2")
c10=pyfits.open('all_10GHz_data_cal.fits')
p10=pyfits.open('all_10ghz_pointing.fits')
p10mag=pyfits.open('all_10ghz_pointing_mag.fits')
ut10=c10['Time'].data['ut']
utstartend=[16.8 ,38.4]
c101= c10['ch1_'].data['t'][(ut10 > utstartend[0]) & (ut10 < utstartend[1])]
az101mag=p10mag[1].data['az1'][(ut10 > utstartend[0]) & (ut10 < utstartend[1])]
#az101=p10[1].data['az1'][(ut10 > utstartend[0]) & (ut10 < utstartend[1])]
ps101mag=util.phasebin(3600,az101mag*np.pi/180.-np.pi,c101)
#ps101=util.phasebin(3600,az101*np.pi/180.-np.pi,c101)
#need to remove zeros:
ps101magbad=where(ps101mag == 0)
ps101mag[ps101magbad]=ps101mag[ps101magbad[0]-1,ps101magbad[1]]
ps101magbad=where(ps101mag == 0)
ps101mag[ps101magbad]=ps101mag[ps101magbad[0]-1,ps101magbad[1]]
ps101magbad=where(ps101mag == 0)
ps101mag[ps101magbad]=ps101mag[ps101magbad[0],ps101magbad[1]-1]
ps101magbad=where(ps101mag == 0)
ps101mag[ps101magbad]=ps101mag[ps101magbad[0],ps101magbad[1]-1]
ps101magbad=where(ps101mag == 0)
ps101mag[ps101magbad]=ps101mag[ps101magbad[0],ps101magbad[1]-1]
ps=copy(ps101mag)
#now set up to get cross correlations one row at at time
nrows=ps.shape[1]-1
#saves combined file as well as one with phasebin subtracted from each channel to remove
# an azimuth template (first order subtraction of sidelobe effects))
import sys
sys.path.append('/smbmount/labuse/software_github_repos/cofe-python-analysis-tools/utils_meinhold/')
sys.path.append('/smbmount/labuse/software_github_repos/cofe-python-analysis-tools/utils_zonca/')
sys.path.append('/usr/local/lib/python2.7/site-packages/')
import cPickle
from glob import glob
import matplotlib
matplotlib.use('agg')
import cofe_util as cu
import realtime as rt
import numpy as np
yrmodays=['20130926','20130927','20130928','20130929','20130930','20131001','20131002','20131003']
for yrmoday in yrmodays:
#yrmoday=sys.argv[1]
cdata=rt.getdatanow(yrmoday)
chans=[]
for i in range(6):
chans.append('ch'+str(i))
for chan in chans:
for mode in ['T','Q','U']:
psmap=cu.phasebin(360,cdata['az'],cdata['sci_data'][chan][mode],degrees=True)
template=np.mean(psmap,axis=1)
alltemplate=np.interp(cdata['az'],np.arange(360),template)
cdata['sci_data'][chan][mode]=cdata['sci_data'][chan][mode]-alltemplate
cPickle.dump(cdata,open('combined_data_'+yrmoday+'_sub.pkl','w'),protocol=-1)
import pyfits
import cPickle
import cofe_util as util
get_ipython().magic(u"cd /cofe/flight_data/Level1/1.2")
get_ipython().magic(u"autocall 2")
c10 = pyfits.open('all_10GHz_data_cal.fits')
p10 = pyfits.open('all_10ghz_pointing.fits')
p10mag = pyfits.open('all_10ghz_pointing_mag.fits')
ut10 = c10['Time'].data['ut']
utstartend = [16.8, 38.4]
c101 = c10['ch1_'].data['t'][(ut10 > utstartend[0]) & (ut10 < utstartend[1])]
az101mag = p10mag[1].data['az1'][(ut10 > utstartend[0])
& (ut10 < utstartend[1])]
#az101=p10[1].data['az1'][(ut10 > utstartend[0]) & (ut10 < utstartend[1])]
ps101mag = util.phasebin(3600, az101mag * np.pi / 180. - np.pi, c101)
#ps101=util.phasebin(3600,az101*np.pi/180.-np.pi,c101)
#need to remove zeros:
ps101magbad = where(ps101mag == 0)
ps101mag[ps101magbad] = ps101mag[ps101magbad[0] - 1, ps101magbad[1]]
ps101magbad = where(ps101mag == 0)
ps101mag[ps101magbad] = ps101mag[ps101magbad[0] - 1, ps101magbad[1]]
ps101magbad = where(ps101mag == 0)
ps101mag[ps101magbad] = ps101mag[ps101magbad[0], ps101magbad[1] - 1]
ps101magbad = where(ps101mag == 0)
ps101mag[ps101magbad] = ps101mag[ps101magbad[0], ps101magbad[1] - 1]
ps101magbad = where(ps101mag == 0)
ps101mag[ps101magbad] = ps101mag[ps101magbad[0], ps101magbad[1] - 1]
ps = copy(ps101mag)
#now set up to get cross correlations one row at at time
nrows = ps.shape[1] - 1
lags = []
