คำถามติดแท็ก rgeos

ข้อผิดพลาด 'TopologyException: Input geom 1 ไม่ถูกต้อง' ซึ่งเกิดขึ้นจากรูปหลายเหลี่ยมที่ไม่ถูกต้องซึ่งมีการกล่าวถึงกันอย่างกว้างขวาง อย่างไรก็ตามฉันไม่พบวิธีแก้ไขปัญหาที่สะดวกบนเว็บที่อาศัยการทำงานของ R เพียงอย่างเดียว ยกตัวอย่างเช่นฉันมีการจัดการเพื่อสร้างวัตถุ 'SpatialPolygons จากการส่งออกของmap("state", ...)ต่อไปนี้เป็นคำตอบที่ดี Josh โอไบรอันที่นี่ 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 …

24 r  polygon  rgeos  self-intersection 

ฉันมีรายชื่อของบัฟเฟอร์เชิงพื้นที่ (30,000 บัฟเฟอร์) ที่ฉันสร้างขึ้นด้วยฟังก์ชั่น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]] …

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 …

10 r  ggplot2  rgeos