Coğrafi Bilgi Sistemleri

ArcGIS 9.3'te, çokgenlerle "maskelemek" istediğim bir raster dosyası var, çokgenlerle çakışan tüm raster hücrelerine tek bir değer veriyor. Bunu nasıl yapabilirim? İşte siyah çokgenlerin benim çokgen katmanım olduğu ve diğerlerinin raster tabanlı olduğu bir ekran görüntüsü.

14 arcgis-desktop  raster  polygon  arcgis-9.3  masking 

İki rasterim var: bir dijital yükseklik modeli (DEM) ve bir Landsat bandı. Her ikisinin ızgara çözünürlüğü 30 m'dir, ancak ızgaralar çakışmaz (resme bakın: üstteki karanlık tarama DEM'dir; altındaki ışık tarama Landsat bandıdır). Şebekenin DEM şebekesi ile mükemmel bir şekilde hizalanması için Landsat bandında bir dönüşüm yapmak istiyorum. Ortam ayarlarında bir …

14 arcgis-desktop  raster  alignment 

QGIS'de, belirli bir katmandaki çokgenlerin ve çizgilerin köşe sayısını saymanın ve bu değerleri içeren ayrı bir alan oluşturmanın bir yolu var mı? Saha hesap makinesinde bir yerde bir fonksiyon olduğundan şüpheleniyorum ama bulamıyorum.

14 qgis  vertices  count 

Verileri çizmek için bir veya daha fazla QGIS katmanına çok sayıda çoklu çizgi ve / veya çokgen almak istiyorum. Her birinin yalnızca basit nitelikleri olacaktır (örneğin ad). Bilmem gerek: 2. adımı etkinleştirmek için gereken CSV dosyasının biçimi. QGIS'de çok çizgi veya çokgen katmanı oluşturmak için gereken komut veya eklenti. Bilgi …

14 qgis  csv  import 

Python özniteliğine göre bir seçim tamamlamak çalışıyorum ama bir öznitelik bir listede olup olmadığını sorgusuna göre. Böyle bir sorgu en basit haliyle şöyle olmalıdır: qry = " \"OBJECTID\" in oid_list" arcpy.SelectLayersByAttribute_management(inft, "NEW_SELECTION", qry) ancak bu yaklaşım geçersiz bir ifade hatası döndürür. Geçmişte, bu tür bir sorgu için daha karmaşık sytax …

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

Python'da GDAL ile düzenli olarak kendi GeoTIFF rasterlerimi oluşturuyorum, örneğin: 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 ancak sonuç ArcCatalog / ArcGIS ile …

14 python  gdal  geotiff-tiff  statistics 

ArcGIS 10 , bir coğrafi veritabanında saklanan sınıfları öne çıkarmak için bireysel özelliklere ekler ekleme yeni özelliğine sahiptir . Bunun bir kullanımı, belirli bir konuma birden fazla resim eklemek olacaktır, yangın hidrantları popüler bir örnek gibi görünmektedir (sanırım internette bir köpek olup olmadığınızı bilmek önemlidir). Veritabanında, ekli bir özellik sınıfının …

14 arcgis-desktop  python  convert  esri-geodatabase  attachments 

PostGIS'te mesafeleri hesaplamakla ilgili basit bir sorum var. İki geometri arasındaki mesafeyi bulmak istiyorum. Bu sid: 4269'u metre olarak kullanıyorum, şu anda yaptığım şey şu: ST_Distance((a.geom,b.geom)) FROM ...ama derece olarak sonuç alıyorum. Coğrafya ile çalışmam gerektiğini düşünüyorum ama bir coğrafyayı nasıl bir coğrafyaya atabilirim? (A.geom :: coğrafya) ile denedim ama …

14 postgis  coordinate-system  postgresql  distance 

Bazı şekil dosyaları, kendisiyle ilişkilendirilmiş bir .prj dosyasına sahiptir ve .prj dosyası, WKT biçimindeki şekil dosyasının projeksiyon bilgilerini içerir. Bazen WKT'yi proj4 dizgisine dönüştürmem ve bazen de geri dönüştürmem gerekiyor. Bunu yapmak için hazır bir kütüphane var mı?

14 shapefile  coordinate-system  wkt  proj 

Ben kullanıyorum plotrixiçinde RABD devlet düzeyinde harita çizmek için. floating.pieHer duruma pasta grafikler koymak için mükemmel bir işlevi vardır . plotrixHer durumda çubuk grafik görüntülemek için paket içinde benzer işlevler olup olmadığını merak ediyorum ? (Belgelere bir göz attım ve çubuk grafiklerle uğraşan işlevler bu olasılığa sahip görünmüyor, ancak sadece …

14 r  statistics  chart 

Bir haritanın üstünde bir katmanım olan basit bir harita uygulaması oluşturdum. Leaflet'te katmanın zaten mevcut olup olmadığını nasıl kontrol edebilirim?

14 leaflet 

Leaflet'teki işaretçi simgesimin alt kısmında bir açılır pencere açmak istiyorum. kodum: 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 …

14 leaflet  popup 

Milyonlarca GPS noktası için bir ülke (ve bazen devlet) bulmak için aşağıdaki kodu kullanıyorum. Kod şu anda nokta başına yaklaşık bir saniye sürüyor, bu inanılmaz derecede yavaş. Şekil dosyası 6 MB'dir. Geopandaların uzaysal birleşimler için rtrees kullandıklarını ve onları inanılmaz derecede verimli hale getirdiğini okudum, ancak bu burada işe yaramıyor …

14 python  spatial-join  geopandas  rtree 

Küresel su kütleleri (1 bit değerleri 0 ve 1) ile büyük bir raster dosyası (129600 x 64800 piksel) var ve okyanus ve iç su kıyıları çıkarmaya çalışın. Rasterden çoklu çizgiye dönüştürmek için ArcGIS ve QGIS ile denedim, ancak zaman alır. Bu görev için daha iyi / daha hızlı bir yol …

14 python  r  performance 

PostgreSQL'e bir shp dosyası almak istiyorum. Önce sql dosyası oluşturup PostgreSQL'i çalıştırıyorum. Sql dosyası oluşturmak için bu komutu windows cmd'de çalıştırıyorum: shp2pgsql -s 4326 worldCountries.shp worldcountries postgres > worldcountries.sql ve sonra çalıştırın: psql -d postgres -U postgres -p 4321 -f worldcountries.sql ama sonuç: psql: worldcountries.sql: 21: HATA: işlev addGeometrycolumn (bilinmeyen, …

14 postgis  shapefile  postgresql  import