タグ付けされた質問 「rgeos」

無効なポリゴンジオメトリから発生する「TopologyException：Input geom 1 is invalid」自己交差エラーは広く議論されています。しかし、R機能だけに依存する便利な解決策をWeb上で見つけられませんでした。 たとえば、私はここでmap("state", ...) Josh O'Brienの素敵な答えを次の出力から「SpatialPolygons」オブジェクトを作成することができました。 library(maps) library(maptools) map_states = map("state", fill = TRUE, plot = FALSE) IDs = sapply(strsplit(map_states$names, ":"), "[[", 1) spydf_states = map2SpatialPolygons(map_states, IDs = IDs, proj4string = CRS("+init=epsg:4326")) plot(spydf_states) この広く適用されているデータセットの問題は、以下のポイントで自己交差が発生することです。 rgeos::gIsValid(spydf_states) [1] FALSE Warning message: In RGEOSUnaryPredFunc(spgeom, byid, "rgeos_isvalid") : Self-intersection at or …

24 r  polygon  rgeos  self-intersection 

関数で構築した空間バッファー（30000バッファー）のリストがありますlapply。 buff.pts <- lapply(1:nrow(pts.prj), FUN=function(l){ buff <- gBuffer(pts.prj[l,], width=1000) ## 1km return(buff) })) > head(buff.pts) [[1]] class : SpatialPolygons features : 1 extent : 307941.8, 311941.8, 4994518, 4998518 (xmin, xmax, ymin, ymax) coord. ref. : +proj=tmerc +lat_0=0 +lon_0=-73.5 +k=0.9999 +x_0=304800 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs +towgs84=0,0,0 [[2]] class : SpatialPolygons …

16 polygon  r  sp  rgeos 

ggplot2を使用して空間データをプロットするのに問題があります。spplotを使用してプロットするとマップが正常に見えるため、要塞化段階でティアリングが発生すると想定しています。 コードは次のとおりです。 #install the packages library(rgdal) library(mapproj) library(raster) library(rgeos) library(ggplot2) library(plyr) if (!require(gpclib)) install.packages("gpclib", type="source") gpclibPermit() setwd("C:/Users/My Documents") #read in laa to regional mapping #must aggregate to higher level regions as data is provided at this higher level laa_region_mapping <- read.csv("laa_region.csv", header = TRUE) #read in LAA polygons laa_polygons <- …

10 r  ggplot2  rgeos 

で測地線バッファリングの理解、Esriのジオプロセシング開発チームは、ユークリッドと測地線バッファリングを区別します。「投影されたフィーチャクラスで実行されるユークリッドバッファリングは、誤解を招く技術的に不正確なバッファを生成する可能性があります。ただし、測地線バッファは、投影座標系によって導入される歪みの影響を受けないため、常に地理的に正確な結果を生成します。」 ポイントグローバルデータセットを使用する必要があり、座標は投影されません（+proj=longlat +ellps=WGS84 +datum=WGS84）。メートル単位で幅が指定されている場合、Rに測地線バッファーを作成する関数はありますか？パッケージgBufferから承知しておりrgeosます。この関数は、使用される空間オブジェクト（例）の単位でバッファーを作成するため、座標を投影して、目的のX kmのバッファーを作成できるようにする必要があります。投影してから、gBuffer実際にユークリッドバッファーを作成する手段を適用します。以下は、私の懸念を示すためのコードです。 require(rgeos) require(sp) require(plotKML) # Generate a random grid-points for a (almost) global bounding box b.box <- as(raster::extent(120, -120, -60, 60), "SpatialPolygons") proj4string(b.box) <- "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs" set.seed(2017) pts <- sp::spsample(b.box, n=100, type="regular") plot(pts@coords) # Project to Mollweide to be able to apply buffer with …

9 coordinate-system  r  buffer  rgeos  geodesic