combinecbass

0001 clear
0002 indir = '/home/cdickins/cbass_analysis/fits_files/'
0003 % file{1}='atest80-512.fits'
0004 % file{2}='atest81-512.fits'
0005 % file{3}='atest84-512.fits'
0006 % file{4}='atest86-512.fits'
0007 % file{5}='atest87-512.fits'
0008 % file{6}='atest88-512.fits'
0009 % file{7}='atest89-512.fits'
0010 % file{8}='atest90-512.fits'
0011 % file{9}='atest91-512.fits'
0012 % file{10}='atest92-512.fits'
0013 % file{11}='atest93-512.fits'
0014 % file{12}='atest95-512.fits'
0015 % file{13}='atest97-512.fits'
0016 % file{14}='atest98-512.fits'
0017 % file{15}='atest99-512.fits'
0018 
0019 
0020 % numscans=24;
0021 % nside = 256
0022 % file{1}='atest80-256.fits'
0023 % file{2}='atest81-256.fits'
0024 % file{3}='atest84-256.fits'
0025 % file{4}='atest86-256.fits'
0026 % file{5}='atest87-256.fits'
0027 % file{6}='atest88-256.fits'
0028 % file{7}='atest89-256.fits'
0029 % file{8}='atest90-256.fits'
0030 % file{9}='atest91-256.fits'
0031 % file{10}='atest92-256.fits'
0032 % file{11}='atest93-256.fits'
0033 % file{12}='atest95-256.fits'
0034 % file{13}='atest97-256.fits'
0035 % file{14}='atest98-256.fits'
0036 % file{15}='atest99-256.fits'
0037 % file{16}='atest101-256.fits'
0038 % file{17}='atest102-256.fits'
0039 % file{18}='atest103-256.fits'
0040 % file{19}='atest104-256.fits'
0041 % file{20}='atest105-256.fits'
0042 % file{21}='atest106-256.fits'
0043 % file{22}='atest107-256.fits'
0044 % file{23}='atest108-256.fits'
0045 % file{24}='atest109-256.fits'
0046 % file{25}='atest110-256.fits'
0047 % file{26}='atest111-256.fits'
0048 % file{27}='atest112-256.fits'
0049 % file{28}='atest113-256.fits'
0050 % file{29}='atest114-256.fits'
0051 
0052 numscans=29;
0053 nside = 128
0054 file{1}='atest80-128.fits'
0055 file{2}='atest81-128.fits'
0056 file{3}='atest84-128.fits'
0057 file{4}='atest86-128.fits'
0058 file{5}='atest87-128.fits'
0059 file{6}='atest88-128.fits'
0060 file{7}='atest89-128.fits'
0061 file{8}='atest90-128.fits'
0062 file{9}='atest91-128.fits'
0063 file{10}='atest92-128.fits'
0064 file{11}='atest93-128.fits'
0065 file{12}='atest97-128.fits'
0066 file{13}='atest98-128.fits'
0067 file{14}='atest99-128.fits'
0068 file{15}='atest99-128.fits'
0069 file{16}='atest101-128.fits'
0070 file{17}='atest102-128.fits'
0071 file{18}='atest103-128.fits'
0072 file{19}='atest104-128.fits'
0073 file{20}='atest105-128.fits'
0074 file{21}='atest106-128.fits'
0075 file{22}='atest107-128.fits'
0076 file{23}='atest108-128.fits'
0077 file{24}='atest110-128.fits'
0078 file{25}='atest111-128.fits'
0079 file{26}='atest112-128.fits'
0080 file{27}='atest112-128.fits'
0081 file{28}='atest113-128.fits'
0082 file{29}='atest114-128.fits'
0083 file{30}='atest118-128.fits'
0084 file{31}='atest119-128.fits'
0085 file{32}='atest120-128.fits'
0086 
0087 
0088 unix(['rm combine-',num2str(nside),'.fits']);
0089 
0090 
0091 sum=[];
0092 f=[];
0093 g=[];
0094 
0095 h = gov.fnal.eag.healpix.PixTools;
0096 npix= h.Nside2Npix(nside);
0097 sum=zeros(npix,7);
0098 
0099 for k=1:numscans
0100     for i=1:11
0101         infile = strcat(indir,file{k})
0102         f(:,i) = cbass_read_bin_col(infile,2,i); 
0103     end
0104     %zero the additions for certain conditions
0105     indI=find(f(:,1)>=8); %get rid of RFI!
0106     f(indI,6)=0;
0107     indI=find(f(:,6)<=4); %get rid of low sampled pixels
0108     f(indI,6)=0;
0109    % indI=find((f(:,1)./f(:,7))<=50); %Signal to noise ratio
0110    % f(indI,6)=0;
0111     for i=1:5
0112         g(:,i) = f(:,i).*f(:,6);
0113     end
0114     g(:,6)=f(:,6);
0115     g(:,7)=f(:,6).*f(:,7); %weighted sigma
0116     indI = find(isnan(g(:,7)));
0117     g(indI,7)=0;
0118     sum=sum+g;
0119 end
0120 %
0121 % for i=1:6
0122 %     f(:,i) = cbass_read_bin_col(file{2},2,i);
0123 % end
0124 %
0125 % for i=1:5
0126 %     g(:,i) = f(:,i).*f(:,6);
0127 % end
0128 % g(:,6)=f(:,6);
0129 % sum = sum+g;
0130 %
0131 % for i=1:6
0132 %     f(:,i) = cbass_read_bin_col(file{3},2,i);
0133 % end
0134 %
0135 % for i=1:5
0136 %     g(:,i) = f(:,i).*f(:,6);
0137 % end
0138 % g(:,6)=f(:,6);
0139 % sum = sum+g;
0140 %
0141 % for i=1:6
0142 %     f(:,i) = cbass_read_bin_col(file{4},2,i);
0143 % end
0144 %
0145 % for i=1:5
0146 %     g(:,i) = f(:,i).*f(:,6);
0147 % end
0148 % g(:,6)=f(:,6);
0149 % sum = sum+g;
0150 
0151 
0152 for i=1:5
0153     meansum(:,i)=sum(:,i)./sum(:,6);
0154 end
0155 meansum(:,6)=sum(:,6);
0156 meansum(:,7) =sum(:,7)./sum(:,6); %total fractional uncertainty/number of measurements->fractional uncertainty
0157 %meansum(:,8) = meansum(:,7).*meansum(:,1); %uncertainty of the pixel using fractional uncertainty*measurement
0158 %indA = find(meansum(:,1)./meansum(:,7))<=0.01); %find the low SNR pixels
0159 indA = find(meansum(:,6)<=50); %find undersampled points
0160 meansum(indA,1:5)=NaN;
0161 meansum(:,1)=meansum(:,1)-5;
0162 
0163 %meansum(:,1)=sum(:,7)./sum(:,6);
0164 %meansum(indA,:)=0; %and then zero them
0165 %meansum(find(isnan(meansum))) = -1;
0166 
0167  cbass_write_image((['combine-',num2str(nside),'.fits']),meansum,nside,npix,'test','test2')
0168