I just installed MatLab r2012b onto a Dell Precision T1650 (Windows 7 Pro)and find I am unable to get the following code to print the .tif. The exact same code works on a r2011b version of MatLab on a Dell Vostro 460 (also Windows 7 Pro). I'm using the Mapping Toolbox 3.7 and the Vostro has v. 3.4.
Any ideas on what I'm doing wrong?
filename='P:\RDData\2013\01\09\ProcessRD2\20130109RLIV_r08RLIV5_ReportsByLatLon_TBL.csv' ;clipboarddata=dlmread(filename,',',2,1) ;m=1:(6.000000*60/1);n=1:(8.000000*60/1);data=clipboarddata(m,n);d=find(data==0);data(d)=NaN;[x,y,s]=find(data);latout=30.000000-(x*1/60);lonout=-98.000000+(y*1/60);xmin=min(latout);ymin=-98.000000;xmax=30.000000;ymax=max(lonout);A=imread('m:/ChartResevoir/LiverpoolGulf.tif','pixelregion',{[1 1125],[1 1500]});B=flipdim(A,1);R=makerefmat(-98.000000,30.000000-(6.000000),0.005333,(0.005333));refvec=refmat2vec(R,size(B));xres=1+60*6.000000/1;yres=1+60*8.000000/1;xv=linspace(xmin,xmax,xres);yv=linspace(ymin,ymax,yres);[Xinterp,Yinterp]=meshgrid(xv,yv);Zinterp=griddata(latout,lonout,s,Xinterp,Yinterp,'nearest');h=figure;axesm('eqdcylin','mlinelocation',[1],'plinelocation',[1],'mlabelparallel',[30.000000-6.000000],'plabelmeridian',[-98.000000]);geoshow(B,refvec);setm(gca,'Grid','on','glinestyle',':','maplatlimit',[(30.000000-6.000000) (30.000000+0)],'maplonlimit',[(-98.000000-0) (-98.000000+8.000000)],'meridianlabel','on','parallellabel','on','fontweight','bold');setm(gca,'plinelimit',[(-98.000000-0) (-98.000000+8.000000)],'mlinelimit',[(30.000000-6.000000) (30.000000+0)]);setm(gca,'mlabellocation',[2],'plabellocation',[2]);setm(gca,'mlabelround',[0],'plabelround',[0]);set(gcf,'inverthardcopy','off');hold on;
Best Answer