Геопространственные растровые провайдеры данных обычно данные о пакете как смежные мозаики. Например, SRTM Пусто Заполненные данные о вертикальном изменении упакован в 1 степень 1 мозаикой степени. Чтобы анализировать распространение данных через несколько мозаик, такой, вычисляя видимость угла обзора, вы должны мозаика мозаики в один растр.
Процессы для объединения растров ячеек и растров регистрации точек отличаются. Для получения информации о дифференцирующихся растрах, составленных из ячеек и отправляющих точек, смотрите Пространственно Ссылочные Импортированные Растры.
Прежде mosaicking мозаики, необходимо гарантировать, что мозаики смежны и что контуры мозаики выравниваются. В противном случае на элементы mosaicked растра можно пространственно сослаться к неправильным местоположениям.
В этом примере показано, как импортировать два смежных растровых файла данных, сделанные из ячеек, мозаика данные в один растр, и отобразить mosaicked растр на карте.
Файлы используются в этом примере, concord_ortho_e.tif
и concord_ortho_w.tif
, содержите смежные изображения восток - запад с concord_ortho_e.tif
на восток concord_ortho_w.tif
. Степени ячейки и мировые y-пределы изображений идентичны.
Импортируйте изображения с помощью readgeoraster
функция. Выходные параметры eastA
и westA
массивы, которые содержат данные изображения и выходные параметры westR
и eastR
MapCellsReference
объекты, которые содержат пространственную информацию о ссылке.
[eastA,eastR] = readgeoraster('concord_ortho_e.tif'); [westA,westR] = readgeoraster('concord_ortho_w.tif');
Создайте mosaicked растр путем объединения массивов. Изображения смежны и составлены из ячеек, таким образом, восточный контур westA
выравнивается с западным контуром eastA
.
mosaicA = [westA eastA];
Пространственно сошлитесь на mosaicked растр путем создания объекта растровой привязки. Можно создать ссылочный объект для спроектированного растра ячеек при помощи maprefcells
функция.
Задайте мир x-и y-пределы растра. X-пределы mosaicked растра являются минимальным x-пределом западного растра и максимальным x-пределом восточного растра. Мировые y-пределы mosaicked растра совпадают с y-пределами импортированных растров. Создайте ссылочный объект.
xlimits = [westR.XWorldLimits(1) eastR.XWorldLimits(2)]; ylimits = westR.YWorldLimits; mosaicR = maprefcells(xlimits,ylimits,size(mosaicA));
Столбцы массивов, импортированных с помощью readgeoraster
функциональный запуск с севера. Поэтому установите ColumnsStartFrom
свойство ссылочного объекта к 'north'
.
mosaicR.ColumnsStartFrom = 'north';
Отобразите изображение mosaicked на карте при помощи mapshow
функция.
mapshow(mosaicA,mosaicR,'DisplayType','image')
В этом примере показано, как импортировать два смежных растровых файла данных, сделанные из регистрации точек, мозаика данные в один растр, и отобразить mosaicked растр на карте.
Файлы используются в этом примере, n39_w106_3arc_v2.dt1
и n40_w106_3arc_v2.dt1
, содержите между севером и югом смежные сетки вертикального изменения с n39_w106_3arc_v2.dt1
на юг n40_w106_3arc_v2.dt1
. Северный предел широты сетки в n39_w106_3arc_v2.dt1
совпадает с южным пределом широты сетки в n40_w106_3arc_v2.dt1
. Интервал регистрации точек и пределов долготы сеток идентичен.
Импортируйте сетки с помощью readgeoraster
функция. Выходные параметры southZ
и northZ
массивы типа double
это содержит данные о вертикальном изменении и выходные параметры southR
и northR
GeographicPostingsReference
объекты, которые содержат пространственную информацию о ссылке.
[southZ,southR] = readgeoraster('n39_w106_3arc_v2.dt1','OutputType','double'); [northZ,northR] = readgeoraster('n40_w106_3arc_v2.dt1','OutputType','double');
Создайте mosaicked растр путем объединения массивов. Контуры регистрации растров точки составлены из наиболее удаленных точек регистрации и столбцов массивов, импортированных с помощью readgeoraster
функциональный запуск с севера. Поэтому самая южная строка northZ
и самая северная строка southZ
то же самое. Чтобы избежать дублирующейся строки в mosaicked растре, удалите самую южную строку northZ
прежде, чем объединить массивы.
northZ(end,:) = []; mosaicZ = [northZ; southZ];
Пространственно сошлитесь на mosaicked растр путем создания объекта растровой привязки. Можно создать ссылочный объект для географического растра регистрации точек при помощи georefpostings
функция.
Задайте пределы широты и долготы растра. Пределы широты являются минимальной широтой южного растра и максимальной широтой северного растра. Пределы долготы mosaicked растра совпадают с пределами долготы импортированных растров. Создайте ссылочный объект.
latlim = [southR.LatitudeLimits(1) northR.LatitudeLimits(2)]; lonlim = southR.LongitudeLimits; mosaicR = georefpostings(latlim,lonlim,size(mosaicZ));
Установите ColumnsStartFrom
и GeographicCRS
свойства ссылочного объекта так, чтобы они совпадали со свойствами импортированных ссылочных объектов.
mosaicR.ColumnsStartFrom = southR.ColumnsStartFrom; mosaicR.GeographicCRS = southR.GeographicCRS;
Отобразите mosaicked растр на карте. Создайте карту оси путем определения пределов широты и долготы данных. Затем отобразите данные как поверхность при помощи geoshow
функция. Примените палитру, подходящую для данных о вертикальном изменении при помощи demcmap
функция.
usamap(mosaicR.LatitudeLimits,mosaicR.LongitudeLimits) geoshow(mosaicZ,mosaicR,'DisplayType','surface') demcmap(mosaicZ)
Данные о вертикальном изменении, используемые в этом примере, от американской Геологической службы.