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

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

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

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

Сначала считайте данные типы поверхности как массив, карту камер ссылки объект и палитру. Элементы 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 функция. The ismember функция возвращает логический 1 (true), если растр содержит индикатор недостающих данных.

ismember(m,A)

ans = logical
   1

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

A = standardizeMissing(A,m);