Извлечение карты WMS с сервера
[ обеспечивает доступ к Интернету для визуализации и извлечения растровой карты с сервера картографического Веб-сервиса (WMS). The A,R] =
wmsread(layer)ServerURL свойство WMSLayer объект, layer, задает сервер. Если layer имеет более одного элемента, затем сервер накладывает каждый последующий слой на верхнюю часть основы (первого) слоя, образуя одно изображение. Сервер отображает несколько слоев, только если все слои имеют один и тот же ServerURL значение.
Сервер WMS возвращает растровую карту, цветное или полутоновое изображение, в выходных данных A. Второй выход, R, является объектом растровой привязки, который связывает A в EPSG:4326 географическую систему координат. Строки A выровнены с параллелями, с четной дискретизацией по долготе. Аналогично, столбцы A выровнены с меридианами, с четной дискретизацией в широте.
Географические пределы A охватывать полную широту и степень долготы layer. The wmsread функция выбирает больший пространственный размер A для соответствия его большей географической размерности. Больший пространственный размер фиксируется на значении 512. Другими словами, принимая выход RGB, A является 512-by-N-by-3 если ширина широты превышает длину долготы и N-by-512-by-3 в противном случае. В обоих случаях N <= 512. The wmsread наборы функций N к целому значению, которое обеспечивает максимально близкое приближение к равным размерам камер по широте и долготе. Карта охватывает весь объем, поддерживаемый для layer.
[ задаёт пары значение параметров, которые изменяют запрос к серверу. Можно сокращать имена параметров, которые не чувствительны к регистру.A,R] =
wmsread(layer,Name,Value,...)
[ использует входной параметр A,R] =
wmsread(mapRequestURL)mapRequestURL для определения запроса к серверу. The mapRequestURL содержит serverURL WMS с дополнительными параметрами WMS. Эти параметры WMS включают BBOX, GetMap и EPSG:4326 или CRS:84 ключевое слово. Получите mapRequestURL из выхода wmsread, а RequestURL свойство WMSMapRequest объект или поиск в Интернете.
[ возвращает A,R,mapRequestURL]
= wmsread(...)GetMap WMS запросить URL в вектор символов mapRequestURL. Можно вставить mapRequestURL в браузер, чтобы сделать запрос на сервер, который затем возвращает растровую карту. Браузер открывает возвращенную карту, если ее тип мима понятен, или сохраняет растровую карту на диск.
Поиск в базе данных WMS слоев, содержащих строку «NASA». Уточните поиск, чтобы найти слои, содержащие Blue Marble: Next Generation, указав строку поиска как «bluemarbleng».
nasa = wmsfind('NASA','SearchField','serverurl'); layer = nasa.refine('bluemarbleng', ... 'SearchFields','layername', ... 'MatchType', 'exact');
Прочтите первый слой и отобразите карту.
[A,R] = wmsread(layer(1)); axesm globe axis off geoshow(A,R) title('Blue Marble')
Прочтите и отобразите ортоизображение северной части моста Золотые ворота в Сан-Франциско, Калифорния, используя сервер USGS National Map Seamless.
Сначала задайте пределы широты и долготы моста Золотые ворота.
latlim = [37.78 37.84]; lonlim = [-122.53 -122.40];
Затем найдите слой ортоимагерии высокого разрешения 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); 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'})
Считывайте и отображайте глобальный ежемесячный композит температуры поверхности моря на основе данных с датчика AMSR-E на борту спутника Aqua.
Сначала найдите в базе данных WMS слои, содержащие слово «coastwatch». Уточните поиск, чтобы найти слои с датчика AMSR-E.
coastwatch = wmsfind('coastwatch','SearchField','serverurl'); layers = refine(coastwatch,'erdAAsstamday','Searchfield','serverurl');
Ознакомьтесь и отобразите композит температуры поверхности моря за 16 апреля 2010 года. Включать береговую линию, ландмаск и национальные слои.
time = '2010-04-16T00:00:00Z'; [A,R] = wmsread(layers(end:-1:1),'Time',time); 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 поддерживает чтение данных в версиях WMS 1.0.0, 1.1.1 и 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 значения параметров, если они допустимы.
WebMapServer | wmsfind | wmsinfo | WMSLayer | wmsupdate