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

บางครั้งฉันเห็นsrsName="EPSG:4326"ในเอกสาร GML (เช่นhttps://postgis.net/docs/ST_GeomFromGML.htmlถึงแม้ว่าแหล่งข้อมูลเฉพาะของฉันจะไม่เปิดเผยต่อสาธารณะ) เจตนา (ยกเว้นบางทีสำหรับลำดับแกนที่ไม่สำคัญที่นี่) มีความชัดเจนและอาจทำงานร่วมกันข้ามระบบต่าง ๆ ได้ แต่ฉันสงสัยว่าเป็นไปตามข้อกำหนดจริงหรือไม่ ฉันได้ดูนโยบายการตั้งชื่อ OGCและOGC URNแล้วและพวกเขาใช้urn:ogc:def:crs:EPSG:\d{4,5}แบบฟอร์ม (เช่นจะเป็นurn:ogc:def:crs:EPSG:4326) โปรไฟล์ GML Point (OGC 05-029r4) ใช้srsName="urn:ogc:def:crs:EPSG:6.6:4326"ในตัวอย่างและต้องการให้เป็น "รู้จักกันดี" ฉันไม่เห็นอะไรในนั้นที่จริงไม่รวมใด ๆsrsNameแม้ว่า ข้อกำหนดเฉพาะของsrsName="EPSG:abcd"แบบฟอร์ม (โดยเฉพาะsrsName="EPSG:4326") สอดคล้องกับองค์ประกอบ GML 3.1 หรือ 3.2 (โดยเฉพาะอย่างยิ่งสำหรับโปรไฟล์โปรไฟล์) หรือไม่

11 gml  ogc  specification 

ฉันมีฟังก์ชั่นที่สร้างแผงเซลล์แสงอาทิตย์ Photovolatic ที่แสดงเป็นรูปหลายเหลี่ยม โดยพื้นฐานแล้วมันจะสร้างกริดสี่เหลี่ยมที่ผู้ใช้สามารถระบุพารามิเตอร์ต่อไปนี้: ความยาว ความกว้าง ระยะทางแนวนอน ระยะทางแนวตั้ง รหัสนี้เป็นไปตามปลั๊กอินFeatureGridCreatorแต่เน้นเฉพาะในรูปหลายเหลี่ยม ใช้ได้ดีกับส่วนใหญ่โดยเฉพาะเมื่อสร้างรูปหลายเหลี่ยมที่มีขนาดใหญ่ (เช่นความยาวและความกว้าง 10 ม. ระยะทางแนวนอนและแนวตั้ง 10 ม.) แต่ฉันสังเกตเห็นสองประเด็น: เมื่อระบุรูปหลายเหลี่ยมสำหรับขนาดที่น้อยกว่า 2 เมตรสำหรับทั้งความยาวและความกว้างจะไม่สร้างรูปหลายเหลี่ยม เมื่อระบุรูปหลายเหลี่ยมที่มีขนาดต่างกัน (เช่นความยาว 5 ม. และความกว้าง 7 ม.) ขนาดนั้นไม่เหมือนกันเมื่อวัดด้วยเครื่องมือวัดเส้น สำหรับมิติเหล่านี้ความยาวและความกว้างถูกแสดงเป็น 4m และ 6m ตามลำดับ CRS ที่ใช้สำหรับการฉายภาพและเลเยอร์คือEPSG: 27700แม้ว่าฉันจะไม่คิดว่ามันจะเป็นปัญหา ดังนั้นใครบ้างมีความคิดใด ๆ ที่อาจทำให้เกิดปัญหาเหล่านี้ ฉันยังเปิดให้คำแนะนำเกี่ยวกับวิธีการปรับปรุงรหัสหรือแทนที่ด้วยทางเลือกที่ดีกว่า นี่คือรหัสที่สามารถทำซ้ำได้ในPython Consoleต้องเลือกเลเยอร์รูปหลายเหลี่ยมด้วย CRS ที่เกี่ยวข้องก่อนเรียกใช้ฟังก์ชัน: from PyQt4.QtCore import QVariant from math …

11 coordinate-system  pyqgis  polygon-creation 

ฉันมีโครงการ QGIS ซึ่งมีนักแต่งเพลง แน่นอนว่าถ้าฉันเข้าสู่ GUI ฉันสามารถส่งออกเป็น pdf และอื่น ๆ สำหรับสคริปต์ทำของฉันแม้ว่าฉันต้องการทำสิ่งนี้ผ่านสคริปต์ python จากการอ่านออนไลน์นี่คือความพยายามของฉัน: #!/usr/bin/env python import sys from qgis.gui import QgsMapCanvas, QgsLayerTreeMapCanvasBridge from PyQt4.QtCore import QFileInfo from PyQt4.QtXml import QDomDocument from qgis.core import * qgs = QgsApplication(sys.argv, True) QgsApplication.setPrefixPath("/usr", True) qgs.initQgis() def make_pdf(): canvas = QgsMapCanvas() bridge = QgsLayerTreeMapCanvasBridge( QgsProject.instance().layerTreeRoot(), canvas) …

11 qgis  python  pyqgis  pdf 

ทุกครั้งที่ฉันสร้างสคริปต์เพื่อจัดการเวกเตอร์ (จุดเส้นและรูปหลายเหลี่ยม) ในหลามฉันพบข้อสงสัยนี้ ความแตกต่างระหว่างกล่องขอบ (bbox), ซองจดหมาย, ขอบเขต, ขอบเขตคืออะไร?

11 geometry  extents  gis-principle  terminology  bbox 

ฉันต้องการเร่งความเร็วกระบวนการที่สอดคล้องกับคำสั่งนี้: gdalwarp -t_srs EPSG:4326 -overwrite input.ntf output.tif มีวิธีใช้การประมวลผลแบบขนานใน GDAL หรือไม่ กึ่งนอกหัวข้อ: ถ้าไม่คุณแนะนำวิธีแก้ปัญหาที่ไม่ใช่ GDAL เพื่อเร่งกระบวนการ gdalwarp เหมือน? นี่คือหน้าเว็บที่ฉันได้ดู: เป็นไปได้ไหมที่จะทำการประมวลผลแบบขนานใน GDAL และ QGIS? [gdal-dev] การประมวลผลแรสเตอร์ของ GDAL: การคำนวณแบบขนาน แก้ไข: ฉันถามคำถามนี้เพราะฉันคิดว่าฉันเห็น CPU% อยู่ต่ำกว่า 100% เมื่อประมวลผล gdalwarp อย่างไรก็ตามในการค้นหาครั้งที่สองก็มาถึง 555.5%

11 gdal  parallel-processing 

ฉันกำลังพยายามเขียนสคริปต์ที่จะบันทึกการเรนเดอร์ของเลเยอร์หลาย ๆ อันโดยใช้โปรแกรมแต่งแผนที่ ปัญหาที่ฉันพบคือสคริปต์บันทึกไว้ก่อน qgis แสดงผลเลเยอร์ทั้งหมดเสร็จสิ้นแล้ว จากคำตอบอื่น ๆ อีกหลายข้อ ( 1 , 2 , 3 ) ฉันพยายามใช้iface.mapCanvas.mapCanvasRefreshed.connect()และวางการบันทึกภาพไว้ในฟังก์ชั่น แต่ฉันยังคงพบปัญหาเดียวกัน - ภาพไม่รวมเลเยอร์ทั้งหมด รหัสที่ฉันใช้รวมถึงรูปภาพของหน้าต่างหลักและการเรนเดอร์หน้าตามีลักษณะดังรายการด้านล่าง ฉันสังเกตว่าถ้าฉันเปิดหน้าต่างคอนโซลและไม่แสดงข้อคิดเห็นทั้งสามprint layerListบรรทัดโปรแกรมจะรอการเรนเดอร์ให้เสร็จก่อนบันทึกภาพ ฉันไม่แน่ใจว่านี่เป็นเพราะเวลาในการประมวลผลที่เพิ่มขึ้นหรือหากมีการเปลี่ยนแปลงวิธีการทำงานของโปรแกรม ฉันจะใช้สิ่งนี้อย่างถูกต้องได้อย่างไรเพื่อให้เลเยอร์ทั้งหมดรวมอยู่ในภาพ from qgis.core import * from qgis.utils import * from qgis.gui import * from PyQt4.QtCore import * from PyQt4.QtGui import * import os.path ##StackExchange Version=name ##Map_Save_Folder=folder …

11 pyqgis  rendering  qgis-2.18 

เมื่อฉันเปิด QGIS แผงเลเยอร์จะใช้ค่าเริ่มต้นเป็นความกว้างที่ไม่แสดงปุ่ม "ลบเลเยอร์ / กลุ่ม" ดังนั้นฉันต้องปรับขนาดพาเนลเพื่อแสดงปุ่มทุกครั้งที่ฉันเปิด QGIS เป็นไปได้ไหมที่จะตั้งค่าความกว้างเริ่มต้นสำหรับพาเนลนี้ดังนั้นฉันไม่จำเป็นต้องปรับขนาดหรือไม่ ฉันใช้ QGIS 2.18.5 บน Windows 10

11 qgis  pyqgis  qgis-2.18  table-of-contents  panel 

ฉันต้องการสร้างรหัสอัตโนมัติในขณะที่สร้างคุณสมบัติใหม่ ดังนั้นทุกครั้งที่ฉันสร้างจุดใหม่ ID ของมันจะเพิ่มขึ้น 1 เช่น: Feature 1: id 1 Feature 2: id 2 Feature 3: id 3 ฉันจะโจมตีสิ่งนี้ได้อย่างไร

11 qgis  point  features  qgis-custom-function 

ฉันเพิ่มองค์ประกอบสองสเกลลงในองค์ประกอบแผนที่ของฉัน สำหรับการประสานกันของเลย์เอาต์ฉันต้องการเห็บและข้อความฉลากด้านล่างแถบขนาด มาตรฐานคือ: สิ่งที่ฉันต้องการคือสิ่งนี้ (เปลี่ยนภาพด้วยตนเอง): มีตัวเลือกใน QGIS สำหรับสิ่งนี้หรือไม่?

11 qgis  print-composer  scale  scale-bar 

ในขณะนี้ฉันพยายามใช้ r.mapcalc จากกล่องเครื่องมือประมวลผล QGIS อย่างไรก็ตามอัลกอริทึมไม่ได้ให้ผลลัพธ์ใด ๆ แก่ฉัน ฉันลองหลายชั้นแบบแรสเตอร์และนิพจน์ที่แตกต่างกัน แต่อัลกอริทึมไม่เคยส่งผลอะไรกับ "ไดเรกทอรีผลลัพธ์" ด้านล่างมีรูปภาพที่แสดงตัวอย่างง่ายๆของการตั้งค่าอัลกอริทึม บางทีฉันอาจจะพลาดบางสิ่งบางอย่าง แต่ฉันไม่พบสิ่งที่คล้ายกันทางออนไลน์

11 qgis  raster  raster-calculator  qgis-processing  qgis-grass-plugin 

ฉันมีไฟล์ในการฉายภาพ British National Grid: Geometry: 3D Polygon Feature Count: 5378 Extent: (9247.520209, 14785.170099) - (638149.173223, 1217788.569952) Layer SRS WKT: PROJCS["British_National_Grid", GEOGCS["GCS_airy", DATUM["OSGB_1936", SPHEROID["Airy_1830",6377563.396,299.3249646]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",49], PARAMETER["central_meridian",-2], PARAMETER["scale_factor",0.9996012717], PARAMETER["false_easting",400000], PARAMETER["false_northing",-100000], UNIT["Meter",1]] cat: Integer (9.0) ฉันสามารถใช้ GDAL / OGR เพื่อให้ได้พื้นที่รวมของรูปหลายเหลี่ยมทั้งหมดใน shapefile เป็นเฮกตาร์หรือไม่ ฉันสงสัยว่ามันเป็นไปได้ด้วย-sqlหรือเปล่า: ogrinfo -sql "SELECT SUM(ST_Area(geom::geography)) FROM mytable" myshapefile.shp …

11 gdal  ogr  area  ogrinfo 

ฉันจำเป็นต้องรู้หรือไม่ว่าทางวิทยาศาสตร์นั้นมีความถูกต้องในการหลอมรวม / สะท้อนแสงผลิตภัณฑ์พื้นผิวของ Landsat 8 ที่มีแถบแบนด์ของวงนั้นหรือไม่? Landsat รายละเอียดของผลิตภัณฑ์สะท้อนสามารถพบได้ที่นี่ จะต้องมีการกล่าวถึงว่าหนึ่งต้องสั่งซื้อผลิตภัณฑ์สะท้อนแสงพื้นผิวแยกต่างหากเพื่อให้ได้ผลิตภัณฑ์นี้ ผลิตภัณฑ์นี้มีเพียง 7 แบนด์ (30m) ไม่ใช่ IR และ Pan band ดังนั้นอีกครั้งคำถามของฉันคือสามารถใช้ฟิวส์ 7 แบนด์ (30 ม.) ของผลิตภัณฑ์สะท้อนแสงพื้นผิวกับแพนปกติ (ไม่ใช่สะท้อนพื้นผิว) (15m) ฉันต้องการใช้ภาพ pansharpened นี้สำหรับการแบ่งส่วนและทำตามการทำแผนที่ที่ดินดังนั้นฉันจำเป็นต้องรู้ว่ามีการปฏิบัติใด ๆ ที่จัดตั้งขึ้นของการเหลาแพนประเภทนี้ในสถาบันการศึกษาที่มีการอ้างอิงถ้าใช่โปรดอ้างอิง

11 remote-sensing  digital-image-processing  landsat-8  pansharpening 

ฉันมีข้อมูลจำนวนมากในรูปแบบของ CSV ที่มีพิกัด UTM XและYและค่าระดับความสูงเป็นZข้อมูล ฉันจำเป็นต้องแปลงข้อมูลเหล่านี้เป็น DEM เป็น GeoTIFF เพื่อการวิเคราะห์เพิ่มเติม ในกรณีนี้จำนวนมากหมายถึง 16 เมตร เส้นที่มีในจุดหนึ่งX, YและZต่อบรรทัด คะแนนมีการกระจายอย่างเท่าเทียมกันดังนั้นจึงไม่จำเป็นต้องมีการแก้ไข แต่ละจุดจะต้องถูกแปลงเป็นเซลล์แรสเตอร์ ข้อมูลต้นฉบับมาโดยไม่มีตัวคั่นพร้อมความกว้างคอลัมน์คงที่ ฉันคิดแล้วออกวิธีการแปลงไฟล์ไวยากรณ์ที่จะใช้คั่นแทนความกว้างคงที่และกำจัดตัวละครทุกพื้นที่โดยใช้โปรแกรมแก้ไขข้อความกระแสsed จากที่นี่ที่ปกติขั้นตอนการทำงานของฉันจะเป็นที่จะนำเข้าข้อมูลลงใน ArcGIS โดยการสร้างระดับคุณลักษณะจากX, YและZข้อมูลและเป็นขั้นตอนที่สองแปลง shapefile จุดเป็น GeoTIFF ใช้ชี้ไปที่ Rasterเครื่องมือ อย่างไรก็ตามไฟล์ที่ฉันมีในขณะนี้มีขนาดใหญ่เกินไปสำหรับกระบวนการนี้ แทนที่จะเป็นกระบวนการทำงานที่อธิบายไว้ข้างต้นฉันกำลังมองหาทางเลือกที่มีประสิทธิภาพและค้นพบ GDAL อย่างไรก็ตามในgdal_translateรูปแบบที่สนับสนุนที่ใกล้เคียงที่สุดที่ฉันสามารถหาได้ในรายการประเภทไฟล์ที่รองรับคือตาราง ASCII แต่ไม่มี XYZ คั่นด้วยเครื่องหมายจุลภาค ปัญหาอีกประการหนึ่งคือฉันมีพิกัด UTMในขณะที่ตัวอย่างส่วนใหญ่ดูเหมือนจะใช้พิกัดองศาแบบทศนิยม อย่างไรก็ตามฉันต้องอยู่ในระบบ UTM (หรืออย่างน้อย GeoTIFF ของฉันจะต้องอยู่ในระบบพิกัด UTM) ดังนั้นฉันกำลังมองหาวิธีแปลง CSV XYZ เป็น GeoTIFF โดยใช้ …

11 gdal  geotiff-tiff  csv  gdal-translate  xyz 

สคริปต์พื้นฐานเพื่อวนซ้ำแบบซ้ำในโฟลเดอร์ย่อยและรวมรูปร่างไฟล์ทั้งหมดไว้ในไฟล์เดียวคือ: #!/bin/bash consolidated_file="./consolidated.shp" for i in $(find . -name '*.shp'); do if [ ! -f "$consolidated_file" ]; then # first file - create the consolidated output file ogr2ogr -f "ESRI Shapefile" $consolidated_file $i else # update the output file with new file content ogr2ogr -f "ESRI Shapefile" -update -append $consolidated_file …

11 gdal  ogr2ogr  merge  bash  mergeshapes 

ฉันกำลังสร้างเครื่องมือที่ดาวน์โหลดภาพจาก "Sentinel-2 on AWS" ฉันมีปัญหาในการแปลพิกัด Long / Lat เป็นกระเบื้อง Sentinel ที่เฉพาะเจาะจง ฉันรู้ว่า Sentinel 2 ใช้ MGRS และฉันพยายามที่จะแปลง Long / Lat เป็น MGRS และด้วยความแม่นยำต่ำฉันจะได้รับรหัสประจำตัว S2A และใช้งานได้ แต่ไม่เสมอไป ตัวอย่างเช่น Long / Lat 34.665,31.625 แก้ไขเป็น "36R XA" ใน MGRS แต่ไม่มี Sentinel ไทล์ 36RXA อะไรคือวิธีที่ถูกต้องในการพิจารณาไทล์เนลโดยใช้พิกัด Long / Lat

11 remote-sensing  latitude-longitude  sentinel-2  mgrs  amazon-web-services