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

ฉันต้องการแบบฝึกหัด / หนังสือ / สื่อการฝึกอบรมอื่น ๆ ซึ่งครอบคลุมถึงวิธีการใช้spplotฟังก์ชั่นในspแพ็คเกจใน R อย่างถูกต้อง เช่นเดียวกับฟังก์ชั่นเหล่านี้ดูเหมือนว่ามีประสิทธิภาพมาก แต่เอกสารดูเหมือนจะไม่เป็นประโยชน์เว้นแต่คุณจะรู้ว่าคุณต้องการทำอะไร ฉันพบตัวอย่างหน้าต่างๆ แต่พวกเขายังไม่ชัดเจน มีแหล่งข้อมูลใดบ้างที่จะนำคุณผ่านวิธีใช้spplotในการสร้างแผนที่ที่มีประโยชน์ (ตัวอย่างเช่นข้อมูลแรสเตอร์รวมถึงการซ้อนทับของรูปหลายเหลี่ยมบนจุดด้วยเกล็ดสีที่แตกต่างกันตำนานลูกศรเหนือ ฯลฯ ) ในแบบทีละขั้นตอน ? ฉันควรจะชอบหนังสือ UseR! -style ที่เน้นไปที่ข้อมูลแรสเตอร์และการสร้างแผนที่ใน R แต่หนังสือเชิงพื้นที่เพียงเล่มเดียวของประเภทนั้นดูเหมือนว่าจะเน้นเรื่องภูมิศาสตร์และข้อมูลจุดเป็นหลัก ความคิดใด ๆ

11 cartography  r 

ฉันพยายามทำแผนที่ระดับน้ำใต้ดินด้วย Kriging สามัญใน R แต่มักเกิดจากการขาดความจำ ไม่สามารถจัดสรรเวกเตอร์ที่มีขนาด 146.5 Mb การจัดสรร (windows xp 32 บิต) ฉันควรเพิ่มขนาดหน่วยความจำและวิธีการหรือฉันควรลดขนาดกริดและวิธีการ (ชนิดของ newbie ใน R) ทำตามตัวอย่างนี้http://casoilresource.lawr.ucdavis.edu/drupal/node/438

11 grass  r  error  kriging 

นี่เป็นกระบวนการที่ใช้เวลาเพียงไม่กี่วินาทีในซอฟต์แวร์ GIS ความพยายามของฉันใน R ใช้หน่วยความจำจำนวนมากแล้วล้มเหลว มีอะไรผิดปกติในรหัสของฉันหรือนี่เป็นเพียงสิ่งที่ R ไม่สามารถทำได้? ฉันได้อ่าน R สามารถทำงานใน Grass ฉันสามารถใช้ฟังก์ชัน Grass จากภายใน R ได้หรือไม่ library(raster) # I have many environmental rasters in this format new_r <- raster(ncol=615, nrow=626, xmn=-156.2, xmx=-154.8, ymn=18.89, ymx=20.30) res(new_r) <- 0.00225 projection(new_r) <- "+proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs +towgs84=0,0,0" R> new_r ### not too …

11 raster  grass  r 

ฉันมีสองคลาสที่แบ่งปัน CRS เดียวกัน (Latitute และลองจิจูด): bolognaQuartieriMap: SpatialPolygonDataFrameข้อมูลที่มีข้อมูลของเมือง crashPoints: SpatialPointsDataFrameข้อมูลที่มีอุบัติเหตุ พวกเขาวางแผนโดยใช้: plot(bolognaQuartieriMap) title("Crash per quartiere") plot(crashPoints, col="red",add=TRUE) สิ่งที่ฉันต้องการคือการได้รับจำนวนของจุด ( crashPoints) bolognaQuartieriMapในแต่ละรูปหลายเหลี่ยมที่เป็น ฉันแนะนำให้ใช้over()แต่ฉันไม่ประสบความสำเร็จ

11 r  sp 

ฉันพยายามเขียนหลายเลเยอร์ไปยัง 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 

ฉันมีข้อมูลบางอย่างเช่นนี้: longitude latitude speed direction 1 6.10722222 46.23639 4 360 2 6.95416667 43.54694 4 360 3 7.21472222 43.66556 13 330 4 4.01666667 48.32167 7 290 5 2.30833333 43.21611 14 290 6 2.48305556 44.40806 13 320 7 5.21500000 43.43694 19 330 8 4.92361111 43.52278 32 320 9 5.10805556 43.60306 26 330 …

11 r  leaflet  labeling  cartography  line 

ฉันต้องการสร้างแรสเตอร์ที่มีความละเอียด 25 เมตร× 25 เมตรซึ่งแต่ละเซลล์มีระยะห่างจากแนวชายฝั่งที่ใกล้ที่สุดซึ่งคำนวณจากจุดศูนย์กลางของเซลล์ การทำเช่นนี้ทั้งหมดที่ฉันต้องเป็นshapefile ของแนวชายฝั่งของนิวซีแลนด์ ฉันได้ลองทำตามบทช่วยสอนของ Dominic Royeเพื่อใช้ใน R ซึ่งใช้งานได้ดี มันปรับลงไปที่ความละเอียดประมาณ 1 km × 1 km แต่ถ้าฉันพยายามที่จะใช้ RAM ที่สูงกว่านี้ก็ต้องใช้งานได้ดีกว่าบนพีซีของฉัน (ต้องใช้ RAM 70 ~ GB) หรืออื่น ๆ ที่ฉันมีการเข้าถึงด้วย ในการบอกว่าฉันคิดว่านี่เป็นข้อ จำกัด ของ R และฉันสงสัยว่า QGIS อาจมีวิธีการคำนวณที่มีประสิทธิภาพมากกว่าในการสร้างแรสเตอร์นี้ แต่ฉันใหม่สำหรับมันและฉันไม่สามารถเข้าใจได้ว่าจะทำอย่างไร ฉันได้ลองติดตามการสร้างแรสเตอร์ด้วยระยะทางถึงฟีเจอร์โดยใช้ QGIS แล้วหรือยัง เพื่อสร้างใน QGIS แต่มันกลับข้อผิดพลาดนี้: _core.QgsProcessingException: ไม่สามารถโหลดเลเยอร์แหล่งที่มาสำหรับ INPUT: C: /..../ Coastline / …

10 qgis  shapefile  r  rasterization 

ขณะนี้ฉันกำลังรวบรวมรายการตัวอย่างจากฐานข้อมูลของพิพิธภัณฑ์ประวัติศาสตร์ธรรมชาติต่างๆเพื่อการวิจัยเพิ่มเติม อย่างไรก็ตามปัญหาที่ทราบกันดีซึ่งเกี่ยวข้องกับข้อมูลประวัติส่วนใหญ่คือการขาดละติจูดและลองจิจูดที่เหมาะสมซึ่งจะช่วยป้องกันไม่ให้มีการใช้ข้อมูลนั้น มีวิธีที่จะเอาชนะข้อมูลนั้นได้เช่นการวาดบัฟเฟอร์รอบ ๆ ภูมิภาคและจัดหาช่วงของความไม่แน่นอนที่เกี่ยวข้องกับสถานที่นั้น ตัวอย่างเช่นฟังก์ชั่น - biogeomancer จากแพ็คเกจ 'อวกาศ' ใน R จะดำเนินการกระบวนการอ้างอิงทางภูมิศาสตร์โดยอัตโนมัติหากมีคำอธิบายเกี่ยวกับใจสองสามอย่างเช่น "2 ไมล์ทางตะวันตกของ XYZ" ดูเอกสารที่นี่ อย่างไรก็ตามข้อกังวลหลักของฉันคือการใช้โปรโตคอลดังกล่าวสำหรับภูมิภาคที่มีขนาดใหญ่ถึง 200 ตารางกิโลเมตร มีวิธีหนึ่งที่สามารถเอาชนะปัญหานั้นได้หรือไม่? ฉันชอบที่จะใช้ข้อมูลพิพิธภัณฑ์มากมายหากฉันสามารถจัดการกับความไม่แน่นอนที่เกี่ยวข้องกับที่ตั้ง ตัวอย่างของตัวอย่างในชุดข้อมูลของฉันแสดงอยู่ด้านล่าง โปรดทราบว่าหลายคนมาพร้อมกับกล่าวถึงระดับความสูง แต่บันทึกส่วนใหญ่นั้นคลุมเครือมาก แก้ไข ในส่วนความเห็นคุณหนึ่งคนพูดถึงจุดประสงค์ของคำถามนี้และสิ่งที่ฉันมักจะได้รับจากสิ่งเดียวกัน 1. ฉันสนใจว่าจะลดรัศมีของความไม่แน่นอนจากขอบเขตเหลี่ยมที่กว้างจริง ๆ ให้เหลือรัศมีความไม่แน่นอนน้อยลงได้อย่างไร (ถ้าเป็นไปได้) 2. ข้อมูลนี้จะช่วยให้ฉันทำการวิเคราะห์เชิงพื้นที่ในอนาคตเช่นการสร้างแบบจำลองการกระจายสายพันธุ์ / การสร้างแบบจำลองการครอบครองตัวอย่างเช่น

10 r  coordinates  buffer  georeferencing 

ฉันกำลังพยายามใช้ ArcGIS 10.2 เพื่อสร้างบัฟเฟอร์ของจุดตามพื้นที่ที่กำหนดไว้ล่วงหน้า (เช่น 400 ตารางกิโลเมตร) ยิ่งไปกว่านั้นบัฟเฟอร์ของบางจุดอยู่ใกล้กับแนวชายฝั่งที่ต้องการให้บัฟเฟอร์ถูกตัดที่แนวชายฝั่งและยังคงมีพื้นที่เดียวกับที่อยู่ในบก (400 ตารางกิโลเมตร) ไม่มีใครรู้ว่าสิ่งนี้สามารถทำได้ด้วยตัวสร้างแบบจำลองหรือ Arcpy? ฉันมีทักษะ จำกัด กับ Arcpy และ R แต่ยินดีที่จะใช้สคริปต์บางอย่างเพื่อแก้ปัญหานี้ โปรดดูภาพด้านล่างแสดงภาพกราฟิกของสิ่งที่ฉันพยายามที่จะบรรลุ

10 arcpy  arcgis-10.2  r  buffer 

มีวิธีการกำจัด "slivers" ขนาดเล็กระหว่างรูปหลายเหลี่ยมที่ใช้Rหรือไม่? วิธีการแก้ปัญหาที่ต้องการจะสร้างใหม่SpatialPolygonsDataFrameที่ขอบเขตร่วมกันระหว่างรูปหลายเหลี่ยมเป็นเรื่องบังเอิญ ฉันสนใจโซลูชันที่ใช้Rเป็นพิเศษมากกว่า ArcMap หรือ QGIS ฉันก็อยากจะได้ยินคำอธิบายว่าทำไมช่องว่างเหล่านี้อยู่ในตอนแรก นี่คือตัวอย่างของข้อมูลเชิงพื้นที่ที่ฉันทำงานด้วย: library(rgdal) library(sp) library(tigris) library(magrittr) library(leaflet) library(gplots) # This project will use WGS 84 projected coordinate system crs_proj <- CRS("+init=epsg:4326") # These are the FIPS codes of the specific block groups in my study area sel <- c("530330079005", "530330079001", "530330079004", "530330085002", …

10 polygon  leaflet  r  tiger  slivers 

ฉันทำ semivariogram ใน R โดยใช้ package gstat, variogram()function ฉันต้องการที่จะดูว่ามี autocorrelation เชิงพื้นที่ในแบบจำลองของฉัน (ความอุดมสมบูรณ์ของสปีชีส์เป็นหน้าที่ของที่อยู่อาศัยในพื้นที่ต่าง ๆ ห่างกันเพียงไม่กี่กม. ถึง 900 กม. โดยใช้ glmm) หน่วยของฉันอยู่ในกม. และดังนั้นการตีความของฉันก็คือช่วงนั้นอยู่ที่มากกว่า 100 กม. จนกว่าการเปลี่ยนแปลงเชิงพื้นที่สัมพันธ์จะไม่เป็นปัญหาอีกต่อไป ฉันสงสัยว่าบางคนสามารถอธิบายได้ว่าทำไมนักเก็ตดูจึงสูง นี่หมายความว่าแม้ในสถานที่ใกล้เคียงกัน แต่ก็มีความแตกต่างค่อนข้างสูงใช่ไหม? หรือรูปคลื่นที่แปรปรวนแบบนี้หมายความว่าฉันควรจะปรับจำนวนความล่าช้าและระยะทางที่ช้ากว่าของฉันจนกว่าฉันจะได้รูปร่างที่มากกว่าปกติหรือไม่? เพื่อที่จะตรวจสอบเพิ่มเติมอีกเล็กน้อยฉันยังใช้ฟังก์ชั่นvariog()ในแพคเกจ geoR และใช้breaks=seq(0,100,10)เพื่อลองดูที่ระยะทางที่ใกล้กว่า (ใช้จุดเดียวกันและส่วนที่เหลือแบบเดียวกัน) สิ่งนี้บ่งชี้ว่าจุดที่ใกล้ที่สุดนั้นแตกต่างกันมากขึ้นซึ่งก็ไม่สมเหตุสมผล บางทีนี่อาจบ่งบอกว่าไม่มีความสัมพันธ์เชิงพื้นที่อัตโนมัติและแบบจำลองของฉันได้ทำสิ่งนี้แล้ว ฉันพบแหล่งข้อมูลที่ยอดเยี่ยมนี้"Geostats ที่ไม่มีน้ำตา"และในหน้า 51 มีคำแนะนำที่ดีเกี่ยวกับรูปแปรที่เหมาะสม ตามคำแนะนำนี้คนแรกของฉันดูเหมือนจะมีช่วงที่ถูกต้อง ดังนั้นสิ่งนี้กลับไปที่คำถามแรก - ฉันจะตีความสิ่งนี้ได้อย่างไร

10 r  interpolation  spatial-statistics  autocorrelation  variogram 

ฉันต้องตรวจสอบการสังเกตการณ์ของนกในระยะเวลานานขึ้นเพื่อดูรายการที่ซ้ำกัน / ทับซ้อนกัน ผู้สังเกตการณ์จากจุดต่าง ๆ (A, B, C) ทำการสังเกตการณ์และทำเครื่องหมายไว้บนแผนที่กระดาษ เส้นเหล่านั้นซึ่งนำมาสู่คุณสมบัติเส้นที่มีข้อมูลเพิ่มเติมสำหรับสปีชีส์จุดสังเกตและช่วงเวลาที่เห็น โดยปกติผู้สังเกตการณ์จะสื่อสารกันทางโทรศัพท์ในขณะที่สังเกต แต่บางครั้งพวกเขาก็ลืมดังนั้นฉันจึงได้เส้นที่ซ้ำกัน ฉันลดข้อมูลลงไปเป็นเส้นเหล่านั้นซึ่งแตะวงกลมดังนั้นฉันไม่จำเป็นต้องทำการวิเคราะห์เชิงพื้นที่ แต่เพียงเปรียบเทียบช่วงเวลาสำหรับแต่ละสปีชีส์และมั่นใจได้ว่าเป็นบุคคลเดียวกันที่พบโดยการเปรียบเทียบ . ตอนนี้ฉันกำลังมองหาวิธีใน R เพื่อระบุรายการเหล่านั้นซึ่ง: จะทำในวันเดียวกันกับช่วงเวลาที่ทับซ้อนกัน และมันอยู่ที่ไหนสายพันธุ์เดียวกัน และซึ่งทำจากจุดสังเกตที่แตกต่างกัน (A หรือ B หรือ C หรือ ... )) ในตัวอย่างนี้ฉันพบว่ามีรายการที่ซ้ำกันของบุคคลเดียวกัน จุดสังเกตแตกต่างกัน (A <-> B) สปีชีส์เหมือนกัน (Sst) และช่วงเวลาของการเริ่มต้นและสิ้นสุดของเวลาทับซ้อนกัน ตอนนี้ฉันจะสร้างเขตข้อมูลใหม่ "ทำซ้ำ" ใน data.frame ของฉันโดยให้ทั้งสองแถวมีรหัสทั่วไปเพื่อที่จะสามารถส่งออกและตัดสินใจในภายหลังว่าจะทำอย่างไร ฉันค้นหาวิธีแก้ปัญหาที่มีอยู่มากมาย แต่ไม่พบว่าเกี่ยวข้องกับความจริงที่ว่าฉันต้องเซตย่อยกระบวนการสำหรับสปีชีส์ (โดยเฉพาะอย่างยิ่งที่ไม่มีลูป) และต้องเปรียบเทียบแถวสำหรับจุดสังเกต 2 + x ข้อมูลบางส่วนที่จะเล่นด้วย: …

10 r  time  comparison  duplication 

ฉันกำลังพยายามจัดแนวแรสเตอร์สองกริดในอาร์เมื่อจัดตำแหน่งแล้วฉันต้องการที่จะเพิ่มพวกเขาเข้าด้วยกัน ฉันได้ลองตรวจสอบว่าการทำแบบstackนี้จะได้ผลหรือไม่: grid_snap <- stack(habi_sdw, Pop_sdw) และฉันได้รับข้อผิดพลาดต่อไปนี้: ข้อผิดพลาดใน comparRaster (x): ขอบเขตที่ต่างกัน ตะแกรงแรสเตอร์มีคุณสมบัติดังต่อไปนี้: show(habi_sdw) # class : RasterLayer # dimensions : 9187, 9717, 89270079 (nrow, ncol, ncell) # resolution : 0.00892857, 0.00892857 (x, y) # extent : -28.83706, 57.92186, -36.02464, 46.00214 (xmin, xmax, ymin, ymax) # coord. ref. : +proj=longlat +ellps=WGS84 …

10 raster  r  clip  snapping  resolution 

ฉันได้ให้ภาพ Geotiff และข้อมูล Lidar ที่สอดคล้องกัน (x, y, z) ในพิกัด UTM ฉันต้องรวมข้อมูล Lidar เข้ากับค่า RGB จากภาพ นั่นหมายความว่าในตอนท้ายฉันต้องพล็อต (3D) แต่ละจุดของสีคลาวด์ LiDAR ที่เข้ารหัสด้วยค่า RGB ที่สอดคล้องกันจากภาพ Geotiff ฉันแปลงข้อมูล Lidar เป็น shapefile โดยใช้ QGIS ฉันควรทำอย่างไรต่อไป ใน R ฉันลองใช้plot3Dฟังก์ชั่น แต่มันไม่ทำงาน ฉันกำลังแนบdoc text , shapefile และอิมเมจ tif แก้ไข: ฉันได้ทำโปรแกรมต่อไปนี้ตามที่แสดงด้านล่าง: require(raster) require(maptools) # to take shape files #require(car) …

10 raster  r  3d  lidar  3d-model 

เรามีLIDAR DEM 1 เมตรจากเมือง สามารถดาวน์โหลดเซ็ตย่อยขนาดเล็กได้จากลิงค์นี้ : ภาพหน้าจอนี้แสดง DEM แบบดิบพร้อมจานสีเทา (เข็มขัดที่เข้มกว่าคือถนนและสี่เหลี่ยมสีเทาและขาวเป็นบล็อค): สิ่งนี้สอดคล้องกับสถานที่ในเมืองซานโตโดมิงโกซึ่งสามารถเห็นได้ในภาพหน้าจอของ Google นี้: โดยเฉลี่ยแล้วบล็อกจะ "สูงขึ้น" โดยประมาณ 2 เมตรจากถนนซึ่งไม่เป็นความจริง เราต้องการให้ DEM สะอาดเพื่อสร้างเครือข่ายสตรีมและดัชนีความเปียกชื้นของภูมิประเทศ (TWI) ด้วย DEM ที่ให้มา (เราไม่ได้มีวงดนตรีต้นฉบับจากเลเซอร์สแกนเนอร์) ที่เครือข่าย hydrgraphic ดูเหมือนจะเป็นไปตามรูปแบบสี่เหลี่ยมและ TWI ผลในรูปแบบบล็อก ภาพเหล่านี้แสดงผลลัพธ์: นี่คือผลลัพธ์เครือข่ายสตรีมที่สร้างขึ้นด้วยr.watershedในGrass GIS: และนี่คือผลลัพธ์ TWI สร้างโดยSAGA: เราได้ลองขั้นตอนบางอย่างเพื่อแก้ไขความไม่ถูกต้องนี้โดยไม่สำเร็จ: 1) เครื่องมือ denoising เราใช้r.denoiseเครื่องมือในGrass GISแต่มีปัญหาบางอย่างกับการติดตั้งโมดูล เรารันมันอีกครั้งโดยใช้เชลล์ใน Windows และได้รับข้อความหน่วยความจำไม่เพียงพอ 2) ตัวกรอง เราวิ่งประเภทที่แตกต่างกันของตัวกรอง (low-pass …

10 qgis  r  dem  lidar  topographic-wetness-index