Можно создать карты одной и той же географической области в разное время и просмотреть их как фильм. В течение семи дней, прочитайте и отобразите ежедневный набор визуальных изображений из сцен NASA's Moderate Resolution Imaging Spectrororadiometer (MODIS), захваченных в декабре 2010 года.
Поиск по базе данных WMS для слоя MODIS.
neo = wmsfind('neowms*nasa', 'SearchField', 'serverurl'); modis = neo.refine('true*color*terra*modis'); modis = wmsupdate(modis);
Создайте WebMapServer объект.
server = WebMapServer(modis.ServerURL);
Создайте WMSMapRequest объект.
mapRequest = WMSMapRequest(modis, server);
The Extent поле содержит информацию о способе извлечения отдельных систем координат. Вы можете запросить один день, так как объем определяется днем ('/P1D'). Обратите внимание, что для декабря 2010 года системы координат для 8 и 31 декабря недоступны.
modis.Details.Dimension.Extent
Создайте массив, указывающий первые семь дней.
days = 1:7;
Установите значение startTime по 01 декабря 2010 года и использовать серийный номер даты.
time = '2010-12-01';
startTime = datenum(time);
Откройте окно рисунка с осями, соответствующими области, заданной modis слой.
hFig = figure('Color', 'white'); worldmap(mapRequest.Latlim, mapRequest.Lonlim);
Сохраните каждую систему координат в видео файла.
videoFilename = 'modis_dec.avi';
writer = VideoWriter(videoFilename);
writer.FrameRate = 1;
writer.Quality = 100;
writer.open;
Извлечение карты modis слой для каждого требуемого дня. Установите Time свойство по номеру дня. При получении данных с сервера используйте try/catch оператор, игнорирующий данные, не найденные на сервере, или ошибку, выданную сервером. Задайте startTime на один день меньше для правильной индексации.
startTime = startTime - 1; for k = days try mapRequest.Time = startTime + k; timeStr = datestr(mapRequest.Time); dailyImage = server.getMap(mapRequest.RequestURL); geoshow(dailyImage, mapRequest.RasterReference); title({mapRequest.Layer.LayerTitle, timeStr}, ... 'Interpreter', 'none', 'FontWeight', 'bold') shg frame = getframe(hFig); writer.writeVideo(frame); catch e fprintf(['Server error: %s.\n', ... 'Ignoring frame number %d on day %s.\n'], ... e.message, k, timeStr) end drawnow shg end writer.close
Читайте во всех видеокадрах.
v = VideoReader(videoFilename);
vidFrames = read(v);
numFrames = get(v, 'NumberOfFrames');
Создайте MATLAB® структура фильма из видеокадров.
frames = struct('cdata', [], 'colormap', []); frames(numFrames) = frames(1); for k = 1 : numFrames frames(k).cdata = vidFrames(:,:,:,k); frames(k).colormap = []; end
Воспроизведение фильма один раз на системе координат видео.
movie(hFig, frames, 1, v.FrameRate)
Прочитайте и отобразите анимацию о ледяном шельфе Ларсена, испытывающем драматический крах в период с 31 января по 7 марта 2002 года.
Найдите в базе данных WMS фразу «Ледяной шельф Ларсена».
iceLayer = wmsfind('Larsen Ice Shelf');Попробуйте первый слой.
Создайте WebMapServer объект.
server = WebMapServer(iceLayer(1).ServerURL);
Используйте WebMapServer.updateLayers метод для синхронизации слоя с исходным сервером WMS. Получите самые последние данные и заполните Abstract, CoordRefSysCodes, и Details поля.
iceLayer = server.updateLayers(iceLayer(1));
Просмотрите аннотацию.
fprintf('%s\n', iceLayer(1).Abstract)Создайте WMSMapRequest объект.
request = WMSMapRequest(iceLayer(1), server);
Поскольку вы обновили свой слой, Details теперь поле имеет содержимое. Нажмите Details в редакторе Переменного MATLAB. Затем нажмите Dimension. Имя размерности 'time'. Нажмите Extent. The Extent поле предоставляет доступные значения для размерности, в этом случае времени. Сохраните эту информацию, введя в командной строке следующее:
extent = [',' iceLayer.Details.Dimension.Extent, ','];
Вычислим количество необходимых систем координат. (The extent содержит запятую перед первой системой координат и после последней системы координат. Чтобы получить количество систем координат, вычитайте 1.)
frameIndex = strfind(extent, ',');
numFrames = numel(frameIndex) - 1;Откройте окно рисунка и настройте оси карты с соответствующими географическими пределами.
h = figure; worldmap(request.Latlim, request.Lonlim)
Установите свойства осей карты. MLineLocation устанавливает интервал между отображаемыми меридианами сетки. MLabelParallel определяет параллель, в которой отображаются метки.
setm(gca,'MLineLocation', 1, 'MLabelLocation', 1, ... 'MLabelParallel',-67.5, 'LabelRotation', 'off');
Инициализируйте значение animated в 0.
animated(1,1,1,numFrames) = 0;
Отображение изображения ледяного шельфа Ларсен в разные дни.
for k=1:numFrames request.Time = extent(frameIndex(k)+1:frameIndex(k+1)-1); iceImage = server.getMap(request.RequestURL); geoshow(iceImage, request.RasterReference) title(request.Time, 'Interpreter', 'none') drawnow shg frame = getframe(h); if k == 1 [animated, cmap] = rgb2ind(frame.cdata, 256, 'nodither'); else animated(:,:,1,k) = rgb2ind(frame.cdata, cmap, 'nodither'); end pause(2) end
Сохраните и затем просмотрите анимированный GIF файла.
filename = 'wmsanimated.gif'; imwrite(animated, cmap, filename, 'DelayTime', 1.5, ... 'LoopCount', inf); web(filename)
Снимок из анимации шельфа льда Ларсена
Предоставлено NASA/Goddard Space Flight Center Scientific Visualization Studio
Отображение изображений радара следующего поколения (NEXRAD) для США с помощью данных сервера веб-карты Iowa Environmental Mesonet (IEM). Сервер хранит слои, покрывающие последние 50 минут до настоящего времени, за шаги 5 минут. Считывайте и отображайте объединенные слои.
Найти слои в базе данных WMS, которые включают 'mesonet' и 'nexrad' в своих ServerURL поля.
mesonet = wmsfind('mesonet*nexrad', 'SearchField', 'serverurl');
Основа NEXRAD отражает ток ('nexrad-n0r') измеряет интенсивность осадков. Уточните поиск, включив только слои с этой фразой в одно из полей поиска.
nexrad = mesonet.refine('nexrad-n0r', 'SearchField', 'any');
Удалите слои 900913, потому что они предназначены для наложения Maps™ Google. Также удалите слой WMST, поскольку он содержит данные за разное время.
layers_900913 = nexrad.refine('900913', 'SearchField', ... 'layername'); layer_wmst = nexrad.refine('wmst', 'SearchField', 'layername'); rmLayerNames = {layers_900913.LayerName layer_wmst.LayerName}; index = ismember({nexrad.LayerName}, rmLayerNames); nexrad = nexrad(~index);
Обновите свои nexrad слой для заполнения всех полей и получения самых последних данных.
nexrad = wmsupdate(nexrad, 'AllowMultipleServers', true); 'conus' представляет собой контерминированные 48 состояния США (все, кроме Гавайев и Аляски). Используйте usamap функция для создания осей для контерминированных состояний. Читайте в nexrad слои.
region = 'conus'; figure usamap(region) mstruct = gcm; latlim = mstruct.maplatlimit; lonlim = mstruct.maplonlimit; [A, R] = wmsread(nexrad, 'Latlim', latlim, 'Lonlim', lonlim);
Отображение карты объединенных слоев NEXRAD. Наложение на пограничные многоугольники Соединенных Штатов.
geoshow(A, R); geoshow('usastatehi.shp', 'FaceColor', 'none'); title({'NEXRAD Radar Map', 'Merged Layers'});
Цикл через последовательность временных радиолокационных наблюдений.
hfig = figure; usamap(region) hstates = geoshow('usastatehi.shp', 'FaceColor', 'none'); numFrames = numel(nexrad); frames = struct('cdata', [], 'colormap', []); frames(numFrames) = frames; hmap = []; frameIndex = 0; for k = numFrames:-1:1 frameIndex = frameIndex + 1; delete(hmap) [A, R] = wmsread(nexrad(k), 'Latlim', latlim, 'Lonlim', lonlim); hmap = geoshow(A, R); uistack(hstates,'top') title(nexrad(k).LayerName) drawnow frames(frameIndex) = getframe(hfig); end
Создайте массив для записи в виде анимированного GIF файла.
animated(1,1,1,numFrames) = 0; for k=1:numFrames if k == 1 [animated, cmap] = rgb2ind(frames(k).cdata, 256, 'nodither'); else animated(:,:,1,k) = ... rgb2ind(frames(k).cdata, cmap, 'nodither'); end end
Просмотрите анимированный GIF файла.
filename = 'wmsnexrad.gif'; imwrite(animated, cmap, filename, 'DelayTime', 1.5, ... 'LoopCount', inf); web(filename)