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

ฉันได้ติดตั้งสองแพ็คเกจต่อไปนี้: library(ggmap) library(maps) ฉันพบรหัสออนไลน์เพื่อสร้างแผนที่โลก (ดูด้านล่าง) ฉันจะซูมพล็อตในระดับประเทศได้อย่างไร ยกตัวอย่างเช่นเดนมาร์ก ggplot()+ borders("world", colour="gray50", fill="gray50")

14 r  ggplot2  ggmap 

ฉันมีคำถามเกี่ยวกับการรวมตัวเชิงพื้นที่ใน R สิ่งที่ฉันพยายามทำคือการรวมชุดข้อมูลแบบจุดไปยังกริด ฉันไม่แน่ใจอย่างไรก็ตามวิธีการทำเช่นนี้เพราะฉันมีประสบการณ์น้อยกับสิ่งเหล่านี้ ฉันหวังว่าทุกคนในคุณอาจมีแนวทางที่เป็นประโยชน์ / แนวทางแก้ไขที่เป็นไปได้ จุดได้เปรียบของฉันคือชุดข้อมูลที่มีข้อมูลอ้างอิงทางภูมิศาสตร์สำหรับเหตุการณ์ความขัดแย้งในแอฟริกา (ดู www.acleddata.com) คะแนนจะถูกอ้างอิงทางภูมิศาสตร์ด้วยพิกัดละติจูด / ลองจิจูดและมีข้อมูลเกี่ยวกับประเภทและเวลาของเหตุการณ์ สิ่งที่ฉันต้องการทำคือรวมคะแนนเหล่านี้เข้ากับกริด 1x1 องศา ดังนั้นกริดเซลล์ควรมีข้อมูลของจุดข้อมูลถ้ามีเหตุการณ์เกิดขึ้นภายในกริดเซลล์นั้น ผลิตภัณฑ์สุดท้ายนี้ควรเป็น data frame หรือบางสิ่งที่ฉันสามารถส่งออกไปยัง csv-file ได้เนื่องจากข้อมูลนั้นมีวัตถุประสงค์เพื่อใช้ในชุด data panel สำหรับการวิเคราะห์ทางสถิติ จนถึงตอนนี้ฉันโหลดและลงจุดข้อมูลและ shapefile โดยใช้รหัสด้านล่าง ฉันเชื่อว่าฉันควรใช้ฟังก์ชันoverจากแพ็คเกจspเพื่อรวม แต่ฉันไม่ทราบวิธี หวังว่าหนึ่งในคุณสามารถช่วย รหัสฉันใช้เพื่อให้ห่างไกลสามารถพบได้ที่นี่กับผลภาพที่สอดคล้องกันมากกว่ามี ข้อเสนอแนะสำหรับการทำเช่นนี้ใน QGIS ก็ยินดีต้อนรับเช่นกัน

14 raster  r  point  aggregation 

ให้วัตถุเชิงพื้นที่ใน R ฉันจะคลิปองค์ประกอบทั้งหมดของมันให้อยู่ในกล่องที่มีขอบเขตได้อย่างไร มีสองสิ่งที่ฉันต้องการจะทำ (นึกคิดฉันจะรู้วิธีที่จะทำทั้งสองอย่าง แต่ทั้งสองเป็นวิธีที่ได้รับการยอมรับในการแก้ไขปัญหาปัจจุบันของฉัน - การ จำกัด ไฟล์รูปหลายเหลี่ยมรูปร่างเพื่อทวีปสหรัฐอเมริกา) ปล่อยแต่ละองค์ประกอบไม่เต็มภายในกล่องขอบเขต ดูเหมือนว่าbbox()<-จะเป็นวิธีตรรกะ แต่ไม่มีวิธีการดังกล่าวอยู่ ทำคลิปการดำเนินงานที่แท้จริงเช่นว่าไม่ใช่องค์ประกอบเล็ก (เช่นรูปหลายเหลี่ยมเส้น) ถูกตัดออกไปในขอบเขต sp::bboxไม่มีวิธีการกำหนดดังนั้นวิธีเดียวที่ฉันคิดจะใช้overหรือgContains/ gCrossesร่วมกับวัตถุ SpatialPolygons ที่มีกล่องที่มีพิกัดของกล่องขอบเขตใหม่ จากนั้นเมื่อตัดวัตถุรูปหลายเหลี่ยมคุณจะต้องคิดออกว่ามีรูปกากบาทอยู่ที่ไหนและตัดกันและเปลี่ยนพิกัดของรูปหลายเหลี่ยมเหล่านั้นเพื่อไม่ให้เกินกล่อง หรืออะไรทำนองgIntersectionนั้น แต่แน่นอนว่ามีวิธีที่ง่ายกว่า? ในขณะที่ฉันรู้ว่ามีปัญหามากมายเกี่ยวกับกล่องที่มีขอบเขตและการซ้อนทับเชิงพื้นที่ไปยังรูปหลายเหลี่ยมที่กำหนดภูมิภาคที่น่าสนใจมักเป็นที่นิยมมากกว่าในหลาย ๆ สถานการณ์กล่องที่ทำงานได้ดีและเรียบง่ายกว่า

14 r  clip  extents 

ฉันใช้plotrixในRการพล็อตของแผนที่ระดับรัฐของสหรัฐ มันมีฟังก์ชั่นที่ยอดเยี่ยมfloating.pieในการวางแผนภูมิวงกลมในแต่ละรัฐ ฉันสงสัยว่ามีฟังก์ชั่นที่คล้ายกันภายในplotrixแพ็คเกจเพื่อแสดงแผนภูมิแท่งในแต่ละรัฐหรือไม่? (ฉันลองดูที่เอกสารประกอบและฟังก์ชั่นที่จัดการกับแผนภูมิแท่งดูเหมือนจะไม่มีความเป็นไปได้ แต่ฉันแค่อยากจะแน่ใจ) ฉันชอบทำงานในแพ็คเกจพล็อตริกซ์ถ้าเป็นไปได้ แต่รู้สึกอิสระที่จะตั้งชื่อ แพคเกจอื่น ๆ ตัวอย่างเช่นฉันสนใจที่จะสร้างแผนที่คล้ายกับแผนที่นี้ (แต่สำหรับสหรัฐอเมริกา): สำหรับแผนที่อเมริกาของฉันจะมีแผนภูมิแท่ง 50 แผนภูมิแต่ละแผนภูมิ ฉันได้รับแผนที่นี้จาก/programming/20465070/barplots-on-a-mapแต่ดูเหมือนว่า ggsubplot ไม่ทำงานในเวอร์ชัน R ของฉัน (คล้ายกับที่คนอื่นพูดในโพสต์ )

14 r  statistics  chart 

ฉันมีไฟล์แรสเตอร์ขนาดใหญ่ (129,600 คูณ 64,800 พิกเซล) ที่มีแหล่งน้ำทั่วโลก (1 บิตค่า 0 และ 1) และพยายามที่จะแยกมหาสมุทรและชายฝั่งน้ำจืด ฉันได้ลองกับ ArcGIS และ QGIS เพื่อแปลงจาก raster เป็น polyline แต่ใช้เวลานาน มีใครรู้วิธีที่ดีกว่า / เร็วกว่า (Python หรือ R) หรือเป็นเครื่องมือที่ดีกว่าสำหรับงานนี้หรือไม่? ปรับปรุง R: rasterToContour อาจเร็วและแม่นยำ แต่ถ้าคุณมีชุดข้อมูลขนาดใหญ่มากเช่นฉัน (8,398,080,000 พิกเซล) คุณต้องมี RAM จำนวนมาก (มากกว่า 16GB) หรือบังคับ R เพื่อทำการประมวลผลเพิ่มเติมบนฮาร์ดไดรฟ์และมัน จะใช้เวลานานด้วย Python / GDAL: gdal_poligonize สร้างรูปหลายเหลี่ยมแทนโพลีน อัปเดต …

14 python  r  performance 

ฉันได้สร้างวัตถุชนิดSpatialPointsDataFrameโดยใช้spแพคเกจใน R อย่างไรก็ตามฉันสับสนเกี่ยวกับตัว@, $, . and []ดำเนินการและเมื่อใดจึงใช้เพื่อเข้าถึงคุณสมบัติที่แตกต่างกันของวัตถุของฉัน นี่คือตัวอย่างรหัสของฉัน: library(sp) library(rgdal) #creating a SpatialPointsDataFrame with sample points in UTM x <- c(15.2, 15.3, 15.4, 15.5, 15.7) y <- c(50.4, 50.2, 50.3, 50.1, 50.4) v1 <- c(1.0, 2.0, 3.0, 4.0, 5.0) v2 <- c("a","b","b","c","a") attributes <- as.data.frame(cbind(v1,v2)) xy <- cbind(x,y) locationsDD <- …

14 r  syntax 

ดูเหมือนว่าตัวเลือกการประมวลผลใน QGIS 3.0 จะไม่มี R Scripts อีกต่อไป (ไม่มีอยู่ในรายการพร้อมกับ SAGA, GRASS และ GDAL) นี่เป็นการเปลี่ยนแปลงใน QGIS รุ่นนี้หรือฉันมีปัญหากับการติดตั้งหรือไม่?

13 qgis  r  qgis-3 

ฉันพยายามสร้างภาพแรสเตอร์จากฐานข้อมูลจุดที่เว้นระยะไม่สม่ำเสมอ ข้อมูลดูเหมือนว่า - > head(s100_ras) x y z 1 267573.9 2633781 213.29545 2 262224.4 2633781 69.78261 3 263742.7 2633781 51.21951 4 259328.4 2633781 301.98413 5 264109.8 2633781 141.72414 6 255094.8 2633781 88.90244 ฉันต้องการค่า 'z' เหล่านี้ภายในตาข่ายที่ฉันสร้างขึ้น # Create a fine mesh grid my_mesh=expand.grid(seq(min(s100_ras$Y),max(s100_ras$Y),l=100), seq(min(s100_ras$X),max(s100_ras$X),l=100)) ฉันต้องการให้ค่า z ถูกกำหนดเป็น 'NA' สำหรับจุดตาข่ายที่อยู่นอกจุดข้อมูล คะแนนเหนือตาข่ายมีลักษณะดังนี้: https://drive.google.com/file/d/0B6GUNg-8d30vYzlwTkhvaHBFTnc/edit?usp=sharingเมื่อฉันลงจุด …

13 raster  r  vector-grid 

ฉันมีรูปหลายเหลี่ยมสองรูปนี้: library(sp); library(rgeos); library(maptools) coords1 <- matrix(c(-1.841960, -1.823464, -1.838623, -1.841960, 55.663696, 55.659178, 55.650841, 55.663696), ncol=2) coords2 <- matrix(c(-1.822606, -1.816790, -1.832712, -1.822606, 55.657887, 55.646806, 55.650679, 55.657887), ncol=2) p1 <- Polygon(coords1) p2 <- Polygon(coords2) p1 <- Polygons(list(p1), ID = "p1") p2 <- Polygons(list(p2), ID = "p2") myPolys <- SpatialPolygons(list(p1, p2)) spdf1 = …

13 r  polygon  vector  clip  intersection 

คำถามของฉันเกี่ยวกับการแยกค่าจาก rasters ที่ตำแหน่งของจุด ด้วยฟังก์ชั่นการแยกมันง่ายมากและฟังก์ชั่นให้ฉันดาต้าเฟรมที่มีค่าของตัวแปรทั้งหมดในจุด ฉันต้องการที่จะมีใน dataframe พิกัดของแต่ละจุด ฉันจะทำให้เกิดขึ้นได้อย่างไร เป็นไปได้ไหมที่จะบอกกับ R ว่าเมื่อทำการแยกค่าจากแรสเตอร์จะเพิ่มคอลัมน์ของตำแหน่งที่ตั้ง นี่คือขั้นตอนของฉัน: presencias=read.table("c:/SDM_R/presencias/P_lentiscus_pres.csv",header=TRUE,sep=";") lista_variables <-list.files(path="Variables_modelizacion/solo_ascii",pattern='*.asc',full.names=TRUE) variables <- stack(lista_variables) variables_presencia<-extract(variables, presencias) ผลลัพธ์เป็นดังนี้: > bio1 bio12 bio18 bio2 bio4 > 90 875 165 95 4886 > 115 1085 158 83 4075 > 135 1153 153 67 3402 > 85 1026 137 99 …

13 raster  r 

ฉันมีคำถามเกี่ยวกับการพล็อตแผนที่เป็นนักร้อง สิ่งที่แนะนำเกี่ยวกับการพล็อตแผนที่ด้วยตัวแบ่งที่กำหนดเอง? และทำให้ดูน่าสนใจขึ้นมาหน่อย ในสคริปต์ของฉันฉันใช้คำสั่ง plot และ spplot เพื่อแสดงข้อมูลบางอย่าง แต่ฉันไม่พอใจกับผลลัพธ์ทั้งหมด มันง่ายที่จะดูว่าทำไม ใช้spplotฉันได้รับทางลาดสีที่ทำซ้ำสีสำหรับค่าที่แตกต่างซึ่งเป็นบิตแปลก การใช้พล็อตที่มีตัวแบ่งที่กำหนดเองฉันไม่สามารถทำให้ช่วงเวลาที่ดีขึ้นสำหรับค่าเช่นว่าค่าศูนย์เป็นสีขาวบนแผนที่ ฉันคิดว่าสัดส่วนที่เบ้ระหว่างคำอธิบายแผนภูมิและแผนที่อาจเป็นผลมาจาก Rstudio ดังนั้นฉันสงสัยว่ามีใครบางคนมีคำแนะนำที่ดีกว่าสำหรับการวางแผนแผนที่ choropleth ใน R โดยใช้ตัวแบ่งที่กำหนดเองหรือไม่? แก้ไข 20-02-2013 ตามที่แนะนำฉันได้ใช้ ggplot2 เพื่อสร้างแผนที่ choropleth ฉันใช้รหัสด้านล่าง ฉันพบปัญหาสองข้อ หนึ่งคือหลังจากเสริมข้อมูลไม่ได้ข้อมูลทั้งหมดจาก shapefile จะถูกโอนไปยังกรอบข้อมูล ดังนั้นฉันต้องเปลี่ยนชื่อตัวแปร "id" เพื่อที่ฉันจะได้รวมข้อมูล ปัญหาเล็กน้อย ปัญหาที่ใหญ่กว่าคือแม้ว่าข้อเท็จจริงที่ว่าข้อมูลใน dataframe นั้นถูกต้อง ggplot จะแสดงแผนที่ที่ข้อมูลไม่ถูกต้อง ค่าสีตามคำอธิบายไม่สอดคล้องกับค่าของตัวแปรในข้อมูล ฉันมองอะไรบางอย่างที่นี่ # Plotting polygon shapefiles (try-out) # Load the packages …

13 r  choropleth 

ฉันกำลังวิเคราะห์เส้นทางที่เป็นไปได้ที่ผู้สัญจรใช้ในการเดินทางตามข้อมูลการสำรวจสำมะโนประชากรการกระจายความหนาแน่นของงานและเครือข่ายการขนส่ง ในขณะที่ฉันใช้ pgRouting และ pgRouting add-on สำหรับ QGIS (อธิบายไว้ที่นี่ ) สำหรับการวิเคราะห์ของฉัน: อย่างไรก็ตามฉันต้องการขั้นตอนทั้งหมดที่ฉันทำเพื่อทำซ้ำแบบเต็ม (ดูบทความที่ยอดเยี่ยมในที่นี่ ) จากนั้นฉันต้องการที่จะทำการวิเคราะห์ของฉันจากบรรทัดคำสั่งและบันทึกไฟล์สคริปต์สำหรับการอ้างอิงในอนาคต แพ็คเกจ R-geo ตอบสนองความต้องการนี้ได้ดีเนื่องจากมีอิสระน้ำหนักเบาและเสถียร: ดาวน์โหลดไฟล์สคริปต์ขนาดเล็กและคุณควรจะสามารถทำทุกสิ่งที่ฉันทำใหม่ได้ (ดูที่นี่สำหรับตัวอย่างที่ไม่เกี่ยวข้องกับ Rpubs ) แพ็คเกจ R ที่ฉันได้ระบุไว้สำหรับการทำเช่นนี้คือ: gdistanceตีพิมพ์เมื่อธันวาคม 2555 ดูมีแนวโน้ม แต่ไม่มีสัญญาณของความเข้ากันได้กับเครือข่ายการขนส่ง osm e1071ดูเหมือนจะค่อนข้างเป็นนามธรรมและอาจเกินความเป็นจริงสำหรับสิ่งที่ฉันต้องการจะทำ ฟังก์ชั่นpairdist.lppในแพคเกจ spatstat ดูมีแนวโน้ม แต่มีเอกสารเล็กน้อย ก่อนที่จะดำน้ำให้คิดว่ามันคุ้มค่าที่จะถามคนที่อาจต้องเผชิญกับปัญหาที่คล้ายคลึงกัน เป็นไปได้ไหมที่จะทำการวิเคราะห์เส้นทางของฉันใน R หรือจะเป็นการดีที่สุดที่จะเรียนรู้ PostGIS อย่างถูกต้องและทำทุกอย่างใน pgRouting?

13 pgrouting  r  routing  osm2pgrouting 

ฉันต้องการฉายแผนที่นี้ในการฉายภาพของ Robinson: library(ggmap) world <- map_data("world") ggplot() + geom_path(data = world, aes(long, lat, group = group)) และฉันต้องการเปลี่ยนการฉายเป็น "Robinson" (ทำตามคำแนะนำจากคำตอบของคำถามก่อนหน้านี้: การคาดการณ์ใดที่แผนที่ภูมิอากาศโลกจาก Wikipedia ใช้? ฉันมีปัญหาในการหาการใช้งานเริ่มต้นของการฉายภาพนี้ฉันทำงานต่อไปนี้สำหรับการใช้proj4ห้องสมุด: library(proj4) robinson <- project(cbind(world$long, world$lat), proj = "+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs") ฉันได้ลองหลายวิธีแล้วรวมถึง: # using ggmap::get.map() get_map("world", projection = mapprojection(robinson)) # using ggplot2::coord_map coord_map(projection …

13 coordinate-system  r  proj  ggmap 

มีวิธีที่ตรงไปตรงมาในการสุ่มตัวอย่างแรสเตอร์เพื่อให้ผลลัพธ์ของกระบวนการเป็นแรสเตอร์หรือไม่? ฉันใช้ตัวอย่างที่ฉันพบในr-sig-geo รายการและฉันได้ลองใช้sampleRandomฟังก์ชันในrasterแพ็คเกจด้วย วิธีการทั้งสองนี้สร้างผลลัพธ์ที่ฉันไม่แน่ใจว่าจะเปลี่ยนเป็นแรสเตอร์ได้อย่างไร ฉันไม่สามารถค้นหาวิธีการหลังจากค้นหาชุดค่าผสมหลายชุดของ "SpatialPointsDataFrame raster" library(raster) # read in raster rasterSource <- 'landsat.TIF' r <- raster(rasterSource) # convert to spatial points data frame r.spgrd<-as(r,"SpatialPointsDataFrame") # elminate NA values r.spgrd = r.spgrd[!is.na(r.spgrd[[1]]),] # sample points selectedPoints = sample(1:length(r.spgrd[[1]]), 1000) r.sampled = r.spgrd[selectedPoints,] # try to make spgrd into a …

13 raster  r 

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

13 raster  r  rgdal