Exponenta Banner

usgs24kdem

Считайте минуту USGS 7.5 (30 метров или 10 метров) Цифровые Модели Повышения

свернуть все на странице

Синтаксис

[lat,lon,Z] = usgs24kdem
[lat,lon,Z] = usgs24kdem(filename)
[lat,lon,Z] = usgs24kdem(filename,samplefactor)
[lat,lon,Z] = usgs24kdem(filename,samplefactor,latlim,lonlim)
[lat,lon,Z] = usgs24kdem(filename,samplefactor,latlim,lonlim,gsize)
[lat,lon,Z,header,profile] = usgs24kdem(___)

Описание

[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 является необработанными данными о профиле, из которых создается определенная геолокацию сетка данных.

Примеры

свернуть все

Скрипт Open Live Script

Этот пример показывает, как считать 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

Документация Mapping Toolbox
Поддержка