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

ฉันมีคุณสมบัติที่มีอยู่ในฐานข้อมูลทางภูมิศาสตร์ที่มีขนาดใหญ่กว่า 2GB เป็นรูปร่างไฟล์ที่ส่งออก ฉันต้องการเรียกใช้ฟังก์ชันแยกใน R เพื่อแอตทริบิวต์รูปหลายเหลี่ยมด้วยข้อมูลจากไฟล์แรสเตอร์ การส่งออกคุณสมบัติเป็นตารางไม่ใช่วิธีแก้ปัญหา ฉันจะอ่านคลาสคุณลักษณะที่อยู่ในฐานข้อมูลไฟล์ Esri ได้อย่างไร

27 r  features  file-geodatabase  feature-class  rgdal 

ฉันพยายามที่จะอ่านตารางโดยตรงจากแฟ้มของ Geodatabase ESRI เข้าไปในอาร์ไฟล์ตัวอย่างเช่นข้อมูลที่สามารถดาวน์โหลดได้ที่นี่ ฐานข้อมูลมีคลาสคุณลักษณะจุด (Zone9_2014_01_Broadcast) และตารางที่เชื่อมโยงสองตาราง (Zone9_2014_01_Vessel และ Zone9_2014_01_Voyage) คุณสามารถอ่าน shapefile ใน R โดยใช้readOGRจากrgeosแพ็คเกจ: library(rgeos) library(downloader) download("https://coast.noaa.gov/htdata/CMSP/AISDataHandler/2014/01/Zone9_2014_01.zip", dest="Zone9_2014_01.zip", mode="wb") unzip("Zone9_2014_01.zip", exdir = ".") # Not Run (loads large point file) # broadcast <- readOGR(dsn = "Zone9_2014_01.gdb", layer = "Zone9_2014_01_Broadcast") ทั้งสองตารางที่เชื่อมโยงยังแสดงให้เห็นเมื่อคุณใช้หรือogrListLayers ogrInfoอย่างไรก็ตามogrInfoให้คำเตือน: ข้อความเตือน: ใน ogrInfo ("Zone9_2014_01.gdb", layer = "Zone9_2014_01_Vessel"): ogrInfo: …

21 r  file-geodatabase  table  rgdal 

Zillow มีชุดรูปร่างที่แตกต่างกันสำหรับละแวกใกล้เคียงของเมืองสำคัญ ๆ ในสหรัฐอเมริกา ฉันต้องการตรวจสอบว่ามีอาคารบางแห่งอยู่ในละแวกใกล้เคียงโดยใช้ R: library(rgeos) library(sp) library(rgdal) df <- data.frame(Latitude =c(47.591351, 47.62212,47.595152), Longitude = c(-122.332271,-122.353985,-122.331639), names = c("Safeco Field", "Key Arena", "Century Link")) coordinates(df) <- ~ Latitude + Longitude wa.map <- readOGR("ZillowNeighborhoods-WA.shp", layer="ZillowNeighborhoods-WA") sodo <- wa.map[wa.map$CITY == "Seattle" & wa.map$NAME == "Industrial District", ] ฉันสามารถลงจุดได้โดยไม่มีปัญหาใด ๆ plot(sodo) points(df$Latitude …

19 r  rgdal  zillow 

ฉันมีรูปหลายเหลี่ยมสองอัน หนึ่งประกอบด้วยเขตข้อมูล (X, Y, Z) และอื่น ๆ ประกอบด้วยชนิดของดิน (A, B, C, D) ฉันต้องการที่จะรู้ว่าสิ่งที่พื้นที่ของทุกเขตข้อมูลมีชนิดของดิน ฉันพยายามต่อไปนี้: library(rgdal) library(rgeos) Field<-readOGR("./","Field") Soil<-readOGR("./","Soil") Results<-gIntersects(Soil,Field,byid=TRUE) rownames(Results)<-Field@data$FieldName colnames(Results)<-Soil@data$SoilType > Results A B C D Z TRUE FALSE FALSE FALSE Y FALSE TRUE TRUE FALSE X TRUE TRUE TRUE TRUE และได้ผลลัพธ์ที่ดีโดยบอกว่าสนามใดมีดินประเภทใด อย่างไรก็ตามฉันจะหาพื้นที่แทนได้อย่างไร

19 shapefile  polygon  r  intersection  rgdal 

ฉันต้องการสร้างรูปหลายเหลี่ยมสองรูป หนึ่งในนั้นคือรูปสี่เหลี่ยมผืนผ้าของแรสเตอร์ของฉัน ฉันรู้ว่าฉันสามารถได้รับขอบเขตของการใช้ raster r <- raster("band5.tif") e <- extent(r) plot(e)ให้ฉันนี้ ฉันจะสร้างหนึ่งรูปหลายเหลี่ยมของขอบเขตของแรสเตอร์ที่แสดงด้านล่างได้อย่างไร นี่คือลักษณะแรสเตอร์ของฉัน:

13 raster  r  rgdal 

ฉันใช้gSimplify(แพ็คเกจ rgeos) เพื่อทำให้รูปทรงเรขาคณิตของ shapefile ง่ายขึ้น funcion ทำงานได้ดี แต่ตอนนี้ฉันไม่สามารถเขียนผลลัพธ์ในรูปแบบไฟล์ใหม่ได้ ฉันลองวิธีบางอย่าง: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') ฉันเข้าใจ obj ต้องเป็น SpatialPointsDataFrame, SpatialLinesDataFrame หรือ SpatialPolygonsDataFrame และด้วย: writePolyShape(simplyshape, file) ฉันเข้าใจ: ข้อผิดพลาด: คือ (x, "SpatialPolygonsDataFrame") ไม่ใช่ TRUE

12 shapefile  r  simplify  rgdal 

ฉันมี และNE lat lng SW lat lngเป้าหมายของฉันง่ายมากฉันใช้ RStudio และฉันต้องการสร้างกล่องกระโดดจากสองด้านบนและด้านล่างจากนั้นสร้างรูปหลายเหลี่ยมและฉันต้องการเพิ่มจุดสุ่มลงในรูปหลายเหลี่ยมสี่เหลี่ยมจัตุรัส ฉันมี NELat / Lng และ SW-Lat / Lng ในไฟล์ csv ที่ฉันสามารถนำเข้าสู่ data.frame ฉันแค่ต้องการคำแนะนำเล็กน้อยในการเริ่มต้นฉันยังใหม่กับภาษา R มันแตกต่างจากที่ฉันคุ้นเคย จนถึงตอนนี้ฉันมีสิ่งนี้: coords = cbind(78.46801, 19.53407) coordsmax = cbind(78.83157, 19.74557 ) sp = SpatialPoints(coordsmax) sp2 = SpatialPoints(coords) r1 = rbind(coords, coordsmax[1, ]) # join P1 = Polygon(r1) …

11 r  rgdal 

ฉันพยายามเขียนหลายเลเยอร์ไปยัง GeoPackage เดียวกันใน R แต่ฉันได้รับข้อผิดพลาดCreation of output file failed. ฉันพยายามค้นหาเอกสารการอ่านและเขียนไปยังไฟล์. gpkg ด้วย RGDAL เช่นเพื่อค้นหาว่าwriteOGR()รองรับหลายเลเยอร์จริงหรือไม่ เป็นไปได้ไหมถ้าเป็นเช่นนั้นจะทำอย่างไร? ตัวอย่างการทำงานที่น้อยที่สุด: library(sp) library(maptools) library(rgdal) data(wrld_simpl) norway <- wrld_simpl[wrld_simpl$NAME == "Norway", ] sweden <- wrld_simpl[wrld_simpl$NAME == "Sweden", ] file <- tempfile("scandinavia", fileext = c(".gpkg")) writeOGR(norway, dsn = file, layer = "norway", driver = "GPKG") writeOGR(sweden, dsn …

11 gdal  r  rgdal  geopackage  sf 

ฉันมีไฟล์รูปร่าง ESRI ที่มีขนาดใหญ่มาก (> 8000 polygons, .shp ไฟล์> 32MB) ซึ่งฉันต้องการอ่านเป็น R บน Windows 7 ฉันคุ้นเคยกับ rgdal และ readOGR () อย่างไรก็ตามด้วยขนาดของ shapefile การอ่าน shapefile ทั้งหมดในหน่วยความจำจึงเป็นกระบวนการที่ช้ามาก นอกจากนี้ฉันสนใจเพียงไม่กี่โพลิกอนคลัสเตอร์ (<100) จาก 8,000 ที่มีอยู่ในรูปร่างไฟล์ กระบวนการที่ฉันต้องทำคือทำซ้ำอย่างมากนั่นคือต้องโหลดรูปร่างไฟล์หลายครั้งดังนั้นจึงช้ามาก ดังนั้นฉันสงสัยว่ามันเป็นไปได้ที่จะอ่านเฉพาะในส่วนของ shapefile เพื่อเร่งกระบวนการ ดูเหมือนว่าไม่น่าจะโหลดข้อมูลจำนวนมหาศาลเช่นนั้นทิ้งส่วนใหญ่

10 shapefile  r  rgdal 

ฉันกำลังแปลงเวกเตอร์เป็นแรสเตอร์ใน R อย่างไรก็ตามกระบวนการนี้ยาวเกินไป มีความเป็นไปได้ไหมที่จะนำสคริปต์ไปใช้ในการประมวลผลแบบมัลติเธรดหรือ GPU เพื่อให้เร็วขึ้น? สคริปต์ของฉันเป็นเวกเตอร์ที่แรสเตอร์ r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do raster res(r.raster) = 10 #definindo o tamanho do pixel setor.r = rasterize(setor, r.raster, 'dens_imov') r.raster คลาส: RasterLayer ขนาด: 9636, 11476, 110582736 (nrow, ncol, ncell) ความละเอียด: 10, 10 (x, y) ขอบเขต: 505755, 620515, 8555432, …

9 raster  r  rasterization  rgdal  parallel