Считайте минуту USGS 7.5 (30 метров или 10 метров) Цифровые Модели Вертикального изменения
[ читает USGS 1:24,000 файл цифровой карты вертикального изменения (DEM) в стандартном формате. Вы выбираете файл в интерактивном режиме. lat,lon,Z]
= usgs24kdemusgs24kdem читает целый файл, подпроизведенный фактором 5, возвращая определенную геолокацию сетку данных с массивом широты, lat, массив долготы, lon, и массив вертикального изменения, Z. Горизонтальные модули в градусах, вертикальные модули могут варьироваться. 1:24,000 серии демократов хранят данные как сетку вертикальных изменений, расположенных с интервалами или в 10 или на расстоянии в 30 метров. Число точек в файле меняется в зависимости от географического положения.
[ читает подмножество данных DEM из lat,lon,Z]
= usgs24kdem(filename,samplefactor)filename, где samplefactor скалярное целое число, которое задает демонстрационную частоту.
[ читает подмножество данных о вертикальном изменении из lat,lon,Z]
= usgs24kdem(filename,samplefactor,latlim,lonlim)filename заданный двухэлементными векторами latlim и lonlim. Вы задаете пределы широты и долготы в градусах. Элементы в векторах должны быть в порядке возрастания. Данные могут расширить вне требуемой области.
[ также возвращает содержимое заголовка и необработанные профили файла DEM. lat,lon,Z,header,profile]
= usgs24kdem(___)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
Американская Геологическая служба создала ряд цифровых моделей вертикального изменения на основе их статьи 1:24,000 карты шкалы. Интервал сетки для этих моделей вертикальных изменений составляет или 10 метров или 30 метров на Универсальной Поперечной Меркаторской сетке. Каждый файл покрывает 7.5-минутный четырехугольник. Карта и ряд данных доступны для большой части совпадающих Соединенных Штатов, Гавайев и Пуэрто-Рико. Данные были выпущены в нескольких форматах. Эта функция считывает данные в “стандартном” формате файла.
Эта функция читает файлы USGS DEM, хранившие в проекции UTM. Функциональные непроекты сетка назад к широте и долготе. Используйте usgsdem для данных, хранимых в географических сетках.
Число точек в файле меняется в зависимости от географического положения. В отличие от продуктов USGS DEM, которые используют сетку равного угла, демократы сетки проекции UTM не могут просто быть конкатенированы, чтобы покрыть большие области. Между демократами могут быть разрывы данных.
Можно получить файлы данных из американской Геологической службы и от коммерческих поставщиков. Другие агентства сделали некоторые доступные данные местности онлайн. Смотрите Находят Картографические данные Онлайн. Файлы DEM являются ASCII-файлами и могут быть переданы как текст. Заканчивающее линию преобразование не обязательно требуется.