чтение и графическое изображение esri файла форм в R

Я испытываю затруднения при чтении в .shp (esri файл форм) в R. Я попробовал несколько опций в R и попытался преобразовать файл форм в ArcMap к чему-то, что правильно читает в файле форм, но ничто еще не работало. (В ArcMap я исправил геометрию, преобразованную от единственного до мультиполигона, и т.д. который был, вероятно, не необходим или релевантен),

Это, вероятно, имеет что-то к с тем, что мой файл форм содержит 'регионы' (мультиполигоны) вместо 'полигонов'...

Как я могу считать тот тип файла форм правильно в R для графического изображения? (это похоже на нормальную форму в ArcMap),

В ArcMap файл форм похож на это: http://bfast.r-forge.r-project.org/arcmapshape.jpg

(показывает файл форм с полигонами в других полигонах),

В R это похоже на это: enter image description here

(показывает файл форм, где некоторые полигоны неправильно заполнены),

Я использовал следующий код в 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

Как это может быть решено?

15
задан iPhoneProcessor 29 April 2013 в 18:44
поделиться

1 ответ

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.

11
ответ дан 1 December 2019 в 04:53
поделиться
Другие вопросы по тегам:

Похожие вопросы: