Issue charting onto a tif.

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;

[Walter Roberson]

What happens instead of the tif being printed?

When you say "printed" are you referring to the imread()? If so then is the data accessible in that exact location on the m: drive on the machine that it is not working on? Do you have permissions to read it?

dir('m:/ChartResevoir/*.tif')

[Renee]

It prints out a plain white grid. When I run it on the other machine it prints the .tif with the labels on the side. This is just the grid on a white background.

I have read and write permissions on that drive (M:\).

What it's doing: http://i579.photobucket.com/albums/ss236/Tasare13/Work/WhatItsDoing_zps301fe1b2.png

What it should be doing: http://i579.photobucket.com/albums/ss236/Tasare13/Work/WhatItShouldDo_TomMachine_zps968d0006.png

[Walter Roberson]

You mean "display" rather than "print".

Use the debugger and put a breakpoint in at the statement B=flipdim(A,1); and run. When it stops, please show

size(A)
class(A)

and the same thing on the system it is working on.

You could step one line on the two systems until you find the place that the behavior differs.

[Renee]

Thank you! I'm still quite new to MatLab and don't have the verbage down yet.

I added the breakpoint and Siza(a) and Class(A) on the broken one are: K>> size(A)

      ans =

              1125        1500           3

      K>> class (A)

      ans =

      uint8

The Size(A) and Class(A) on the functioning one are exactly the same.

We stepped through the rest of the code and found that:

geoshow(B,refvec);

displays the chart that I want. However the following:

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');

removes the chart and just gives me the white graphic with the grid. While on the other machine it keeps the chart and overlays the gridlines and numbers.

[Walter Roberson]

On the system that is giving the problem, experiment with using the command

opengl software

[Renee]

Thank you. That seems to have fixed it!

Answer 1

[Walter Roberson]

As above, fixed by switching to "opengl software"