mosaic_map_plot

 [moss,mosl]=mosaic_map_plot(ad,ads,o,p,beam,mrk)

 Make and plot mosaic maps for short and long base

0001 function [moss,mosl]=mosaic_map_plot(ad,ads,o,p,beam,mrk)
0002 % [moss,mosl]=mosaic_map_plot(ad,ads,o,p,beam,mrk)
0003 %
0004 % Make and plot mosaic maps for short and long base
0005 
0006 if(~exist('mrk'))
0007   mrk=[];
0008 end
0009 
0010 % make mosaic for uv range 0-2000
0011 moss=mosaic_map2(ad,ads,p,beam,[0,2000]);
0012 % make mosaic for uv range 2000-inf
0013 mosl=mosaic_map2(ad,ads,p,beam,[2000,Inf]);
0014 
0015 % NaN all map points where the noise is >2 times the
0016 % minimum noise
0017 moss.signal(moss.noise>2*min(moss.noise(:)))=NaN;
0018 mosl.signal(mosl.noise>2*min(mosl.noise(:)))=NaN;
0019 
0020 moss.sigma(moss.noise>2*min(moss.noise(:)))=NaN;
0021 mosl.sigma(mosl.noise>2*min(mosl.noise(:)))=NaN;
0022 
0023 setwinsize(gcf,900,700); clf
0024 
0025 subplot(2,3,1)
0026 imagesc(ad.t_val_deg,ad.t_val_deg,moss.signal);
0027 skyplot(' '); colorbar('horiz')
0028 mark_field_cen(p,mrk);
0029 title('Jy/beam short baseline'); ylabel('deg on sky')
0030 zoom(0.7)
0031 subplot(2,3,2)
0032 imagesc(ad.t_val_deg,ad.t_val_deg,moss.sigma);
0033 skyplot(' '); colorbar('horiz')
0034 mark_field_cen(p,mrk);
0035 title('Significance short baseline'); ylabel('deg on sky')
0036 zoom(0.7)
0037 subplot(2,3,3)
0038 [bc,n]=hfill(moss.sigma,100); hplot(bc,n,'L');
0039 hold on; plot(bc,gauss([max(n),0,1],bc),'r'); hold off;
0040 title('Pixel val. dist. short baseline')
0041 xlabel('sigma'); ylabel('count per bin')
0042 
0043 subplot(2,3,4)
0044 imagesc(ad.t_val_deg,ad.t_val_deg,mosl.signal);
0045 skyplot(' '); colorbar('horiz')
0046 mark_field_cen(p,mrk);
0047 title('Jy/beam long baseline'); ylabel('deg on sky')
0048 zoom(1.0)
0049 subplot(2,3,5)
0050 imagesc(ad.t_val_deg,ad.t_val_deg,mosl.sigma);
0051 skyplot(' '); colorbar('horiz')
0052 mark_field_cen(p,mrk);
0053 title('Significance long baseline'); ylabel('deg on sky')
0054 zoom(1.0)
0055 subplot(2,3,6)
0056 [bc,n]=hfill(mosl.sigma,100); hplot(bc,n,'L');
0057 hold on; plot(bc,gauss([max(n),0,1],bc),'r'); hold off;
0058 title('Pixel val. dist. long baseline')
0059 xlabel('sigma'); ylabel('count per bin')
0060 
0061 return
0062 
0063 %%%%%%%%%%%%%%%%%%%%%%%%%%
0064 function mark_field_cen(p,mrk)
0065 
0066 for i=1:length(p.px)
0067   text(p.px(i),p.py(i),num2str(p.n(i)),'HorizontalAlignment','center','VerticalAlignment','middle','Color','red');
0068 end
0069 
0070 if(~isempty(mrk))
0071   for i=1:length(mrk.x)
0072     hold on; plot(mrk.x(i),mrk.y(i),mrk.sym{i}); hold off
0073   end
0074 end
0075 
0076 return
0077