Извлечь «цвет» из RasterStack по точке в R, игнорируя черный текст - PullRequest
0 голосов
/ 17 июня 2020

Я хотел бы извлечь видимый цвет из этой карты торгового потока из Танзании здесь для каждой точки в 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 и указать один цвет. Во-вторых, я хочу сначала сгладить цвет текста и небольших линий дороги, потому что некоторые точки будут попадать на текст или дороги, но я хочу знать общий цвет под текстом. Я тоже не знаю, как это сделать.

...