0001 ad=calc_ad(256*2/60,256); 0002 0003 [dm.x,dm.y,dm.dmap]=fitsget('../make_dmap_test/map_test.fits'); 0004 dm.dmap=fliplr(dm.dmap); 0005 figure(1) 0006 imagesc(ad.t_val_deg,ad.t_val_deg,dm.dmap); skyplot; colorbar 0007 0008 load('../make_dmap_test/map_test'); 0009 0010 figure(2) 0011 tic; [dmap,mapvar]=make_dmap(ad,d.u,d.v,d.vis,d.var); toc 0012 imagesc(ad.t_val_deg,ad.t_val_deg,dmap); skyplot; colorbar 0013 0014 figure(3) 0015 tic; [dmap2,mapvar2]=difmap(ad,d.u,d.v,d.vis,d.var); toc 0016 imagesc(ad.t_val_deg,ad.t_val_deg,dmap2); skyplot; colorbar 0017 0018 [sqrt(mapvar) sqrt(mapvar2)] 0019 0020 figure(4) 0021 imagesc(ad.t_val_deg,ad.t_val_deg,dm.dmap-dmap2); skyplot; colorbar