Геопространственные растровые провайдеры данных обычно данные о пакете как смежные мозаики. Например, 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)
Данные о вертикальном изменении, используемые в этом примере, от американской Геологической службы.