lire et tracer un fichier de forme d'esri en R

Je vais avoir des difficultés de lecture dans une .shp (esri shape file) dans R. j'ai essayé plusieurs options dans R, et a essayé de convertir le fichier de forme dans ArcMap pour quelque chose qui lit correctement dans le fichier de formes, mais rien n'a fonctionné encore. (Dans ArcMap, j'ai corrigé la géométrie, converti à partir d'unique à multipolygon, etc, qui n'était probablement pas nécessaire ou pertinent)

Il a probablement quelque chose avec le fait que ma forme de fichier contient des "régions" (multi-polygones) au lieu de " polygones...

Comment puis-je lire ce type de fichier de forme correctement dans la R pour le dessin? (il ressemble à une forme normale dans ArcMap)

Dans ArcMap la forme de fichier ressemble à ceci:
http://bfast.r-forge.r-project.org/arcmapshape.jpg

(montre un fichier de forme avec les polygones dans d'autres polygones)

Dans R, il ressemble à ceci:
lire et tracer un fichier de forme d'esri en R

(montre un fichier de forme où certains polygones sont mal remplis)

J'ai utilisé le code suivant dans 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)

Les deux plot() commandes de donner le même résultat R. Aucune erreur ne se produit. seulement le message suivant

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

Comment cela peut-il être résolu?

source d'informationauteur Janvb

  1. 11

    Jan,

    Votre fichier a probablement été lu correctement. Il est plus probable que le tracé n'était pas correct. Essayez spplot() au lieu de plot() :

    spplot(shp, col.regions="gray", col="blue")

    Si votre fichier a plus d'une colonne (ce qui est probable), ajouter zcol= 1 (ou n'importe quelle autre colonne), de sorte que vous obtenez un seul panneau

    spplot(shp, zcol=1, col.regions="gray", col="blue")

    Si vous voulez vraiment utiliser de la parcelle, vous devez définir une couleur pour l'arrière-plan, parce que sinon il va imprimer les trous transparent. Donc, mettre pbg="blanc". Essayez:

    plot(shp, col="gray", border="blue", axes=TRUE, pbg="white")

    Par la, la liste je pense que Dirk a été fait allusion à R-sig-Geo et vous pouvez y trouver beaucoup de réponses et de poser beaucoup de sig-questions connexes.