Время ускорения для построения пространственно-полигональных данных очень длинное в R
Я делаю заговор в R некоторых SpatialPolygonDataFrames, но это занимает вечность. Мне понадобилось три часа, чтобы сохранить в png график одного из этих слоев, который имеет 314 элементов и весов (как сказал RStudio) 35 Мб. Я не знаю, много это или нет, но для его отображения в программном обеспечении ГИС, таком как QGIS или ArcMap, требуется не более доли секунды.
Это имеет настолько малый смысл, что мне интересно, если это так, как это должно быть, или что-то не так в моей настройке. Это код, который я использую, где APpat это шейп-файл защищенных областей
> library(rgdal)
> library(rasterVis)
> APpat <- readOGR("./PAT", "AP_PAT_edited1")
> PAT.alt <- raster("./altitude") # I plot the polygons over this raster, which is very light to plot on its own
# This is how APpat looks like
> APpat
class : SpatialPolygonsDataFrame
features : 314
extent : -92.58878, -57.4542, -22.9007, 13.39352 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0
variables : 6
names : id, category, name, year, country, IUCN
min values : 1, Area Nacional de Recreacion, Abra de Rio Frio, 1936, Bolivia, I
max values : SN09, Via Parque, Yurubi, 2015, Venezuela, VI
# And then save the plot to png
> png(paste0("./ProtectedAreas.png"), width=1300, heigh=1300, res=300)
> levelplot(PAT.alt, margin=FALSE,
at=seq(0,6800, by=2000),
main=paste0("Protected areas, ", years[y]),
par.settings=mytheme,
xlab=NULL, ylab=NULL,
contour = TRUE, col='dark green',
colorkey=NULL) +
> layer(sp.polygons(myAP, lwd=0.8, col='black', fill="#fc9272"))
> par(cex=.7)
> dev.off()
На сюжет ушло 3 часа 7 минут. Я попытался использовать нормальный график вместо уровня, с эквивалентными результатами. То же самое, если я строю не растр, а полигоны самостоятельно.
> sessionInfo()
R version 3.3.2 (2016-10-31)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: macOS Sierra 10.12.4
Я установил XQuartz. Я нахожу похожие времена, когда просто строю график в RStudio и сохраняю график в файл PDF. Спасибо за вашу помощь!
РЕДАКТИРОВАТЬ: это слой шейп-файла, который я пытаюсь построить: https://drive.google.com/open?id=0B0XRH7FN95-RSzFRckhXNjVWYlk
РЕДАКТИРОВАТЬ 2: Это не обобщаемая проблема, кажется, происходит только со мной. Мы смогли нанести рисунок на разные машины за ~20 секунд. Тем не менее, если у кого-то есть идея, что может быть не так с моим компьютером, это будет высоко ценится!