Объедините отличающиеся сетки путем преобразования обычной сетки в определенную геолокацию сетку данных

В этом примере показано, как объединить сетку данных о вертикальном изменении и атрибут (цвет) сетка данных, которые покрывают ту же область, но с координатной сеткой по-другому. Пример драпирует наклонные данные из обычной сетки данных сверху данных о вертикальном изменении из определенной геолокацию сетки данных. Пример использует определенную геолокацию сетку данных в качестве источника для поверхностных вертикальных изменений и преобразовывает обычную сетку данных в наклонные значения, которые затем произведены, чтобы соответствовать определенной геолокацию сетке данных (создающий набор наклонных значений для ромбовидной сетки) и нанесены цветную маркировку для поверхностного отображения. Этот подход работает с любыми отличающимися сетками, несмотря на то, что эти два набора данных в этом примере на самом деле имеют тот же источник (определенная геолокацию сетка выводит из topo60c набор данных).

Загрузите определенную геолокацию сетку данных от mapmtx файл и обычная сетка данных от topo60c файл. mapmtx файл на самом деле содержит две области, но этот пример только использует ромбовидный фрагмент, lt1, lg1, и map1, сосредоточенный на Ближнем Востоке.

load mapmtx lt1 lg1 map1 
load topo60c

Вычислите поверхностный аспект, наклон и градиенты для topo60c. Этот пример только использует наклоны на последующих шагах.

[aspect,slope,gradN,gradE] = gradientm(topo60c,topo60cR);

Используйте geointerp функция, чтобы интерполировать наклонные значения к определенной геолокацию сетке, заданной lt1 и lg1 . Выход 50 50 сетка вертикальных изменений, совпадающих с покрытием map1 переменная.

slope1 = geointerp(slope,topo60cR,lt1,lg1);

Настройте фигуру с проекцией Миллера и используйте surfm отобразить наклонные данные. Задайте z - значения для поверхности явным образом как map1 данные, которые являются вертикальным изменением ландшафта. Карта в основном изображает крутые утесы, которые представляют горы (Гималаи на северо-востоке), и континентальные шельфы и канавки.

figure 
axesm miller
surfm(lt1,lg1,slope1,map1)

Окраска изображает крутизну наклона. Измените палитру, чтобы сделать самый крутой пурпурный наклонов, более нежные наклоны темно-синий, и плоские области голубой:

colormap cool

Используйте view получить юго-восточную перспективу поверхности с низкой точки зрения. В 3-D вы сразу видите топографию, а также наклон.

view(20,30)
daspectm('meter',200)

Использование рендеринга значения по умолчанию фасетировало штриховку (никакая сплайн-интерполяция). Представьте поверхность снова, на этот раз делая его блестящим со штриховкой Гура и подсветкой с востока (значение по умолчанию camlight световые сигналы появляются из-за правого плеча средства просмотра).

material shiny
camlight
lighting Gouraud

Наконец, удалите пробел и повторно представьте фигуру в перспективном режиме.

axis tight
ax = gca;
ax.Projection = 'perspective';