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

ฉันจำเป็นต้องเปิด shapefile จาก ArcMap ใน R เพื่อใช้สำหรับการวิเคราะห์ทางภูมิศาสตร์เพิ่มเติม ฉันแปลงมันเป็นไฟล์ข้อความ ASCII แต่ใน R มันได้รับการยอมรับว่าเป็น data.frame ฟังก์ชันพิกัดไม่ทำงานทันทีที่ x และ y ได้รับการยอมรับว่าไม่ใช่ตัวเลข คุณช่วยจัดการกับมันได้ไหม?

64 shapefile  r 

Rกำลังกลายเป็นเครื่องมือที่แข็งแกร่งสำหรับการจัดการและวิเคราะห์ข้อมูลเชิงพื้นที่ ฉันเรียนรู้สิ่งที่มีประโยชน์ผ่านคำถามเช่น นี้ที่ SO และคิดว่าอาจมีประโยชน์ในการมีบางสิ่งที่คล้ายคลึงกัน แต่มีการ "เชิงพื้นที่" มากกว่า คุณสามารถแบ่งปันเคล็ดลับและเทคนิค R เชิงพื้นที่ที่คุณเห็นว่ามีประโยชน์หรือไม่

44 tips-and-tricks  r 

ไม่มีใครรู้บทเรียนที่ดีในการจัดการข้อมูลเชิงพื้นที่ใน R หรือไม่? ฉันรู้ว่ามีแพคเกจบางอย่างเช่น rgdal และ maptools แต่ฉันไม่ได้พบบทเรียนที่มีประโยชน์เพื่อติดต่อกับฟังก์ชั่นพื้นฐาน (อ่าน, พล็อต, เปลี่ยนสัญลักษณ์, เข้าร่วมข้อมูล ฯลฯ ) และ / หรือเพื่อดูภาพรวมที่กว้างขึ้น อะไรคือความสามารถแบบเฉียบพลันของ R สำหรับการวิเคราะห์เชิงพื้นที่และการจัดการข้อมูลเชิงพื้นที่ ยินดีต้อนรับสู่ recomendations สำหรับหนังสือเอกสารบล็อก ฯลฯ

43 r  spatial-statistics  references 

มี sf-native (เช่น "ถูกต้อง") วิธีการแยกทุกอย่างยกเว้นคอลัมน์รูปทรงเรขาคณิตจากวัตถุคุณสมบัติที่เรียบง่าย? วิธีนี้ใช้ได้ผล df <- dplyr::select(as.data.frame(sf), -geometry) แต่select( , -geometry)ขั้นตอนนั้นไม่จำเป็น นอกจากนี้ยังไม่ได้ลบคุณลักษณะทางเรขาคณิต

43 r  sf 

มีความเป็นไปได้ค่อนข้างดีในการวิเคราะห์ข้อมูลเชิงพื้นที่ใน R และในบริบทของโครงการปัจจุบันของฉันฉันต้องการใช้ R บ่อยขึ้นในการทำเช่นนั้น จนถึงตอนนี้ฉันพล็อตแผนที่ของฉันด้วยแพ็คเกจ ggplot2 ที่นำเครื่องมือที่ใช้งานได้มากมายมาพล็อตและสำรวจข้อมูล แผนที่ของฉันยังดูไม่ค่อยดีเท่าแผนที่ที่ใช้ใน ArcGIS ดังนั้นเพื่อวัตถุประสงค์ในการเผยแพร่ฉันสงสัยว่ามีแบบฝึกหัดที่ดีหนังสือเคล็ดลับการปฏิบัติแพคเกจอื่น ๆ ที่ฉันสามารถใช้ในการทดสอบเล็กน้อยและทำให้แผนที่ของฉันดูดีขึ้นหรือไม่

41 r  cartography 

ฉันใช้เวลาเล็กน้อยในการหาคำตอบสำหรับคำถามนี้ ไม่ชัดเจนในทันทีจากการค้นหาของ Googleดังนั้นคิดว่าอาจมีประโยชน์ในการโพสต์คำตอบที่นี่ นอกจากนี้ยังมีคำถามเพิ่มเติมเกี่ยวกับรูปหลายเหลี่ยมไม่ต่อเนื่องกัน คำตอบง่าย ๆ ทันที: ใช้คำสั่ง: centroids <- getSpPPolygonsLabptSlots(polys) (ซึ่งพบได้ในคำอธิบายคลาสของคลาสข้อมูล SpatialPolygonsDataFrame R สำหรับแพ็คเกจอวกาศที่ครอบคลุมใน R, sp ) ดูเหมือนว่าจะทำสิ่งเดียวกัน cents <- SpatialPointsDataFrame(coords=cents, data=sids@data, proj4string=CRS("+proj=longlat +ellps=clrk66")) ในรหัสต่อไปนี้ซึ่งควรทำซ้ำในการติดตั้ง R ใด ๆ (ลองใช้!) #Rcentroids install.packages("GISTools") library(GISTools) sids <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1], proj4string=CRS("+proj=longlat +ellps=clrk66")) class(sids) plot(sids) writeSpatialShape(sids, "sids") cents <- coordinates(sids) cents <- SpatialPointsDataFrame(coords=cents, data=sids@data, …

41 r  centroids 

ฉันทำงานกับไฟล์. kml ขนาดใหญ่ (มากถึง 10 Gb) และต้องการวิธีที่มีประสิทธิภาพในการอ่านไฟล์ให้เป็นไฟล์ R จนกระทั่งตอนนี้ฉันได้แปลงไฟล์เหล่านี้ให้เป็นไฟล์รูปร่างผ่าน QGIS แล้วกลับไปเป็น R ด้วย readShapePoly และ readOGR โดยวิธีการคือ 1,000 เร็วกว่าเดิม) ฉันชอบที่จะลดขั้นตอนการเป็นตัวกลาง QGIS เนื่องจากยุ่งยากและช้าลง จะอ่านไฟล์. kml โดยตรงได้อย่างไร ฉันเห็นนี้ยังสามารถทำได้ด้วย readOGR น่าเสียดายที่ผมไม่สามารถดูวิธีการใช้ตัวอย่างที่ทำงาน (หลังจากการเตรียมความยาวของไฟล์ .kml: xx <- readOGR(paste(td, "cities.kml", sep="/"), "cities")) ดูเหมือนว่า "เมือง" ที่นี่เป็นชื่อของวัตถุเชิงพื้นที่ Roger Bivand ยอมรับว่า "วิธีค้นพบชื่อนี้ไม่ชัดเจนเนื่องจากไดรเวอร์ KML ใน OGR ต้องการให้เข้าถึงไฟล์ได้สิ่งหนึ่งที่เป็นไปได้คือ: system(paste("ogrinfo", paste(td, "cities.kml", …

41 kml  r 

ฉันมีชุดข้อมูลของค่าเหนือกริดกม. ในทวีปอเมริกาคอลัมน์คือ "ละติจูด", "ลองจิจูด" และ "การสังเกต" เช่น: "lat" "lon" "yield" 25.567 -120.347 3.6 25.832 -120.400 2.6 26.097 -120.454 3.4 26.363 -120.508 3.1 26.630 -120.562 4.4 หรือในฐานะกรอบข้อมูล R: mydata <- structure(list(lat = c(25.567, 25.832, 26.097, 26.363, 26.63), lon = c(-120.347, -120.4, -120.454, -120.508, -120.562), yield = c(3.6, 2.6, 3.4, 3.1, 4.4)), …

40 raster  convert  r 

ฉันใช้PL/Rฟังก์ชั่นและPostGISสร้างรูปหลายเหลี่ยม voronoi รอบจุดต่างๆ ฟังก์ชั่นที่ผมใช้ถูกกำหนดไว้ที่นี่ เมื่อฉันใช้ฟังก์ชั่นนี้บนชุดข้อมูลเฉพาะฉันได้รับข้อความแสดงข้อผิดพลาดต่อไปนี้: Error : ERROR: R interpreter expression evaluation error DETAIL: Error in pg.spi.exec(sprintf("SELECT %3$s AS id, st_intersection('SRID='||st_srid(%2$s)||';%4$s'::text,'%5$s') AS polygon FROM %1$s WHERE st_intersects(%2$s::text,'SRID='||st_srid(%2$s)||';%4$s');", :error in SQL statement : Error performing intersection: TopologyException: found non-noded intersection between LINESTRING (571304 310990, 568465 264611) and LINESTRING (568465 264611, 594406 …

38 postgis  polygon  r  voronoi-thiessen 

ฉันมีรูปร่างหลายแบบใน CRS ที่แตกต่างกัน (ส่วนใหญ่เป็น WGS84 lat / lon) ที่ฉันต้องการแปลงเป็นรูปแบบทั่วไป (น่าจะเป็น Albers Equal Area Conic แต่ฉันอาจขอความช่วยเหลือในการเลือกคำถามอื่นเมื่อปัญหาของฉันดีขึ้น -defined) ฉันใช้เวลาสองสามเดือนในการทำสถิติเชิงพื้นที่ใน R แต่เมื่อ 5 ปีที่แล้ว สำหรับชีวิตของฉันฉันจำไม่ได้ว่าจะเปลี่ยนspวัตถุอย่างไร (เช่นSpatialPolygonsDataFrame) จากการฉายภาพหนึ่งไปอีกภาพหนึ่ง รหัสตัวอย่าง: P4S.latlon <- CRS("+proj=longlat +datum=WGS84") hrr.shp <- readShapePoly("HRR_Bdry"), verbose=TRUE, proj4string=P4S.latlon) # Shapefile available at # http://www.dartmouthatlas.org/downloads/geography/hrr_bdry.zip # but you must rename all the filenames to have …

35 r  coordinate-system 

ฉันมีจุดข้อมูลมากมายพร้อมละติจูดและลองจิจูด ฉันต้องการใช้ R เพื่อจัดกลุ่มตามระยะทาง ฉันลองดูที่หน้านี้แล้วลองใช้แพ็คเกจclustTool แต่ฉันไม่แน่ใจว่าฟังก์ชัน clust ใน clustTool พิจารณาจุดข้อมูล (lat, lon) เป็นข้อมูลเชิงพื้นที่และใช้สูตรที่เหมาะสมเพื่อคำนวณระยะทางระหว่างพวกเขาหรือไม่ ฉันหมายถึงฉันไม่สามารถเห็นความแตกต่างระหว่างข้อมูลเชิงพื้นที่กับข้อมูลลำดับได้ ฉันเชื่อว่าการคำนวณระยะทางระหว่างจุดสองจุดบนแผนที่ (พื้นที่) และตัวเลขสองปกติจะแตกต่างกัน (ไม่ใช่เหรอ?) จะเกิดอะไรขึ้นถ้าฉันต้องการพิจารณาพารามิเตอร์ที่สามในการจัดกลุ่มของฉัน เช่นถ้าฉันมี (lat, lon) และพารามิเตอร์อื่น ๆ ระยะทางคำนวณอย่างไร ปัญหาอื่น ๆ ที่ฉันมีกับ clustTool คือมันถูกออกแบบโดยคำนึงถึง GUI ฉันไม่รู้ว่าจะข้ามค่าใช้จ่าย GUI ในห้องสมุดไปได้อย่างไรเพราะฉันไม่ต้องการ ฉันมีตัวเลือกอะไรใน R สำหรับการวิเคราะห์กลุ่มของข้อมูลเชิงพื้นที่?

34 r  distance  software-recommendations  spatial-statistics  clustering 

ฉันกำลังสร้างแผนที่สำหรับภาคตะวันออกเฉียงเหนือของสหรัฐอเมริกาพื้นหลังของแผนที่จำเป็นต้องเป็นแผนที่ระดับความสูงหรือแผนที่อุณหภูมิเฉลี่ยรายปี ฉันมีแรสเตอร์สองตัวจาก Worldclim.org ซึ่งให้ตัวแปรเหล่านี้แก่ฉัน แต่ฉันต้องคลิปพวกมันให้มากเท่าที่ฉันสนใจข้อเสนอแนะใด ๆ เกี่ยวกับวิธีการทำเช่นนี้ นี่คือสิ่งที่ฉันมี: #load libraries library (sp) library (rgdal) library (raster) library (maps) library (mapproj) #load data state<- data (stateMapEnv) elevation<-raster("alt.bil") meantemp<-raster ("bio_1.asc") #build the raw map nestates<- c("maine", "vermont", "massachusetts", "new hampshire" ,"connecticut", "rhode island","new york","pennsylvania", "new jersey", "maryland", "delaware", "virginia", "west virginia") map(database="state", …

32 raster  r  clip  maps 

วิธีหนึ่งสามารถใช้ Rเพื่อ แยก Shapefileใน 200 เมตรสี่เหลี่ยม / รูปหลายเหลี่ยมย่อย พล็อตกริดนี้ (รวมถึงหมายเลข ID สำหรับแต่ละช่อง) เหนือแผนที่ดั้งเดิมด้านล่างและ ประเมินว่าพิกัดทางภูมิศาสตร์เฉพาะของรูปสี่เหลี่ยมจัตุรัสตั้งอยู่ที่ไหน ฉันเป็นผู้เริ่มต้นใน GIS และนี่อาจเป็นคำถามพื้นฐาน แต่ฉันไม่พบบทเรียนเกี่ยวกับวิธีการทำเช่นนี้ใน R สิ่งที่ฉันทำไปแล้วคือการโหลดรูปร่างไฟล์ของ NYC และวางแผนพิกัดทางภูมิศาสตร์ที่เป็นแบบอย่าง ฉันกำลังมองหาตัวอย่าง (รหัส R) วิธีการนี้ด้วยข้อมูลด้านล่าง # Load packages library(maptools) # Download shapefile for NYC # OLD URL (no longer working) # shpurl <- "http://www.nyc.gov/html/dcp/download/bytes/nybb_13a.zip" shpurl <- "https://www1.nyc.gov/assets/planning/download/zip/data-maps/open-data/nybb_13a.zip" tmp <- …

32 r  vector-grid  point-in-polygon 

ฉันอยู่หลังคำอธิบาย / ระเบียบวิธีสำหรับการดำเนินการวิเคราะห์องค์ประกอบหลักทางภูมิศาสตร์ถ่วงน้ำหนัก(GWPCA) ฉันมีความสุขที่ใช้ Python สำหรับส่วนใดส่วนหนึ่งของนี้และฉันคิดว่า SPSS หรือ R ถูกใช้เพื่อเรียกใช้ PCA กับตัวแปรถ่วงน้ำหนักทางภูมิศาสตร์ ชุดข้อมูลของฉันประกอบด้วยตัวแปรอิสระ 30 ตัวที่วัดได้ทั่วระบบสำรวจสำมะโนประชากร ~ 550 (เรขาคณิตเวกเตอร์) ฉันรู้ว่านี่เป็นคำถามที่โหลด แต่เมื่อฉันค้นหาและค้นหาดูเหมือนว่าจะไม่มีทางออกใด ๆ สิ่งที่ฉันเจอคือสมการทางคณิตศาสตร์ที่อธิบายองค์ประกอบพื้นฐานของ GWPCA (และ GWR) สิ่งที่ฉันเป็นหลังจากนั้นถูกนำไปใช้มากขึ้นในแง่ที่ว่าฉันกำลังมองหาสิ่งที่ขั้นตอนสำคัญที่ฉันต้องประสบความสำเร็จในการรับจากข้อมูลดิบกับผลลัพธ์ GWPCA ฉันต้องการขยายในส่วนแรกของการแก้ไขนี้เนื่องจากความคิดเห็นที่ได้รับด้านล่าง เพื่อที่อยู่พอล ... ฉันอ้างอิงความสนใจของฉันใน GWPCA จากเอกสารต่อไปนี้: Lloyd, CD, (2010) การวิเคราะห์ลักษณะประชากรโดยใช้การวิเคราะห์องค์ประกอบหลักทางภูมิศาสตร์: กรณีศึกษาไอร์แลนด์เหนือในปี 2544 คอมพิวเตอร์สิ่งแวดล้อมและระบบเมือง 34 (5), หน้า 339-399 สำหรับผู้ที่ไม่สามารถเข้าถึงวรรณกรรมฉันได้แนบภาพหน้าจอของส่วนเฉพาะที่อธิบายคณิตศาสตร์ด้านล่าง: และเพื่อที่อยู่ whuber ... โดยไม่ต้องลงรายละเอียด (การรักษาความลับ) …

32 r  spatial-statistics 

นี่คือการติดตามเพื่อโพสต์ก่อนหน้านี้: เครื่องเรียนรู้ขั้นตอนวิธีสำหรับการจำแนกประเภทที่ดินปก ดูเหมือนว่าวิธีการจำแนกแบบสุ่มฟอเรสต์ (RF) กำลังได้รับแรงผลักดันมากมายในโลกการรับรู้จากระยะไกล ฉันสนใจ RF โดยเฉพาะเนื่องจากจุดแข็งหลายประการ: วิธีการแบบไม่มีพารามิเตอร์ที่เหมาะสมกับข้อมูลการสำรวจระยะไกล ความถูกต้องของการจำแนกประเภทที่รายงานสูง มีการรายงานความสำคัญของตัวแปร ด้วยจุดแข็งเหล่านี้ฉันต้องการทำการจำแนกดินแดนป่าไม้แบบสุ่มโดยใช้ภาพความละเอียดสูง 4 แบนด์ มีวัสดุจำนวนมากและงานวิจัยโน้มน้าวข้อได้เปรียบของป่าสุ่ม แต่มีข้อมูลน้อยมากเกี่ยวกับวิธีการวิเคราะห์การจำแนกประเภท ฉันคุ้นเคยกับการถดถอย RF โดยใช้ R และต้องการใช้สภาพแวดล้อมนี้เพื่อเรียกใช้อัลกอริทึมการจำแนก RF ฉันจะรวบรวมประมวลผลและป้อนข้อมูลการฝึกอบรม (เช่นจากภาพถ่ายทางอากาศ CIR ความละเอียดสูง) ลงในอัลกอริทึมRandom Forestโดยใช้ R ได้อย่างไร คำแนะนำขั้นตอนที่ชาญฉลาดเกี่ยวกับวิธีการผลิตแรสเตอร์ครอบคลุมพื้นที่ประเภทจะได้รับการชื่นชมอย่างมาก

32 r  land-classification  random-forest  machine-learning