«python» etiketlenmiş sorular

Bazı enterpolasyonlar yapmak istediğim bir rasterim var. İşte buradayım: from osgeo import gdal from numpy import array # Read raster source = gdal.Open('my_raster.tif') nx, ny = source.RasterXSize, source.RasterYSize gt = source.GetGeoTransform() band_array = source.GetRasterBand(1).ReadAsArray() # Close raster source = None # Compute mid-point grid spacings ax = array([gt[0] + ix*gt[1] …

12 python  algorithm  interpolation 

Python API ile QGIS'de katman görünürlüğünü nasıl değiştirebilirim?

12 python  qgis 

Birden çok poligonun kesişimini almak istiyorum. Python'un shapelypaketini kullanarak, intersectionfonksiyonu kullanarak iki çokgenin kesişimini bulabilirim . Birden çok poligonun kesişimini elde etmek için benzer etkili bir fonksiyon var mı? İşte ne demek istediğimi anlamak için bir kod snippet'i: from shapely.geometry import Point coord1 = ( 0,0 ) point1 = Point(coord1) …

12 python  polygon  intersection  polygon-creation  shapely 

Bir Ubuntu 14.4 çalıştırıyorum ve bir python sanal ortamı kurdum. Sonra birkaç python kütüphanesi (gdal, lxml, psycopg2 vb.) Kurdum. Ayrıca sistemimde QGIS için modül var. Python yazdığımda bu modülü şu şekilde aktarabilirim: import qgis.core Ama aynı şeyi virtualenv üzerinden yapmaya çalıştığımda bir hata alıyorum: modül bulunamadı. Virtualenv'i oluştururken şu seçenekle …

12 qgis  python 

Düz İskelet algoritmasının uygulanmasını sağlayan bir Python paketi var mı? Açık kaynak kodlu (C ++) proje CGAL'ın bir uygulama içerdiğinin farkındayım, ancak cgal-bağlamaları bu CGAL paketini içermiyor gibi görünüyor . Her durumda, ihtiyaçlarıma göre değiştirebileceğim / genişletebileceğim saf bir Python uygulamasını tercih ederim. Delikli çokgenleri işleyebilen bir uygulama tercih edilirken, …

12 python  polygon  geometry  skeleton  computational-geometry 

Bir özelliğin bir özelliğini güncellemek istiyorum. Ancak, özellik nesnesini kullanarak güncelleştirmeyi bilmiyorum. Güncellemek için veri sağlayıcı kullanmak zorundayım. layers = QgsMapLayerRegistry.instance().mapLayersByName('my_line') layer = layers[0] dp = layer.dataProvider() it = dp.getFeatures() for i in range(0, dp.featureCount()): feat = it.next() attrs = { 2 : 30 } layer.dataProvider().changeAttributeValues({ feat.id() : attrs }) …

12 qgis  python  pyqgis  features 

Python Polygonize işlevini kullanarak sorun yaşıyorum. Bunun için yemek kitabı örneğini burada bulabilirsiniz . Kodumun ilgili kısmı: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile = "polygonized" driver = ogr.GetDriverByName("ESRI Shapefile") if os.path.exists(outShapefile+".shp"): driver.DeleteDataSource(outShapefile+".shp") outDatasource = driver.CreateDataSource(outShapefile+ ".shp") outLayer = outDatasource.CreateLayer("polygonized", srs=None) gdal.Polygonize( band, None, outLayer, -1, [], …

12 python  gdal 

Fiona'yı kullanarak bir şekil dosyası okuyorum: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] Bu bana koordinatların listesini içeren bir dikte verir ve 'MultiPolygon' yazın. Ancak, gerekli olduğu için Shapely MultiPolygon oluşturmak için koordinat verilerinin nasıl kullanılacağından emin değilim a sequence of exterior ring …

12 python  shapely  fiona 

GDAL Python bağlamaları kullanarak gdalbuildvrt yardımcı programıyla aynı görevi gerçekleştirmenin bir yolu var mı? Şimdiye kadar bunu tek bir veri kümesi vrt oluşturmak ve xml el ile düzenleme dışında herhangi bir yol bulamadım. Birden fazla rasterden (aslında bir mozaik gerçekleştirerek) bir vrt oluşturmak istiyorum. Saf Python kullanmak mümkün mü? Diğer …

12 python  raster  gdal  vrt  gdalbuildvrt 

[EDIT: Bunun çözümü, şekil dosyalarını okumak için OGR'yi kullanmaktı. Geographika örneğine bakın.] ESRI şekil dosyasında, Çokgenler ve Çokgenler arasında bir ayrım yoktur. Ayrıca, iç delikler ile dış halkalar arasında belirli bir ayrım yoktur (belirli bir çokgenin "elle tutulması" yanında). Bir şekil dosyasını okuduktan sonra, halkaları tanımlayan bir koordinat dizileri listem …

12 python  shapefile  polygon  shapely 

ülke alt bölümleri bir shapefile var ve her bölüm için lat ve lon sınırlayıcı bir dizi ayıklamak istiyorum .. mümkün mü?

12 python  shapefile  ogr 

Birden çok çizginin iki farklı coğrafi veri çerçevesi için bir çokgenin kesiştiği yeri bulmaya çalışıyorum: from shapely.geometry import Polygon, LineString import geopandas as gpd polygon = Polygon([(0, 0), (1, 0), (1, 1), (0, 1), (0, 0)]) line1 = LineString([(0.5, 0.5), (0.7, 0.7)]) line2 = LineString([(0.9, 0.9), (0.2, 0.6)]) poly_gdf = …

12 python  intersection  shapely  geopandas 

Python paketinin hangi sürücüleri fionavar? Kılavuzu kontrol ettiğimde [...] and the possible formats are enumerated in the fiona.drivers list. Ancak python yazdığımda yazıyor from fiona import drivers print drivers > <function drivers at 0x108763050> Buna nasıl "bakabilirim"?

12 python  fiona 

Raster katmanları NumPy ile işlemek için dizilere dönüştüren bir araç üzerinde çalışıyorum ve ideal olarak hepsini ihraç etmeden bir .gdb içinde gelen rasterlerle çalışabilmek istiyorum (özellikle bu ArcGIS veya ArcPy ile ilgileniyorsa) ). Bu OpenFileGDB sürücüsü ile mümkün mü? Bu sürücüye söyleyebileceğimden, raster katmanlarını vektör katmanlarıyla aynı şekilde ele alıyor …

12 python  gdal  file-geodatabase 

Gdal / OGR kullanarak bir shp, kml veya PostGIS çizgi katmanını düğümler ve segmentler ağına dönüştürmenin bir yolu vardır, böylece networkx gibi paketlerde kullanılabilir mi?

12 python  topology  ogr  networkx