У меня есть простой растр, который проходит через функцию fillContour в растровом пакете, чтобы получить заполненный контурный контур.
Я пытаюсь наложить контур многоугольника поверх поверх этого заполненного контура plot, но возникают проблемы.
чтение заметок в fill.contour - добавление линий к закрашенному контурному графику должно быть сделано в разделе plot.axis, но безуспешно.
Итак пока у меня есть что-то вроде этого:
library(raster)
library(rgdal)
test_raster <- raster("test_raster.grd")
test_Shapefile <- readOGR("test_shapefile.shp")
test_filled_contour <- filledContour(test, maxpixels=100000, plot.axes = {
axis(1)
axis(2)
map(test_Shapefile, add=T)
})
И растр, и шейп-файл используют одну и ту же систему координат.
Любая помощь приветствуется, я могу предоставить данные, если это необходимо.