Pytania otagowane jako python

Python to język programowania z interpretacją open source wykorzystywany w wielu programach GIS.

4
Dwuliniowa interpolacja danych punktowych na rastrze w Pythonie?
Mam raster, z którym chciałbym wykonać interpolację punktową. Oto gdzie jestem: 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] …
3
Skuteczne uzyskiwanie przecięcia wielu wielokątów w Pythonie
Chciałbym uzyskać przecięcie wielu wielokątów. Korzystając z shapelypakietu Pythona , mogę znaleźć przecięcie dwóch wielokątów za pomocą intersectionfunkcji. Czy istnieje podobna skuteczna funkcja uzyskiwania przecięcia wielu wielokątów? Oto fragment kodu, aby zrozumieć, co mam na myśli: from shapely.geometry import Point coord1 = ( 0,0 ) point1 = Point(coord1) circle1 = …
1
Jak uzyskać dostęp do modułu QGIS z mojego virtualenv
Korzystam z systemu Ubuntu 14.4 i zainstalowałem środowisko wirtualne Python. Potem zainstalowałem kilka bibliotek Pythona (gdal, lxml, psycopg2 itp.) Również w moim systemie mam moduł QGIS. Kiedy więc wpisuję python, mogę zaimportować ten moduł jako: import qgis.core Ale kiedy próbuję zrobić to samo za pośrednictwem virtualenv, pojawia się błąd: moduł …
12 qgis  python 
2
Jak zmienić wartość atrybutu za pomocą QgsFeature w PyQGIS?
Chciałbym zaktualizować jeden atrybut funkcji. Nie wiem jednak, jak go zaktualizować za pomocą obiektu funkcji. Muszę użyć dostawcy danych, aby go zaktualizować. 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({ …
12 qgis  python  pyqgis  features 
1
GDAL poligonizuje w pythonie tworząc pusty wielokąt?
Mam problem z użyciem funkcji Polygonize w Pythonie. Przykład książki kucharskiej na ten temat można znaleźć tutaj . Odpowiednia część mojego kodu to: 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) …
12 python  gdal 
2
Tworzenie zgrabnych wielokątów z kształtnego wielokąta
Czytam plik kształtu za pomocą Fiony: with fiona.open('data/boroughs/boroughs_n.shp') as source: mpolys = [p for p in source] candidate = polys[0]['geometry'] daje mi to dyktando zawierające listę list współrzędnych i wpisz „MultiPolygon”. Nie jestem jednak pewien, jak użyć danych współrzędnych, aby utworzyć Shapely MultiPolygon, ponieważ wymaga a sequence of exterior ring …
12 python  shapely  fiona 
3
Odpowiednik gdalbuildvrt w języku Python
Czy istnieje sposób wykonania tego samego zadania, co narzędzie gdalbuildvrt przy użyciu powiązań GDAL Python? Do tej pory nie znalazłem innego sposobu niż utworzenie vrt pojedynczego zestawu danych i ręczna edycja pliku XML. Chciałbym utworzyć vrt z wielu rastrów (zasadniczo wykonujących mozaikę). Czy jest to możliwe przy użyciu czystego Pythona? …
4
przy użyciu foremnego: tłumaczenia między wielokątami i wielokątami
[EDYCJA: rozwiązaniem tego było po prostu użycie OGR do odczytu plików kształtów. Zobacz przykład geographika.] W pliku kształtu ESRI nie ma rozróżnienia między wielokątami i wielokątami. Co więcej, nie ma wyraźnego rozróżnienia między wewnętrznymi otworami a zewnętrznymi pierścieniami (poza „przydatnością” danego wielokąta). Po przeczytaniu pliku kształtu mam listę sekwencji współrzędnych …
1
Kierowcy Fiony
Jakie sterowniki ma pakiet Python fiona? Kiedy sprawdzam instrukcję, to [...] and the possible formats are enumerated in the fiona.drivers list. jednak napisałem w Pythonie from fiona import drivers print drivers > <function drivers at 0x108763050> Jak mogę na to spojrzeć?
12 python  fiona 
2
Uzyskujesz dostęp do rastrów geobazy danych (.gdb) za pomocą GDAL?
Pracuję nad narzędziem, które konwertuje warstwy rastrowe na tablice do przetwarzania za pomocą NumPy, i idealnie chciałbym móc pracować z rastrami dostarczanymi w formacie .gdb bez eksportowania ich wszystkich (szczególnie jeśli wymaga to zaangażowania ArcGIS lub ArcPy ). Czy jest to możliwe dzięki sterownikowi OpenFileGDB? Z tego, co mogę powiedzieć, …
2
Python nie może odczytać pliku kształtu
Próbuję otworzyć plik kształtu za pomocą notatnika Jupyter, ale nie mogę go otworzyć. Plik shapefile znajduje się w tym samym folderze co mój skoroszyt, a Python nie może go zidentyfikować, nawet jeśli użyję pełnej ścieżki. import shapefile as shp import matplotlib.pyplot as plt sf = shp.Reader(r'C:\Users\Public\Documents\1-11-99n.shp') ShapefileException: Unable to open …
11 python  shapefile 
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.