I need to plot weather data (LAT/LON specified) above Denmark using a satellite picture. I chose Blue Marble as it has worked perfectly for something like a stereographic with Greenland in focus, and it can be properly set up for longitude and latitude data. However a mercator projection for Denmark needs a higher resolution image, thus the image resolution will change (WIDTH= and HEIGHT= in DenmarkBlueMarbleURL). This however increases the runtime rather exponentially.
The following code is the problematic part from a larger project,
disp(' - Denmark')
DenmarkBlueMarbleURL = 'http://neowms.sci.gsfc.nasa.gov/wms/wms?SERVICE=WMS&LAYERS=BlueMarbleNG&EXCEPTIONS=application/vnd.ogc.se_xml&FORMAT=image/jpeg&TRANSPARENT=FALSE&HEIGHT=6000&BGCOLOR=0xFFFFFF&REQUEST=GetMap&WIDTH=7200&BBOX=-180.0,-90.0,180.0,90.0&STYLES=&SRS=EPSG:4326&VERSION=1.1.1';
[A, R] = wmsread(DenmarkBlueMarbleURL);
FigDenmark = figure('Position', [100, 100, 1440, 900]);
axesm('mercator', 'MapLatLimit', [52 70], 'MapLonLimit', [00 30])
axis off;
geoshow(A, R);
pcolorm(YM, XM, AuroraData, 'FaceAlpha', .5); shading flat;
colormap(AuroraColourmap);
caxis([0 100])
ColourBarAurora = colorbar; ylabel(ColourBarAurora, 'Probability in percentage of visible aurora');
title(['Aurora forecast for ' datestr(CurrentTime, 'yyyy/mm/dd HHMM')]);
gridm;
drawnow
saveas(FigDenmark, 'CurrentDataPlot', 'tif');
It takes me about five minutes to run this section of the code. I believe the problem is geoshow(A,R), since everything goes fine until then, so this is the interesting part, ignore the AuroraData-dependent code. Also, if you change the latitude and longitude (i.e. to only fetch a map of Denmark) in the DenmarkBlueMarbleURL (BBOX=) it seems to take longer as well, so I always just fetch the whole map.
In short: Is there a faster way to use a blue marble plot to plot data on top?
