usgs24kdem

Считайте минуту USGS 7.5 (30 метров или 10 метров) Цифровые Модели Вертикального изменения

Описание

[lat,lon,Z] = usgs24kdem читает USGS 1:24,000 файл цифровой карты вертикального изменения (DEM) в стандартном формате. Вы выбираете файл в интерактивном режиме. usgs24kdem читает целый файл, подпроизведенный фактором 5, возвращая определенную геолокацию сетку данных с массивом широты, lat, массив долготы, lon, и массив вертикального изменения, Z. Горизонтальные модули в градусах, вертикальные модули могут варьироваться. 1:24,000 серии демократов хранят данные как сетку вертикальных изменений, расположенных с интервалами или в 10 или на расстоянии в 30 метров. Число точек в файле меняется в зависимости от географического положения.

[lat,lon,Z] = usgs24kdem(filename) читает USGS DEM, заданный filename и возвращает результат как определенную геолокацию сетку данных.

пример

[lat,lon,Z] = usgs24kdem(filename,samplefactor) читает подмножество данных DEM из filename, где samplefactor скалярное целое число, которое задает демонстрационную частоту.

[lat,lon,Z] = usgs24kdem(filename,samplefactor,latlim,lonlim) читает подмножество данных о вертикальном изменении из filename заданный двухэлементными векторами latlim и lonlim. Вы задаете пределы широты и долготы в градусах. Элементы в векторах должны быть в порядке возрастания. Данные могут расширить вне требуемой области.

[lat,lon,Z] = usgs24kdem(filename,samplefactor,latlim,lonlim,gsize) задает размер координатной сетки в gsize. gsize двухэлементный вектор, задающий количество строк и столбцов в широте, и долгота скоординировала сетку.

[lat,lon,Z,header,profile] = usgs24kdem(___) также возвращает содержимое заголовка и необработанные профили файла DEM. header структура содержит описания данных из заголовка файла. profile структура является необработанными данными о профиле, из которых создается определенная геолокацию сетка данных.

Примеры

свернуть все

В этом примере показано, как считать USGS 24K Цифровой файл Модели Вертикального изменения.

Разархивируйте USGS 24K файл DEM. Тулбокс включает файл DEM sanfranciscos.dem.gz .

filenames = gunzip('sanfranciscos.dem.gz', tempdir); 
demFilename = filenames{1};

Считайте любую точку 1:24,000 файл DEM.

[lat,lon,Z,header,profile] = usgs24kdem(demFilename,2);

Удалите временный gunzipped файл.

delete(demFilename);

Когда никакие отрицательные вертикальные изменения не существуют, перемещают все точки на уровне моря к-1, чтобы окрасить их в синий.

Z(Z==0) = -1;

Вычислите пределы широты и долготы для DEM.

latlim = [min(lat(:)) max(lat(:))]
latlim = 1×2

   37.6249   37.7504

lonlim = [min(lon(:)) max(lon(:))]
lonlim = 1×2

 -122.5008 -122.3740

Отобразите значения DEM.

figure
usamap(latlim, lonlim)
geoshow(lat, lon, Z, 'DisplayType','surface')
demcmap(Z)
daspectm('m',1)

Исследуйте метаданные в заголовке.

header
header = struct with fields:
                    Quadranglename: 'SAN FRANCISCO SOUTH, CA  BIG BASIN DEM  '
                       TextualInfo: 'WMC                             CTOG    '
                            Filler: ''
                       ProcessCode: ''
                           Filler2: ''
                SectionalIndicator: ''
                      MCoriginCode: ''
                      DEMlevelCode: 2
              ElevationPatternCode: 'regular'
    PlanimetricReferenceSystemCode: 'UTM'
                              Zone: 10
              ProjectionParameters: [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0]
                   HorizontalUnits: 'meters'
                    ElevationUnits: 'feet'
               NsidesToBoundingBox: 4
                       BoundingBox: [1x8 double]
                  MinMaxElevations: [0 1314]
                     RotationAngle: 0
                      AccuracyCode: 'accuracy information in record C'
                    XYZresolutions: [30 30 1]
                         NrowsCols: [1 371]
                    MaxPcontourInt: NaN
                SourceMaxCintUnits: 'unknown'
                   SmallestPrimary: NaN
                SourceMinCintUnits: 'unknown'
                    DataSourceDate: NaN
                   DataInspRevDate: NaN
                       InspRevFlag: ''
                DataValidationFlag: 'unknown'
                   SuspectVoidFlag: 'unknown'
                     VerticalDatum: 'unknown'
                   HorizontalDatum: 'unknown'
                       DataEdition: NaN
                       PercentVoid: NaN

Входные параметры

свернуть все

Имя файла, содержащего цифровую карту вертикального изменения, заданную как скаляр строки или символьный массив.

Типы данных: char | string

Фактор выборки данных, заданный как скалярное целое число. Например, если samplefactor равно 1, usgs24kdem считывает данные в его полном разрешении, то есть, каждом пикселе. Если вы задаете samplefactor значение n это больше 1, usgs24kdem читает каждый nточка th.

Типы данных: double

Пределы желаемых данных, заданных как двухэлементный вектор, в градусах. Пределы должны быть в порядке возрастания. Данные могут расширить вне требуемой области.

Типы данных: double

Пределы желаемых данных, заданных как двухэлементный вектор, в градусах.

Типы данных: double

Размер координатной сетки, заданный как двухэлементный вектор. gsize задает количество строк и столбцов в широте, и долгота скоординировала сетку. Если не использовано, usgs24kdem возвращает координатную сетку тот же размер как определенная геолокацию сетка данных. Чтобы задать скоординированный размер сетки, не задавая географические пределы, используйте пустые матрицы в latlim и lonlim.

Типы данных: double

Выходные аргументы

свернуть все

Массив широты, возвращенный как матрица класса double.

Массив долготы, возвращенный как матрица класса double.

Массив вертикального изменения, возвращенный как матрица класса double.

Описания данных из заголовка файла, возвращенного как struct.

Необработанные данные о профиле, из которых создается определенная геолокацию сетка данных, возвратились как struct.

Советы

  • Американская Геологическая служба создала ряд цифровых моделей вертикального изменения на основе их статьи 1:24,000 карты шкалы. Интервал сетки для этих моделей вертикальных изменений составляет или 10 метров или 30 метров на Универсальной Поперечной Меркаторской сетке. Каждый файл покрывает 7.5-минутный четырехугольник. Карта и ряд данных доступны для большой части совпадающих Соединенных Штатов, Гавайев и Пуэрто-Рико. Данные были выпущены в нескольких форматах. Эта функция считывает данные в “стандартном” формате файла.

  • Эта функция читает файлы USGS DEM, хранившие в проекции UTM. Функциональные непроекты сетка назад к широте и долготе. Используйте usgsdem для данных, хранимых в географических сетках.

  • Число точек в файле меняется в зависимости от географического положения. В отличие от продуктов USGS DEM, которые используют сетку равного угла, демократы сетки проекции UTM не могут просто быть конкатенированы, чтобы покрыть большие области. Между демократами могут быть разрывы данных.

  • Можно получить файлы данных из американской Геологической службы и от коммерческих поставщиков. Другие агентства сделали некоторые доступные данные местности онлайн. Смотрите Находят Картографические данные Онлайн. Файлы DEM являются ASCII-файлами и могут быть переданы как текст. Заканчивающее линию преобразование не обязательно требуется.

Смотрите также

| | | | | |

Представлено до R2006a

Для просмотра документации необходимо авторизоваться на сайте