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

ฉันสังเกตเห็นว่าการคำนวณ heatmap ของ qgis เวอร์ชั่นใหม่นั้นได้รับแผนที่แรสเตอร์โดยไม่มีระบบพิกัดใด ๆ คุณมีปัญหาเดียวกันหรือไม่?

9 qgis  coordinate-system 

ฉันมีปัญหาในการค้นหา ฉันมีเครือข่ายของ linestrings แต่ละอันมีค่าในคอลัมน์ n_type นี่อาจเป็นหนึ่งในตัวเลือกมากมาย ฉันต้องการสร้างตารางใหม่ที่จัดกลุ่ม linestrings ใด ๆ ที่เป็นประเภทเดียวกันและสร้างบรรทัดต่อเนื่อง ก่อน: หลังจาก: นี่คือสิ่งที่ฉันได้รับจนถึง มันส่งคืนผลลัพธ์ แต่ไม่สมเหตุสมผล - ประเภทไม่ตรงกันและส่งคืนคุณสมบัติมากเกินไป โปรดทราบว่าฉันได้กำหนด "ต่อเนื่อง" เป็นบรรทัดใด ๆ ที่อยู่ในระยะ 5 ฟุตจากเพื่อนบ้านและพบกันที่มุมต่ำกว่า 30 องศา WITH RECURSIVE all_links (i, pk_uid, n_type, geom) AS ( SELECT 1 AS i, pk_uid, n_type, geom FROM network WHERE n_type != 'none' UNION …

9 postgis 

ฉันมีแผนที่ OpenLayers 3.2.0 ซึ่งมีบางแหล่งเวกเตอร์ ( ol.source.Vector) และเลเยอร์เวกเตอร์ที่เกี่ยวข้อง ( ol.layer.Vector) เมื่อเพิ่มคุณสมบัติ ( ol.Feature) ลงในแหล่งที่มาของเวกเตอร์คุณสมบัติเหล่านั้นจะได้รับdataคุณสมบัติที่ตั้งค่าเป็นวัตถุจาวาสคริปต์ที่คุณสมบัติดังกล่าวแสดง TypeScript มีดังนี้ ... vectorSource.addFeature(new ol.Feature({ geometry: /* ... */, data: vectorData, })); เลเยอร์เวกเตอร์นั้นมีฟังก์ชั่นสไตล์ที่อ่านdataคุณสมบัติและดึงสไตล์ของมัน: vectorLayer = new ol.layer.Vector({ source: vectorSource, renderBuffer: /* ... */, style: function (feature: ol.Feature, resolution: any) { var data = </* TypeScript Type */>feature.get('data'); …

9 openlayers-2 

ฉันเพิ่งเริ่มทำงานกับ Python 3.4.1 ซึ่งมาพร้อมกับ ArcGIS Pro 1.0.0 บนแล็ปท็อปของฉันซึ่งมี Python 2.7.8 ซึ่งมาพร้อมกับ ArcGIS 10.3 สำหรับเดสก์ท็อปที่ติดตั้งและฉันต้องการลดโอกาสในการทำงาน Python เวอร์ชันหนึ่งเมื่อฉันต้องการและจำเป็นต้องใช้อีกเวอร์ชันหนึ่ง ใน Windows Explorer ตอนนี้ฉันมีโฟลเดอร์ Python สองโฟลเดอร์: C: \ Python27 และ C: \ Python34 ซึ่งคาดว่าจะได้ ในเมนูเริ่มของ Windows ภายใต้โปรแกรมทั้งหมด> ArcGIS มีโฟลเดอร์ Python 2.7 ที่มี IDLE (Python GUI) ที่ฉันมักใช้เพื่อเริ่ม IDLE ดังนั้นการเริ่มต้นจากที่นี่ควรจะดี ฉันทดสอบและเริ่ม 2.7.8 ฉันคาดว่าจะเห็นโฟลเดอร์ Python 3.4 ถัดจาก …

9 arcpy  arcgis-pro  idle  python-3 

ฉันกำลังพยายามนำเข้าไฟล์ CSV ลงใน QGIS แต่ดูเหมือนว่าหน้าจอพิมพ์ของฉันด้านล่าง (โปรแกรมเป็นภาษาสวีเดน แต่ในหน้าต่างแสดงตัวอย่างคุณสามารถดูว่า QGIS อ่านข้อมูลของฉันไม่ถูกต้อง) ฉันไม่สามารถโหลดข้อมูลได้เนื่องจากเหตุผลบางอย่าง QGIS อ่านคอลัมน์ตาราง CSV ของฉันเป็นแถวเดียว เกิดอะไรขึ้น? ฉันต้องการให้ QGIS อ่านข้อมูลของฉันเป็นตัวเลขดังนั้นฉันจึงสร้างไฟล์ csvt เพิ่มเติมด้วยรหัส "String", "String", "Integer" มีอะไรอีกบ้างที่ฉันต้องทำเพื่อให้ QGIS อ่านคอลัมน์ที่สามเป็นตัวเลข? ถึงกระนั้นหลังจากทำอย่างนั้นมันก็แบบนี้ ...

9 qgis  csv  import  mac 

ฉันมีไฟล์ GeoJSON ชื่อ mygeojson.json และฉันต้องการเพิ่มเป็นเลเยอร์ใน OpenLayers 3 ด้านบนของเลเยอร์ openstreetmap จนถึงตอนนี้ฉันสามารถแสดงโลก openstreetmap รวมถึงการซูมเป็นต้น แต่ด้วยเหตุผลบางอย่างที่ฉันไม่สามารถรับ mygeojson.json ได้ geojson มีรูปหลายเหลี่ยมและมีลักษณะดังนี้: { "type": "FeatureCollection", "crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:OGC:1.3:CRS84" } }, "features": [ { "type": "Feature", "properties": { "DN": 2 }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 13.559093915055664, 52.545214330050563 …

9 javascript  openstreetmap  geojson  openlayers 

ฉันเลื่อนภาพ Landsat หนึ่งภาพบน ERDAS โดยใช้ทรัพยากรผสานความละเอียดพร้อมกับการแปลงรูป Brovey และเทคนิคเพื่อนบ้านที่ใกล้ที่สุด แต่มีคนบอกฉันว่าในแง่ของการจำแนกการใช้ที่ดินมันเป็นความผิดพลาดที่ทำเนื่องจากการปรับความคมชัดของภาพสร้างพิกเซลปลอมที่จะเพิ่มความผิดพลาด มันเป็นเรื่องจริงเหรอ?

9 landsat  erdas-imagine  land-classification 

ฉันต้องการโหลดไฟล์. adm ลงในอาร์ข้อมูลมาจากหน้านี้: http://www.fao.org/geonetwork/srv/en/metadata.show?id=14057 ฉันลองใช้รหัสต่อไปนี้ที่ฉันพบหลังจากการวิจัยในอินเทอร์เน็ต ปัญหาคือว่าในคลาส RasterLayer ฉันได้รับค่า negativ ที่ไม่ควรมี ฉันไม่รู้ว่าทำไมสิ่งนี้จึงเกิดขึ้นดังนั้นหวังว่าใครบางคนจะสามารถช่วยฉันได้! รหัส: library(rgdal) library(RColorBrewer) dpath<- path... x <- new("GDALReadOnlyDataset", dpath) getDriver(x) getDriverLongName(getDriver(x)) xx<-asSGDF_GROD(x) r <- raster(xx) ผลลัพธ์สำหรับ 'r' คือ: คลาส r: RasterLayer ขนาด: 2160, 4320, 9331200 (nrow, ncol, ncell) ความละเอียด: 0.08333333, 0.08333333 (x, y) ขอบเขต: -180, 180, -90, 90 (xmin, …

9 raster  data  r  esri-grid-format 

ก่อนอื่นฉันเป็นมือใหม่ใน GIS ฉันมีคลาวด์จุด lidar กับ GPS และข้อมูล telemetry ที่เกี่ยวข้อง จากนั้นฉันเปลี่ยนพิกัดของทุกจุดให้เป็นระบบพิกัด ECEF และในที่สุดฉันก็เปลี่ยนพิกัดเชิงมุมเป็นระบบพิกัดคาร์ทีเซียนโดยใช้การประมาณ Local Tangent Plane (LTP) (การประมาณนี้ใช้พิกัด GPS ของจุดแรกเป็นจุดกำเนิดของระบบอ้างอิงจากนั้นอ้างอิงจุดที่เหลือตามลำดับแรก) ฉันกำลังเขียนไฟล์ LAS โดยใช้ liblas และฉันไม่รู้ว่าจะตั้งค่าข้อมูลอ้างอิงทางภูมิศาสตร์เป็นไฟล์ LAS ได้อย่างไร ฉันกำลังเขียนพิกัดของทุกจุดในการประมาณ LTP ดังนั้นฉันจึงมี 2 คำถาม: มันถูกต้องหรือไม่ที่จะเขียนพิกัดจุดที่เกี่ยวข้องกับการประมาณ LTP หรือฉันควรจะเขียนพิกัดในหน่วย lat / lon (gps)? หากสิ่งนี้ถูกต้องฉันจะตั้งค่าการอ้างอิงทางภูมิศาสตร์ให้สอดคล้องกับ LTP ได้อย่างไร หาก '1' ไม่ถูกต้องและฉันต้องเขียนพิกัดในหน่วย lat / lon ฉันจะตั้งค่าการอ้างอิงทางภูมิศาสตร์ได้อย่างไร การเขียน VLR …

9 coordinate-system  georeferencing  lidar  latitude-longitude  liblas 

มีวิธีใดบ้างที่บอกให้ QGIS ใช้การติดตั้ง Python ที่แตกต่างจากค่าเริ่มต้นระบบทั่วทั้งระบบ? เช่นใช้ virtualenv Python บาง/myvirtualenv/bin/pythonตัวแทน/usr/bin/python? แบ็คกราวน์: ฉันกำลังพัฒนาปลั๊กอินที่กำหนดเองซึ่งใช้โมดูลบางตัวที่ติดตั้งใน virtualenv เท่านั้น แต่ไม่ใช่ทั้งระบบ virtualenv ได้รับอนุญาตให้เข้าถึงแพ็คเกจทั่วทั้งระบบดังนั้น qgis core libs ควรยังคงสามารถเข้าถึงได้เมื่อใช้ virtualenv นั้น ทางออกเดียวที่ฉันเห็นในตอนนี้คือการติดตั้งโมดูลที่จำเป็นใน Python ทั้งระบบ อย่างไรก็ตามฉันต้องการเก็บไว้ใน virtualenv เท่านั้นถ้าเป็นไปได้

9 qgis  qgis-plugins  pyqgis  qgis-2 

ฉันมีโครงการ QGIS พร้อมเลเยอร์ตามแบบสอบถาม PostGIS SQL ที่ใช้วันที่ปัจจุบัน เนื่องจากข้อมูลในฐานข้อมูลมีการเปลี่ยนแปลงเคียวรีเหล่านี้จึงส่งคืนข้อมูลที่แตกต่างกันทุกวัน เป็นไปได้ไหมที่จะเปลี่ยนชื่อของเลเยอร์ในเลเยอร์พาเนลแบบไดนามิกดังนั้นจึงเป็นวันที่ปัจจุบัน (เช่นชื่อเลเยอร์เช่นTraffic on 24.01.2015ที่เปลี่ยนทุกวัน) ข้อมูลสำหรับคำอธิบายภาพที่ควรจะนำมาจากแอตทริบิวต์ของคุณลักษณะใด ๆ ในชั้น - traffic_dateพวกเขาทุกคนมีค่าเท่ากันสำหรับเขตข้อมูล

9 qgis  pyqgis  layers  query-layer  dynamic-layer 

ฉันต้องการนำเข้าแรสเตอร์โมเสคหลาย ๆ อัน (แต่ละโมเสคประกอบด้วย 7 ชั้น) จากโฟลเดอร์ที่ระบุลงในอาร์ของฉันจากนั้นเข้าถึงพวกมันเป็นแรสเตอร์แบบหลายเลเยอร์เดียว ดังนั้นฉันจึง: # read all mosaics named "mos....img" in R raster_data<-list.files(path=getwd(), pattern="mos.*.img$") # read files as rasters s <- stack(raster_data) # check my imported rasters p.ex. raster n°8 from "s" raster stack s[[8]] และแรสเตอร์ของฉัน [[8]] มีเพียงชั้นเดียวเท่านั้นดังนั้นจึงไม่มีการนำเข้าโมเสคทั้งหมด! nlayers(s[[8]]) [[1]] ถ้าฉันอ่านแต่ละโมเสกแยกกันมันใช้งานได้: # read 1 mosaic (composed …

9 raster  r  remote-sensing  landsat 

ไม่มีใครทราบรายการความแตกต่างที่แม่นยำระหว่างสตริงรูปแบบการฉาย ESRI WKT และ OGC WKT หรือไม่ ฉันรู้ว่ามีเครื่องมือหลายอย่างที่ช่วยแปลงจาก ESRI WKT เป็น OGC WKT รวมถึงอรรถประโยชน์ของ GDAL และบริการเว็บไซต์ต่างๆ แต่คำถามของฉันไม่ใช่แบบที่ใช้งานได้จริงฉันเพียงต้องการเข้าใจความแตกต่างของการจัดรูปแบบ / ไวยากรณ์ที่บริการเหล่านี้ใช้งาน คำถาม Stackexchange ก่อนหน้านี้ได้พูดคุยเกี่ยวกับความแตกต่างในตัวอย่างเฉพาะเท่านั้นหรือเกี่ยวกับเครื่องมือและบริการที่มีให้เท่านั้น แม้ว่าคุณจะเพิ่งทราบความแตกต่างมันก็คงจะดีถ้าคุณเพิ่งโพสต์มัน จากประสบการณ์ของฉันเองควรมีความแตกต่างเล็กน้อย ความแตกต่างที่ฉันรู้คือ: รายการข้อความส่วนใหญ่ในการกำหนด esri ใช้ขีดล่างที่ ogc ใช้พื้นที่ ข้อความที่กำหนด datum ใน esri wkt เหมือนกับ ogc wkt ยกเว้นว่าเริ่มต้นด้วย "D_" บางครั้งตัวระบุข้อความสำหรับ PROJCS ที่กำหนดไว้ล่วงหน้า, PROJECTION, GEOGCS และ DATUM จะถูกเขียนต่างกัน (เช่น "NAD83" …

9 coordinate-system  wkt  ogc 

ฉันต้องการแสดงภาพซ้อนทับเชิงพื้นที่เพื่อระบุรูปหลายเหลี่ยมที่มีจุดตก อย่างไรก็ตามก่อนอื่นฉันต้องการบัฟเฟอร์และยุบรูปหลายเหลี่ยมดังกล่าวที่จุดตกที่ใดก็ได้ภายในรูปหลายเหลี่ยมผสาน (แต่ไม่ได้อยู่ในหลุมใด ๆ ) จะถูกติดป้ายคล้ายกันโดยขั้นตอนการซ้อนทับ น่าเสียดายที่กระบวนการบัฟเฟอร์และ / หรือการละลายที่ใช้อยู่นั้นเป็นการลดSpatialPolygonsวัตถุให้เป็นรูปหลายเหลี่ยม การใช้gBufferจะสร้างหลายรูปหลายเหลี่ยมเมื่อbyid=FALSEแต่ไม่สามารถรวมรูปหลายเหลี่ยมที่ทับซ้อนกันได้เมื่อbyid=TRUEใด ในกรณีหลังนี้การละลายรูปหลายเหลี่ยมด้วยการgUnaryUnionสร้างรูปหลายเหลี่ยมอีกครั้ง การวางซ้อนSpatialPointsกับผลลัพธ์หลายเหลี่ยมหลายเหลี่ยมเหล่านี้ในทุกจุดที่มีการรายงานว่าอยู่ในรูปหลายเหลี่ยม 1 นี่คือตัวอย่างของเล่นที่เกี่ยวข้องกับคะแนนบัฟเฟอร์ที่ฉันต้องการทำซ้อนทับคะแนนกับรูปหลายเหลี่ยมบัฟเฟอร์: library(sp) library(rgeos) pts <- SpatialPoints(cbind(c(1, 1, 2, 3), c(1, 2, 1.5, 2.5))) plot(gBuffer(pts, width=0.6), lwd=2) points(pts, pch=20, cex=2) text(coordinates(pts), labels=seq_len(length(pts)), pos=4, font=2) และผลลัพธ์ของการวางซ้อนบางอย่าง ... ด้วยbyid=FALSE: b <- gBuffer(pts, width=0.6) over(pts, b) # [1] 1 1 1 …

9 r  sp 

PyQGIS ตำราอธิบายถึงวิธีการตั้งค่าดัชนีเชิงพื้นที่ แต่เพียงอธิบายถึงครึ่งหนึ่งของการใช้งาน: สร้างดัชนีเชิงพื้นที่ - รหัสต่อไปนี้สร้างดัชนีที่ว่างเปล่า index = QgsSpatialIndex() เพิ่มคุณสมบัติในการทำดัชนี - ดัชนีใช้วัตถุ QgsFeature และเพิ่มลงในโครงสร้างข้อมูลภายใน คุณสามารถสร้างวัตถุด้วยตนเองหรือใช้จากการเรียกไปยังคุณลักษณะถัดไปของผู้ให้บริการ () index.insertFeature(feat) เมื่อดัชนีเชิงพื้นที่เต็มไปด้วยค่าบางอย่างคุณสามารถทำแบบสอบถามได้ # returns array of feature IDs of five nearest features nearest = index.nearestNeighbor(QgsPoint(25.4, 12.7), 5) ขั้นตอนที่มีประสิทธิภาพมากที่สุดในการรับคุณลักษณะจริงที่เป็นของรหัสคุณลักษณะที่ส่งคืนคืออะไร

9 qgis  pyqgis  spatial-index