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

ฉันได้อ่านใน OpenLayers แล้ววิธีโหลด KML, GeoJSON ฯลฯ ไฟล์เวกเตอร์ใน OpenLayers แต่ฉันจะเชื่อมต่อกับข้อมูล PostGIS ได้อย่างไร มิฉะนั้นฉันควรอัพโหลดข้อมูลของฉันไปที่ PostGIS แทนที่จะวางข้อมูลของฉันโดยตรงโดยใช้ WebServices ฉันไม่มีอะไรเลยหรือ และฉันต้องการที่จะโหลดข้อมูลเวคเตอร์นี้เพื่อที่เราจะได้เห็นแอททริบิวต์โดยคลิกที่ฟีเจอร์ ฉันต้องการ PostGIS เพราะในอนาคตฉันต้องการเปิดใช้งานแบบสอบถาม ฉันคิดว่าเป็นไปได้โดยการรักษาฐานข้อมูลเช่น PostGIS แต่วิธีการใช้ OpenLayers? ฉันได้อ่านคู่มือผู้เริ่มต้นของ OpenLayers แล้ว แต่ไม่ได้มีการกล่าวถึง กรุณาช่วย!!

10 openlayers-2  postgis  vector 

ฉันมีชุดข้อมูลแรสเตอร์ความละเอียด 100 ม. ขนาดใหญ่ (ขนาดต่ำกว่า 48 US) ที่มีตารางแอตทริบิวต์แรสเตอร์ (RAT) ข้อมูลถูกสร้างโดยคนอื่น หนึ่งในแอตทริบิวต์ใน RAT มีค่าข้อมูลที่ฉันสนใจในขณะที่ค่าเซลล์ raster จริงนั้นไม่มีความหมาย ฉันพยายามหาวิธีกำหนดค่าของแอตทริบิวต์ที่ฉันสนใจเกี่ยวกับค่าของเซลล์แรสเตอร์จริง แต่ไม่สามารถหาวิธีการทำสิ่งนี้ได้ เหตุผลที่ฉันต้องการทำเช่นนี้คือฉันมีเครื่องมือ Python แบบกำหนดเองที่เกี่ยวข้องอย่างเป็นธรรมซึ่งคนอื่นพัฒนาขึ้นมาเพื่อทำการคำนวณที่ฉันต้องการ แต่น่าเสียดายที่มันไม่มีความสามารถในการใช้แอตทริบิวต์ RAT ใช้เฉพาะค่าที่กำหนดให้กับแรสเตอร์ของเซลล์ซึ่งอย่างที่ฉันบอกในตัวอย่างนี้ไม่มีความหมาย ฉันลองใช้วงเวียนในการรันแรสเตอร์เพื่อชี้เครื่องมือเพื่อสร้างคลาสฟีเจอร์แบบชี้ด้วยฟิลด์เป้าหมายที่เป็นแอททริบิวที่ฉันสนใจโดยคิดว่าฉันสามารถสร้างแรสเตอร์ใหม่จากคลาสฟีเจอร์จุดได้ อย่างไรก็ตาม raster to point tool ไม่สามารถระบุได้ว่าฉันมีหน่วยความจำชั่วคราวไม่เพียงพอ (ฉันใช้เครื่อง Alienware ที่มี RAM ขนาด 16 GB และไดรฟ์ข้อมูล 1 TB) ฉันมีตัวเลือกหรือไม่? ฉันใช้ ArcGIS 10 แต่ถ้าฉันสามารถทำสิ่งที่ฉันอธิบายได้อย่างง่ายดายในแพ็คเกจอื่นฉันยินดีที่จะลองวิธีอื่น (ฉันไม่ได้แต่งงานกับ ArcGIS)

10 arcgis-desktop  arcgis-10.0  spatial-analyst  fields-attributes 

ฉันต้องการฉายโครงการ QGIS อีกครั้ง แต่ไม่มีตัวเลือกดังกล่าว เมื่อค้นหาคำตอบฉันเห็นogr2ogrคำแนะนำสำหรับงานนี้โดยฉายแต่ละเลเยอร์เวกเตอร์ใหม่ (shapefile ในตัวอย่างนี้) อย่างไรก็ตามฉันไม่ได้รับผลลัพธ์ที่คาดหวัง: C:\temp>ogrinfo -al -so misc.shp INFO: Open of `misc.shp' using driver `ESRI Shapefile' successful. Layer name: misc Geometry: Line String Feature Count: 10 Extent: (21.267388, 42.015857) - (21.270225, 42.017470) Layer SRS WKT: PROJCS["MGI_Balkans_zone_7_deprecated", GEOGCS["GCS_MGI", DATUM["Militar_Geographische_Institute", SPHEROID["Bessel_1841",6377397.155,299.1528128]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",21], PARAMETER["scale_factor",0.9999], PARAMETER["false_easting",7500000], …

10 shapefile  coordinate-system  ogr2ogr 

ฉันเป็นนักสำรวจและใช้ QGIS ในการดูแลรักษาบันทึกการทดสอบดินในโครงการโยธาและได้รับการจัดให้มีพื้นที่ทดสอบอย่างต่อเนื่องหรือจำนวนที่ต้องการ เหล่านี้มักจะมากกว่าในรูปทรงที่ผิดปกติและฉันสงสัย มีปลั๊กอิน / โซลูชันเพื่อตัดเลเยอร์ที่มีรูปหลายเหลี่ยมที่ผิดปกติเป็นพื้นที่เท่ากันหรือจำนวนล็อตที่มีพื้นที่ที่ระบุผลลัพธ์ที่ได้คือจำนวนล็อตมากที่สุดตามพื้นที่ที่ให้ไว้หรือไม่

10 qgis  qgis-plugins  area  splitting 

ฉันกำลังใช้ Leaflet รวมกับคุณสมบัติ GeoJSON มีวิธีติดป้ายกำกับคุณลักษณะ GeoJSON (ในกรณีนี้คือรูปหลายเหลี่ยม) หรือไม่ ควรได้รับฉลากจาก feature.properties.name นี่คือรหัสของฉันที่ฉันคิดว่าฉันสามารถแทรกป้ายกำกับ: function style(feature) { return { weight: 2, opacity: 1, color: 'white', dashArray: '3', fillOpacity: 0.7, fillColor: getColor(feature.properties.coloring) }; }

10 javascript  web-mapping  labeling  geojson  leaflet 

ฉันต้องการสร้าง vectorLayer QGISเพื่อแสดงจุดทั้งหมดที่เป็นแบบฟอร์มLINESTRINGหรือMULTILINESTRINGเก็บไว้ในPostGISฐานข้อมูล ฉันคิดว่าฉันต้องการที่จะเปลี่ยนทุกจุดของLINESTRING(หรือMULTILINESTRING) POINTเพื่อ คำถามแรก : มีPostGISฟังก์ชั่นที่จะทำเช่นนั้น? คำถามที่สอง : หากฉันต้องการแยกคะแนนที่ซ้ำกันออกจากPOINTตารางผลลัพธ์ฉันจะทำอย่างไร

10 qgis  postgis  point  polyline-creation  linestring 

ฉันมีคอลเลคชั่นสองแบบ: 1 - รหัสไปรษณีย์เขตการปกครองและการเมือง 2 - หมวดหมู่การใช้ที่ดินต่าง ๆ เช่นสวนสาธารณะ & การจัดสรรและคะแนนเช่นต้นไม้ ฉันต้องการจัดทำสถิติเกี่ยวกับพื้นที่ของหมวดหมู่การใช้ที่ดินเหล่านั้นแยกตามรหัสไปรษณีย์ผู้ดูแลและพื้นที่ทางการเมืองและจำนวนคะแนนด้วยวิธีการเดียวกัน ตัวอย่างเช่นมีสวนสาธารณะ 4.1ha ใน SE15, 879 ต้นไม้ใน SE22 และอื่น ๆ ฉันกำลังใช้ QGIS ฉันได้ประมวลผลข้อมูลดิบของฉันเพื่อเตรียมรูปร่างให้พร้อมแล้ว แต่หลังจากการค้นหาที่ไร้ผลมากในเว็บไซต์นี้และเว็บไซต์อื่น ๆ ฉันไม่เห็นวิธีการทำเช่นนี้

10 qgis  intersection  statistics 

clip_fc = arcpy.GetParameter(0)สคริปต์ของฉันทำงานได้ดีกับชั้นคุณลักษณะเดียวเลือกสำหรับ อย่างไรก็ตามเมื่อฉันตั้ง clip_fc เป็น MultiValue ในคุณสมบัติของสคริปต์และจากนั้นเลือกคุณสมบัติหลายชั้นฉันได้รับข้อผิดพลาด "พารามิเตอร์ไม่ถูกต้อง" ทุกคนสามารถแสดงให้ฉันเห็นว่าฉันกำลังจะไปไหนผิดหรือเปล่า? import arcpy import os.path from arcpy import env arcpy.env.workspace = r"C:\GIS Home\project_1\clipshps\Combined_mf.gdb" arcpy.env.overwriteOutput = True clip_fc = arcpy.GetParameter(0) clip_name = str(clip_fc) output_dir = r"C:\GIS Home\project_1\clipshps\outputs_clip" for fds in arcpy.ListDatasets(): for fc in arcpy.ListFeatureClasses('','',fds): out = os.path.join(output_dir, fc + "_" + …

10 arcgis-desktop  arcpy  clip 

ฉันต้องการสร้าง.tiffไฟล์อ้างอิงทางภูมิศาสตร์จาก.pngไฟล์ ดังนั้นฉันแปลง. png เป็น.tiffไฟล์ปกติโดยใช้บรรทัดคำสั่งต่อไปนี้: convert image.png image.tiff ที่นี่ทำงานได้ดีฉันมี.tiffไฟล์ที่ถูกต้องซึ่งแสดงผลได้ดี จากนั้นฉันต้องการเพิ่มข้อมูลการอ้างอิงทางภูมิศาสตร์ลงไป ในการทำเช่นนั้นฉันพบเครื่องมือชื่อtiff2geotiffฉันใช้มันเช่นนี้: tiff2geotiff -4 "+proj=latlong +ellps=WGS84 +datum=WGS84 +no_defs" -c none -n "-122.5575664 38.5818201 -122.5449425 38.5896175" notgeoreferenced.tiff georeferenced.tiff สิ่งนี้ทำให้ฉันมีภาพสีดำซึ่งนำเข้ามาใน QGis แสดงได้ดีในสถานที่ที่ดียกเว้นว่าเป็นสีดำทุกที่ (พิสูจน์ได้ว่ามีการเพิ่มข้อมูลทางภูมิศาสตร์ที่ได้รับการพิสูจน์มาอย่างดี ) ดังนั้นเมื่อฉัน: gdalinfo georeferenced.tiff ฉันได้รับ: Warning 1: TIFFReadDirectory:Bogus "StripByteCounts" field, ignoring and calculating from imagelength Driver: GTiff/GeoTIFF Files: araujo_tmp.tiff Size is …

10 qgis  gdal  php  geotiff-tiff 

เรามีคุณสมบัติจากข้อมูลการสำรวจที่มีข้อมูล 3 มิติบางส่วน ตัวอย่างที่พบบ่อยที่สุดคือ LineString 2D ซึ่งเป็นถนนที่มีข้อมูลระดับความสูงในบางจุดที่มีการสำรวจ ตัวอย่างอื่น ๆ ได้แก่ รูปร่างหลังคา - MultiLineString ที่บางประเด็นสำคัญมีการยกระดับที่กำหนดจากแผนอาคาร แต่ไม่ใช่ทั้งหมด การใช้ PostGIS คุณควรแนะนำให้ใช้แบบจำลองข้อมูลใดในการจัดเก็บข้อมูลประเภทนี้เพื่อให้สามารถเข้าถึงได้มากที่สุดโดยไม่สูญเสียหรือสร้างข้อมูลที่เกี่ยวข้อง

10 postgis  data  3d  model 

ฉันกำลังสร้างซีรี่ส์แผนที่ที่มีรูปหลายเหลี่ยมหลายแผ่นต่อแผ่น มีคำอธิบายประกอบสำหรับรูปหลายเหลี่ยม แต่รูปหลายเหลี่ยมบางส่วนมีขนาดเล็กมากและไม่มีคำอธิบายประกอบสำหรับพวกเขา (คำอธิบายประกอบจะซ่อนรูปหลายเหลี่ยมทั้งหมด) ดังนั้นฉันมีตารางแอตทริบิวต์ที่เก็บป้ายกำกับสำหรับรูปหลายเหลี่ยมเฉพาะเหล่านั้น ฉันจะให้ตารางนั้นแสดงบนแผนที่โดยมีเฉพาะค่าแอททริบิวต่อแผ่นมากกว่าค่า 2000+ ทั้งหมดได้อย่างไร ฉันกำลังมองหาโซลูชันที่เรียบง่ายและมีพลวัต

10 arcgis-desktop  arcgis-10.0  arcmap  attribute-table  legend 

เป็นไปได้ไหมที่จะตั้งค่า CRS โครงการของโครงการใน QGIS ผ่านคำสั่ง Python? ในกรณีของฉันผู้ใช้ควรเปิด QGIS เริ่มปลั๊กอินที่เขียนของฉันและปลั๊กอินนี้จะโหลดข้อมูลเชิงพื้นที่ซึ่งมีรหัส epsg 31467 (Gauss Krueger Zone 3) น่าเสียดายที่ CRS เริ่มต้นของ QGIS คือ WGS84 และฉันก็เลยมีปัญหากับยูนิตและสเกล วิธีแรกของฉันคือการใช้รหัสนี้: my_crs = core.QgsCoordinateReferenceSystem(31467, core.QgsCoordinateReferenceSystem.EpsgCrsId) self.iface.mapCanvas().mapRenderer().setDestinationCrs(my_crs) การเปลี่ยนแปลงนี้ (หลังจากข้อมูลถูกแสดงให้เห็นแล้ว) CRS เป็น Gauss Krueger Zone 3 อย่างน้อยโค้ด espg ที่มุมขวาล่างได้เปลี่ยนไป หน่วยยังคงอยู่ในระดับและขนาดผิดทั้งหมด 1: 5000 ใน Gauss Krügerคือ 1: 528822376 ใน WGS 84 …

10 qgis  python 

ฉันมีปัญหา / คำถามเกี่ยวกับหลามและ QGIS เป็นไปได้หรือไม่ที่จะได้รับข้อมูลของฟีเจอร์หนึ่งที่เลือกใน QGIS โดยใช้ python และถ้าเป็นเช่นนั้น ในตัวอย่างของฉันฉันมีเลเยอร์ซึ่งแสดงพัสดุบางส่วนสำหรับการวางผังเมือง ถ้าฉันเลือกหนึ่งในพัสดุที่ฉันต้องการรับข้อมูลสำหรับพัสดุตรงนี้ ข้อมูลนี้ควรอ่านโดยสคริปต์ไพ ธ อนจากนั้นเพิ่มลงในไฟล์ PDF ที่สร้างขึ้นโดยอัตโนมัติ ขั้นตอนสุดท้ายนั้นง่าย แต่ฉันไม่ได้รับข้อมูลจากคุณสมบัติที่เลือกนี้ เป็นไปได้ที่จะรับข้อมูลจากชั้นที่ใช้งานอยู่ แต่มันเป็นไปได้หรือไม่ที่จะได้รับข้อมูลจากคุณลักษณะที่เลือกหนึ่งของชั้นนี้ หวังว่าคุณจะช่วยฉัน ความนับถือ

10 qgis  python  pyqgis 

ด้วยรูปหลายเหลี่ยมหลายรูปที่ซ้อนทับกันหลายวิธีฉันต้องการส่งออกจากรูปหลายเหลี่ยมเหล่านี้ทั้งหมดที่ไม่ทับซ้อนกับรูปอื่น ๆ ซ้ำแล้วซ้ำอีก ผลิตภัณฑ์จะเป็นคุณสมบัติจำนวนหนึ่งโดยไม่มีการทับซ้อนกันซึ่งเมื่อรวมเข้าด้วยกันจะทำให้เป็นของแท้ ผลิตภัณฑ์นั้นสามารถนำมาใช้เป็นข้อมูลเข้าสู่ Zonal Statistics และมันจะเร็วกว่าการทำซ้ำ Zonal Statistics มากกว่ารูปหลายเหลี่ยมแต่ละรูป ฉันได้ลองรหัสนี้ใน ArcPy ไม่สำเร็จ รหัสสำหรับการทำเช่นนี้มีอยู่แล้ว?

10 arcpy  algorithm  overlapping-features 

ฉันได้เรียนรู้วิธีการเพิ่มปุ่มเครื่องมือลงในปลั๊กอิน - แถบเครื่องมือผ่าน python ตอนนี้ฉันสงสัยว่าจะเพิ่มแถบเครื่องมือที่สมบูรณ์ด้วยปุ่มแถบเครื่องมือผ่าน python ได้อย่างไร ใครสามารถให้โค้ดตัวอย่างได้บ้าง?

10 pyqgis  toolbar