Я испытываю затруднения при чтении в .shp (esri файл форм) в R. Я попробовал несколько опций в R и попытался преобразовать файл форм в ArcMap к чему-то, что правильно читает в файле форм, но ничто еще не работало. (В ArcMap я исправил геометрию, преобразованную от единственного до мультиполигона, и т.д. который был, вероятно, не необходим или релевантен),
Это, вероятно, имеет что-то к с тем, что мой файл форм содержит 'регионы' (мультиполигоны) вместо 'полигонов'...
Как я могу считать тот тип файла форм правильно в R для графического изображения? (это похоже на нормальную форму в ArcMap),
В ArcMap файл форм похож на это: http://bfast.r-forge.r-project.org/arcmapshape.jpg
(показывает файл форм с полигонами в других полигонах),
В R это похоже на это:
(показывает файл форм, где некоторые полигоны неправильно заполнены),
Я использовал следующий код в R:
require(maptools)
require(rgdal)
newproj <- "+proj=utm +zone=55 +south +ellps=GRS80 +units=m"
shape<- readShapeSpatial(pdir, proj4string = CRS(newproj),repair=TRUE,force_ring=T,verbose=TRUE) # without any errors
plot(shape, col="gray",border="blue", axes=TRUE)
# via rgdal
folder <- c("spatial")
lyr <- c("clipped_forest_mga")
shp <- readOGR(dsn=folder,layer=lyr)
plot(shp, col="gray",border="blue", axes=TRUE)
Оба графика () команды дают тот же результат R. Никакие ошибки не происходят. только следующее сообщение
OGR data source with driver: ESRI Shapefile
Source: "P:/Victoria_DSE/BFAST_spatial/vector/PLM_excl_fire03_09_GDA94LL/mgaz94z55/clipped_EG", layer: "clipped_forest_mga"
with 1 features and 4 fields
Feature type: wkbMultiPolygon with 2 dimensions
Как это может быть решено?
http://msdn.microsoft.com/en-us/library/system.windows.forms.application.restart.aspx
-121--3329824-Приложения .Restart () - ваш метод:)
Здесь - еще один ответ StackOverflow, который указывает на пару «watch-out-for» с использованием этого метода.
-121--3329823-Ян,
Вероятно, файл формы прочитан правильно. Более вероятно, что построение графика было неверным. Попробуйте spplot () вместо plot ():
spplot(shp, col.regions="gray", col="blue")
Если файл формы имеет более одного столбца (что вероятно), добавьте zcol = 1 (или любой другой столбец), чтобы получить только одну панель
spplot(shp, zcol=1, col.regions="gray", col="blue")
Если вы действительно хотите использовать plot, необходимо задать цвет для фона, так как в противном случае отверстия будут напечатаны прозрачными. Так что установите pbg = «белый». Попробуйте:
plot(shp, col="gray", border="blue", axes=TRUE, pbg="white")
Кстати, список, о котором я думаю, что Дирк говорил, это R-sig-Geo , и вы можете найти там много ответов и задать много вопросов, связанных с sig.