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

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

2
อัปเดตแถวฐานข้อมูลด้วยคะแนนที่อยู่ในรูปหลายเหลี่ยม
ฉันมี PostGIS / Postgresql DB ที่มีสองตารางอยู่ในนั้น หนึ่งที่มีรูปทรงเรขาคณิตจุดและอื่น ๆ ที่แสดงถึงขอบเขตของประเทศเป็นรูปหลายเหลี่ยม ฉันต้องการเพิ่มชื่อประเทศที่แต่ละจุดตัดกับแต่ละแถวในตารางพอยต์ของฉัน อาจเป็นแบบสอบถามแบบใช้ปรับปรุงที่มีขนาดใหญ่ ฉันคิดว่านี่เป็นไปได้ที่จะใช้ SQL แบบตรง แต่ไม่รู้ว่าจะเริ่มต้นตรงไหน คำแนะนำใด ๆ เกี่ยวกับเรื่องนี้จะได้รับการชื่นชมอย่างมาก ...

3
ทราบแหล่งที่มาของปริมาณน้ำฝนตลอด 24 ชั่วโมงที่ผ่านมาสำหรับสถานที่ในสหรัฐอเมริกาหรือไม่
ฉันกำลังมองหาข้อมูลทางประวัติศาสตร์เกี่ยวกับขนาดของปริมาณน้ำฝนสำหรับสถานที่ที่ระบุ (Boston, Massachusetts USA) ตัวอย่างเช่นในการสร้างฮิสโตแกรมของเหตุการณ์ปริมาณน้ำฝนตามปริมาณการตกตะกอนสำหรับสถานที่ที่กำหนด จนถึงตอนนี้ฉันเพิ่งพบชุดข้อมูลประเภท "3 วันล่าสุด" หรือ "เดือนที่แล้ว" จาก NOAA และอื่น ๆ ฉันต้องการข้อมูลมูลค่าหลายปีถ้าเป็นไปได้ ฉันไม่ได้มองหาไฟล์รูปร่างเช่นเดียวกับข้อมูลประเภทตาราง (ฐานข้อมูล, Excel, CSV และอื่น ๆ ) ไม่จำเป็นต้องมีข้อมูลสำหรับ Lat / Lon ใด ๆ โดยพลการ แต่อาจมีไว้สำหรับ POI ใกล้เคียงเช่นสนามบิน

2
แผนที่สแกนต้องได้รับการอ้างอิงทางภูมิศาสตร์ในการฉายภาพต้นฉบับหรือไม่
ฉันมีแผนที่กระดาษประวัติศาสตร์สำหรับอนุทวีปอินเดียซึ่งถูกสแกน แผนที่อยู่ในการฉายภาพ Albers Equal Area Conical และมีเส้นรุ้งและเส้นแวงบน ฉันพยายามหาตำแหน่งทางภูมิศาสตร์ในการคาดการณ์ของแหล่งที่มา แต่ฉันไม่ได้มีพารามิเตอร์ที่จำเป็นทั้งหมดสำหรับการประมาณการนี้ มีเพียงสองมาตรฐานแนวเท่านั้นที่จะถูกบันทึกไว้บนแผนที่ ฉันไม่ทราบว่า Central Meridian หรือ Latitude ของ Origin หรือ datum ของเรื่องนั้น สิ่งที่ควรทำในสถานการณ์เช่นนี้? เนื่องจากนี่เป็นแผนที่ที่มีระยะทาง 4,000km ถึง 4000km ฉันไม่คิดว่าเป็นไปได้ที่จะได้รับผลการพิมพ์ที่มีความแม่นยำสูง ฉันคิดว่าการอ้างอิงทางภูมิศาสตร์ในระบบพิกัดภูมิศาสตร์ wgs84 นี่เป็นความคิดที่ดีหรือไม่? หรือมีอย่างอื่นที่คุณจะแนะนำ?

1
การปรับคุณสมบัติแรสเตอร์ให้เรียบเป็นรูปหลายเหลี่ยม
ฉันมีข้อมูลแบบแรสเตอร์ที่แสดงถึงพื้นที่ในเมืองซึ่งฉันต้องการติดตามและแปลงเป็นรูปหลายเหลี่ยม การใช้เครื่องมือ Raster to Polygon ของ ArcMap ทำให้ฉันสามารถเข้าใกล้สิ่งนั้นได้ Unsimplified: ย่อ: อย่างไรก็ตามคุณจะเห็นว่ามีเส้นทแยงมุมหลายเส้น (ถนน) ซึ่งสร้างโซ่ของรูปหลายเหลี่ยมที่แบ่งมุม อัพรูปหลายเหลี่ยมของฉันนับอย่างมีนัยสำคัญ - ชุดข้อมูลแรสเตอร์เหล่านี้กำลังถูกประมวลผลในระดับรัฐ มีเครื่องมือ Arc หรือ GeoTools เพื่อทำให้เส้นเหล่านี้เรียบขึ้นและสร้างรูปหลายเหลี่ยมน้อยลงหรือไม่? ไม่เป็นไรหากผลลัพธ์สูญเสียความแม่นยำ ในท้ายที่สุดฉันกำลังมองหาสิ่งที่คล้ายกับสิ่งที่ฉันได้รับจากรูปหลายเหลี่ยมที่เรียบง่ายเพียงเข้าร่วมรูปหลายเหลี่ยมที่อยู่ติดกัน

2
วิธีเพิ่มจุดยอดนิยมไปยัง linestrings ที่มีอยู่
ถ้าฉันมี Linestring(1 2, 1 5, 1 9) และ Point(1 3) มีฟังก์ชั่นใด ๆ ที่สามารถรวมการคืนค่าและการเก็บรักษาจุดของคำสั่งเอาท์พุทดังนั้นจะเป็น: Linestring(1 2, 1 3, 1 5, 1 9)

3
วิธีการเวกเตอร์ผ้าม่านอย่างถูกต้องมากกว่า DEM?
ฉันมี DEM ของพื้นที่ศึกษาของฉันและลำธารในนั้น เมื่อใดก็ตามที่ฉันเพิ่มทั้งสองลงใน ArcScene และพับเวกเตอร์เหนือ DEM เวกเตอร์บางส่วนจะหายไปจากพื้นผิวและสามารถมองเห็นได้ภายใต้พื้นผิว (ดูภาพด้านล่าง) มีการตั้งค่าใด ๆ ที่เฉพาะเจาะจงเพื่อให้เวกเตอร์ทับ DEM อย่างถูกต้องหรือไม่? มีตัวเลือกหรือซอฟต์แวร์อื่นใดที่ให้ผลลัพธ์ที่ดีกว่า (มองเห็น) เมื่อคุณต้องใส่เวกเตอร์เหนือ DEM หรือไม่?

2
สคริปต์ Python ที่ทำงานใน ArcMap เทียบกับสคริปต์ที่ทำงานภายนอก?
ฉันเพิ่งเริ่มทำงานกับสคริปต์ Python เพื่อทำงาน ฉันกำลังสร้างสคริปต์เพื่อทำให้กระบวนการเป็นไปโดยอัตโนมัติ โดยทั่วไปแล้วมันจะถามชื่อลูกค้าของผู้ใช้รับการฉายถ้ามีสร้างไดเรกทอรีบนไดรฟ์ C: สำหรับลูกค้าสร้างฐานข้อมูลไฟล์ทางภูมิศาสตร์เฉพาะสำหรับลูกค้าสร้างชุดข้อมูลที่ต้องการและสร้างคลาสคุณลักษณะเฉพาะสำหรับข้อมูลลูกค้า ในที่สุดมันจะเพิ่มฟิลด์ที่จำเป็นให้กับแต่ละฟีเจอร์คลาสและอาจเป็นอย่างอื่น ฉันเริ่มสิ่งนี้โดยไม่รู้ถึงมารยาทที่เหมาะสมของการเขียนสคริปต์ Python สำหรับ ArcMap แต่สิ่งที่ฉันสร้างขึ้นจนถึงตอนนี้จะทำงานนอก ArcMap ที่ฉันเชื่อเท่านั้น เป็นที่ยอมรับหรือไม่? แทนที่จะได้รับการป้อนข้อมูลของผู้ใช้ผ่าน arcpy.getparamaterastext () ซึ่งฉันเพิ่งค้นพบเกี่ยวกับฉันใช้ raw_input () ไม่เป็นไร มันใช้งานได้ฉันไม่แน่ใจว่านี่เป็นวิธีที่เหมาะสมในการเขียนสคริปต์หรือไม่ นี่คือรหัสที่ฉันมี import sys import arcpy import os #Records name of the client client = raw_input("Enter the name of the client: (letters and underscores only) \n") #Records …
10 arcpy 

1
การแบ่งกลุ่มของ OpenLayers แบ่งเป็นชั้น ๆ
ฉันมีแผนที่ OpenLayers พร้อมการจัดกลุ่มจุด แต่ฉันต้องการใช้ข้อมูลที่ละเอียดยิ่งขึ้นกับข้อมูลที่แสดงในแผนที่ โดยเฉพาะฉันต้องการแบ่งชั้นกลุ่มตามแอตทริบิวต์ในคุณสมบัติ (เราจะเรียกมันว่า 'หมวดหมู่') ดังนั้นสมมติว่าฉันมีห้าหมวดหมู่: อย่างมีประสิทธิภาพฉันต้องจัดกลุ่มจุดเฉพาะกับจุดอื่น ๆ จากหมวดหมู่เดียวกัน ฉันคิดว่านี่สามารถทำได้ด้วยห้าชั้นที่แตกต่างกัน แต่สิ่งนี้กลายเป็นเรื่องยุ่งยากอย่างรวดเร็วเนื่องจากจำนวนหมวดหมู่มีขนาดใหญ่ (ฉันมี 50+ ในขณะนี้) มีใครเห็นอะไรเช่นนี้กับ OpenLayers?

1
กำลังมองหาทางเลือกบนคลาวด์บน Google Fusion Tables ใช่ไหม [ปิด]
ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ปิดเมื่อปีที่แล้ว คำตอบของคำถามนี้เป็นความพยายามของชุมชน แก้ไขคำตอบที่มีอยู่เพื่อปรับปรุงโพสต์นี้ ขณะนี้ไม่ยอมรับคำตอบหรือการโต้ตอบใหม่ ฉันเป็นผู้เริ่มต้นในการทำแผนที่เว็บ ฉันเล่นกับ Fusion Tables แต่ฉันรู้ว่ามีข้อมูลมากเกินกว่าจะแสดงได้ มีทางเลือกอะไรบ้างในตาราง Fusion ฉันมีข้อมูล 9 เลเยอร์ (หนึ่งเลเยอร์คือข้อมูลพัสดุดังนั้นมีมากกว่า 500,000 เรคคอร์ด) และฉันต้องสามารถสืบค้นข้อมูลได้ (เช่นจำนวนพัสดุที่อยู่ในระยะ 150m ของโรงเรียน ฯลฯ )

4
สอนเกี่ยวกับวิธีจำลองแบบน้ำท่วมโดยใช้ข้อมูลจากสถานีวัดตามแม่น้ำใน QGIS?
ใครช่วยกรุณาพาฉันไปสอนที่ฉันสามารถเรียนรู้เพิ่มเติมเกี่ยวกับวิธีการแบบจำลองน้ำท่วมโดยใช้ข้อมูลจากสถานีวัดตามแม่น้ำโดยใช้ QGIS?

2
ป้ายรูปหลายเหลี่ยมซ้ำกันสำหรับแต่ละกระเบื้อง
ฉันใช้ GeoWebCache แบบสแตนด์อโลนเพื่อให้บริการกระเบื้องจาก GeoServer ระยะไกล ปัญหาของฉันคือฉลากรูปหลายเหลี่ยมถูกเพิ่มลงในแผ่นกระเบื้องแต่ละแผ่นที่ให้บริการแทนที่จะแสดงเพียงครั้งเดียวในรูปหลายเหลี่ยม centroid ฉันพบโพสต์ที่กล่าวถึงปัญหา: http://osgeo-org.1560.n6.nabble.com/polygon-label-repeated-for-each-tile-td4995203.html คำตอบแรกที่กล่าวถึงวิธีแก้ปัญหาที่เป็นไปได้: "โดยรวมแล้วฉันแนะนำให้ใช้เครื่องมือการแสดงผลแบบเรียงต่อกัน (GeoWebCache, MapProxy, TileCache) ต่อไปแทนที่จะขอภาพขนาดเล็กจาก GeoServer และให้เอ็นจิ้นการแสดงผลแบบเรียงต่อกันมีการแบ่งกระเบื้องน้อยลงหลังจากนั้นคุณจะส่งคำขอน้อยลง 1 ภาพใหญ่แทนที่จะเป็นภาพเล็กหลายภาพ) ดังนั้นความเร็วในการสร้างไทล์แคชโดยรวมจะเพิ่มขึ้น " ปัญหาคือฉันไม่สามารถหาวิธีทำได้โดยอ้างถึงเอกสาร GeoWebCache และการโพสต์ดังกล่าวไม่ได้อธิบายวิธีการใช้งาน ฉันยังพบโพสต์พร้อมคำตอบที่เชื่อมโยงไปยังเอกสาร "GeoWebCache" Tiled "ของ GeoWebCacheแต่โค้ดของฉันใช้คุณสมบัติที่จำเป็นทั้งหมดแล้วและป้ายกำกับยังปรากฏขึ้นหลายครั้ง: var Layer_1874 = new OpenLayers.Layer.WMS( 'Grundkort', '/wms10.ashx' , { format: 'image/png', srs: 'EPSG:25832', layers: 'ballerupkommune_grundkort_bk', tiled: true, tilesOrigin: '698804,6173460' } , …

2
กำหนดจำนวนรูปหลายเหลี่ยมขั้นต่ำจาก shapefile เพื่อครอบคลุมพื้นที่ที่สนใจ
ฉันมีไฟล์รูปร่างเป็นจำนวนมากที่แสดงถึงพื้นที่ที่น่าสนใจสำหรับการวิเคราะห์ที่จะดำเนินการโดยใช้แหล่งที่มาของภาพถ่ายดาวเทียม (IKONOS, RapidEye ฯลฯ ) น่าเสียดายที่ภาพไม่ได้ใช้ระบบ Pathrow เช่น Landsat เป็นต้นขอบเขตจึงแตกต่างกันอย่างมาก ฉันมีรูปร่างที่ตัดกับ AOI แต่ละตัวซึ่งแสดงถึงขอบเขตของการได้รับภาพที่แตกต่างกันซึ่งทั้งหมดนี้ได้รับการยอมรับแล้ว บางรูปร่างเหล่านี้มีรูปหลายเหลี่ยมมากกว่า 500 รูป ฉันต้องการค้นหาวิธีที่ควรเป็นแบบอัตโนมัติ (Python และ ArcInfo 10 โดยเฉพาะอย่างยิ่ง FOSS ก็เป็นที่ยอมรับเช่นกัน) เพื่อกำหนดรูปหลายเหลี่ยมจำนวนน้อยที่สุดเพื่อครอบคลุมพื้นที่ที่ฉันสนใจแต่ละแห่ง

2
ใช้ pyshp เพื่อแปลงไฟล์. csv เป็น. shp หรือไม่
ฉันพยายามที่จะเข้าใจว่าฉันสามารถใช้โมดูล csv ใน python เพื่อเปิดไฟล์ csv ในโฟลเดอร์เดียวกับสคริปต์ python ได้อย่างไรจากนั้นสร้าง shapefile โดยใช้ shapefile module pyshp ไฟล์ csv มีลักษณะเช่นนี้ แต่สามารถมีเรคคอร์ดสองสามพันแถว: id_nr;date;target;start_lat;start_lon 1;2012-05-21;navpoint 25x;55.123654;13.456954 1;2012-05-23;navpoint 11f;55.143654;12.456954
10 python  shapefile  csv  pyshp 

1
Postgis 2.0 - มุมมอง public.geometry_columns ไม่มีข้อมูลที่ถูกต้องสำหรับมุมมอง
เราเพิ่งปรับปรุงฐานข้อมูลของเราเป็น Postgis 2.0 ในระบบของเราเราใช้ในการสร้างมุมมองสำหรับตารางที่มีคอลัมน์รูปทรงเรขาคณิต หลังจากสร้างมุมมองเรายังปรับปรุงตาราง geometry_columns ด้วยตนเองด้วยข้อมูลที่จำเป็น มุมมองใช้เพื่อสร้างเลเยอร์ Geoserver และฟีดข้อมูลสำหรับ ogr2ogr ปัญหาใน Postgis 2.0 คือตาราง geometry_columns ถูกแทนที่ด้วยมุมมองที่ดึงข้อมูลจากแคตตาล็อกระบบ ในกรณีของเราสำหรับมุมมองที่สร้างขึ้นเราจะไม่ดึงข้อมูลที่ถูกต้องในคอลัมน์ 'coord_dimension', 'srid' และ 'type' ข้อมูล : หลังจากการวิจัยบางอย่างฉันพบว่ามุมมอง geometry_columns กำลังใช้ตาราง pg_constraint เพื่อเติมคอลัมน์เหล่านั้น แต่สำหรับมุมมองนั้นไม่มีข้อ จำกัด ที่กำหนดไว้ ไม่มีใครมีความคิดใด ๆ ว่าเราจะแก้ไขปัญหานี้ได้อย่างไร? เราจะบังคับให้มุมมอง gemetry_columns ดึงข้อมูลที่จำเป็นได้อย่างไร มีวิธีแก้ปัญหาหรือวิธีอื่นในการทำเช่นนี้?

5
จะแบ่งถนน OSM ออกเป็นส่วน ๆ ได้ที่ทางแยกได้อย่างไร
ฉันต้องการสร้างเครือข่ายถนนสำหรับใช้กับ pgRouting โดยใช้ข้อมูล OpenStreetMap ฉันโหลด shapefile จาก GeoFabrik ลงในตาราง Postgres (โดยเปิดใช้งาน PostGIS) อย่างไรก็ตามปัญหาหนึ่งที่ฉันมีคือถนนมักจะไม่จบที่สี่แยกดังนั้นฉันจึงตัดสินใจแยกพวกเขาทุกครั้งที่สี่แยกหรือข้าม ในการระบุการแยกทั้งหมดที่ข้ามถนนหรือทางแยกฉันใช้สิ่งต่อไปนี้SQL(คล้ายกับคำถามก่อนหน้า ): CREATE TABLE split_points as SELECT DISTINCT ST_GeometryN(ST_Intersection(a.geom, b.geom),1) as geom FROM roads as a, roads as b WHERE ST_Touches(a.geom, b.geom) OR ST_Crosses(a.geom, b.geom) AND a.gid != b.gid GROUP BY ST_Intersection(a.geom, b.geom); ตอนนี้ฉันต้องการแยกถนนโดยใช้จุดเหล่านี้ ฉันใช้วิธีการต่อไปนี้: CREATE TABLE …

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