Я хотел бы извлечь видимый цвет из этой карты торгового потока из Танзании здесь для каждой точки в SpatialPointsDataFrame, который я показываю здесь , нанесенный на карту . Следующий код создал этот график:
trade_stack=stack("TzMap.tif")
trade.e <- setExtent(trade_stack, extent(tzborder))
trade.m <- mask(trade.e, tzborder)
plotRGB(trade.m)
plot(tzborder, add=TRUE)
plot(tzshp, add=TRUE,pch=18, cex=1)
Мне в основном нужен индикатор для каждой точки, указывающий, находится ли точка на зеленом, оранжевом или бежевом фоне. Несколько проблем. Во-первых, я не знаю, как извлечь из RasterStack и указать один цвет. Во-вторых, я хочу сначала сгладить цвет текста и небольших линий дороги, потому что некоторые точки будут попадать на текст или дороги, но я хочу знать общий цвет под текстом. Я тоже не знаю, как это сделать.