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

ฉันได้ทำการทดลองกับเทคนิคการทำแผนที่ choropleth ใน R โดยได้รับประโยชน์จาก ggplot2 เพื่อสร้างกราฟที่สวยงามภายในแพ็คเกจการวิเคราะห์ข้อมูลที่ทรงพลัง เมื่อพูดถึงการทำแผนที่ฉันไม่สามารถสร้างผลลัพธ์ที่เทียบเคียงได้กับสิ่งที่ฉันได้รับจาก QGIS เป็นประจำ การใช้ข้อมูลการไหลใน Sheffield (ข้อมูลที่จำลองได้อย่างสมบูรณ์และไฟล์. qgs ที่นี่ ) ทำให้ QGIS สร้างรายการต่อไปนี้ได้อย่างง่ายดาย: สิ่งที่ดีที่สุดที่ฉันสามารถทำได้โดยใช้ R (โดยใช้รหัสที่อธิบายไว้ที่นี่ ) คือแผนที่ทั้งสองนี้แผนที่แรกที่ใช้ GISTools ส่วนที่สองใช้ ggplot2: สำหรับฉันแล้วตัวเลือก ggplot2 ดูน่าสนใจยิ่งกว่านี้หากฉันสามารถแก้ไขปัญหาของบรรทัดที่ผิดพลาด (อาจเป็นปัญหากับคำสั่ง fortify () หรือไม่อ่านในไฟล์รูปร่างโดยใช้ readOGR () อธิบายไว้ที่นี่ ) ดังนั้นคำถามคือ 2-fold: เป็นตัวเลือก ggplot2 ที่ดีที่สุดในการแก้ปัญหาการทำแผนที่ choropleth ใน R และถ้าเป็นเช่นนั้นฉันจะแก้ปัญหาของเส้นสีขาวที่ผิดพลาดได้อย่างไร? รหัสจำลองเพื่อหาสิ่งที่ผมทำคือที่นี่ แก้ไข - …

18 r  thematic-map 

ฉันอยากจะสรุป a SpatialPolygonsDataFrame ใน R เพื่อให้การวางแผนมีประสิทธิภาพมากขึ้น รูปหลายเหลี่ยมที่ฉันส่งผลให้เป็นไฟล์ PDF ขนาดใหญ่ที่โหลดช้าและยากที่จะรวมเข้ากับเอกสาร ฉันพยายามใช้gSimplifyฟังก์ชั่นจากrgeosแพ็คเกจ แต่น่าเสียดายที่มันไม่ได้อนุรักษ์โทโพโลยีของรูปหลายเหลี่ยม มีทางเลือกอื่นที่ฉันสามารถใช้ได้หรือไม่?

18 r  topology  generalization 

ขณะนี้ฉันกำลังใช้สคริปต์ต่อไปนี้เพื่อเพิ่มข้อมูลคุณลักษณะจากตารางไปยังรูปร่างหลายไฟล์: library(rgdal) specieslist <- read.csv("SpeciesList1.txt", header=F) attdata <- read.table("TestAtt.csv", sep = ",", header=T) for (n in 1:dim(specieslist)[1]) { speciesname <- specieslist[n,1] shp <- readOGR("Mesoamerica_modified_polygons", speciesname) shp$ENGL_NAME<-attdata[n,2] writeOGR(shp, "PolygonsV2", speciesname, driver="ESRI Shapefile") } ฉันได้รับคำเตือนต่อไปนี้เมื่อสิ้นสุด: 1: In writeOGR(shp, "PolygonsV2", speciesname, driver = "ESRI Shapefile") : Field names abbreviated for ESRI Shapefile driver …

18 shapefile  gdal  r  ogr  attribute-table 

ด้วย RI ต้องการที่จะสร้างแผนที่ที่มีชีวิตชีวาด้วยrMapsแพคเกจและแผนที่ที่กำหนดเอง ผมพยายามที่จะทำตามตัวอย่างของดิเอโกวัลที่อธิบายไว้ที่นี่และการใช้ข้อมูลแผนที่ที่กำหนดเองสำหรับ Datamaps อย่างไรก็ตาม แผนที่ไม่แสดง แต่คำอธิบายสัญลักษณ์พร้อมคีย์เติมที่เหมาะสมจะปรากฏขึ้นพร้อมแถบภาพเคลื่อนไหวปี ดูเหมือนไฟล์ topoJson จะใช้ได้ library(rCharts);library(rMaps);library(rgdal) map = readOGR("https://raw.githubusercontent.com/gupierre/GHA_vol/master/Data/GHA_adm/GHA_adm1.topo.json", layer="admin_1",stringsAsFactors=F) #Data: SDLOG<-read.table("https://raw.githubusercontent.com/gupierre/GHA_vol/master/Data/sdlog.csv", header=T,stringsAsFactors=F, sep =",") head(SDLOG) year name sdlog type >1 2004 Greater Accra 0.20030687 urban >2 2005 Greater Accra 0.15322712 urban >3 2006 Greater Accra 0.10782677 urban #Map: library(rCharts);library(rMaps) d1 <- ichoropleth(sdlog ~ …

17 r  topojson  choropleth 

ฉันมีSpatialPointsDataFrameข้อมูลเพิ่มเติมบางอย่างด้วย ฉันต้องการที่จะดึงจุดเหล่านั้นเข้าไปในรูปหลายเหลี่ยมและในเวลาเดียวกันให้เก็บรักษาSPDFวัตถุและข้อมูลที่เกี่ยวข้อง จนถึงตอนนี้ฉันมีโชคเล็กน้อยและใช้วิธีจับคู่และรวมกันผ่าน ID ทั่วไป แต่วิธีนี้ใช้ได้เพราะฉันมีข้อมูล gridded กับ IDS แต่ละรายการเท่านั้น นี่คือตัวอย่างย่อ ๆ ฉันกำลังมองหาจุดที่อยู่ภายในจัตุรัสแดง library(sp) set.seed(357) pts <- data.frame(x = rnorm(100), y = rnorm(100), var1 = runif(100), var2 = sample(letters, 100, replace = TRUE)) coordinates(pts) <- ~ x + y class(pts) plot(pts) axis(1); axis(2) ply <- matrix(c(-1,-1, 1,-1, 1,1, -1,1, -1,-1), …

17 r  overlay 

ฉันยังใหม่กับการวิเคราะห์เชิงพื้นที่และขอขอบคุณทิศทางทั่วไปในโครงการที่ฉันพยายามอธิบายไว้ด้านล่าง (ฉันเริ่มจากศูนย์) เป้าหมาย:เพื่อค้นหาสถานที่ที่ดีที่สุดในการติดตั้งทางเท้า 2000 ฟุตในบ้านเกิดของฉันเพื่อเชื่อมต่อครัวเรือนส่วนใหญ่กับย่านธุรกิจกลาง (CBD) ที่ซึ่ง "เชื่อมต่อ" หมายถึงเดิน 1.2 กิโลเมตรจาก CBD ฉันมีรูปร่างที่แสดงโครงสร้างที่มีอยู่ (ครัวเรือน) ถนนและทางเท้า (ติดตั้งแล้ว) นี่คือโซลูชัน / กระบวนการคิดที่เสนอของฉัน: แปลงเครือข่ายทางเท้าในสถานที่เป็นฐานข้อมูลของโหนดที่เชื่อมต่อด้วยน้ำหนัก (เช่นระยะทาง) มีวิธีในการทำสิ่งนี้โดยตรงใน QGIS (หรือโปรแกรมอื่น ๆ ) โดยคลิกที่ทางแยกทั้งหมดหรือไม่? คำนวณจำนวนครัวเรือนที่อยู่ในระยะเดิน 1.2 กม. จากย่านธุรกิจกลาง (เช่นจุด lat-long หรือรูปหลายเหลี่ยม) โดยใช้ความสามารถในการกำหนดเส้นทางของ pgRouting หรืออย่างอื่น นี่จะเป็นกรณีฐาน "ค่าการเข้าถึงในครัวเรือน" การใช้เลเยอร์ถนนเป็นแนวทางให้วางเท้าเพิ่มอีก 2000 ฟุต (พูดในเซ็กเมนต์ 10 ฟุต) ของทางเท้าลงบนเลเยอร์ทางเท้า นี่คือสิ่งที่เทียบเท่ากับการสร้างทางเท้าใหม่จำนวนมากโดยพลการ คำนวณโหนดและตุ้มน้ำหนักใหม่โดยใช้เครือข่ายคนเดินเท้าใหม่ตาม (1) แล้วคำนวณจำนวนครัวเรือนที่อยู่ภายใน …

17 qgis  distance  pgrouting  r  simulation 

โพสต์นี้เป็นวิกิพีเดียชุมชน แก้ไขคำตอบที่มีอยู่เพื่อปรับปรุงโพสต์นี้ ขณะนี้ไม่ยอมรับคำตอบใหม่ ฉันสนใจสถิติทางภูมิศาสตร์อย่างมาก (และอีกหลายคำที่มี 3 t ของ :) ฉันไม่เอนเอียงทางโปรแกรม แต่ได้ติดตามรายการหลายรายการด้วยชุดรูปแบบนี้ หนึ่งในนั้นคือตัวกรอง R บน stackexchange (ฉันได้รับอีเมลสรุปคำถามเป็นครั้งคราว) ฉันพบว่ามีภาษาที่ซับซ้อนอย่างยิ่งที่เกี่ยวข้องกับสถิติโดยทั่วไป ฉันรู้สึกสับสนในการเรียนรู้จากมุมมองของบุคคลทั่วไปความหมายของภาษานี้และวิธีการใช้กับ GIS ในชีวิตประจำวัน กรุณาระบุแหล่งที่มาของกระแสสำหรับความรู้ / ความเข้าใจทางภูมิศาสตร์ นี่คือตัวอย่างของข้อมูลที่เป็นประโยชน์ ... เดินทางเพื่อฝึกอบรม โอกาสออนไลน์ แก้ไข: ต่อคำขอ; ฉันสนุกกับการฝึกอบรมออนไลน์ แต่ได้รับเพิ่มเติมจากการเรียนในห้องเรียนหรือทีละขั้นตอน (หรือ PDF)

17 r  training 

ฉันมีรายชื่อของบัฟเฟอร์เชิงพื้นที่ (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 

coord. ref. : NAฉันมีชั้นด้วย ฉันcoord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0จะกำหนดให้ได้coord. ref. : NAอย่างไร class : RasterLayer dimensions : 127, 146, 18542 (nrow, ncol, ncell) resolution : 0.01, 0.01 (x, y) extent : -122.39, -120.93, 36.72, 37.99 (xmin, xmax, ymin, ymax) coord. ref. : NA data source : C:\2012Tif\2012ASC5min_ppt_spas1306_0001_20121222_0005_UTC.asc.tif …

16 raster  coordinate-system  r 

ฉันต้องการทำคลิป GIS จริงใน R ของรูปหลายเหลี่ยมโดยใช้ชุดรูปหลายเหลี่ยมแบบเขตแดนเดียว แต่ฉันไม่พบฟังก์ชัน R เพื่อทำอย่างถูกต้อง มันควรจะทำงานเหมือนกับclipฟังก์ชันใน ArcMap ของ ESRI ฉันลองใช้overวิธีการในspแพ็คเกจ แต่ดูเหมือนว่าจะใช้งานไม่ได้กับ polys มากกว่า polys หนึ่งข้อเสนอแนะคือการใช้gIntersectionในrgeosแพคเกจเป็นคลิปโดยใช้รหัสต่อไปนี้: #------------------------------------ library(rgeos) library(maptools) #Read layers as SpatialPolygonsDataFrame (both the same Albers projection) Soils_poly = readShapePoly("Soils_polygons") #Note - Has 400 polygons clipper_poly = readShapePoly("clipper_polygon") #Note - Has 1 polygon #Try gintersection as clip …

16 r  polygon  layers  vector  clip 

ฉันค่อนข้างใหม่กับ QGIS และฉันกำลังค้นหาวิธีง่ายๆในการสร้างฐานข้อมูลที่ช่วยให้ฉันใช้ตารางใน QGIS และ R ฉันมีข้อมูลจำนวนมากที่เก็บไว้ใน Excel สเปรดชีตซึ่งจะใช้สำหรับการวิเคราะห์ทางสถิติ ในอาร์ งานของฉันในรายละเอียดเพิ่มเติม: งานแรกของฉันจะทำการแปลงข้อมูลให้เป็นรูปเป็นร่างโดยการเข้าร่วมตารางกับ Shapefile ที่สองฉันต้องการสร้างการเปลี่ยนแปลงเชิงพื้นที่ภายใน QGIS ด้วย shapefiles และใช้ข้อมูลที่สร้างขึ้นใหม่ในการวิเคราะห์ทางสถิติของฉัน ดังนั้นฉันจะต้องเข้าร่วมข้อมูลนี้กับฐานข้อมูลของฉัน ท้ายสุด แต่ไม่ท้ายสุดฉันไม่อยากมั่นใจว่าข้อมูลจะถูกจัดเก็บในทางที่ดีสำหรับการวิเคราะห์ในอนาคตของ QGIS และ R เนื่องจากฐานข้อมูลจะถูกใช้ในบริบทอื่น ๆ บางคนโทรมฉันว่า MS Access เป็นวิธีที่ดีในการทำเช่นนั้น ทุกคนในคุณสามารถให้คำแนะนำกับฉันได้ไหม? อาจจะเกี่ยวกับวรรณกรรมและอะไรทำนองนั้นด้วย

16 qgis  r  spatial-database  database  ms-access 

ฉันมีสองชั้น ชั้นรูปหลายเหลี่ยมที่มีแผ่นกระเบื้องจำนวนมากและชั้นแรสเตอร์ที่มีที่ดิน CORINE 2006 ครอบคลุมหลายประเภทใน colourmap ฉันต้องการได้รูปหลายเหลี่ยมใน shapelayer สำหรับผลรวมของหมวดหมู่ที่ดินปกคลุมของแรสเตอร์เลเยอร์ ตัวอย่างเช่นมีรูปหลายเหลี่ยมที่มี id '2' และฉันต้องการคุณสมบัติเช่นนี้สำหรับรูปหลายเหลี่ยมนี้ (เป็นเปอร์เซ็นต์หรือตารางเมตร): ที่ดินทำกิน: 15% ป่าไม้: 11% ถนน: 2% (... และอีกอย่างหนึ่ง) ฉันพยายามที่จะทำมันในสนามหญ้า, qgis (ไม่มีฟังก์ชั่น), นิยายเกี่ยวกับวีรชน (เพียงผลรวมทุกค่ารวม) r (ผลรวมทั้งหมด) แต่ฉันยังไม่พบวิธีแก้ปัญหา ปลั๊กอินส่วนใหญ่ (สถิติโซนใน qgis) รองรับเลเยอร์แรสเตอร์ 0-1 เท่านั้น v.rast.stats ก็ไม่ได้ช่วยอะไรเช่นกัน เอี่ยมเปิดให้โซลูชั่นที่ดีและชาญฉลาด! บางทีฉันอาจใช้วิธีที่ผิดหรือทำผิดพลาด ใน Arcgis งานนี้ค่อนข้างง่ายถ้าจำถูกต้อง แต่ฉันยังขาดวิธีที่ดีสำหรับผู้ใช้ linux ทุกวัน ฉันใช้ระบบ debian linux และนี่คือสาเหตุที่ฉันสามารถใช้โปรแกรมสำหรับระบบปฏิบัติการนี้เท่านั้น …

16 qgis  raster  polygon  r  statistics 

ฉันกำลังโยกย้ายรหัสจาก sp แพ็คเกจไปยังแพ็คเกจ sf ที่ใหม่กว่า รหัสก่อนหน้าของฉันฉันมีรูปหลายเหลี่ยม SpatialDataFrame (censimentoMap) และ SpatialPointDataFrame (indirizzi.sp) และฉันได้รับรูปหลายเหลี่ยมเซลล์ id ("Cell110") สำหรับแต่ละจุดที่วางไว้ภายในด้วยคำแนะนำด้านล่าง: points.data <- over(indirizzi.sp, censimentoMap[,"Cell110"]) ที่จริงฉันสร้างวัตถุ sf สองอัน: shape_sf <- st_read(dsn = shape_dsn) shape_sf <- st_transform(x=shape_sf, crs=crs_string) และ indirizzi_sf = st_as_sf(df, coords = c("lng", "lat"), crs = crs_string) และฉันกำลังมองหาเอสเอฟเทียบเท่าคำสั่งข้างต้น ... ids<-sapply(st_intersects(x=indirizzi_sf,y=shshape_sfpeCrif), function(z) if (length(z)==0) NA_integer_ else …

16 r  overlay  sp  sf 

ฉันต้องการค้นหาแพคเกจใน R หรือซอร์สโค้ดใน R ที่มีการสร้างรหัสภูมิศาสตร์โดยใช้ Bing, Yahoo, OpenStreetmap คุณมีข้อเสนอแนะใด?

15 r  geocoding  batch  bing-maps 

ฉันกำลังมองหาไฟล์รูปร่างถนนสำหรับออสเตรเลียในวันนี้ ฉันสิ้นสุดขึ้นไปGeofabrik Open Street Map มีข้อมูลใน.osmรูปแบบ แต่พวกเขาไม่ได้มีหนึ่งใน.shpสำหรับออสเตรเลีย ถ้าเป็นไปได้ฉันต้องการแปลงosmไฟล์เป็นshapefilesใช้ R. ฉันได้ค้นหาไปแล้ว แต่ยังไม่พบวิธีแก้ไข มีวิธีการนี้หรือไม่? อีกวิธีหนึ่งเป็นไปได้ไหมที่จะอ่านไฟล์ osm เป็น R และแปลงข้อมูลเป็น data frame เพื่อให้ฉันสามารถวาดรูปโดยใช้ggplot2? ฉันได้ประสบความสำเร็จในภาพต่อไปนี้ของนิวซีแลนด์ใช้สอง shapefiles ด้วยและmaptools ggplot2เป็นการดีที่ฉันต้องการสร้างแผนที่ที่คล้ายกันของออสเตรเลีย

15 shapefile  openstreetmap  r  ggplot2