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

ฉันกำลังพยายามเลือกตามคุณลักษณะใน Python ให้เสร็จสิ้น แต่ขึ้นอยู่กับแบบสอบถามว่ามีแอตทริบิวต์อยู่ในรายการหรือไม่ แบบสอบถามที่ง่ายที่สุดควรเป็นดังนี้: qry = " \"OBJECTID\" in oid_list" arcpy.SelectLayersByAttribute_management(inft, "NEW_SELECTION", qry) แต่วิธีการนั้นจะส่งกลับข้อผิดพลาดการแสดงออกที่ไม่ถูกต้อง ก่อนหน้านี้ฉันต้องใช้ระบบที่ซับซ้อนมากขึ้นสำหรับการสืบค้นประเภทนี้เช่น: sqlQuery2 = "nid in (" + ','.join(["'"+x+"'" for x in delta_list]) +")" แต่การปรับตัวอย่างโค้ดนี้ก็ไม่ได้ผลสำหรับฉันเช่น: "OBJECTID_1 in (" + ','.join(["'"+str(x)+"'" for x in oid_list]) +")" ฉันหายไปนี่อะไร

14 arcpy  arcgis-10.0  sql  select-by-attribute  where-clause 

ฉันสร้าง raster GeoTIFF ของตัวเองเป็นประจำกับ GDAL ใน Python เช่น: from osgeo import gdal from numpy import random data = random.uniform(0, 10, (300, 200)) driver = gdal.GetDriverByName('GTiff') ds = driver.Create('MyRaster.tif', 200, 300) band = ds.GetRasterBand(1) band.WriteArray(data) ds = band = None # save, close อย่างไรก็ตามเมื่อดูผลลัพธ์ด้วย ArcCatalog / ArcGIS มันจะเป็นสีดำหรือเทาเนื่องจากไม่มีสถิติ สิ่งนี้สามารถแก้ไขได้ด้วยการคลิกขวาที่แรสเตอร์และเลือก "คำนวณสถิติ ... …

14 python  gdal  geotiff-tiff  statistics 

ArcGIS 10 มีคุณสมบัติความสามารถใหม่ในการเพิ่มสิ่งที่แนบมากับคุณสมบัติแต่ละอย่างเพื่อจัดเก็บคลาสในฐานข้อมูลภูมิศาสตร์ การใช้สิ่งนี้เพียงอย่างเดียวคือการเพิ่มรูปภาพหลายภาพไปยังสถานที่เฉพาะเจาะจงhydrants ไฟดูเหมือนจะเป็นตัวอย่างที่นิยม (สำคัญที่จะต้องรู้ว่าคุณเป็นสุนัขบนอินเทอร์เน็ตในความต้องการของการหยุดหลุมฉันคิดว่า) ในฐานข้อมูลโครงสร้างของคลาสคุณลักษณะที่เปิดใช้งานสิ่งที่แนบมาคือสามเท่า: คลาสฟีเจอร์เอง ตารางที่เก็บสิ่งที่แนบและเมตาดาต้าที่เกี่ยวข้อง (ชื่อไฟล์, ขนาด, ฯลฯ ) และ คลาสความสัมพันธ์ที่กำหนดความสัมพันธ์แบบ 1 ต่อหลายคนระหว่างทั้งสอง คำถามของฉันคือ: วิธีการส่งออกข้อมูลนี้ - คลาสคุณลักษณะ, สิ่งที่แนบ, ข้อมูลเมตาของสิ่งที่แนบมา - เพื่อการบริโภคนอก Arcgis? * การส่งออก FC นั้นเป็นเรื่องง่ายFeatureClassToFeatureClassทำงานเหมือนกับบน FC ที่ไม่ได้เปิดใช้งาน กำลังส่งออกตารางสารบัญหรือไม่ ไม่มากนัก. TableToTableไปยัง output.dbf แยกเมตะดาต้าของไฟล์ที่แนบมาเท่านั้นไปยังข้อผิดพลาด output.csv ด้วยประเภทเขตข้อมูลที่ไม่ได้รับการสนับสนุนและไปที่ตารางข้อมูลผลลัพธ์ล้มเหลวพร้อมข้อผิดพลาดฟิลด์ทั่วไป นี่ไม่ใช่สิ่งที่น่าแปลกใจเนื่องจากไฟล์ประเภทเหล่านั้นไม่สนับสนุนประเภทข้อมูลไบนารี blob ฉันคาดหวัง / หวังว่าจะพบเครื่องมือที่จะแปลงไบนารีที่แนบมาเป็นรูปแบบดั้งเดิมเช่น แล้วมันเป็นอย่างไร ฉันต้องทำอะไรเพื่อให้ได้ข้อมูลออกมา * แดกดันโครงการของฉันในช่วงเวลาที่วางไข่ Q นี้คือการส่งออกเพื่อใช้ภายใน …

14 arcgis-desktop  python  convert  esri-geodatabase  attachments 

ฉันมีคำถามง่ายๆเกี่ยวกับการคำนวณระยะทางใน PostGIS ฉันอยากได้ระยะห่างระหว่างรูปทรงเรขาคณิตสองอัน ฉันกำลังใช้ sid นี้: 4269 หน่วยเป็นเมตรสิ่งที่ฉันกำลังทำอยู่ตอนนี้คือ: ST_Distance((a.geom,b.geom)) FROM ...แต่ฉันได้ผลลัพธ์เป็นองศา ฉันคิดว่าฉันควรทำงานกับภูมิศาสตร์ แต่ฉันจะขว้าง geom ไปยังภูมิศาสตร์ได้อย่างไร ฉันลองด้วย (a.geom :: ภูมิศาสตร์) แต่ฉันได้รับข้อผิดพลาด ฉันจะทำอย่างไรเพื่อให้ได้ผลลัพธ์เป็นเมตร ขอขอบคุณ

14 postgis  coordinate-system  postgresql  distance 

บาง Shapefiles มีไฟล์. prj ที่เชื่อมโยงกับมันและไฟล์. prj มีข้อมูลการฉายของ Shapefile ในรูปแบบของ WKT บางครั้งฉันต้องแปลง WKT เป็นสตริง proj4 และบางครั้งฉันต้องแปลงกลับ มีห้องสมุดสำเร็จรูปที่จะทำเช่นนี้หรือไม่?

14 shapefile  coordinate-system  wkt  proj 

ฉันใช้plotrixในRการพล็อตของแผนที่ระดับรัฐของสหรัฐ มันมีฟังก์ชั่นที่ยอดเยี่ยมfloating.pieในการวางแผนภูมิวงกลมในแต่ละรัฐ ฉันสงสัยว่ามีฟังก์ชั่นที่คล้ายกันภายในplotrixแพ็คเกจเพื่อแสดงแผนภูมิแท่งในแต่ละรัฐหรือไม่? (ฉันลองดูที่เอกสารประกอบและฟังก์ชั่นที่จัดการกับแผนภูมิแท่งดูเหมือนจะไม่มีความเป็นไปได้ แต่ฉันแค่อยากจะแน่ใจ) ฉันชอบทำงานในแพ็คเกจพล็อตริกซ์ถ้าเป็นไปได้ แต่รู้สึกอิสระที่จะตั้งชื่อ แพคเกจอื่น ๆ ตัวอย่างเช่นฉันสนใจที่จะสร้างแผนที่คล้ายกับแผนที่นี้ (แต่สำหรับสหรัฐอเมริกา): สำหรับแผนที่อเมริกาของฉันจะมีแผนภูมิแท่ง 50 แผนภูมิแต่ละแผนภูมิ ฉันได้รับแผนที่นี้จาก/programming/20465070/barplots-on-a-mapแต่ดูเหมือนว่า ggsubplot ไม่ทำงานในเวอร์ชัน R ของฉัน (คล้ายกับที่คนอื่นพูดในโพสต์ )

14 r  statistics  chart 

ฉันสร้างแอพแผนที่อย่างง่ายที่มีเลเยอร์อยู่ด้านบนของแผนที่ ในแผ่นพับฉันจะตรวจสอบว่ามีเลเยอร์อยู่แล้วได้อย่างไร

14 leaflet 

ฉันต้องการเปิดป๊อปอัพที่ด้านล่างของไอคอนเครื่องหมายของฉันใน Leaflet รหัสของฉัน: var mymap = L.map('mapid').setView([51.505, -0.09], 13); L.tileLayer('https://api.tiles.mapbox.com/v4/{id}/{z}/{x}/{y}.png?access_token=pk.eyJ1IjoibWFwYm94IiwiYSI6ImNpandmbXliNDBjZWd2M2x6bDk3c2ZtOTkifQ._QA7i5Mpkd_m30IGElHziw', { maxZoom: 18, attribution: 'Map data © <a href="http://openstreetmap.org">OpenStreetMap</a> contributors, ' + '<a href="http://creativecommons.org/licenses/by-sa/2.0/">CC-BY-SA</a>, ' + 'Imagery © <a href="http://mapbox.com">Mapbox</a>', id: 'mapbox.streets' }).addTo(mymap); var lati = 51.51; var longi = -0.09; var popupLocation1 = new L.LatLng(lati, longi); var popupContent1 = …

14 leaflet  popup 

ฉันใช้รหัสด้านล่างเพื่อค้นหาประเทศ (และบางครั้งระบุ) สำหรับจุด GPS นับล้าน รหัสในปัจจุบันใช้เวลาประมาณหนึ่งวินาทีต่อจุดซึ่งช้ามากอย่างไม่น่าเชื่อ ไฟล์รูปร่างคือ 6 MB ฉันอ่านว่า geopandas ใช้ rtrees สำหรับการรวมเชิงพื้นที่ทำให้พวกมันมีประสิทธิภาพอย่างเหลือเชื่อ แต่สิ่งนี้ดูเหมือนจะไม่ทำงานที่นี่ ผมทำอะไรผิดหรือเปล่า? ฉันหวังเป็นพันคะแนนต่อวินาทีหรือมากกว่านั้น สามารถดาวน์โหลด shapefile และ csv ได้ที่นี่ (5MB): https://www.dropbox.com/s/gdkxtpqupj0sidm/SpatialJoin.zip?dl=0 import pandas as pd import geopandas as gpd from geopandas import GeoDataFrame, read_file from geopandas.tools import sjoin from shapely.geometry import Point, mapping,shape import time #parameters shapefile="K:/.../Shapefiles/Used/World.shp" …

14 python  spatial-join  geopandas  rtree 

ฉันมีไฟล์แรสเตอร์ขนาดใหญ่ (129,600 คูณ 64,800 พิกเซล) ที่มีแหล่งน้ำทั่วโลก (1 บิตค่า 0 และ 1) และพยายามที่จะแยกมหาสมุทรและชายฝั่งน้ำจืด ฉันได้ลองกับ ArcGIS และ QGIS เพื่อแปลงจาก raster เป็น polyline แต่ใช้เวลานาน มีใครรู้วิธีที่ดีกว่า / เร็วกว่า (Python หรือ R) หรือเป็นเครื่องมือที่ดีกว่าสำหรับงานนี้หรือไม่? ปรับปรุง R: rasterToContour อาจเร็วและแม่นยำ แต่ถ้าคุณมีชุดข้อมูลขนาดใหญ่มากเช่นฉัน (8,398,080,000 พิกเซล) คุณต้องมี RAM จำนวนมาก (มากกว่า 16GB) หรือบังคับ R เพื่อทำการประมวลผลเพิ่มเติมบนฮาร์ดไดรฟ์และมัน จะใช้เวลานานด้วย Python / GDAL: gdal_poligonize สร้างรูปหลายเหลี่ยมแทนโพลีน อัปเดต …

14 python  r  performance 

ฉันต้องการนำเข้าไฟล์ shp ไปยัง postgreSQL ก่อนอื่นฉันสร้างไฟล์ sql แล้วเรียกใช้ PostgreSQL เพื่อสร้างไฟล์ sql ฉันเรียกใช้คำสั่งนี้ใน windows cmd: shp2pgsql -s 4326 worldCountries.shp worldcountries postgres > worldcountries.sql จากนั้นเรียกใช้: psql -d postgres -U postgres -p 4321 -f worldcountries.sql แต่ผลลัพธ์คือ: psql: worldcountries.sql: 21: ข้อผิดพลาด: ฟังก์ชัน addGeometrycolumn (ไม่ทราบ, ไม่ทราบ, ไม่ทราบ, ไม่ทราบ, ไม่ทราบ, จำนวนเต็ม)

14 postgis  shapefile  postgresql  import 

ฉันกำลังมองหาวิธีในการสร้างแผนภาพ voronoi ตามส่วนของเส้นแทนที่จะเป็นจุด ดูตัวอย่างด้านล่าง (นำมาจากการนำเสนอนี้ ) เป็นการดีที่ฉันต้องการบางสิ่งบางอย่างที่ฉันสามารถเขียนสคริปต์จาก Python แต่โซลูชันที่ใช้ ArcGIS หรือที่คล้ายคลึงกันก็เป็นที่ยอมรับได้เช่นกัน ห้องสมุดแห่งเดียวที่ทำสิ่งนี้ที่ฉันได้ค้นพบในตอนนี้คือopenvoronoiซึ่งดูมีแนวโน้ม มีคนอื่นบ้างไหม?

14 line  voronoi-thiessen 

ฉันได้สร้างวัตถุชนิดSpatialPointsDataFrameโดยใช้spแพคเกจใน R อย่างไรก็ตามฉันสับสนเกี่ยวกับตัว@, $, . and []ดำเนินการและเมื่อใดจึงใช้เพื่อเข้าถึงคุณสมบัติที่แตกต่างกันของวัตถุของฉัน นี่คือตัวอย่างรหัสของฉัน: library(sp) library(rgdal) #creating a SpatialPointsDataFrame with sample points in UTM x <- c(15.2, 15.3, 15.4, 15.5, 15.7) y <- c(50.4, 50.2, 50.3, 50.1, 50.4) v1 <- c(1.0, 2.0, 3.0, 4.0, 5.0) v2 <- c("a","b","b","c","a") attributes <- as.data.frame(cbind(v1,v2)) xy <- cbind(x,y) locationsDD <- …

14 r  syntax 

ฉันมี 23 ตารางใน schema ของ PostGIS ซึ่งฉันจำเป็นต้องนับจำนวนจุดยอดของ ตารางเป็นส่วนผสมของเส้นและรูปหลายเหลี่ยมดังนั้นฉันจึงตระหนักว่าต้องใช้ ST_NPoints (geom) ดังนั้นฉันจึงเรียกใช้แบบสอบถามต่อไปนี้ SELECT count(ST_NPoints(geom) FROM lines; ผลลัพธ์ / การนับเท่ากับจำนวนของคุณสมบัติในตารางนั้นและไม่ใช่จำนวนยอดรวมของคุณลักษณะทั้งหมดในตารางนั้น ฉันจะต้องหายไปบางอย่าง แต่ไม่สามารถคิดออก (ต้องเช้าวันจันทร์;)

14 postgis 

โปรแกรม ArcGIS ของแพลตฟอร์ม ESRI พร้อมใช้งานบนUbuntuหรือไม่? ฉันรักโอเพ่นซอร์สและถ้าฉันสามารถลองและย้ายไปยังโอเพ่นซอร์สระบบปฏิบัติการนี้ได้

14 arcgis-desktop  arcgis-server  ubuntu  arcgis-platform