Наложение растра
-
- Новоприбывший
- Сообщения: 2
- Зарегистрирован: 03 мар 2017, 22:19
- Репутация: 0
Наложение растра
Здравствуйте, имеется 2 растра со значениями(допустим высота и уклон поверхности), как можно их объединить, чтобы можно было раскрасить растр одновременно по двум параметрам, например уклон от белого до красного, а высоту от белого до черного?
- rhot
- Гуру
- Сообщения: 1727
- Зарегистрирован: 25 янв 2011, 17:50
- Репутация: 194
- Ваше звание: доктор
- Откуда: Архангельск
Re: Наложение растра
1) Растры перегнать и объединить в одну таблицу.
2) Использовать пакет colorplaner для ggplot в R.
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='Где-то в Альпийских горах...')
Последний раз редактировалось rhot 28 фев 2018, 14:29, всего редактировалось 2 раза.
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________
-
- Гуру
- Сообщения: 534
- Зарегистрирован: 30 ноя 2006, 13:31
- Репутация: 116
- Откуда: Moscow
Re: Наложение растра
можно по аналогии с этой статьей
http://gis-lab.info/qa/gdal-hillshade-colormap.html
можно Qgis
открыть оба продукта, выставить визуализацию - Одноканальное псевдоцветное - задать нужные палитры и диапазоны значений, выставить прозрачность.
Получится нечто вроде этого
http://gis-lab.info/qa/gdal-hillshade-colormap.html
можно Qgis
открыть оба продукта, выставить визуализацию - Одноканальное псевдоцветное - задать нужные палитры и диапазоны значений, выставить прозрачность.
Получится нечто вроде этого
- rhot
- Гуру
- Сообщения: 1727
- Зарегистрирован: 25 янв 2011, 17:50
- Репутация: 194
- Ваше звание: доктор
- Откуда: Архангельск
Re: Наложение растра
Добавил решение в R.
___________(¯`·.¸(¯`·.¸ Scientia potentia est _/ {SILVA}:::{FOSS}:::{GIS} \_ Знание сила ¸.·´¯)¸.·´¯)___________
Кто сейчас на конференции
Сейчас этот форум просматривают: Ahrefs [Bot] и 10 гостей