Наложение растра

Не знаете, где задать вопрос? Задавайте здесь.
Ответить
dimen4ik
Новоприбывший
Сообщения: 2
Зарегистрирован: 03 мар 2017, 22:19
Репутация: 0

Наложение растра

Сообщение dimen4ik » 28 фев 2018, 11:20

Здравствуйте, имеется 2 растра со значениями(допустим высота и уклон поверхности), как можно их объединить, чтобы можно было раскрасить растр одновременно по двум параметрам, например уклон от белого до красного, а высоту от белого до черного?

Аватара пользователя
rhot
Гуру
Сообщения: 1727
Зарегистрирован: 25 янв 2011, 17:50
Репутация: 194
Ваше звание: доктор
Откуда: Архангельск

Re: Наложение растра

Сообщение rhot » 28 фев 2018, 13:01

1) Растры перегнать и объединить в одну таблицу.
2) Использовать пакет colorplaner для ggplot в R.

Код: Выделить всё

library(sp)
library(raster)
library(ggplot2)
library(colorplaner)

dem <- getData('SRTM', lon=5, lat=45)

elev <- crop(dem, extent(6, 6.5, 44.5, 44.75))
slope <- terrain(elev, unit='degrees')

br <- brick(elev, slope)
names(br) <- c('elevation', 'slope')

spdf.br <- as(br, "SpatialPixelsDataFrame")
df.br <- as.data.frame(spdf.br)

ggplot(df.br, aes(x, y, fill=elevation, fill2=slope))+
  geom_tile()+
  scale_x_continuous(breaks=seq(6.1, 6.4, .1))+
  scale_y_continuous(breaks=seq(44.5, 44.7, .1))+
  scale_fill_colorplane(name='',
                        color_projection = interpolate_projection,
                        zero_color = "white",
                        horizontal_color = "black",
                        vertical_color = "red",
                        axis_title='высота',
                        axis_title_y='уклон')+
  theme(panel.backgroun=element_blank())+
  guides(fill=guide_colorplane(
    label.y.theme=theme_gray(),
    label.theme=element_text(angle=90, vjust=0.5)))+
  labs(x='Долгота', y='Широта', title='Где-то в Альпийских горах...')
gislab_upload.png
gislab_upload.png (190.51 КБ) 3242 просмотра
Последний раз редактировалось rhot 28 фев 2018, 14:29, всего редактировалось 2 раза.
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________

alexandr cherepanov
Гуру
Сообщения: 534
Зарегистрирован: 30 ноя 2006, 13:31
Репутация: 116
Откуда: Moscow

Re: Наложение растра

Сообщение alexandr cherepanov » 28 фев 2018, 13:38

можно по аналогии с этой статьей
http://gis-lab.info/qa/gdal-hillshade-colormap.html

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

Аватара пользователя
rhot
Гуру
Сообщения: 1727
Зарегистрирован: 25 янв 2011, 17:50
Репутация: 194
Ваше звание: доктор
Откуда: Архангельск

Re: Наложение растра

Сообщение rhot » 28 фев 2018, 14:54

Добавил решение в R.
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________

Ответить

Вернуться в «Я новичок!»

Кто сейчас на конференции

Сейчас этот форум просматривают: Ahrefs [Bot] и 10 гостей