ระบบสารสนเทศทางภูมิศาสตร์

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

41 r  cartography 

ดูเหมือนว่าฉันจะชอบคำถามง่ายๆ (และอาจเป็น) แต่ฉันก็ไม่สามารถหาตัวอย่างที่ให้คำตอบได้ ใช้ PostGIS ฉันแค่ต้องการเลือกจุดที่อยู่นอกรูปหลายเหลี่ยม ท้ายที่สุดนี่คือสิ่งที่ตรงกันข้ามของ ST_Intersects เท่าที่ฉันจะเห็น ตัวอย่าง: ฉันมีเลเยอร์ภาษีและชั้นที่อยู่ ฉันคิดว่าฉันควรใช้ ST_Intersects แต่ฉันจะบอกให้ทำการเลือกแบบย้อนกลับได้อย่างไร ฉันคิดว่าอาจเพิ่มคำสั่งNOTด้านหน้าโค้ดด้านล่าง แต่นั่นไม่ได้ผล CREATE table t_intersect AS SELECT hp.gid, hp.st_address, hp.city, hp.st_num, hp.the_geom FROM public.parcel as par, public.housepoints as hp WHERE ST_Intersects(hp.the_geom,par.the_geom);

41 postgis  geoprocessing  sql  intersection 

ฉันใช้เวลาเล็กน้อยในการหาคำตอบสำหรับคำถามนี้ ไม่ชัดเจนในทันทีจากการค้นหาของ 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 

หากฉันมีการแจกแจงแบบแผนสองจุดในพื้นที่ทางภูมิศาสตร์เดียวกันฉันจะเปรียบเทียบการแจกแจงสองแบบนี้ได้อย่างไร สมมติว่าฉันมีหลายจุดในพื้นที่เล็ก ๆ ดังนั้นการแสดงแผนที่พินนั้นไม่เป็นไปตามปกติ

41 point  spatial-statistics  visualisation  gis-principle 

สำหรับโครงการที่ฉันกำลังดำเนินการอยู่ฉันต้องเป็นตัวแทนเครือข่ายที่ไหลผ่านโซน 1,000 โซนด้วยทิศทาง มีข้อเสนอแนะใด ๆ เกี่ยวกับการแสดงข้อมูลในแบบที่ดีและไม่ท่วมท้นหรือไม่? ตัวอย่างใด ๆ ในการแสดงกระแสเครือข่ายในเขตเมืองใหญ่ ข้อเสนอแนะสีและโครงร่างความหนา? "ไม่ต้อง dos" ใด ๆ เมื่อแสดงกระแสเครือข่าย

41 network  transportation  flow-map 

ฉันใช้ปลั๊กอินเลเยอร์ QGIS ของ Google เพื่อแปลงพัสดุเป็นรูปแบบไฟล์ หลังจากนั้นฉันต้องการนำเข้ามาที่ PostGIS ฉันรู้ว่า Google แผนที่ใช้การฉายภาพ 'google mercator' พิเศษ รหัสที่ถูกต้องของการฉายภาพนี้คืออะไร: 900913 หรือ 3857 หรือสองรหัสนี้อธิบายสมการการฉายภาพที่เหมือนกันหรือไม่

41 google-maps  coordinate-system  epsg 

ฉันทำงานกับไฟล์. 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 

ในความเป็นจริงเมื่อ Sinnott เผยแพร่สูตร haversine ความแม่นยำในการคำนวณมี จำกัด ทุกวันนี้ JavaScript (และคอมพิวเตอร์และภาษาที่ทันสมัยที่สุด) ใช้เลขทศนิยมแบบ 64 บิตของ IEEE 754 ซึ่งให้ความแม่นยำ 15 ตัวเลขที่สำคัญ ด้วยความแม่นยำนี้กฎหมายทรงกลมเรียบง่าย ของสูตรโคไซน์ ( cos c = cos a cos b + sin a sin b cos C) ให้ผลลัพธ์ที่มีสภาพดีไปจนถึงระยะทางที่มีขนาดเล็กเพียงประมาณ 1 เมตร ในมุมมองนี้มันอาจจะคุ้มค่าในสถานการณ์ส่วนใหญ่โดยใช้กฎที่เรียบง่ายของโคไซน์หรือสูตรรูปไข่ทรงรีที่แม่นยำยิ่งขึ้นเพื่อใช้กับแฮเวอร์ซีน! (โปรดจำไว้ว่าข้อ จำกัด ด้านความแม่นยำของแบบจำลองทรงกลมด้านล่าง) ที่มา: http://www.movable-type.co.uk/scripts/latlong.html ทำไมกฎโคไซน์จึงเป็นที่นิยมมากกว่า? หมายเหตุ:ข้อความที่ยกมาได้รับการปรับปรุงโดยผู้เขียนดังกล่าวข้างต้น

41 distance  algorithm  spherical-geometry  haversine 

เป็นไปได้ไหมที่จะใช้ gdal API โทรgdal_translateจากรหัส Python? ฉันไม่ได้หมายถึงเพียงแค่เรียกใช้ gdal_translate.exe จากระบบไฟล์ แต่จะเรียกมันว่าอย่างใดในรหัสดังนั้นฉันไม่จำเป็นต้องรู้ไดเรกทอรีที่แน่นอนว่าไฟล์ gdal_translate นั้นอยู่ในอะไร

40 python  gdal 

ฉันมีรูปหลายเหลี่ยมหลายรูปตั้งอยู่นอกชายฝั่งรอบสหราชอาณาจักร ฉันกำลังพยายามคำนวณพื้นที่ของรูปหลายเหลี่ยมเหล่านี้ในหน่วยตารางกิโลเมตรโดยไม่แยกพวกมันด้วยเขต UTM เป็นไปได้ไหมที่จะทำสิ่งนี้ให้ถูกต้อง? ฉันใช้ World Mercator แต่มีการบิดเบือนขนาดใหญ่ในระดับนี้ ถ้ารูปหลายเหลี่ยมถูก จำกัด ที่ชายฝั่งตะวันออกบอกว่าฉันจะใช้ UTM31N อย่างไรก็ตามในขณะที่พวกเขาอยู่ทั่วน่านน้ำสหราชอาณาจักรฉันไม่แน่ใจว่าจะใช้อะไร มีระบบพิกัดที่ครอบคลุมพื้นที่นี้หรือไม่? กริดแห่งชาติของอังกฤษไม่ได้ขยายออกไปไกลพอที่จะใช้งานได้

40 coordinate-system  utm  united-kingdom 

ฉันสงสัยว่ามันเป็นไปได้ที่จะดูเนื้อหาของ shapefile โดยใช้ Python โดยไม่ต้องมีและสิทธิ์การใช้งาน ArcMap สถานการณ์คือคุณสามารถสร้าง shapefiles จากหลาย ๆ แอปพลิเคชั่นไม่เพียง แต่จากซอฟต์แวร์ ESRI ฉันต้องการสร้างสคริปต์ Python ที่ตรวจสอบการอ้างอิงเชิงพื้นที่ประเภทสถานที่ชื่อคุณลักษณะและคำจำกัดความและเนื้อหาของเขตข้อมูลในรูปแบบไฟล์และเปรียบเทียบกับชุดของค่าที่ยอมรับได้ ฉันต้องการให้สคริปต์นี้ทำงานแม้ว่าองค์กรจะไม่มีใบอนุญาต ESRI ในการทำสิ่งนี้คุณต้องใช้ ArcPy หรือคุณสามารถขุดลงในไฟล์รูปร่างโดยไม่ต้องใช้ ArcPy ได้หรือไม่?

40 python  shapefile 

ฉันได้รับไฟล์ TIFF จากใครบางคน แต่ฉันไม่แน่ใจว่ามีการอ้างอิงทางภูมิศาสตร์หรือไม่ มีวิธีง่าย ๆ ในการตรวจสอบเรื่องนี้หรือไม่?

40 geotiff-tiff  georeferencing 

วันนี้ฉันได้อ่านมากเกี่ยวกับอุปกรณ์ Android และ iphone ที่ใช้สำหรับการรวบรวมข้อมูลเชิงพื้นที่ในฟิลด์และดูเหมือนว่าจะเป็นอุปกรณ์วัตถุประสงค์เดียวเช่น Garmin eTrex ที่ได้รับความนิยมอย่างมากเป็นจุดจบ ทำไมต้องซื้ออะไรที่มีหน้าจอเล็ก ๆ ไม่สามารถเพิ่มอินเทอร์เฟซการรวบรวมข้อมูลแบบกำหนดเองแผนที่ฐานต้นทุนบวกและไม่มีการถ่ายภาพเมื่อคุณได้รับทุกสิ่งและอื่น ๆ ฟรีหรือราคาถูกมากด้วยแท็บเล็ต / โทรศัพท์ ใกล้เท่าที่ฉันสามารถบอกได้ถึงประโยชน์ที่แท้จริงเพียงอย่างเดียวจีพีเอสจริงมีมากกว่าอุปกรณ์อื่น ๆ เหล่านี้มีความทนทานและยืดอายุแบตเตอรี่ เว้นแต่จะมีความแม่นยำ ความแม่นยำของสมาร์ทโฟนหรือแท็บเล็ตที่เปิดใช้งาน GPS เปรียบเทียบกับ GPS ระดับ "ของจริง" สำหรับผู้บริโภคอย่างไร สิ่งที่เกี่ยวกับเมื่อกริด? สิ่งที่เกี่ยวกับเวลาเริ่มเย็นเพื่อล็อค?

40 gps  android  accuracy 

ฉันมีไฟล์ CSV ซึ่งมีข้อมูลรหัสไปรษณีย์ในเขตข้อมูล แต่ไม่ใช่เชิงพื้นที่เช่นไม่มีการ Easting และ Northing ฉันมี shapefile พร้อมรหัสไปรษณีย์ในฟิลด์และมีข้อมูลเพิ่มเติม ฉันต้องการเข้าร่วม CSV กับ shapefile ตามแอตทริบิวต์ (เขตข้อมูลรหัสไปรษณีย์) ฉันรู้ว่าฉันสามารถทำได้ใน ArcMap แต่ฉันจะบรรลุสิ่งเดียวกันใน QGIS ได้อย่างไร

40 qgis  csv  attribute-joins 

มีกฎที่ยากและรวดเร็วเกี่ยวกับวิธีการแก้ไขที่เหมาะกับข้อมูลแรสเตอร์แต่ละประเภทหรือไม่?

40 raster  interpolation  gis-principle  resampling