Документация

Ищите Базу данных WMS все слои с 'nasa' в ServerURL поле .

nasa = wmsfind('nasa','SearchField','serverurl');

Используйте WMSLayer.refine метод, чтобы совершенствовать ваш поиск, чтобы включать только те слои с фразой 'bluemarbleng' в LayerName поле . Этот синтаксис создает точный поиск.

layer = nasa.refine('bluemarbleng', 'SearchField','layername', ...
   'MatchType','exact');

Используйте wmsread функция, чтобы получить первый слой Blue Marble Next Generation.

[A,R] = wmsread(layer(1));

wmsread функция возвращает A, растровая карта, на которую географически ссылаются и R, объект растровой привязки, который связывает A к географической системе координат EPSG:4326. Географические пределы A охватите полную степень широты и долготы layer.

Создайте окно фигуры, настройте свои оси карты и отобразите вашу карту.

figure
axesm globe
axis off
geoshow(A,R)
title('Blue Marble: Next Generation')

Слой, используемый в этом примере, является любезностью NASA/JPL-Caltech.

Ищите базу данных WMS слои, которые содержат foundation.gtopo30 в LayerName поле . GTOPO30, цифровая модель вертикального изменения, разработанная Геологической службой США (USGS), имеет интервал сетки по горизонтали 30 секунд дуги.

gtopo30Layer = wmsfind('foundation.gtopo30');

Задайте цвет фона, задав красные, зеленые, и синие уровни.

oceanColor = [0 170 255];

Используйте BackgroundColor и CellSize параметры wmsread функционируйте, чтобы выбрать цвет фона и размер ячейки вашей полученной карты.

cellSize = 0.5;
[A,R] = wmsread(gtopo30Layer,'BackgroundColor',oceanColor, ...
    'CellSize', cellSize);

Создайте окно фигуры и настройте мировую карту оси. Отобразите свою карту с заголовком.

figure
worldmap world
geoshow(A,R)
title({'GTOPO30 Elevation Model',gtopo30Layer.LayerTitle})

Ищите слои от Космического полета Годдарда НАСА Сервер SVS Изображений. Этот сервер содержит слои, которые имеют изображения легенды. Можно сказать, что изображения легенды доступны, потому что слои имеют содержимое в LegendURL поле .

layers = wmsfind('svs.gsfc.nasa.gov','SearchField','serverurl');
serverURL = layers(1).ServerURL;
gsfc = wmsinfo(serverURL);

Найдите слой, содержащий городские температурные подписи, и отобразите краткий обзор:

urban_temperature = gsfc.Layer.refine('urban*temperature');
disp(urban_temperature.Abstract)

Big cities influence the environment around them. For example,
urban areas are typically warmer than their surroundings.
Cities are strikingly visible in computer models that simulate
the Earth's land surface. This visualization shows average 
surface temperature predicted by the Land Information System (LIS)
for a day in June 2001. Only part of the global computation
is shown, focusing on the highly urbanized northeast corridor
in the United States, including the cities of Boston, New York,
Philadelphia, Baltimore, and Washington.

Additional Credit:
NASA GSFC Land Information System (http://lis.gsfc.nasa.gov/)

Считайте и отобразите слой. Карта появляется с различными цветами в различных областях, но без легенды не ясно, что представляют эти цвета.

[A,R] = wmsread(urban_temperature);
figure
usamap(A,R)
geoshow(A,R)
title('Urban Temperature Signatures')
axis off

Исследуйте Details поле urban_temperature слой. Этот слой имеет только одну структуру в Style поле . Style поле определяет, как сервер представляет слой.

urban_temperature.Details

ans = 

  struct with fields:

     MetadataURL: 'http://svs.gsfc.nasa.gov/vis/a000000/a003100/a003152/a003152.fgdc'
      Attributes: [1x1 struct]
     BoundingBox: [1x1 struct]
       Dimension: [1x1 struct]
    ImageFormats: {'image/png'}
     ScaleLimits: [1x1 struct]
           Style: [1x1 struct]
         Version: '1.3.0'

Отобразите Style поле в Командном окне:

urban_temperature.Details.Style

ans = 

        Title: 'Opaque'
         Name: 'opaque'
     Abstract: [1x319 char]
    LegendURL: [1x1 struct]

Каждый Style элемент имеет только один LegendURL. Исследуйте LegendURL:

urban_temperature.Details.Style.LegendURL

ans = 

    OnlineResource: [1x65 char]
            Format: 'image/png'
            Height: 90
             Width: 320

Загрузите легенду URL:

url = urban_temperature.Details.Style.LegendURL.OnlineResource

URL появляется в командном окне:

url =

http://svs.gsfc.nasa.gov/vis/a000000/a003100/a003152/temp_bar.png

Отобразите изображение легенды с помощью image команда и свойства набора, такие, что изображение отображается с непосредственным, разрешением экрана к пикселю. Эта легенда является просто изображением шкалы палитры, не легендой в графике MATLAB^®.

temperatureLegend = webread(url);
figure('Color','white')
axis off image
set(gca,'units','pixels','position',...
[0 0 size(temperatureLegend,2) size(temperatureLegend,1)]);
pos = get(gcf,'position');
set(gcf,'position',...
[pos(1) pos(2) size(temperatureLegend,2) size(temperatureLegend,1)]);
image(temperatureLegend)

Теперь карта имеет больше смысла. Области к красному концу спектра теплее.

Шаги 7-10 демонстрируют, как получить выход из карты, структурируют и добавляют легенду. Путем добавления легенды этим способом вы стараетесь не деформировать текст в изображении легенды. (Текст легенды деформируется, если вы отображаете изображение с geoshow.)

Сначала установите свои пределы широты и долготы совпадать с пределами вашей карты и чтения в файле форм с мировыми городскими данными:

latlim = R.LatitudeLimits;
lonlim = R.LongitudeLimits;
S = shaperead('worldcities','UseGeoCoords',true,...
    'BoundingBox',[lonlim(1) latlim(1);lonlim(2) latlim(2)]);

Определите положение окна текущей фигуры. Варьируйтесь pos(1) и pos(2) 'Position' параметры по мере необходимости на основе разрешения вашего экрана.

colValue = [1 1 1];
dimension = size(A,1)/2;
figure
set(gcf,'Color',[1,1,1])
pos = get(gcf, 'Position');
set(gcf, 'Position', [pos(1) pos(2) dimension dimension])

Отобразите карту и добавьте городские маркеры, государственные границы, меридиан и найдите что-либо подобное меткам, заголовку и Северной стреле:

usamap(A,R)
geoshow(A,R)
geoshow(S, 'MarkerEdgeColor', colValue, 'Color', colValue)

geoshow('usastatehi.shp', 'FaceColor', 'none',...
   'EdgeColor','black')
mlabel('FontWeight','bold')
plabel('FontWeight','bold')
axis off
title('Urban Temperature Signatures', 'FontWeight', 'bold')

for k=1:numel(S)
   textm(S(k).Lat, S(k).Lon, S(k).Name, 'Color', colValue,...
       'FontWeight','bold')
end

lat = 36.249;
lon = -71.173;
northarrow('Facecolor', colValue, 'EdgeColor', colValue,...
    'Latitude', lat, 'Longitude', lon);

Отобразите карту и легенду как одно, объединенное изображение:

f = getframe(gcf);
legendImg = uint8(255*ones(size(temperatureLegend,1),size(f.cdata,2),3));
offset = dimension/2;
halfSize = size(temperatureLegend, 2)/2;
legendImg(:,offset-halfSize:offset+halfSize-1,:) = temperatureLegend;
combined = [f.cdata; legendImg];
figure
pos = get(gcf,'position');
set(gcf,'position',[10 100 size(combined,2) size(combined,1)])
set(gca,'units','normalized','position', ...
   [0 0 1 1]);
image(combined)
axis off image

Другой способ отобразить карту и легенду вместе состоит в том, чтобы записать легенду в карту в заданном местоположении. Чтобы просмотреть изображение, используйте image команда, устанавливая параметры положения, таким образом, что существует непосредственный пиксель к разрешению экрана. (Текст легенды деформируется, если изображение отображено с geoshow.)

A_legend = A;
A_legend(end-size(temperatureLegend,1):end-1,...
   end-size(temperatureLegend,2):end-1,:) = temperatureLegend;
figure
image(A_legend)
axis off image
set(gca,'Units','normalized','position',...
   [0 0 1 1]);
set(gcf,'Position',[10 100 size(A_legend,2) size(A_legend,1)]);
title('Urban Temperature Signatures', 'FontWeight', 'bold')

Объедините карту и легенду в одном файле, и затем опубликуйте его к сети. Сначала запишите изображения в файл:

mkdir('html')
imwrite(A_legend, 'html/wms_legend.png')
imwrite(combined, 'html/combined.png')

Откройте редактор MATLAB и вставку в этом коде:

%%
% <<wms_legend.png>>

%%
% <<combined.png>>

Добавьте любой другой текст, который вы хотите включать в свой опубликованный документ. Затем выберите одну из ячеек и выберите File> Save File and Publish в меню.

Получите географические координаты, которые являются случайными с изображением в файле boston.tif.

filename = 'boston.tif';
info = georasterinfo(filename);
R = info.RasterReference;
proj = R.ProjectedCRS;
[latlim,lonlim] = projinv(proj, ...
                          R.XWorldLimits,R.YWorldLimits);

Преобразуйте географические пределы Меркаторской сети. Для этого сначала создайте projcrs объект путем определения системы координат EPSG:3857. Затем спроектируйте использование координат projfwd. Чтобы получить формирование изображений в этой системе координат, необходимо использовать WMSMapRequest и WebMapServer с тех пор wmsread только запрашивает данные в системе EPSG:4326.

p = projcrs(3857);
[x,y] = projfwd(p,latlim,lonlim);
xlimits = [min(x) max(x)];
ylimits = [min(y) max(y)];

Вычислите значения высоты изображения и ширины для объема выборки 5 метров.

metersPerSample = 5;
imageHeight = round(diff(ylimits)/metersPerSample);
imageWidth  = round(diff(xlimits)/metersPerSample);

Повторно вычислите новые пределы.

yLim = [ylimits(1), ylimits(1) + imageHeight*metersPerSample];
xLim = [xlimits(1), xlimits(1) + imageWidth*metersPerSample];

Найдите Национальную Карту USGS от базы данных WMS и выберите слой Digital Ortho-Quadrangle.

doqLayer = wmsfind('usgsnaipplus','SearchField','serverurl');
doqLayer = wmsupdate(refine(doqLayer,'image'));
doqLayer = refine(doqLayer,'USGSNaip','SearchField','abstract');

Создайте WebMapServer и WMSMapRequest объекты.

server  = WebMapServer(doqLayer.ServerURL);
numberOfAttempts = 50;
attempt = 0;
request = [];
while(isempty(request))
    try
        request = WMSMapRequest(doqLayer, server);
    catch e 
        attempt = attempt + 1;
        fprintf('%d\n', attempt)
        if attempt > numberOfAttempts
            throw(e)
        end
    end
end

Используйте WMSMapRequest свойства изменить различные аспекты вашего запроса карты, такие как карта ограничивает, размер изображения и код системы координат. Установите пределы карты, чтобы покрыть ту же область, как найдено в boston.tif файл.

request.CoordRefSysCode = 'EPSG:3857';
request.ImageHeight = imageHeight;
request.ImageWidth  = imageWidth;
request.XLim = xLim;
request.YLim = yLim;

Запросите карту ортоформирования изображений в сети Меркаторские координаты.

A_PCS = getMap(server, request.RequestURL);
R_PCS = request.RasterReference;

Получите карту для той же области, но в координатах EPSG:4326.

request.CoordRefSysCode = 'EPSG:4326';
request.Latlim = latlim;
request.Lonlim = lonlim;
A_Geo = getMap(server, request.RequestURL);
R_Geo = request.RasterReference;

Читайте в Бостонских названиях места из файла форм и наложите их сверху карт. Преобразуйте координаты функций к сети Меркаторские и географические координаты. Координаты точки в файле форм исчисляются в метрах и Координаты на плоскости состояния Массачусетса, но проекция GeoTIFF задана в ногах обзора.

S = shaperead('boston_placenames');
x = [S.X]*unitsratio('sf','meter');
y = [S.Y]*unitsratio('sf','meter');
names = {S.NAME};
[lat, lon] = projinv(proj, x, y);
[xPCS, yPCS] = projfwd(p, lat, lon);

Проект и отображение ортоформирование изображений, полученное в использовании координат EPSG:4326, геопоказывают.

figure
axesm('mercator','MapLatLimit',latlim,'MapLonLimit',lonlim)
geoshow(A_Geo,R_Geo)
textm(lat, lon, names, 'Color',[0 0 0], ...
    'BackgroundColor',[0.9 0.9 0],'FontSize',6);
axis tight
title({'USGS Digital Ortho-Quadrangle - Boston', ...
    'Geographic Layer'})

Отобразите ортоформирование изображений, полученное в сети Меркаторские координаты.

figure
mapshow(A_PCS, R_PCS);
text(xPCS, yPCS, names, 'Color',[0 0 0], ...
    'BackgroundColor',[0.9 0.9 0],'FontSize',6,'Clipping','on');
axis tight
title({'USGS Digital Ortho-Quadrangle - Boston', 'Web Mercator Layer'})

Отобразите изображение от boston.tif для сравнения.

figure
mapshow(filename)
text(x, y, names, 'Color',[0 0 0], ...
    'BackgroundColor',[0.9 0.9 0],'FontSize',6,'Clipping','on');
axis tight
title({filename, proj.Name})