Анимационные слои данных

Создайте Фильм Карт Terra/MODIS

Можно создать карты той же географической области в разное время и просмотреть их как фильм. Сроком на семь дней считайте и отобразите ежедневный составной объект визуальных изображений от Умеренного Спектрорадиометра Обработки изображений Разрешения НАСА (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. 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. Extent поле вводит доступные значения для размерности в это время случая. Сохраните эту информацию путем ввода следующего в командной строке:

    extent = [',' iceLayer.Details.Dimension.Extent, ','];
    
  7. Вычислите количество необходимых систем координат. (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
    end
  12. Сохраните и затем просмотрите анимированный файл GIF.

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

    Создайте снимки от анимации Шельфового ледника Ларсена

Анимационные радарные наблюдения таймлапса

Отобразите Радар Следующего поколения (NEXRAD) изображения для данных об использовании Соединенных Штатов из веб-сервера карты Айовы экологической мезосети (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. Удалите 900 913 слоев, потому что они предназначаются для наложения Google Maps™. Также удалите слой 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)

Смотрите также

| |

Похожие темы

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