В этом примере показано, как отобразить радиолокационные изображения NEXRAD. Изображения охватывают последние 24 часа, взятые за один час интервалов, для Соединенных Штатов с использованием данных сервера IEM WMS. Используйте слой JPL Daily Planet в качестве фона.
Найдите 'nexrad-n0r-wmst' слой и обновить его.
wmst = wmsfind('nexrad-n0r-wmst', 'SearchField', 'layername'); wmst = wmsupdate(wmst);
Найдите сгенерированный CONUS композит изображений GOES IR и обновите его.
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);