gtopo30
Читайте 30, образуют дугу вторые глобальные цифровые данные о повышении (GTOPO30)
Синтаксис
[Z,refvec] = gtopo30(tilename)
[Z,refvec] = gtopo30(tilename,samplefactor)
[Z,refvec] = gtopo30(tilename,samplefactor,latlim,lonlim)
[Z,refvec] = gtopo30(foldername, ...)
Описание
[Z,refvec] = gtopo30(tilename) читает мозаику GTOPO30, заданную tilename, и возвращает результат как обычную сетку данных. tilename является скаляром строки или вектором символов, который не включает расширение и указывает на мозаику GTOPO30 в текущей папке или на пути MATLAB®. Если tilename будет пуст или не использован, файловый браузер откроется для интерактивного выбора заголовочного файла GTOPO30. Данные возвращены в полном разрешении с пределами широты и долготы, определенными от мозаики GTOPO30. Сетка данных, Z, возвращена как массив повышений. Повышения даны в метрах над средним уровнем моря с помощью WGS84 в качестве горизонтальной данной величины. refvec является связанным вектором ссылки.
[Z,refvec] = gtopo30(tilename,samplefactor) читает подмножество данных о повышении из tilename. samplefactor является скалярным целым числом, которое, когда равный 1 считывает данные в его полном разрешении. Когда samplefactor является целым числом n больше, чем один, каждая энная точка читается. Если samplefactor не использован или пуст, он принимает значение по умолчанию к 1.
[Z,refvec] = gtopo30(tilename,samplefactor,latlim,lonlim) читает подмножество данных о повышении из tilename с помощью широты, и долгота ограничивает latlim и lonlim, заданный в градусах. latlim является двухэлементным вектором формы:
[southern_limit northern_limit]
Аналогично, lonlim имеет форму:
[western_limit eastern_limit]
Если latlim и lonlim не использованы, координатные пределы определяются из файла. Пределы широты и долготы сфотографированы исходящие, чтобы задать самую маленькую прямоугольную сетку ячеек GTOPO30, которая полностью заключает область, заданную входными пределами. Любые ячейки в этой сетке, которые выходят за пределы степени мозаики, заполнены NaN.
[Z,refvec] = gtopo30(foldername, ...) подобно синтаксисам выше за исключением того, что данные GTOPO30 считаны и конкатенированы от нескольких мозаик в GTOPO30 CD-ROM или структуре папок. Вход foldername является скаляром строки или вектором символов с именем папки, которая содержит папки мозаики GTOPO30 или мозаики GTOPO30. В мозаике папки являются несжатыми файлами данных. foldername для CD-ROMs, распределенных USGS, является именем устройства дисковода для компакт-дисков. Как со случаем с одной мозаикой, любыми ячейками в сетке, заданной latlim и lonlim, является NaN, заполненный, если они не покрыты мозаикой в foldername. samplefactor, если не использовано или пустые значения по умолчанию к 1. latlim, если не использовано или пустые значения по умолчанию к [-90 90]. lonlim, если не использовано или пустые значения по умолчанию к [-180 180].
Для получения дополнительной информации при определении местоположения данных GTOPO30 для загрузки по Интернету, смотрите, Находят Картографические данные Онлайн.
Примеры
Отобразите сетку данных и наложите контур Стейтлайна
Чтобы запустить этот пример, необходимо загрузить набор данных GTOPO30. Для получения дополнительной информации при определении местоположения этого набора данных для загрузки по Интернету, смотрите, Находят Картографические данные Онлайн.
Извлеките и отобразите данные о полном разрешении для Массачусетса. Считайте контур полигона Стейтлайна и вычислите граничные пределы.
Massachusetts = shaperead('usastatehi','UseGeoCoords',true, ... 'Selector',{@(name) strcmpi(name,'Massachusetts'),'Name'}); latlim = [min(Massachusetts.Lat(:)) max(Massachusetts.Lat(:))]; lonlim = [min(Massachusetts.Lon(:)) max(Massachusetts.Lon(:))];
Считайте данные GTOPO30 в полном разрешении.
[Z,refvec] = gtopo30('W100N90',1,latlim,lonlim);Отобразите сетку данных и наложите контур Стейтлайна.
figure ax = usamap(Z,refvec); ax.SortMethod = 'ChildOrder'; geoshow(Z,refvec,'DisplayType','surface') demcmap(Z) geoshow(Massachusetts,'DisplayType','polygon',... 'facecolor','none','edgecolor','y')
Извлеките каждую 20-ю точку от набора данных
Чтобы запустить этот пример, необходимо загрузить набор данных GTOPO30. Для получения дополнительной информации при определении местоположения этого набора данных для загрузки по Интернету, смотрите, Находят Картографические данные Онлайн.
Извлеките каждую 20-ю точку от мозаики. Обеспечьте пустое имя файла и выберите файл в интерактивном режиме.
[Z,refvec] = gtopo30([],20);
Извлеките данные, которые колеблются между двумя мозаиками
Чтобы запустить этот пример, необходимо загрузить набор данных GTOPO30. Для получения дополнительной информации при определении местоположения этого набора данных для загрузки по Интернету, смотрите, Находят Картографические данные Онлайн.
Извлеките данные для Таиланда, область, которая колеблется между двумя мозаиками.
latlim = [5.22 20.90]; lonlim = [96.72 106.38]; gtopo30s(latlim,lonlim)
Извлеките каждую пятую точку данных для Таиланда. Код принимает, что данные находятся на CD номер 3, распределенный USGS. Устройство CD является 'F:\'
[Z,refvec] = gtopo30('F:\',5,latlim,lonlim);
worldmap(Z,refvec);
geoshow(Z,refvec,'DisplayType','surface')
demcmap(Z)
Извлеките каждую 10-ю точку из столбца данных
Чтобы запустить этот пример, необходимо загрузить набор данных GTOPO30. Для получения дополнительной информации при определении местоположения этого набора данных для загрузки по Интернету, смотрите, Находят Картографические данные Онлайн.
Извлеките каждую 10-ю точку из столбца данных 5 градусов вокруг главного меридиана. Этот пример принимает, что текущая папка содержит данные GTOPO30.
[Z,refvec] = gtopo30(pwd,10,[],[-5 5]);