Получите карту WMS из сервера
[A,R] =
wmsread(layer)[A,R] =
wmsread(layer,Name,Value,...)[A,R] =
wmsread(mapRequestURL)[A,R,mapRequestURL]
= wmsread(...)[ получает доступ к Интернету, чтобы представить и получить растровую карту из сервера Веб-картографического сервиса (WMS). Свойство A,R] =
wmsread(layer)ServerURL объекта WMSLayer, layer, задает сервер. Если layer имеет больше чем один элемент, то сервер накладывает каждый последующий слой сверх основного (первого) слоя, формируя одно изображение. Сервер представляет несколько слоев, только если все слои совместно используют то же значение ServerURL.
Сервер WMS возвращает растровую карту, или цветное или полутоновое изображение, в выводе A. Второй вывод, R, является растровым объектом ссылки, который связывает A с географической системой координат EPSG:4326. Строки A выравниваются с параллелями с ровной выборкой в долготе. Аналогично, столбцы A выравниваются с меридианами с ровной выборкой в широте.
Географические пределы A охватывают полную степень широты и долготы layer. Функция wmsread выбирает больший пространственный размер A, чтобы совпадать с его большей географической размерностью. Больший пространственный размер фиксируется в значении 512. Другими словами, принимая RGB вывод, A является 512-by-N-by-3, если степень широты превышает степень долготы и N-by-512-by-3 в противном случае. В обоих случаях N <= 512. Функция wmsread устанавливает N на целочисленное значение, которое обеспечивает самое близкое приближение, чтобы равняться размерам ячейки в широте и долготе. Карта охватывает полный объем, поддержанный для layer.
[ задает пары значения параметров, которые изменяют запрос к серверу. Можно сократить названия параметра, которые являются нечувствительными к регистру.A,R] =
wmsread(layer,Name,Value,...)
[ использует входной параметр A,R] =
wmsread(mapRequestURL)mapRequestURL, чтобы задать запрос к серверу. mapRequestURL содержит serverURL WMS с дополнительными параметрами WMS. Эти параметры WMS включают BBOX, GetMap и ключевое слово CRS:84 или EPSG:4326. Получите mapRequestURL из вывода wmsread, свойства RequestURL объекта WMSMapRequest или интернет-поиска.
[ возвращается A,R,mapRequestURL]
= wmsread(...)GetMap WMS запрашивают URL в векторе символов mapRequestURL. Можно вставить mapRequestURL в браузер, чтобы выполнить запрос к серверу, который затем возвращает растровую карту. Браузер открывает возвращенную карту, если ее тип пантомимы понят или сохраняет растровую карту на диск.
Считайте слой NASA и отобразите его.
nasa = wmsfind('nasa', 'SearchField', 'serverurl'); layer = nasa.refine('bluemarbleng', 'SearchField', 'layername', ... 'MatchType', 'exact'); [A, R] = wmsread(layer(1)); figure axesm globe axis off geoshow(A, R) title('Blue Marble')
Считайте и отобразите ортоизображение северного раздела моста Золотые Ворота в Сан-Франциско, Калифорния, с помощью Национальной Карты USGS Бесшовный сервер.
Задайте видимую область.
latlim = [37.78 37.841]; lonlim = [-122.53 -122.4];
Найдите слой ортоформирования изображений с высоким разрешением USGS. Национальная Карта USGS обеспечивает ортоформирование изображений из различных областей Соединенных Штатов. Один метод, чтобы получить слой ортоформирования изображений с высоким разрешением должен получить документ возможностей из сервера. Слой ортоформирования изображений является единственным слоем с этого сервера. Используйте несколько попыток соединиться с сервером в случае, если это занято.
numberOfAttempts = 5; attempt = 0; info = []; serverURL = 'http://basemap.nationalmap.gov/ArcGIS/services/USGSImageryOnly/MapServer/WMSServer?'; while(isempty(info)) try info = wmsinfo(serverURL); orthoLayer = info.Layer(1); catch e attempt = attempt + 1; if attempt > numberOfAttempts throw(e); else fprintf('Attempting to connect to server:\n"%s"\n', serverURL) end end end
Получите изображение и отобразите его в проекции UTM.
imageLength = 1024; [A, R] = wmsread(orthoLayer, 'Latlim', latlim, 'Lonlim', lonlim, ... 'ImageHeight', imageLength, 'ImageWidth', imageLength); figure axesm('utm', 'Zone', utmzone(latlim, lonlim), ... 'MapLatlimit', latlim, 'MapLonlimit', lonlim, ... 'Geoid', wgs84Ellipsoid) geoshow(A,R) axis off title({'San Francisco','Northern Section of Golden Gate Bridge'})
Считайте и отобразите глобальный ежемесячный составной объект температуры поверхности моря на 16 апреля 2010 на основе данных из датчика AMSR-E на борту спутника Воды. Включайте береговую линию, landmask, и слои страны.
coastwatch = wmsfind('coastwatch', 'SearchField', 'serverurl'); layers = coastwatch.refine('erdAAsstamday', ... 'Searchfield','serverurl'); time = '2010-04-16T00:00:00Z'; [A, R] = wmsread(layers(end:-1:1), 'Time', time); figure axesm('pcarree', 'Maplonlimit', [0, 360], ... 'PLabelLocation', 45, 'MLabelLocation', 90, ... 'MLabelParallel', -90, 'MeridianLabel', 'on', ... 'ParallelLabel', 'on'); geoshow(A, R); title({layers(end).LayerTitle, time})
Установите Интернет-соединение, чтобы использовать wmsread. Периодически, сервер WMS недоступен. Получение карты может занять несколько минут. wmsread связывается с сервером с помощью объекта WebMapServer, представляющего сервер WMS. Объект действует как прокси к серверу WMS и находится физически на клиентской стороне. Объект получает карту из сервера. Объект автоматически времена после 60 секунд, если связь не установлена с сервером.
Чтобы задать прокси-сервер, чтобы соединиться с Интернетом, выберите File> Preferences> Web и введите свою информацию прокси. Используйте эту функцию, если у вас есть брандмауэр.
Поддержки wmsread, считывающие данные в версиях 1.0.0, 1.1.1 WMS, и 1.3.0. Для версии 1.3.0 только, состояний спецификации WMS, "EPSG:4326 отсылает к WGS 84 географическую широту, затем долгота. Таким образом, в этом CRS x - ось соответствует широте и y - ось к долготе". Большинство серверов обеспечивает данные этим способом; однако, некоторые серверы соответствуют версии 1.1.1, где x - ось соответствует долготе и y - ось к широте.
wmsread пытается подтвердить, подтверждает ли сервер к спецификации. Это проверяет ограничительную рамку EPSG:4326, и если значения XLim превышают область значений широты, то оси подкачиваются, чтобы соответствовать версии 1.1.1 а не 1.3.0. Если wmsread не обнаруживает, что значения XLim превышают область значений широты, и вы замечаете, что пределы широты и долготы инвертируются, то необходимо подкачать их. Можно или изменить параметры bbox в mapRequestURL или изменить Latlim и значения параметров Lonlim, если допустимый.
WMSLayer | WebMapServer | wmsfind | wmsinfo | wmsupdate