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

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

2
วิธีการแปลงรอยเท้าอาคารด้วยดิจิทัลเป็นมุมฉาก
ด้วยการใช้ QGIS ฉันกำลังพยายามทำแผนที่แรสเตอร์แบบดิจิทัลและเนื่องจากในที่สุดฉันต้องการโหลดเลเยอร์ 2d เหล่านี้ลงในแบบร่างเพื่อเปลี่ยนเป็นแบบจำลอง 3 มิติฉันต้องการให้ทุกอย่างเป็นรูปสี่เหลี่ยมจัตุรัส นี่คือตัวอย่างของแผนที่ที่ฉันทำงานด้วย: เส้นสีเขียวที่คุณเห็นนั้นเป็นเลเยอร์เวกเตอร์ของรางรถไฟที่ฉันได้ทำการแปลงเป็นดิจิทัลโดยใช้รูปทรงของเส้นส่วนที่เหลือเป็นเลเยอร์ตามแรสเตอร์ที่ฉันอ้างอิงทางภูมิศาสตร์ได้อย่างถูกต้องและเป็นพื้นหลังทุกชั้น สิ่งที่ฉันตั้งใจจะทำคือการจับภาพรอยเท้าอาคารตัวอย่างของฉันและการใช้เครื่องมือรูปหลายเหลี่ยมจับภาพปกติสิ่งนี้: ขณะนี้ทำงานได้อย่างเพียงพอสำหรับการสร้างแผนที่เวกเตอร์ 2 มิติ (EG: เพื่อแสดงเป็นแผ่นแผนที่เว็บหรือใช้ใน QGis) มันไม่ดีสำหรับการส่งออกเป็น KML จากนั้นนำเข้าสู่ Sketch-up เนื่องจากทุกอย่างไม่ได้เป็นสี่เหลี่ยมจัตุรัส 100% ยิ่งมันยากมากที่จะได้ของ 100% สแควร์ อย่างที่คุณเห็นจากเครื่องมือของฉันใน QGis ฉันได้ติดตั้งรูปสี่เหลี่ยมดิจิทัลและปลั๊กอินเครื่องมือ CAD แต่สิ่งเหล่านี้ไม่ดีสำหรับสิ่งที่ฉันต้องทำตัวอย่างเช่นถ้าเราใช้เครื่องมือจับรูปหลายเหลี่ยมแบบ orthogonal: ฉันสามารถเริ่มต้นวาดเส้นแรกของฉันที่มุมที่ฉันต้องการ แต่แล้วทันทีที่ฉันกด CTRL เพื่อให้ได้สิ่งต่าง ๆ ในมุมที่ถูกต้องและ / หรือสี่เหลี่ยมจัตุรัสถึงบรรทัดแรกทุกอย่างผิดปกติเพราะปลั๊กอินไม่ได้ นำมุมการหมุนของอาคารมาพิจารณาซึ่งหมายความว่าจุดเพิ่มเติมใด ๆ เป็นจริงโดยเฉพาะขอบเขตทิศเหนือ / ทิศตะวันออก / แนวตั้งและ 90 องศากับแกนเหล่านั้น หากอาคารที่ฉันพยายามตามรอยอยู่บนแกนเหล่านั้นก็ไม่มีปัญหาทุกอย่างสมบูรณ์แบบ แต่ส่วนมากไม่อยู่ในมุมต่าง …
10 qgis  digitizing 

2
ระเบียบวิธีสำหรับการจำแนกประเภท Hansen ของ Global Forest Watch?
ฉันได้ค้นพบเพียงแค่การจัดหมวดหมู่ Hansen's น่าทึ่งเกี่ยวกับการสูญเสียป่าไม้และเรื่องการเจริญเติบโตได้ที่http://www.globalforestwatch.org/ตีพิมพ์ในวิทยาศาสตร์ 2013เมื่อ: แฮนเซน, MC, Potapov พีวีมัวร์ R ., Hancher, M. , Turubanova, SA, & Tyukavina, A. (2013) แผนที่ความละเอียดสูงระดับโลกของการเปลี่ยนแปลงการปกคลุมป่าในศตวรรษที่ 21 วิทยาศาสตร์, 342 (6160) (15 พฤศจิกายน), 850–854 ดอย: DOI: 10.1126 / science.1244693 อย่างไรก็ตามฉันไม่สามารถหาวิธีการที่แน่นอนในการทำซ้ำแผนที่ดังกล่าวในเว็บไซต์นี้ / บนเว็บไซต์เพื่อให้การจำแนกประเภทที่ Hansen ใช้? สิ่งเดียวที่ฉันสามารถหาได้คืออัลกอริทึมการเรียนรู้ภายใต้การดูแลใช้เพื่อระบุต้นไม้ปกคลุม แต่มันค่อนข้างกว้าง ถ้าเป็นไปได้ฉันต้องการใช้วิธีการเดียวกัน (แต่ใช้กับ 90 ปี) ดังนั้นก่อนการจัดประเภทของ Hansen ในพื้นที่ที่ฉันเลือก

2
รับพิกัดและค่าพิกเซลที่สอดคล้องกันจาก GeoTiff โดยใช้ python gdal และบันทึกเป็นอาร์เรย์แบบ numpy
ฉันจะได้รับพิกัดที่ฉายรวมถึงค่าพิกเซลจริงที่พิกัดเหล่านั้นจากไฟล์ GeoTiff จากนั้นบันทึกลงในอาร์เรย์แบบ numpy ได้อย่างไร ฉันมีไฟล์ arsenci020l.tif และพิกัดอยู่ในหน่วยเมตร ด้านล่างคือผลลัพธ์ย่อของ gdalinfo ที่ฉันวิ่งไป ~$ gdalinfo arsenci020l.tif Driver: GTiff/GeoTIFF Files: arsenci020l.tif arsenci020l.tfw Size is 10366, 7273 Coordinate System is: PROJCS["Lambert Azimuthal Equal Area projection with arbitrary plane grid; projection center 100.0 degrees W, 45.0 degrees N", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], AUTHORITY["EPSG","6326"]], …

1
การคำนวณการเดินและขี่จักรยานแผนที่ Isochrone โดยใช้ QGIS & pgRouting?
เพิ่งติดตั้ง QGIS เพราะฉันเชื่อว่าฉันจะสามารถผลิตไอโซนิกเดินและขี่จักรยานได้จากจุดที่กำหนด ฉันทำสิ่งนี้ได้อย่างง่ายดายก่อนที่จะใช้ ArcMap แต่ฉันกำลังดิ้นรนกับ QGIS จริงๆ จากสิ่งที่ฉันได้อ่านฉันควรใช้การทำแผนที่ PGRouting และ OSM2PGrouting (ซึ่งฉันได้ติดตั้งแล้ว) ฉันไม่ได้ใช้ SQL มาก่อนและค้นหาแนวทางในการทำไอโซโทปเหล่านี้ ฉันได้อ่านเกี่ยวกับวิธีการต่างๆโดยละเอียดแล้ว SQL ส่วนใหญ่ที่ฉันกำลังดิ้นรนอยู่

2
ใช้ Python 3.x ใน ArcGIS หรือไม่
ฉันได้อ่าน " Whats ใหม่ใน ArcGIS " สำหรับการอัปเดต ArcGIS 10.3 และเจอกับส่วน Python และ ArcPy มันระบุว่า: ArcGIS 10.3 ได้รับการอัพเกรดให้เป็น Python 2.7.8 ตอนนี้ใช้ได้แล้วฉันไม่มีปัญหาอะไรกับงูหลามรุ่น 2.x ความกังวลของฉันคือ Esri จะผลักดัน ArcGIS ไปสู่ ​​Python 3x ในที่สุด ฉันหวังว่าพวกเขาจะไม่ได้ในอนาคตอันใกล้ แต่ในเวลาเดียวกัน Python 3x จะต้องมีการใช้งานในบางครั้ง ฉันทำการเขียนโปรแกรมในหลามมากมายสำหรับแอปพลิเคชั่น GIS และที่ไม่ใช่ GIS ทุกประเภท ฉันใช้เวอร์ชันปัจจุบันเสมอกับการเปิดตัว ArcGIS ของฉัน หากพวกเขาเปลี่ยนเป็น 3x ฉันจะใช้เวลาหนึ่งเดือนในการอัปเกรด 2x เป็น 3x สำหรับแอปพลิเคชันของฉัน ฉันแค่สงสัยว่าใครมีข้อมูลภายในที่เกี่ยวข้องกับหัวข้อนี้ …

1
วิธีการแสดง Mapbox Vector Tiles จากข้อมูล OSM
ฉันต้องการแสดงภาพเวกเตอร์ของตัวเองในรูปแบบ Mapbox Vector Tile รูปแบบ (MVT) โดยใช้ข้อมูล OpenStreetMap (สำหรับใช้กับ Mapbox GL JS client renderer ที่ร้อนแรงที่สุดของ Mapbox) TileStache + Mapnik + mapnik-vector-tile + PostGIS เป็นวิธีการทำเช่นนี้หรือไม่ ผู้ใช้กำหนดค่า Mapnik ให้แสดงผลแผ่น MVT ได้อย่างไร

1
ตั้งค่าเริ่มต้นสำหรับตารางค่าในเครื่องมือกล่องเครื่องมือหลาม
ฉันสร้างเครื่องมือ Python Toolbox เพื่อจัดลำดับฟิลด์ใหม่และสร้างคลาสฟีเจอร์ใหม่ด้วยฟิลด์ที่เรียงลำดับใหม่ เครื่องมือทำงานได้เป็นอย่างดีและฉันสามารถใช้ตารางค่าเพื่อให้ผู้ใช้จัดเรียงเขตข้อมูลตามลำดับที่พวกเขาเลือกหรือพวกเขาสามารถกรอกค่าอันดับสำหรับแต่ละเขตข้อมูล อย่างไรก็ตามส่วนที่น่ารำคาญของเครื่องมือนี้คือต้องเพิ่มเขตข้อมูลทั้งหมดในตารางค่าทีละครั้งก่อนที่จะจัดลำดับใหม่ ฉันกำลังพยายามตั้งค่านี้เพื่อนำฟิลด์ทั้งหมดไปยังตารางค่าโดยค่าเริ่มต้นและสามารถลบฟิลด์ที่ไม่ต้องการใด ๆ ก่อนที่จะจัดลำดับใหม่ มีใครเคยประสบความสำเร็จทำอะไรแบบนี้มาก่อนหรือไม่ ฉันพยายามทำสิ่งนี้ให้สำเร็จในวิธี UpdateParameters นี่คือรหัสที่ฉันพยายาม: import arcpy import os class Toolbox(object): def __init__(self): """Define the toolbox (the name of the toolbox is the name of the .pyt file).""" self.label = "Reorder Fields" self.alias = "Reorder Fields" # List of tool classes associated …

2
ST_ ระยะทางไม่ได้ใช้ดัชนีสำหรับการสืบค้นเชิงพื้นที่
ฉันไม่สามารถให้ PostGIS 2.1 ทำงานบน PostgreSQL 9.3.5 เพื่อใช้ดัชนีเชิงพื้นที่ได้แม้จะเป็นการสืบค้นที่ง่ายที่สุด ชุดข้อมูลทั้งหมดเป็น 8 ล้านคะแนน (จำนวนประชากรตารางจากที่นี่) ตารางถูกสร้างขึ้นเป็น CREATE TABLE points ( population DOUBLE PRECISION NOT NULL, location GEOGRAPHY(4326, POINT) NOT NULL ) CREATE INDEX points_gix ON points USING GIST(location); การสืบค้นนั้นง่ายมากตามที่ได้รับ SELECT SUM(population) FROM points WHERE ST_Distance( location, ST_GeographyFromText('SRID=4326; POINT(0 0)') ) < 1000 PostgreSQL ใช้การสแกน …

2
การไพพ์ข้อมูลไปที่ ogr2ogr
เป็นไปได้ไหมที่จะส่งข้อมูลโดยตรงไปที่ogr2ogr? ฉันพยายามcurlไฟล์ GeoJSON และไพพ์ไปที่ogr2ogrเพื่อแปลงเป็นรูปแบบอื่น โดยไม่ต้องใช้ไพพ์ฉันสามารถทำสิ่งนี้ได้โดยการเขียนไฟล์ก่อนแล้วแปลงมันโดยใช้คำสั่งต่อไปนี้ แต่ฉันอยากรู้ว่ามันเป็นไปได้ไหมที่จะข้ามขั้นตอนนั้น curl "https://raw.githubusercontent.com/nvkelso/natural-earth-vector/master/geojson/ne_50m_admin_0_countries.geojson" -o "countries.geojson" && ogr2ogr -f "KML" countries.kml countries.geojson ฉันยังพยายามสร้างตัวแปรด้วยเอาต์พุต curl แทนที่จะเป็น piping แต่ogr2ogrไม่สามารถอ่านได้: DATA=$(curl -s "https://raw.githubusercontent.com/nvkelso/natural-earth-vector/master/geojson/ne_50m_admin_0_countries.geojson") && ogr2ogr -f "KML" countries.kml $DATA
10 gdal  ogr2ogr 

2
ฟังก์ชั่นหลักใน QGIS จะใช้ประโยชน์จากมัลติเธรดหรือไม่?
หลังจากอ่านการเปิดใช้งานการประมวลผลแบบมัลติเธรดใน QGIS? ฉันสงสัยว่า QGIS 2.6 จะรวมสิ่งนี้หรือไม่เมื่อใช้กระบวนการต่างๆ ฉันตรวจสอบออนไลน์และฉันพบได้เพียงว่าการเรนเดอร์แบบมัลติเธรดนั้นได้รับการแนะนำให้รู้จักกับ QGIS 2.4 (ฉันคิดว่ามันจะใช้งานได้อีกใน 2.6) ตอนนี้แผนที่สามารถวาดใหม่ได้อย่างรวดเร็วโดยเฉพาะเมื่อจัดการกับชุดข้อมูลขนาดใหญ่ ฉันอ่านว่ามันค่อนข้างยากที่จะรวมฟังก์ชั่นมัลติเธรดเข้ากับคอร์และผู้ใช้ควรใช้โค้ด Python เพื่อให้บรรลุผลเช่นในการปฏิบัติการ Parallelising GIS ใน PyQGIS? . ฉันได้ตรวจสอบคำขอคุณสมบัติแล้วแต่ถูกปิดไป 9 เดือนโดย Tim Sutton ด้วยความคิดเห็นล่าสุด: "ฉันกำลังปิดตัวลงนี้ - Martin Dobias มีการนำไปใช้ในสาขาซึ่งจะถูกรวมเข้าด้วยกัน QGIS 2.2" ฟังก์ชั่นบางอย่างใน QGIS 2.6 ใช้ประโยชน์จากมัลติเธรด (หรือจะเน้นไปที่การเรนเดอร์อีกครั้ง) และถ้าไม่จะเปลี่ยนที่ QGIS 3.0 หรือไม่

1
OpenLayers 3 - วาดเส้น / เส้นทางหลายเส้นตามพิกัด
ฉันพยายามวาดเส้นตามพิกัดที่ให้ (จุดเริ่มต้นและจุดสิ้นสุด) Googled พบตัวอย่างเล็ก ๆ น้อย ๆ แต่ดูเหมือนว่าไม่ใช่พวกมันอาจจะทำงานได้เพราะมันใช้กับ OL2 ดังนั้นนี่จึงเป็นทางเลือกสุดท้ายของฉัน พิกัดตั้งอยู่ในอาร์เรย์เครื่องหมาย <!doctype html> <html lang="en"> <head> <link rel="stylesheet" href="css/ol.css" type="text/css"> <style> .map { height: 100%; width: 100%; } </style> <script src="build/ol.js" type="text/javascript"></script> </head> <body> <div id="map" class="map"></div> <script type="text/javascript"> // inicijalizacija mape var map = new ol.Map({ target: 'map', layers: …
10 openlayers 

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

2
QGIS ระยะทางของส่วนข้ามของโปรไฟล์โดยใช้ Ellipsoid?
ฉันกำลังทำงานกับ QGIS 1.8 และ 2.6 และไม่สามารถใช้เครื่องมือโปรไฟล์ข้ามส่วนใด ๆ เพื่อให้ระยะทางรี ("เหมาะสม") ระยะทางทั้งหมดใช้พิกัดทางภูมิศาสตร์ซึ่งยาว / ไกลเกินไป ตัวอย่างเช่นระยะทางเส้นตรงจาก Vancouver BC ถึง Kelowna ประมาณ 270km ภาพตัดขวางของฉันบอกว่ามันคือ 420km เครื่องมือวัดระยะทาง QGIS ช่วยให้คุณใช้ Ellipsoid (ฉันตั้งค่าของฉันเป็น WGS 84) และให้ฉัน 271km ทุกคนสามารถอธิบายวิธีรับเครื่องมือโปรไฟล์เพื่อให้ระยะทาง "ถูกต้อง" ได้หรือไม่

2
ดาวน์โหลด SRTM ใหม่ 30 เมตร
ฉันอ่านเกี่ยวกับรุ่นใหม่ของรุ่นระดับความสูงดิจิตอลที่ 30 เมตร (SRTM 30 เมตร) แต่ฉันไม่สามารถหาได้จากที่ฉันสามารถดาวน์โหลดได้ แหล่งที่มาบางส่วน: http://glcfapp.glcf.umd.edu:8080/esdi/index.jsp http://asterweb.jpl.nasa.gov/data.asp http://www2.jpl.nasa.gov/srtm/

3
การจัดการเงาของภูเขาเมื่อจำแนก Landsat 8
ฉันพยายามจำแนกบางพื้นที่ทางตอนเหนือของประเทศมองโกเลียโดยใช้ภาพถ่ายดาวเทียม Landsat 8 จากปี 2013 เนื่องจากการบันทึกถูกสร้างขึ้นในฤดูหนาวดวงอาทิตย์จึงต่ำมากในเวลาที่ซื้อกิจการ ดังนั้นจึงมีเงาที่ยาวและมืดมากจากภูเขา ฉันสามารถระบุเงานี้โดยใช้ DEM ตามที่อธิบายไว้ในคำถามนี้: วิธีลบเอฟเฟกต์ภูมิประเทศเช่นเดียวกับเงาจากภาพ Landsat ฉันจะจัดหมวดหมู่ภายใต้การดูแลสำหรับพื้นที่เงาที่ถูกคลิปได้อย่างไร เป็นไปได้หรือไม่ที่จะยกระดับพื้นที่นี้ ฉันลองอัตราส่วนวงบ้าง แต่ไม่รู้ว่าอะไรดีที่สุดสำหรับงานของฉัน ในภาพนี้คุณจะเห็นว่ามีบางพื้นที่ที่ปลูกในเงามืด แต่ไม่สามารถจัดประเภทได้

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