I'm having some problems with contourm and geoshow.
This code produces the intended output, which is a contour map overlaid on top of a map of North America.
if true
lats = 35:54;
lons = -89:-70;
mval = peaks(20);
H1 = worldmap([25 55],[-90,-60]);
setm(gca,'mapprojection','mercator')
coast = load('coast');
states = shaperead('usastatehi', 'UseGeoCoords', true);
geoshow(coast.lat,coast.long);
contourm(lats,lons,mval,'linewidth',2);
end
However, if change to the uncommented geoshow line, which fills in the coastline data to make the land green, the land mass appears to be plotted over the contour plots:
if true
lats = 35:54;
lons = -89:-70;
mval = peaks(20);
H1 = worldmap([25 55],[-90,-60]);
setm(gca,'mapprojection','mercator')
coast = load('coast');
states = shaperead('usastatehi', 'UseGeoCoords', true);
geoshow(coast.lat,coast.long,'DisplayType', 'polygon', 'FaceColor', [.45 .60 .30],'edgecolor','black')
contourm(lats,lons,mval,'linewidth',2);
end
I am running version 2012b on a windows 7 machine. When I try this in the linux system (R2012a), the second set of code produces the expected output. Any thoughts?
