exponenta event banner

readgeoraster

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

Описание

пример

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

пример

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

пример

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

Примеры

свернуть все

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

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

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

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

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

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

Figure contains an axes. The axes contains an object of type image.

Данные, использованные в этом примере, предоставлены Национальным управлением по исследованию океанов и атмосферы (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 функция возвращает логический 1 (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', 'single', 'double', 'int16', 'int32', 'int64', 'uint8', 'uint16', 'uint32', 'uint64', или '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-by-N или M-by-N-by-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 пуст. Пространственную привязку можно указать как 'geographic' или 'planar' с использованием 'CoordinateSystemType' пара имя-значение.

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

Подробнее

свернуть все

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

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

Формат файлаРасширение

GeoTIFF

.tif или .tiff

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

.adf

Сетка Esri ASCII

.asc или .grd

Esri GridFloat

.flt

DTED

.dt0, .dt1, или .dt2

SDTS

.DDF

USGS DEM

.dem

Мэпперы ER

.ers

ENVI

.dat

ERDAS IMAGINE

.img

Числовая сетка отображения по вертикали

.grd

Классифицированная сетка вертикального отображения

.grc

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

На форматы файлов можно ссылаться с использованием различных имен. Например, формат Esri GridFloat также может называться Esri .hdr Маркированные или ITT ESRI .hdr Растр RAW. Формат двоичной сетки Esri может также упоминаться как двоичная сетка ArcGrid, двоичная сетка Esri ArcGIS или сетка Esri ArcInfo.

Совет

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

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

Представлен в R2020a