readgeoraster

Считайте геопространственный растровый файл данных

Описание

пример

[A,R] = readgeoraster(filename) создает массив путем чтения географических или спроектированных растровых данных из файла. Выходной аргумент R содержит пространственную информацию о ссылке для массива. Поддерживаемые форматы файлов включают Бинарную Сетку Esri, Esri GridFloat, GeoTIFF и DTED. Для полного списка поддерживаемых форматов смотрите Поддерживаемые форматы и Расширения.

пример

[A,R] = readgeoraster(___,Name,Value) задает опции с помощью одного или нескольких Name,Value парные аргументы.

пример

[___,cmap] = readgeoraster(___) также возвращает палитру A.

Примеры

свернуть все

Считайте изображение GeoTIFF Бостона как массив и объект ссылки ячеек карты. Массив имеет размер 2881 4 481 3 и задает красные, зеленые, и синие компоненты изображения. Отобразите изображение с помощью mapshow функция.

[A,R] = readgeoraster('boston.tif');
mapshow(A,R)

Данные, используемые в этом примере, включают материал, защищенный авторским правом GeoEye, все права защищены.

Считайте и отобразите классификацию типов поверхности Оаху, Гавайи.

Во-первых, считайте данные о типе поверхности как массив, объект ссылки ячеек карты и палитру. Элементы A индексируйте в палитру. Каждая строка палитры задает красные, зеленые, и синие компоненты одного цвета. Затем отобразите данные о типе поверхности.

[A,R,cmap] = readgeoraster('oahu_landcover.img');
mapshow(A,cmap,R)

Данные, используемые в этом примере, являются любезностью Национального управления океанических и атмосферных исследований (NOAA).

Считайте и отобразите данные о вертикальном изменении для области вокруг Южного Пика Валуна в Колорадо.

Во-первых, считайте данные о вертикальном изменении как массив и географический объект ссылки регистраций. Отобразить данные как поверхность, geoshow функция требует данных типа double или single. В этом случае сохраните точность путем определения выходного типа как 'double'.

[A,R] = readgeoraster('n39_w106_3arc_v2.dt1','OutputType','double');

Создайте карту. Во-первых, создайте оси карты путем определения пределов широты и долготы данных. Затем отобразите данные как поверхность с помощью geoshow функция. Примените палитру, подходящую для данных о вертикальном изменении с помощью demcmap функция.

latlim = R.LatitudeLimits;
lonlim = R.LongitudeLimits;
usamap(latlim,lonlim)
geoshow(A,R,'DisplayType','surface')
demcmap(A)

Данные о вертикальном изменении, используемые в этом примере, являются любезностью Геологической службы США.

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

Импортируйте растровые данные и ссылочный объект с помощью readgeoraster функция. Найдите индикатор недостающих данных с помощью georasterinfo функция.

[A,R] = readgeoraster('MtWashington-ft.grd');
info = georasterinfo('MtWashington-ft.grd');
m = info.MissingDataIndicator
m = -32766

Проверьте, что растровые данные содержат недостающие данные с помощью ismember функция. ismember функция возвращает логическую единицу (true) если растр содержит индикатор недостающих данных.

ismember(m,A)
ans = logical
   1

Замените недостающие данные на NaN значения с помощью standardizeMissing функция.

A = standardizeMissing(A,m);

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

свернуть все

Имя файла, чтобы читать в виде вектора символов или строкового скаляра. Форма filename зависит от местоположения вашего файла.

  • Если файл находится в вашей текущей папке или в папке на пути MATLAB®, то задайте имя файла, такого как 'myFile.dem'.

  • Если файл не находится в текущей папке или в папке на пути MATLAB, то задайте полное имя или относительный путь, такой как 'C:\myfolder\myFile.tif' или 'dataDir\myFile.dat'.

Для списка поддерживаемых форматов файлов смотрите Поддерживаемые форматы и Расширения.

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

Аргументы в виде пар имя-значение

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример: 'OutputType','double','Bands',1:2

Тип данных для AВ виде разделенной запятой пары, состоящей из 'OutputType' и вектор символов или строковый скаляр, содержащий одно из этих значений: 'native'единственный'double'int16int32int64uint8uint16uint32uint64, или 'logical'.

Значение по умолчанию для 'OutputType' 'native', который возвращает A использование собственного типа данных встраивается в filename. Используя тип данных кроме 'native' может привести к потере точности.

Пример: 'OutputType','double'

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

Полосы, чтобы читать в виде разделенной запятой пары, состоящей из 'Bands' и 'all', положительное целое число или вектор положительных целых чисел. Например, если вы задаете значение 3, readgeoraster читает третью полосу в файле. Полосы возвращены в заданном порядке.

Значение по умолчанию для 'Bands' 'all', где readgeoraster чтения все полосы в файле.

Пример: 'Bands',3

Тип системы координат для RВ виде одного из этих значений:

  • 'auto' – Возвращает R когда объект растровой привязки определяется содержимым файла.

  • 'geographic' – Возвращает R как географический объект ссылки ячеек или регистраций.

  • 'planar' – Возвращает R как карта ячейки или объект ссылки регистраций.

Задайте тип системы координат, когда ваши данные не будут содержать информацию о проекции.

Пример: 'CoordinateSystemType','geographic'

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

свернуть все

Изображение, на которое геоссылаются, или сетка данных, возвращенная как M на n или M N P числовым массивом.

По умолчанию, тип данных A совпадает с собственным типом данных, встроенным в filename. Задайте тип данных с помощью 'OutputType' пара "имя-значение".

Независимо от того, как данные закодированы, первая строка A представляет самые северные данные и последнюю строку A представляет самые южные данные.

Пространственная ссылка для A, возвращенный как GeographicCellsReference объект, GeographicPostingsReference объект, MapCellsReference объект или MapPostingsReference объект. Значение R зависит от данных в filename:

  • Если данные в filename ссылается к географической системе координат, затем R GeographicCellsReference объект или GeographicPostingsReference объект.

  • Если данные в filename ссылается к спроектированной системе координат, затем R MapCellsReference объект или MapPostingsReference объект.

Если файл не содержит достаточно информации, чтобы определить, спроектированы ли данные или географические, то R MapCellsReference или MapPostingsReference объект. Если файл не содержит допустимой пространственной справочной информации, то R isempty. Можно задать пространственную ссылку как 'geographic' или 'planar' использование 'CoordinateSystemType' пара "имя-значение".

Палитра сопоставлена с индексируемым изображением, возвращенным как n-3 числовая матрица со значениями в области значений [0,1]. Каждая строка cmap трехэлементный триплет RGB, который задает красные, зеленые, и синие компоненты одного цвета в палитре. Значение cmap пусто если A индексируемое изображение.

Больше о

свернуть все

Поддерживаемые форматы и расширения

readgeoraster и georasterinfo функции поддерживают эти форматы файлов и расширения. В некоторых случаях можно считать поддерживаемые форматы файлов с помощью расширений кроме тех перечисленных.

FileFormat Расширение

GeoTIFF

.tif или .tiff

Двоичная сетка Esri

.adf

Сетка ASCII Esri

.asc или .grd

Esri GridFloat

.flt

DTED

.dt0, .dt1, или .dt2

SDTS

.DDF

DEM USGS

.dem

Картопостроитель ER ERS

.ers

ENVI

.dat

ЭРДАС ИМАГИНЯ

.img

Некоторые форматы файлов состоят из файла данных и нескольких вспомогательных файлов. Например, файлы Esri GridFloat могут иметь заголовочные файлы поддержки (.hdr). Когда вы читаете файл данных со вспомогательными файлами с помощью readgeoraster или georasterinfo, задайте расширение файла данных.

Форматы файлов могут быть отнесены в использование различных имен. Например, формат Esri GridFloat может также упоминаться как Esri .hdr Помеченный или ITT ESRI .hdr НЕОБРАБОТАННЫЙ Растр. Бинарный формат Сетки Esri может также упоминаться как Двоичный файл ArcGrid, Двоичная Сетка Esri ArcGIS или Сетка Esri ArcInfo.

Советы

  • Некоторые функции требуют входных параметров типа single или double, такой как geoshow функция для того, чтобы отобразить поверхности. Использовать выход readgeoraster с этими функциями задайте выходной тип как 'single' или 'double' использование 'OutputType' пара "имя-значение".

  • Независимо от формата файла, массив, возвращенный readgeoraster имеет столбцы, запускающиеся с севера и ColumnsStartFrom свойство ссылочного объекта имеет значение 'north'.

Введенный в R2020a

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