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]);