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

ถามตอบสำหรับนักทำแผนที่นักภูมิศาสตร์และผู้เชี่ยวชาญด้าน GIS

2
วิธีการรับพิกัดของจุดในแบบโต้ตอบได้ QGIS
ฉันมีรูปร่างไฟล์ที่มีเลเยอร์พิมพ์จุดเดียว มันอยู่ใน Google CRS (EPSG: 900913) เมื่อฉันใช้เครื่องมือ 'ระบุคุณสมบัติ' มันให้ข้อมูลมากมายเกี่ยวกับจุดที่ฉันคลิก แต่ฉันไม่เห็นพิกัด มีวิธีหา lat / lng ของจุดแบบโต้ตอบหรือไม่?
12 qgis 

2
การลบแถวใน shapefile ใน R
ฉันได้นำเข้า shapefile ลงใน R และเข้าร่วมกับตาราง Shapefile ของฉันมีรหัสสำรวจสำมะโนประชากรทั้งหมดในขณะที่ตารางของฉันมีรหัสสำรวจสำมะโนประชากรที่เลือกไว้เท่านั้น ตอนนี้ฉันกำลังพยายามลบแถวทั้งหมดไม่ได้รับการจับคู่ นี่คือชุดข้อมูลของฉันดูเหมือน (ฉันพยายามลบแถวทั้งหมดด้วย NA ดังนั้นทั้งสองจะต้องถูกลบออก): CTUID Cluster Average 5350007.01 1 124.53 5350007.02 1 234.87 5350010.01 4 110.11 5350010.02 5 187.68 5350001 NA NA 5350002 NA NA ฉันพยายามใช้รหัสบรรทัดนี้: shape2[!(rowSums(is.na(shape2))==NCOL(shape2)),] ซึ่งทำให้ฉันข้อผิดพลาดนี้: Error in rowSums(is.na(shape2)) : 'x' must be an array of at least two dimensions …

1
เส้นวงกลมขนาดใหญ่ในการฉายภาพ Equirectangular
เพียงเพื่อตรวจสอบว่าฉันอยู่ในเส้นทางที่ถูกต้อง: วงกลมใหญ่ทั้งหมดอยู่ในทรงกลมและในการฉายภาพทางอ้อม (เช่นละติจูด, ลองจิจูดคู่): เส้นเมอริเดียน (เช่นไปขั้วโลก) ของแบบฟอร์ม tan latitude = sin360(longitude + rotation) * amplitude + offset (ด้วยข้อ จำกัด เพิ่มเติมเกี่ยวกับการรวมกันของออฟเซ็ต / แอมพลิจูด - เห็นได้ชัดว่าเส้นทางวงกลมที่ยิ่งใหญ่ทั้งหมดที่มีแอมพลิเคชั่น 0 ยังมีออฟเซต 0 - เส้นศูนย์สูตรด้วย) หรือมีเส้นทางวงกลมขนาดใหญ่ที่ไม่เหมาะกับรูปแบบนี้ (อีกครั้งเฉพาะในระบบพิกัดละติจูด - ลองจิจูดไม่ใช่บนแผนที่อื่น ๆ ) หมายเหตุ: ฉันได้เพิ่มtanข้างต้นหลังจากโพสต์คำถามในการตอบสนองต่อการตอบที่ดีเยี่ยม whubers ปรากฎว่าoffsetเป็น 0 เสมอ

3
ตัวอย่างอุตสาหกรรมที่ใช้ ArcPy ใน Python สำหรับการประมวลผลทางภูมิศาสตร์?
คำตอบของคำถามนี้เป็นความพยายามของชุมชน แก้ไขคำตอบที่มีอยู่เพื่อปรับปรุงโพสต์นี้ ขณะนี้ไม่ยอมรับคำตอบหรือการโต้ตอบใหม่ เมื่อเร็ว ๆ นี้นักเรียนบางคนถูกถามว่าประเภทของการดำเนินการ GIS Analyst / Developer มักจะพยายามทำให้เป็นอัตโนมัติโดยใช้ Python โดยใช้กระบวนการทางภูมิศาสตร์ในแพ็คเกจเว็บไซต์ ArcGIS และ ArcPy ฉันคิดว่ามันอาจมีประโยชน์ที่จะรู้เมื่อค้นหาแบบฝึกหัดบางอย่างเพื่อให้แน่ใจว่าสิ่งที่คุณสร้างขึ้นมีความเกี่ยวข้องกับอุตสาหกรรมและอาจถูกนำมาใช้ซ้ำในภายหลังในที่ทำงาน คำตอบที่ง่ายที่สุดคือ "อ่านความช่วยเหลือของ Esri และดูตัวอย่าง" แต่ฉันกำลังมองหาสถานการณ์ที่เฉพาะเจาะจงซึ่งเป็นเรื่องธรรมดาที่จะนำไปใช้ ดังนั้นการแบ่งปันเวิร์กโฟลว์จะเกิดขึ้นในฐานะ "เราได้รับไฟล์. zip พร้อมกับ shapefiles เราใช้ Python เพื่อแกะมันออกมาฉายพวกมันทั้งหมดไปยังระบบพิกัด X โหลดลงในฐานข้อมูลทางภูมิศาสตร์ของ ArcSDE ให้สิทธิ์ผู้ใช้เข้าถึงข้อมูลเหล่านี้" โปรดรู้สึกสะดวกสบายที่จะให้คำอธิบายที่กระชับของเวิร์กโฟลว์โดยไม่ต้องมีรายละเอียดมากนัก
12 arcpy  references 

3
อ่านไฟล์ Geodatabase โดยใช้ R?
ไฟล์ฐานข้อมูลภูมิศาสตร์ (fgdb) รวมถึงตารางฐานข้อมูลไฟล์จำนวนมาก เท่าที่ผมรู้ว่าพวกเขามีอยู่เป็นdbfไฟล์ Database.gdbแต่อยู่ภายใน ใน ArcCatalog เส้นทางของไฟล์มีลักษณะคล้ายC:\Users\...\Database.gdb\Stats_AAกัน จะอ่านdbfไฟล์เหล่านี้ทั้งหมดได้อย่างไรR(ซอฟต์แวร์ทางสถิติ) เส้นทางที่ถูกต้องในการจัดหาคืออะไร? ฟังก์ชั่นที่ใช้คือread.dbf(ในแพ็คเกจต่างประเทศ) สายพันธุ์ของ test<-read.dbf(file="C:/Users..Database.gdb/Stats_AA.dbf") และ test<-read.dbf(file="C:/Users..Database/Stats_AA.dbf") ไม่ทำงาน "รูปแบบ" ที่ถูกต้องของชื่อไฟล์ที่จะใช้คืออะไรหรือฉันจำเป็นต้องส่งออกตารางฐานข้อมูลไฟล์ทั้งหมดไปยังรูปแบบหรือตำแหน่งอื่นหรือไม่?

5
GML, KML, GeoJSON - การเรนเดอร์ 3109 รูปหลายเหลี่ยม?
ฉันกำลังทำงานกับ Geoserver ซึ่งให้บริการเขตลดระดับ 48 เขตให้แก่ openlayers (3109 รูปหลายเหลี่ยม - จุดยอดอีกมากมาย) มณฑลถูกโหลดลงในฐานข้อมูล postgis ฉันอยากรู้เกี่ยวกับประสบการณ์ของนักพัฒนาซอฟต์แวร์เมื่อพยายามผลักดันจำนวนจุดยอดนั้นไปยังลูกค้า รูปแบบ WFS ใดที่คุณได้รับผลลัพธ์ที่ดีที่สุด มีการใช้การปรับแต่ง Geoserver เพิ่มเติมหรือไม่ ฉันรู้ว่า WMS แบบเรียงต่อกันจะเร็วขึ้น แต่ฉันต้องการอนุญาตให้มีการเปลี่ยนแปลงแบบไดนามิกในแผนที่นักร้องโดยใช้ openLayers เช่น ผู้ใช้ส่งแบบฟอร์มสคริปต์ Python ถูกเรียกใช้และถังขยะข้อมูลใหม่จะถูกส่งคืนเพื่อ openlayers เพื่อโหลด div แผนที่ใหม่ ฉันต้องการลองแบบเต็มความละเอียดก่อนที่จะลดความซับซ้อนของรูปหลายเหลี่ยมใน openlayers
12 geoserver  kml  wfs  geojson  gml 

2
บัฟเฟอร์ในหน่วยเมตร / กม. ​​โดยใช้ชั้น WGS84 โดยใช้ QGIS
ฉันใช้ QGIS 1.7.4-Wroclaw และทำงานกับเลเยอร์ SHP ใน WGS84 ฉันต้องการบัฟเฟอร์จำนวนกิโลเมตรจากเลเยอร์นี้ด้วยเครื่องมือ ftools buffer ฉันเข้าใจว่าเครื่องมือบัฟเฟอร์ใช้หน่วยชั้นเสมอซึ่งสำหรับ WGS84 เป็นองศาทศนิยม เราทุกคนรู้ดีว่าองศาไม่ได้แปลงเป็นเมตรอย่างสม่ำเสมอดังนั้นฉันควรทำบัฟเฟอร์อย่างไร จำเป็นหรือไม่ที่ฉันจะต้องเปลี่ยนรูปร่างของไฟล์ให้เป็น CRS อื่นที่ใช้หน่วยวัดเป็นเกณฑ์? ถ้าเป็นเช่นนั้นฉันจะเลือกได้อย่างไร
12 qgis  buffer  units 

5
การจับคู่แผนที่เวลาจริงโดยใช้ ArcGIS Desktop?
ฉันกำลังทำงานในโครงการที่ต้องการจับคู่แผนที่ตามเวลาจริง ฉันกำลังใช้ ArcGIS อยู่ แต่ก็ไม่แน่ใจว่ามันจะใช้งานได้หรือไม่ นี่คือสิ่งที่เกิดขึ้นกับระบบของฉัน: โทรศัพท์มือถือส่งสัญญาณ GPS แก้ไขทุก ๆ ห้าวินาทีไปยังเซิร์ฟเวอร์ของฉันและฉันต้องทำการจับคู่แผนที่แบบเรียลไทม์เพื่อสแน็ปอินจุดหนึ่งไปยังลิงก์ใดลิงก์หนึ่งของแผนที่และเก็บคะแนนที่ตรงกัน ArcGIS มีเครื่องมือ 'ใกล้' เพื่อหักล้างจุดเชื่อมโยงที่ใกล้ที่สุด อย่างไรก็ตามอินพุตไปยังเครื่องมือต้องเป็นรูปร่างไฟล์สำหรับจุด ฉันสงสัยว่าฉันสามารถจัดการกับข้อมูลเรียลไทม์ที่ฉันได้รับใน ArcGIS ได้อย่างไร การสร้างไฟล์รูปร่างสำหรับแต่ละจุด? มีข้อเสนอแนะอื่น ๆ เพื่อให้งานนี้สำเร็จหรือไม่

3
วิธีจัดศูนย์กลางแผนที่ใน Print Composer ให้อยู่กึ่งกลางเช่นเดียวกับในหน้าต่างแผนที่
ฉันไม่สามารถคิดออกว่าฉันจะอยู่ตรงกลางแผนที่ในพิมพ์นักแต่งเพลงชั้นวัตถุ / ศูนย์ / แผนที่ในลักษณะเดียวกับในหน้าต่างแผนที่ในQGIS 1.8 สมมติว่าฉันมีเลเยอร์ที่ฉันต้องการอยู่ตรงกลางของแผนที่ ในหน้าต่างแผนที่มันค่อนข้างง่าย - ฉันแค่คลิกที่ " Zoom to layer scope " น่าเสียดายที่เมื่อฉันสร้างกรอบแผนที่ในผู้แต่งแผนที่ชั้นจะอยู่กึ่งกลางที่ความกว้าง (ดี) แต่ไม่อยู่กึ่งกลางที่ความสูง ถ้าฉันใช้ปุ่มตั้งค่าเพื่อแม็พขอบเขตของผืนผ้าใบในคุณสมบัตินักแต่งเพลง / รายการ / ขอบเขตการพิมพ์จะอยู่กึ่งกลางความกว้างอีกครั้ง ฉันคิดว่าในทั้งสองกรณีมันจะดีกว่าถ้าจะให้อยู่ตรงกลาง มีวิธีใดที่จะให้ศูนย์กลางของแผนที่ใน Print Composer ไปยังศูนย์กลางเดียวกันกับแผนที่ในหน้าต่างแผนที่หรือไม่?

3
จะกำหนดวันที่ของภาพถ่าย Bing สำหรับพื้นที่ได้อย่างไร
เป็นไปได้หรือไม่ที่จะตัดสินว่าปัจจุบันภาพของ Bing ในบางพื้นที่เป็นอย่างไร ฉันตรวจสอบบล็อกของชุมชนและค้นหาผ่านการอัปเดตรายเดือนแล้ว แต่ฉันไม่พบพื้นที่ศึกษาของฉัน ฉันพบโพสต์ที่พวกเขากล่าวถึงการอัปเดตภาพสำหรับทั้งประเทศ แต่ไม่มีการเจาะจง (เป็นเรื่องจริงทั่วทั้งประเทศหรือไม่หรือแค่ส่วนใหญ่?) ในเครดิตในLayout Viewใน ArcMap 10.1 ผมสังเกตเห็นว่าในขณะที่ไมโครซอฟท์คอร์ปอเรชั่น2012, 2010สินเชื่ออื่น นี่บ่งบอกว่าภาพนั้นมาจากปีนั้นหรือไม่?

3
ชุดข้อมูล Gdal.ReadAsArray () ขัดข้อง Python
ฉันใช้ Python 2.6.5 (32 บิต) กับ Numpy 1.3 และ Gdal 1.9.1 ติดตั้งบน Windows 7 64 บิต ฉันพยายามอ่านชุดข้อมูลแรสเตอร์ 800 MB Imagine (.img) ในอาร์เรย์ Numpy เพื่อทำพีชคณิตแบบแรสเตอร์ แต่ทันทีที่ฉันเรียกใช้รหัสต่อไปนี้ Python.exe จึงล้มเหลว from osgeo import gdal g = gdal.Open(r'path\to\dataset', gdal.GA_Readonly) b = g.GetRasterBand(1) data = b.ReadAsArray() Python.exe ขัดข้องขณะb.ReadAsArray()โทร ฉันได้ทำการค้นหาโดย Google และพบบทความจาก Gdal 1.6 ที่กล่าวถึงปัญหานี้กับ Windows …
12 gdal  numpy  python-2.6  array 

1
Arcgis ใช้ดัชนีเชิงพื้นที่แบบใด
หัวข้อการวิจัยของอาจารย์คนแรกของฉันเกี่ยวกับดัชนีเชิงพื้นที่ ฉันเรียนรู้เกี่ยวกับ R-Trees และ R + -tree, R * -trees และหัวข้ออื่น ๆ ที่คล้ายคลึงกัน แต่มีความคล้ายคลึงกันมาก ข้อมูลนี้อาจไม่ได้รับการเผยแพร่ แต่ ESRI ใช้ดัชนีเชิงพื้นที่ชนิดใด คำถามประเภทนี้สามารถเปลี่ยนเป็นคำถามรองได้อย่างไร มีบางสถานการณ์ที่ดัชนีประเภทอื่นจะเพิ่มประสิทธิภาพ ในกรณีของ Oracle สามารถใช้ประเภทข้อมูลเชิงพื้นที่ ST_Geometry ซึ่งสามารถจัดทำดัชนีเป็น R-tree หรือ quad tree หรือทั้งสองอย่าง: http://docs.oracle.com/html/A88805_01/sdo_intr.htm ฐานข้อมูลอื่น ๆ ใช้กลไกการจัดเก็บแบบไบนารีที่ไม่ใช่เชิงพื้นที่ตามค่าเริ่มต้นและดูเหมือนจะไม่ใช้การจัดทำดัชนีเชิงพื้นที่ที่มีอยู่ ใน SQL 2005 และประเภทข้อมูลเชิงพื้นที่ที่ต่ำกว่าไม่มีอยู่ แต่ ESRI ยังคงสามารถใช้ฐานข้อมูล SQL ได้ ตามกรณีศึกษา SQL 2012 นี้: http://www.microsoft.com/casestudies/Case_Study_Detail.aspx?CaseStudyID=710000000547มันใช้ประเภทประสิทธิภาพเชิงพื้นที่อัตโนมัติ แต่ไม่ได้ระบุประเภท R-Tree …

2
ตัวเลือกการแก้ไข / การวิเคราะห์สมการของ PostGIS
ฉันกำลังมองหาเวิร์กโฟลว์โอเพนซอร์ซเพื่อทำให้การสร้าง DEM เป็นแบบอัตโนมัติ เรามีชุดของเว็บไซต์ที่อยู่ในชุดข้อมูล LIDAR เราต้องการสร้าง DEM เฉพาะไซต์สำหรับแต่ละไซต์และเรากำลังมองหากระบวนการอัตโนมัติ จนถึงตอนนี้เรามีระบบอัตโนมัติ: การโหลดข้อมูลลงใน PostGIS (LIDAR และจุดไซต์) การสร้างขอบเขตไซต์ (การรวมกันของ st_buffer และ st_envelope) ตอนนี้เรากำลังมองหาตัวเลือกในการแก้ไขชุดย่อยข้อมูลในแต่ละไซต์และส่งออกไปยังไฟล์กริด Surfer7 ขณะนี้เรากำลังค้นหาข้อมูลในฐานข้อมูล PostGIS โดยใช้ QuantumGIS ส่งออกไปยังไฟล์ csv จากนั้นนำเข้าและกรอข้อมูลด้วยตนเองใน Surfer7 หวังว่าเราจะทำสิ่งนี้โดยอัตโนมัติเช่นกัน ดังนั้นสำหรับกรณีการใช้งานนี้เราต้องการที่จะวิเคราะห์ข้อมูล lidar ของเราเป็น DEM กรณีที่เป็นไปได้อื่น ๆ ที่เราสามารถนึกถึงเกี่ยวข้องกับวิธีการแก้ไขที่แตกต่างกัน - ดังนั้นหากมีตัวเลือกสำหรับระยะทางผกผันและ kriging - เราสนใจ! นี่เป็นแบบฝึกหัดการเรียนรู้สำหรับเรา - เรากำลังทำงานในขั้นตอนของทารก!


3
GDAL และ Python: จะรับพิกัดสำหรับเซลล์ทั้งหมดที่มีค่าเฉพาะได้อย่างไร
ฉันมีArc / Info Binary Grid --- โดยเฉพาะแรสเตอร์การไหลของ ArcGIS --- และฉันต้องการระบุเซลล์ทั้งหมดที่มีค่าเฉพาะ (หรืออยู่ในช่วงของค่า) ในที่สุดฉันต้องการไฟล์รูปร่างของคะแนนที่แสดงถึงเซลล์เหล่านี้ ฉันสามารถใช้ QGIS เพื่อเปิด hdr.adf และรับผลลัพธ์นี้เวิร์กโฟลว์คือ: QGIS> เมนู Raster> Raster Calculator (ทำเครื่องหมายคะแนนทั้งหมดด้วยมูลค่าเป้าหมาย) QGIS> เมนู Raster> รูปหลายเหลี่ยม QGIS> เมนูเวกเตอร์> เมนูย่อยเรขาคณิต> รูปหลายเหลี่ยม centroids แก้ไขเซนทรอยด์เพื่อลบโพลีเซนทรอยด์ที่ไม่ต้องการ (เหล่านั้น = 0) วิธีการนี้ "ทำงานได้" แต่มันไม่ได้ดึงดูดความสนใจฉันเพราะมันสร้างไฟล์ 2 ไฟล์ที่ฉันต้องลบแล้วฉันต้องลบบันทึกที่ไม่ต้องการออกจาก shapefile ของ centroids (นั่นคือ = 0) คำถามที่มีอยู่วิธีการเรื่องนี้ แต่มันเหมาะสำหรับ ArcGIS …
12 python  raster  gdal 

โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.