Получение карты WMS с сервера
[ осуществляет доступ к Интернету для визуализации и извлечения растровой карты с сервера службы веб-карт (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 содержит WMS serverURL с дополнительными параметрами WMS. Эти параметры WMS включают в себя BBOX, GetMap и EPSG:4326 или CRS:84 ключевое слово. Получить mapRequestURL из выходных данных wmsread, RequestURL свойство WMSMapRequest объект или поиск в Интернете.
[ возвращает WMS A,R,mapRequestURL] = wmsread(...)GetMap URL запроса в символьном векторе mapRequestURL. Можно вставить mapRequestURL в браузер, чтобы сделать запрос на сервер, который затем возвращает растровую карту. Обозреватель открывает возвращаемую карту, если ее тип mime понятен, или сохраняет растровую карту на диск.
Поиск в базе данных 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 слоев, содержащих слово «береговые часы». Уточните поиск, чтобы найти слои от датчика 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 секунд.
Чтобы указать прокси-сервер для подключения к Интернету, выберите «Файл» > «Установки» > «Веб» и введите сведения о прокси-сервере. Эта функция используется при наличии брандмауэра.
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