Поставщики геопространственных растровых данных обычно упаковывают данные в виде смежных фрагментов. Например, данные отметок, заполненных полостями SRTM, упаковываются в плитки 1 градус за 1 градус. Для анализа данных, распределенных по нескольким плиткам, например, при расчете видимости линии визирования, необходимо создать мозаику плиток в одном растре.
Процессы объединения растров ячеек и растров точек проводки различны. Сведения о дифференциации растров, состоящих из ячеек, и точках разноски см. в разделе Пространственно ссылочные импортированные растры.
Перед захватом мозаичных плиток необходимо убедиться в том, что плитки расположены рядом и границы плиток выровнены. В противном случае элементы мозаичного растра могут быть пространственно привязаны к неправильным расположениям.
В этом примере показано, как импортировать два соседних файла растровых данных из ячеек, мозаику данных в один растр и отобразить мозаику растра на карте.
Файлы, используемые в этом примере, 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');
Создание мозаичного растра путем объединения массивов. Изображения смежны и состоят из клеток, поэтому восточная граница westA выравнивается по западной границе eastA.
mosaicA = [westA eastA];
Пространственная ссылка на мозаичный растр путем создания объекта ссылки на растр. Можно создать ссылочный объект для проецируемого растра ячеек с помощью команды maprefcells функция.
Укажите мировой предел x и y растра. X-пределы мозаичного растра являются минимальным x-пределом западного растра и максимальным x-пределом восточного растра. Мировые пределы y мозаичного растра такие же, как и пределы y импортированных растров. Создайте ссылочный объект.
xlimits = [westR.XWorldLimits(1) eastR.XWorldLimits(2)]; ylimits = westR.YWorldLimits; mosaicR = maprefcells(xlimits,ylimits,size(mosaicA));
Столбцы массивов, импортированных с помощью readgeoraster функция начинается с севера. Поэтому установите ColumnsStartFrom свойство объекта ссылки на 'north'.
mosaicR.ColumnsStartFrom = 'north';Отображение мозаичного изображения на карте с помощью mapshow функция.
mapshow(mosaicA,mosaicR,'DisplayType','image')
В этом примере показано, как импортировать два соседних файла растровых данных, состоящих из точек разноски, мозаичные данные в один растр и отобразить мозаичный растр на карте.
Файлы, используемые в этом примере, 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');
Создание мозаичного растра путем объединения массивов. Границы растров точек проводки состоят из самых внешних точек проводки и столбцов массивов, импортированных с помощью readgeoraster функция начинается с севера. Следовательно, самый южный ряд northZ и самый северный ряд из southZ одинаковые. Чтобы избежать дублирования строки в мозаичном растре, удалите самый южный ряд из northZ перед объединением массивов.
northZ(end,:) = []; mosaicZ = [northZ; southZ];
Пространственная ссылка на мозаичный растр путем создания объекта ссылки на растр. Можно создать ссылочный объект для географического растра точек проводки с помощью georefpostings функция.
Задайте пределы широты и долготы растра. Пределы широты - это минимальная широта южного растра и максимальная широта северного растра. Границы долготы мозаичного растра совпадают с границами долготы импортированных растров. Создайте ссылочный объект.
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;
Отображение мозаичного растра на карте. Создайте оси карты, указав пределы широты и долготы данных. Затем отобразите данные как поверхность с помощью geoshow функция. Примените карту цветов, соответствующую данным отметок, используя demcmap функция.
usamap(mosaicR.LatitudeLimits,mosaicR.LongitudeLimits) geoshow(mosaicZ,mosaicR,'DisplayType','surface') demcmap(mosaicZ)
Данные о высотных отметках, используемые в этом примере, получены из Геологической службы США.