В этом примере показано отображение радиолокационных изображений NEXRAD. Изображения охватывают последние 24 часов, дискретизированные с интервалом в один час, для США с использованием данных с сервера IEM WMS. В качестве фона используйте слой JPL Daily Planet.
Найти 'nexrad-n0r-wmst' слой и обновите его.
wmst = wmsfind('nexrad-n0r-wmst', 'SearchField', 'layername'); wmst = wmsupdate(wmst);
Найдите сгенерированный композит CONUS из IR-изображений GOES и обновите его.
goes = wmsfind('goes*conus*ir', 'SearchField', 'layername'); goes = wmsupdate(goes);
Создайте фигуру с требуемой географической протяженностью.
hfig = figure;
region = 'conus';
usamap(region)
mstruct = gcm;
latlim = mstruct.maplatlimit;
lonlim = mstruct.maplonlimit;
Считывайте слой GOES как фоновое изображение.
cellsize = .1; [backdrop, R] = wmsread(goes, 'ImageFormat', 'image/png', ... 'Latlim', latlim, 'Lonlim', lonlim, 'Cellsize', cellsize);
Вычислите текущее время минус 24 часов и настройте кадры для хранения данных из getframe.
now_m24 = datestr(now-1); hour_m24 = [now_m24(1:end-5) '00:00']; hour = datenum(hour_m24); hmap = []; numFrames = 24; frames = struct('cdata', [], 'colormap', []); frames(numFrames) = frames;
За каждый час получайте почасовые данные карты NEXRAD и объединяйте их с копией фона. Из-за того, как этот веб-сервер обрабатывает PNG формат, результирующие данные карты имеют изображение с классом double. Таким образом, вы должны преобразовать его в uint8 перед объединением.
borders = geoshow('usastatehi.shp', 'FaceColor', 'none'); black = [0,0,0]; threshold = 0; for k=1:numFrames time = datestr(hour); [A, R] = wmsread(wmst, 'Latlim', latlim, 'Lonlim', lonlim, ... 'Time', time, 'CellSize', cellsize, ... 'BackgroundColor', black, 'ImageFormat', 'image/png'); delete(hmap) index = any(A > threshold, 3); combination = backdrop; index = cat(3,index,index,index); combination(index) = uint8(255*A(index)); hmap = geoshow(combination, R); uistack(borders,'top') title({wmst.LayerName, time}) drawnow frames(k) = getframe(hfig); hour = hour + 1/24; end
Просмотр цикла ролика.
numTimes = 10; fps = 1.5; movie(hfig, frames, numTimes, fps);