Анимируйте слои данных

Создание фильма Terra/MODIS Maps

Можно создать карты одной и той же географической области в разное время и просмотреть их как фильм. В течение семи дней, прочитайте и отобразите ежедневный набор визуальных изображений из сцен NASA's Moderate Resolution Imaging Spectrororadiometer (MODIS), захваченных в декабре 2010 года.

  1. Поиск по базе данных WMS для слоя MODIS.

    neo = wmsfind('neowms*nasa', 'SearchField', 'serverurl'); 
modis = neo.refine('true*color*terra*modis'); 
modis = wmsupdate(modis);

  2. Создайте WebMapServer объект.

    server = WebMapServer(modis.ServerURL);

  3. Создайте WMSMapRequest объект.

    mapRequest = WMSMapRequest(modis, server);

  4. The Extent поле содержит информацию о способе извлечения отдельных систем координат. Вы можете запросить один день, так как объем определяется днем ('/P1D'). Обратите внимание, что для декабря 2010 года системы координат для 8 и 31 декабря недоступны.

    modis.Details.Dimension.Extent

  5. Создайте массив, указывающий первые семь дней.

    days = 1:7;

  6. Установите значение startTime по 01 декабря 2010 года и использовать серийный номер даты.

    time = '2010-12-01';
startTime = datenum(time);

  7. Откройте окно рисунка с осями, соответствующими области, заданной modis слой.

    hFig = figure('Color', 'white');
worldmap(mapRequest.Latlim, mapRequest.Lonlim);

  8. Сохраните каждую систему координат в видео файла.

    videoFilename = 'modis_dec.avi';
writer = VideoWriter(videoFilename);
writer.FrameRate = 1;
writer.Quality = 100;
writer.open;

  9. Извлечение карты 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

  10. Читайте во всех видеокадрах.

    v = VideoReader(videoFilename);
vidFrames = read(v);
numFrames = get(v, 'NumberOfFrames');

  11. Создайте MATLAB® структура фильма из видеокадров.

    frames = struct('cdata', [], 'colormap', []);
frames(numFrames) = frames(1);
for k = 1 : numFrames
    frames(k).cdata = vidFrames(:,:,:,k);
    frames(k).colormap = [];
end

  12. Воспроизведение фильма один раз на системе координат видео.

    movie(hFig, frames, 1, v.FrameRate)

Создание анимированного GIF Файла из карт WMS

Прочитайте и отобразите анимацию о ледяном шельфе Ларсена, испытывающем драматический крах в период с 31 января по 7 марта 2002 года.

  1. Найдите в базе данных WMS фразу «Ледяной шельф Ларсена».

    iceLayer = wmsfind('Larsen Ice Shelf');

    Попробуйте первый слой.

  2. Создайте WebMapServer объект.

    server = WebMapServer(iceLayer(1).ServerURL);

  3. Используйте WebMapServer.updateLayers метод для синхронизации слоя с исходным сервером WMS. Получите самые последние данные и заполните Abstract, CoordRefSysCodes, и Details поля.

    iceLayer = server.updateLayers(iceLayer(1));

  4. Просмотрите аннотацию.

    fprintf('%s\n', iceLayer(1).Abstract)

  5. Создайте WMSMapRequest объект.

    request = WMSMapRequest(iceLayer(1), server);

  6. Поскольку вы обновили свой слой, Details теперь поле имеет содержимое. Нажмите Details в редакторе Переменного MATLAB. Затем нажмите Dimension. Имя размерности 'time'. Нажмите Extent. The Extent поле предоставляет доступные значения для размерности, в этом случае времени. Сохраните эту информацию, введя в командной строке следующее:

    extent = [',' iceLayer.Details.Dimension.Extent, ','];

  7. Вычислим количество необходимых систем координат. (The extent содержит запятую перед первой системой координат и после последней системы координат. Чтобы получить количество систем координат, вычитайте 1.)

    frameIndex = strfind(extent, ',');
numFrames = numel(frameIndex) - 1;

  8. Откройте окно рисунка и настройте оси карты с соответствующими географическими пределами.

    h = figure;
worldmap(request.Latlim, request.Lonlim)

  9. Установите свойства осей карты. MLineLocation устанавливает интервал между отображаемыми меридианами сетки. MLabelParallel определяет параллель, в которой отображаются метки.

    setm(gca,'MLineLocation', 1, 'MLabelLocation', 1, ...
   'MLabelParallel',-67.5, 'LabelRotation', 'off');

  10. Инициализируйте значение animated в 0.

    animated(1,1,1,numFrames) = 0;

  11. Отображение изображения ледяного шельфа Ларсен в разные дни.

    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

  12. Сохраните и затем просмотрите анимированный 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 минут. Считывайте и отображайте объединенные слои.

  1. Найти слои в базе данных WMS, которые включают 'mesonet' и 'nexrad' в своих ServerURL поля.

    mesonet = wmsfind('mesonet*nexrad', 'SearchField', 'serverurl');

  2. Основа NEXRAD отражает ток ('nexrad-n0r') измеряет интенсивность осадков. Уточните поиск, включив только слои с этой фразой в одно из полей поиска.

    nexrad = mesonet.refine('nexrad-n0r', 'SearchField', 'any');

  3. Удалите слои 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);

  4. Обновите свои nexrad слой для заполнения всех полей и получения самых последних данных.

    nexrad = wmsupdate(nexrad, 'AllowMultipleServers', true);

  5. '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);

  6. Отображение карты объединенных слоев NEXRAD. Наложение на пограничные многоугольники Соединенных Штатов.

    geoshow(A, R);
geoshow('usastatehi.shp', 'FaceColor', 'none');
title({'NEXRAD Radar Map', 'Merged Layers'});

  7. Цикл через последовательность временных радиолокационных наблюдений.

    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

  8. Создайте массив для записи в виде анимированного 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

  9. Просмотрите анимированный GIF файла.

    filename = 'wmsnexrad.gif';
imwrite(animated, cmap, filename, 'DelayTime', 1.5, ...
   'LoopCount', inf);
web(filename)

См. также

| |

Похожие темы

Mapping Toolbox

Поддержка

Для просмотра документации необходимо авторизоваться на сайте