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

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

4
ปัญหาเกี่ยวกับค่า NA เมื่ออ่านไฟล์. DEM พร้อมแพ็คเกจ R 'raster' ใน Windows
ฉันพยายามอ่านไฟล์แรสเตอร์ในรูปแบบ. DEM บน windows โดยใช้แพ็คเกจ 'raster' ใน R ฉันพบปัญหากับค่า NA เมื่อโหลดข้อมูลลงใน R ใน Windows 7 แต่ฉันไม่มีปัญหาใน Mac ที่มี OSX Lion บน windows ค่า NA ดูเหมือนจะอ่านไม่ถูกต้อง คำถามคือทำไมสิ่งนี้เกิดขึ้น? ไฟล์แรสเตอร์ที่ใช้ถูกดาวน์โหลดจาก USGS ด้วยรหัส R ต่อไปนี้: download.file('http://edcftp.cr.usgs.gov/pub/data/gtopo30/global/e020n90.tar.gz', 'e020n90.tar.gz') untar('e020n90.tar.gz') จากนั้นฉันก็อ่านแรสเตอร์เป็น R โดยใช้แพ็คเกจ 'raster' ใน OSX Lion และ R64 เวอร์ชัน 2.13.1 ค่า NA จะได้รับการยอมรับ: > onMac …
10 raster  r  windows  osx 

3
ข้อมูล Shapefiles และ OpenStreetMap ไม่ทับซ้อนกันใน QGIS?
ดังนั้นฉันค่อนข้างใหม่กับ GIS และฉันกำลังทำงานในโครงการทำแผนที่โดยใช้ข้อมูลสำหรับแวนคูเวอร์แคนาดา ฉันต้องการรวมข้อมูล OpenStreetMap เข้ากับข้อมูลที่จัดทำโดยเมืองแวนคูเวอร์ แต่ข้อมูลเมืองไม่สอดคล้องกับข้อมูล OSM เมื่อฉันจับคู่กับข้อมูลเหล่านั้นใน QGIS ทำตามคำแนะนำของคำตอบของแผนที่ที่ฉาย แต่ควรไม่ทับซ้อนกันใน ArcMap หรือไม่? ฉันตรวจสอบแหล่งข้อมูลทั้งสองจากแหล่งที่สาม (แผนที่โลก) ข้อมูล OSM นั้นถูกต้อง - บริติชโคลัมเบียอยู่ในจุดที่เหมาะสม - แต่ข้อมูลของเมืองนั้นใส่แวนคูเวอร์ทางตอนใต้ของฝรั่งเศส แหล่งข้อมูลของฉันคือ: bikeways และห้องสมุดจากเมืองแวนคูเวอร์ (ทั้งคู่เป็นไฟล์ SHP) ข้อมูลเหล่านี้เป็นข้อมูลที่ไม่สอดคล้องกับไฟล์อื่น แผนที่บริติชโคลัมเบียจากแผนที่ถนนเปิด ฉันใส่ข้อมูลลงในฐานข้อมูลโดยใช้ osm2pgsql แผนที่โลกที่เป็นส่วนหนึ่งของโครงการ Open Street Style ประมาณการ: ข้อมูลเมืองจะถูกฉายใน UTM zone10, NAD83 (CSRS) ฉันเลือกสิ่งนี้ในหน้าจอคุณสมบัติของเลเยอร์ของ QGIS สำหรับเลเยอร์เหล่านั้น สตริงการฉาย:+proj=utm +zone=5 +ellps=GRS80 +datum=NAD83 +units=m …

1
วิธีการแสดงข้อมูลแรสเตอร์แบบไดนามิก (ไลดาร์) บนแผนที่เว็บ?
ฉันมีข้อมูล LiDAR คุณภาพสูง ฉันต้องการประมวลผลแบบไดนามิกและอัปเดตฟิลด์ (ไม่ใช่ฟิลด์ระดับความสูง) ขึ้นอยู่กับอินพุตและระดับความสูงของฉัน ณ จุดนั้น เมื่อเสร็จแล้วฉันต้องแสดงข้อมูลแรสเตอร์นี้บน Web Map (Bing Map, Google Map ฯลฯ ) ตามเวลาจริงที่สุด นี่คือสิ่งที่ฉันสามารถทำจนถึงตอนนี้: ฉันได้มากับสคริปต์บางอย่างโดยใช้ ArcGIS Python Lib ที่ประมวลผลคลิป lidar ของฉันเพิ่มฟิลด์ปรับปรุงฟิลด์ตามการป้อนข้อมูลมาดังนั้นตอนนี้ฉันมี คะแนนแรสเตอร์แบ่งออกเป็น 3 ประเภททุก ๆ 5 นาที ถ้าฉันดูแรสเตอร์นี้ใน arcGIS ฉันสามารถดูแรสเตอร์สี 3 สีที่ซ้อนกันบนแผนที่ฐาน แต่ตอนนี้ฉันต้องแสดงคลิปแรสเตอร์เหล่านี้บนแผนที่เว็บ ครั้งแรกของฉันแม้ว่าจะเป็นรูปหลายเหลี่ยมคลิปจัดเหล่านี้และแสดงรูปหลายเหลี่ยมบนแผนที่ JavaScript แต่การสร้างรูปหลายเหลี่ยมนั้นใช้เวลานาน ตอนนี้ตัวเลือกอื่นที่ฉันมีคือการใช้ความสามารถในการเรนเดอร์ Raster ของ GeoServer เพื่อแสดงภาพ จากนั้นใช้ WMS นี้ที่ฉันสร้างเพื่อซ้อนทับบน Bing …

2
จะละลายรูปหลายเหลี่ยมจาก shapefile โดยใช้เครื่องมือโอเพนซอร์สได้อย่างไร
ในโครงการปัจจุบันของฉันฉันเห็นภาพข้อมูลบนแผนที่แสดงค่าบางอย่างสำหรับแต่ละเขตการปกครองของประเทศ ฉันได้รับไฟล์รูปร่างของประเทศนี้ซึ่งฉันแปลงเป็น GeoJSON เพื่อใช้กับ PolyMaps ตอนนี้ฉันมีความต้องการในการจัดกลุ่มหลาย ๆ พื้นที่ให้ใหญ่ขึ้นเนื่องจากไม่มีข้อมูลสำหรับชิ้นส่วนเล็ก ๆ บางส่วน แต่สามารถรวบรวมข้อมูลสำหรับชิ้นส่วนขนาดใหญ่ได้ กลุ่มเหล่านั้นถูกกำหนดไว้ล่วงหน้า (เช่น A + B + C = D) ดังนั้นฉันคิดว่ามันจะง่ายกว่าในการสร้างรูปร่างไฟล์คลัสเตอร์และแปลงเป็น GeoJSON หลังจากนั้น แต่ฉันจะสร้างคลัสเตอร์ดังกล่าวได้อย่างไรโดยอัตโนมัติและมีเครื่องมือโอเพนซอร์ส? แน่นอนว่าวิธีอื่น ๆ ก็เหมาะสมเช่นกัน

2
วิธีการแปลง Distance, Azimuth, Dip ถึง XYZ?
ฉันมีสเปรดชีต Excel พร้อมข้อมูลส่วนหัวและแบบสำรวจ ข้อมูลส่วนหัวประกอบด้วย ID ของหลุมและพิกัดตำแหน่งและข้อมูลการสำรวจมีการสำรวจ downhole ที่เกี่ยวข้องกับค่า Distance, Azimuth และ Dip เนื่องจากฉันรู้ตำแหน่งของรูและระดับความสูงของพื้นผิวฉันจึงสามารถแปลงตารางสำรวจเป็นพิกัด XYZ ได้เช่นกัน ใครบ้างมีฟังก์ชั่น / ขั้นตอน / ตัวอย่าง? (VB & ArcObjects) ข้อมูลส่วนหัว: ข้อมูลการสำรวจ:

3
วิธีการจัดกลุ่มคะแนนตามความหนาแน่นของเคอร์เนล
ฉันมีชุดข้อมูลขนาดใหญ่ที่มี 36k คะแนนซึ่งแสดงถึงการใช้ประโยชน์ที่ดินเพื่อการพาณิชย์แต่ละแห่งมีเขตข้อมูลที่มีวิดีโอสแควร์ ฉันได้ทำการวิเคราะห์ความหนาแน่นของเคอร์เนลในชุดข้อมูลนี้สร้างภาพแรสเตอร์ที่แสดงความหนาแน่นของวิดีโอสแควร์เชิงพาณิชย์ในพื้นที่รถไฟใต้ดินทั้งหมด ฉันต้องแบ่งแรสเตอร์นี้ออกเป็นภูมิภาคที่สอดคล้องกับ maxima ท้องถิ่นซึ่งฉันเรียกว่า "ศูนย์กลาง" ฉันได้กำหนดที่ตั้งของศูนย์เรียบร้อยแล้วและตอนนี้ฉันต้องทำหนึ่งในสองสิ่งต่อไปนี้: ใช้เครื่องมือการทำคลัสเตอร์แบบจุดเช่น "การแบ่งพาร์ติชันรอบยาเม็ด" เพื่อจัดกลุ่มจุดให้เป็นกระจุกรอบ ๆ ศูนย์ที่ฉันระบุ ปัญหาของวิธีนี้ก็คือมันมีความเข้มข้นของการคำนวณและมากยิ่งขึ้นดังนั้นหากฉันพยายามใช้เมทริกซ์ที่แตกต่างกันเพื่อให้ได้คะแนนตามขนาด อย่างใดแบ่งความหนาแน่นเคอร์เนลแรสเตอร์ (ซึ่งประมาณคล้ายกับแรสเตอร์ภูมิประเทศ) เป็น "เนินเขา" แต่ละรอบแต่ละศูนย์ แต่ฉันไม่สามารถคิดถึงเครื่องมือใด ๆ สำหรับการทำเช่นนี้ ปัญหานี้ทำให้ฉันล้มเหลวในขณะที่และฉันหวังว่าฉันจะสามารถทำวิธีการจัดกลุ่มใน R แต่มันใช้เวลานานและฉันหมดเวลา ไม่มีใครรู้วิธีการง่ายๆในการแบ่งแรสเตอร์ความหนาแน่นเป็นละแวกใกล้เคียงของความเข้มหรือสำหรับการจัดกลุ่มชุดข้อมูลขนาดใหญ่ได้อย่างรวดเร็ว?

2
วิธีติดฉลากเซลล์แรสเตอร์ใน QGIS?
ฉันพยายามแสดงไฟล์แรสเตอร์ความละเอียดขนาดเล็กในเครื่องมือ QGIS Print Manager ฉันต้องการแสดงป้ายกำกับขนาดเล็กที่มีค่าแต่ละเซลล์ (ตำแหน่ง - ประมาณกึ่งกลางของแต่ละเซลล์) ฉันจะทำมันได้อย่างไร มีปลั๊กอินใดบ้างที่สามารถช่วยให้ฉันทำงานนี้ได้สำเร็จ ความคิดของฉันคือการใช้ gdal2xyz โหลดไฟล์ดังกล่าวโดยใช้ปลั๊กอิน CSV แล้วบันทึกเป็นไฟล์เวกเตอร์ น่าเสียดายที่มันใช้งานได้ดีมากเพราะฉันมีไฟล์แรสเตอร์มากมาย อีกแนวคิดหนึ่งคือการใช้ gdal2xyz.py แล้ว ogr2ogr เพื่อบันทึกเป็นไฟล์เวกเตอร์ มีวิธีการชุดกระบวนการนี้หรือไม่ มีโปรแกรมโอเพนซอร์ซใดที่ฉันสามารถทำได้ง่ายกว่านี้อีกไหม?
10 qgis  raster  gdal  ogr 

3
วงกลมจริงปรากฏในระบบพิกัดใด
ฉันได้รับการพัฒนาใน Google Maps สำหรับชั่วขณะและกำลังเปลี่ยนเป็น OpenLayers แต่ติดอยู่ที่จุดหนึ่ง: ฉันกำลังวาดเวกเตอร์ที่ส่วนใหญ่เป็นวงกลม แต่ยังเกี่ยวข้องกับเรขาคณิตที่กำหนดเองดังนั้นฉันจึงไม่สามารถใช้ createRegularPolygon ฉันกำลังสร้างอาร์เรย์ของคะแนนและส่งไปยังวัตถุ OpenLayers.Layer.Vector วิธีนี้ใช้งานได้ดี แต่ 'วงกลม' ที่ได้นั้นเป็นเพียงรอบที่อยู่ใกล้เส้นศูนย์สูตรไม่ใช่ที่ละติจูดที่สูงกว่า พวกมันถูกบีบอัดในปริมาณที่เท่ากันกับแผนที่ฐานมีละติจูดที่สูงกว่าซึ่งดูเหมือนจะสมเหตุสมผล ถ้าฉันใช้ Google เป็นเลเยอร์พื้นฐานแล้วพวกเขาจะมีรูปร่างที่ถูกต้องในทุกละติจูด แต่ฉันกำลังสร้างโซลูชันมือถือออฟไลน์เพื่อที่จะไม่ทำเช่นนั้น ฉันจะใช้อะไรกับเลเยอร์ฐาน (หรือพารามิเตอร์ที่จะตั้งค่า) เพื่อให้เรขาคณิตของฉันยังคงอยู่ในทุกละติจูดโดยไม่ต้องใช้ Google หรือแผนที่ออนไลน์

2
Mapnik สามารถแสดงผลข้อมูลเชิงพื้นที่จาก SQL Server ได้หรือไม่
นี่เป็นคำถามง่าย ๆ - ฉันแค่อยากรู้ว่าใครสามารถยืนยันได้ว่าพวกเขาใช้ Mapnik ในการเรนเดอร์ไทล์โดยตรงจากข้อมูลที่เก็บไว้ใน SQL Server (และถ้าเป็นเช่นนั้นคุณทำอะไรเพื่อให้มันใช้งานได้) โปรดอย่าตอบกลับจนกว่าคุณจะได้รับการยืนยันเป็นการส่วนตัวว่าเป็นไปได้ ... ไม่ใช่แค่ว่ามันควรจะเป็นไปได้ ฉันเคยเห็นผู้คนแนะนำ Mapnik สำหรับการแสดงข้อมูลจาก SQL Server ในเธรดเช่นนี้: เอ็นจินการเรนเดอร์แผนที่ที่ดีสำหรับ SqlGeometry / SqlGeography นอกจากนี้เว็บไซต์ mapnik ระบุว่า "รองรับรูปแบบ OGR ทั้งหมด": http://mapnik.org/faq/#data อย่างไรก็ตามฉันไม่สามารถใช้การเชื่อมต่อกับฐานข้อมูล SQL Server ในพื้นที่ของฉันเป็นแหล่งข้อมูลสำหรับเลเยอร์ Mapnik ได้ เป็นการทดสอบฉันใช้ OGR VRT ทดสอบตรงไปตรงมาซึ่งดำเนินการคำสั่ง SQL บน SQL Server เพื่อเลือกรูปหลายเหลี่ยมแบบอินไลน์ภูมิศาสตร์เดียวกับข้อความที่รู้จักกันดีดังต่อไปนี้: <OGRVRTDataSource> <OGRVRTLayer name="test"> <SrcDataSource>MSSQL:server=zangief\DENALICTP3;database=TempDB;</SrcDataSource> <SrcSQL>SELECT geography::STPolyFromText('POLYGON((0 …

3
กรอก Donut Polygons โดยใช้ ArcGIS สำหรับเดสก์ท็อปหรือไม่
ฉันมีชุดรูปหลายเหลี่ยมที่แสดงถึงแหล่งน้ำ ในลักษณะเดียวกันนี้มีรูปหลายเหลี่ยมที่เป็นเกาะในแหล่งน้ำ ดังนั้นสำหรับทะเลสาบที่กำหนดมีรูปหลายเหลี่ยมขนาดเล็กกว่าที่เป็นตัวแทนของเกาะในลักษณะเดียวกัน อย่างไรก็ตามรูปหลายเหลี่ยมน้ำที่เกิดขึ้นจริงสำหรับทะเลสาบนั้นมี "หลุม" ที่มีเกาะอยู่ ดังนั้นสำหรับทะเลสาบที่มีเกาะหนึ่งเกาะอยู่ตรงกลางรูปหลายเหลี่ยมน้ำจะดูเหมือนโดนัท ตอนนี้ฉันต้องการสร้างชั้นของรูปหลายเหลี่ยมน้ำที่เต็มไปด้วยหลุมดังนั้นเมื่อฉันต้องการดูทะเลสาบทั้งหมดในแผนที่ฉันสามารถเห็นรูปหลายเหลี่ยมน้ำที่ไม่มีรู หมายเหตุ: เกาะอาจมีแหล่งน้ำอยู่ภายในตัวเอง อาจมีแหล่งน้ำซ้อนกันภายในเกาะซึ่งเป็นส่วนหนึ่งของแหล่งน้ำที่เป็นส่วนหนึ่งของเกาะอื่น ชนิดของการลงทะเบียนน้ำ / เกาะ :-) ฉันจะทำสิ่งนี้ใน ArcGIS for Desktop ได้อย่างไร

2
คุณลักษณะ DXF เป็นค่า z
เราได้รับ DXF ที่มีแอตทริบิวต์ Elevation แต่เราต้องการสร้าง 3D DXF ฉันรู้ว่าฉันสามารถสร้าง 3D dxf จาก shapefile โดยใช้: ogr2ogr -f "DXF" out.dxf in.shp -zfield Elev แต่ไดรเวอร์ ogr ดูเหมือนจะไม่รับแอตทริบิวต์เป็นเขตข้อมูลใน dxf ไม่มีใครรู้ว่าสิ่งนี้เป็นไปได้ ฉันต้องการโซลูชันฟรีหากเป็นไปได้ เรามี FME แต่ฉันกำลังมองหาตัวเลือกอื่น ๆ
10 3d  ogr  fme  dxf 

1
ประสิทธิภาพที่เพิ่มขึ้นจากดัชนี GIST สำหรับจุดในรูปหลายเหลี่ยมแบบสอบถาม
ฉันมีสองตาราง: สถานที่ตั้ง (id, region_id, the_geom) และภูมิภาค (id, the_geom) สำหรับจุดที่ตั้งแต่ละแห่งฉันต้องการระบุภูมิภาคที่อยู่ใน: UPDATE locations SET region_id = (SELECT id FROM regions WHERE ST_Within(locations.the_geom,regions.the_geom) ); มันสมเหตุสมผลหรือไม่ที่จะสร้างดัชนี GIST ที่จุดตำแหน่ง? ฉันจะสร้างดัชนีในรูปหลายเหลี่ยมภูมิภาค แต่ฉันไม่แน่ใจเกี่ยวกับประเด็นนี้ มันจะเพิ่มความเร็วของแบบสอบถามหรือไม่

4
SHAPE_Length คำนึงถึงค่าระดับความสูงของบัญชีหรือไม่
ที่ไม่Shape_Lengthแอตทริบิวต์คำนึงถึงระดับความสูง? ตัวอย่างเช่นถ้าฉันมี 3 จุด: A, BและC A.X = 10 A.Y = 10 A.Z = 10 B.X = 10 B.Y = 10 B.Z = 50 C.X = 0 C.Y = 0 C.Z = 10 แล้วฉันจะสร้างสองบรรทัด: และL1L2 L1.fromPoint = A L1.toPoint = C L2.fromPoint = B L2.toPoint = C สิ่งที่ฉันสังเกตในผลก็คือว่าSHAPE_Lengthแอตทริบิวต์จะเหมือนกันสำหรับทั้งสองและL1 L2ฉันได้ตรวจสอบแล้วว่าคะแนนและ polylines …

2
การหาระยะห่างระหว่างจุดตามแนวแม่น้ำ (การอ้างอิงเชิงเส้น)
ฉันใหม่กับระบบสารสนเทศภูมิศาสตร์และหลง ฉันมีเครือข่ายของแม่น้ำซึ่งประกอบด้วย linestrings ขนาดเล็ก (polylines) ใน Shapefile ที่แตกต่างกันฉันมีที่ตั้งของสถานีหลายแห่งในแม่น้ำ ฉันต้องหาระยะทางตามแม่น้ำจากแต่ละสถานีจนถึงจุดสิ้นสุดของแม่น้ำ ปัญหาคือว่าฉันไม่มีแม่น้ำที่สมบูรณ์หนึ่งแห่ง และไม่มีคุณสมบัติทั่วไปที่ฉันสามารถใช้ในการรวมหรือรวม linestrings นอกจากนี้บางสถานีอยู่บนแคว ดังนั้นสิ่งที่สำคัญที่สุดคือฉันต้องสร้างเครือข่ายแม่น้ำขนาดใหญ่รวมถึงสาขาและทำการอ้างอิงสายการบิน แต่ฉันไม่รู้จะทำอย่างไร ความช่วยเหลือใด ๆ ที่ชื่นชมอย่างลึกซึ้ง โปรดอธิบายให้ชัดเจนที่สุด โปรดแสดงความคิดเห็นหากคุณต้องการข้อมูลเพิ่มเติม

4
เป็นการแนะนำที่ดีสำหรับ INSPIRE?
ทุกคนสามารถแนะนำให้ฉันรู้จักกับคำแนะนำ INSPIRE ได้ดีหรือไม่? ฉันรู้ว่ามันเป็น WMS 1.3 พร้อมกับเสียงระฆังและนกหวีดที่เพิ่มเข้ามา แต่ตอนนี้ฉันกลับมาที่ยุโรปฉันน่าจะรู้ว่ามันคืออะไร

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