Pytania otagowane jako python

Zacznę tworzyć programy do obsługi danych geograficznych z serwera PostGIS i będę używać GDAL. Zgodnie z twoim doświadczeniem, z najlepszym językiem do programowania? Mogę programować w Perlu i Pythonie, ale chciałbym znać język z większą ilością kodu do ponownego użycia, przykładami, bibliotekami itp.

20 python  postgis  gdal 

Rozwijam się w Pythonie i używam GDAL z OSGEO do manipulacji i interakcji z rastrami i plikami kształtów. Chcę wziąć plik kształtu, który ma cechy punktowe i interpolować go do rastra powierzchniowego. Obecnie używam metody „RasterizeLayer”, która wypala wartość z funkcji punktowej do rastra (która jest ustawiona dla wszystkich wartości …

20 python  raster  gdal  interpolation 

Muszę użyć Pythona, aby utworzyć dużą odległość od okręgu - zarówno liczbę, a najlepiej jakąś „krzywą”, której mogę użyć do narysowania mapy po stronie klienta. Nie dbam o format krzywej - czy to WKT, czy o zestaw par współrzędnych - ale chcę tylko pobrać dane. Jakie są tam narzędzia? Czego …

20 python  distance  polyline-creation  great-circle 

Mam 2 ramki geodata: import geopandas as gpd from shapely.geometry import Point gpd1 = gpd.GeoDataFrame([['John',1,Point(1,1)],['Smith',1,Point(2,2)],['Soap',1,Point(0,2)]],columns=['Name','ID','geometry']) gpd2 = gpd.GeoDataFrame([['Work',Point(0,1.1)],['Shops',Point(2.5,2)],['Home',Point(1,1.1)]],columns=['Place','geometry']) i chcę znaleźć nazwę najbliższego punktu w gpd2 dla każdego wiersza w gpd1: desired_output = Name ID geometry Nearest 0 John 1 POINT (1 1) Home 1 Smith 1 POINT (2 2) …

20 python  distance  geopandas  pandas 

Zamknięte . To pytanie musi być bardziej skoncentrowane . Obecnie nie przyjmuje odpowiedzi. Chcesz poprawić to pytanie? Zaktualizuj pytanie, aby skupiało się na jednym problemie, edytując ten post . Zamknięte 3 lata temu . Jestem zainteresowany opracowaniem aplikacji do mapowania stron internetowych. Stosunkowo dobrze znam Python, dlatego szukam frameworka Python. …

19 python  web-mapping  software-recommendations 

Szukam algorytmu, rozwiązania na wysokim poziomie, a nawet biblioteki, która może pomóc mi ustalić, czy dwa wielokąty przecinają się w Pythonie. Mam wierzchołki dwóch wielokątów (są to wielokąty jednoczęściowe bez otworów) w dwóch różnych układach. Wielokąty są 2D (tzn. Tylko współrzędne X i Y) Chciałbym utworzyć funkcję, która zwróci wartość …

19 python  algorithm  intersection 

Jak korzystać z bibliotek stron trzecich na wtyczkach QGIS w systemie Windows? Opracowałem wtyczkę, która używa rasterio i numpy dla klienta, ale ma on problemy z instalacją rasterio i numpy. W rzeczywistości rasterio i numpy zostały zainstalowane w jego głównym systemie Python (C: \ Python27), ale potrzebuję QGIS Python, aby …

18 python  pyqgis  installation  windows 

Chcę obliczyć statystyki ogniskowe dla każdej komórki rastra, w granicach określonych kryteriów. Tło - Mam trzy binarne rastry, z których każdy reprezentuje jeden interesujący rodzaj roślinności. Chciałbym obliczyć procent pokrycia każdego rodzaju roślinności w obrębie (np.) 20 km ^ 2 dowolnej komórki w moim obszarze badań (suma / całkowita liczba …

18 python  raster-calculator  spatial-statistics  shapely  fiona 

Próbuję napisać samodzielny skrypt Pythona w systemie Windows 7, który używa qgis.core i różnych innych bibliotek qgis. Znalazłem tutaj stare pytanie na ten temat: Pisanie samodzielnych skryptów Python przy użyciu PyQGIS? i postępował zgodnie ze swoją radą, korzystając z przewodnika tutaj: http://docs.qgis.org/testing/en/docs/pyqgis_developer_cookbook/intro.html#python-applications Kod sugeruje, że jest: set PYTHONPATH=c:\qgispath\python śledzony przez: …

18 qgis  python  pyqgis  windows 

W Pythonie, używając GDAL, wyodrębniłem projekcję rastra jako ciąg WKT w następujący sposób: wkt = dataset.GetProjection() # wkt is 'PROJCS["GDA_1994_Transverse_Mercator",GEOGCS["GDA_1994",DATUM["GDA_1994",SPHEROID["GRS_1980",6378137,298.2572221010002],TOWGS84[0,0,0,0,0,0,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",117],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",10000000],UNIT["Meter",1]]' Za pomocą ciągu WKT mogę utworzyć instancję SpatialReference w następujący sposób: src = osr.SpatialReference() src.ImportFromWkt(wkt) To jest łatwe. Mogę srcdość łatwo wyodrębnić kilka parametrów projekcji, takich jak strefa UTM itp. …

18 python  gdal  coordinate-system  python-2.7 

Mam bezzałogowy statek powietrzny (UAV, zwany także dronem) i chciałbym go użyć do stworzenia cyfrowego modelu elewacji i ortomozy na mniejszym obszarze (10 km2). Moja platforma będzie w stanie wytwarzać nakładające się zdjęcia na wysokości około 120 metrów. Wydaje się, że kilka produktów komercyjnych zapewnia kompletne rozwiązania; na przykład. Pix4D …

18 python  dem  unmanned-aerial-vehicle 

Zastanawiam się, czy istnieje biblioteka GIS w języku open source Python, która ma interfejsy API do obsługi połączeń WMS / WFS z innego serwera GIS (np. GeoServer), a następnie zapisuje dane odpowiedzi (mapa bazowa WMS i warstwa WFS) jako obrazy. jakieś rekomendacje? dzięki za wszelkie dane wejściowe! AKTUALIZACJA : to, …

18 python  open-source-gis  wms  wfs 

Ten post jest Wiki społeczności . Edytuj istniejące odpowiedzi, aby poprawić ten post. Obecnie nie przyjmuje nowych odpowiedzi. Zajmuję się tworzeniem kodu w Pythonie i chcę używać powiązań Pythona w GDAL / OGR. Jakie minimum muszę zainstalować, aby móc korzystać z powiązań pythonowych GDAL i OGR? Jaki jest najprostszy sposób …

18 python  gdal  ogr 

Mam plik kształtu z wielokątami. I mam globalny plik rastrowy. Chcę nałożyć wieloboki pliku kształtu na siatkę rastrową i obliczyć średnią wartość rastra dla każdego wielokąta. Jak mogę to zrobić za pomocą GDAL, zapisując wyniki w pliku shapefile?

18 python  raster  shapefile  gdal  spatial-statistics 

Mam trochę problemów z interfejsem API Python OGR. Staram się uzyskać wszystkie współrzędne każdego wierzchołka zewnętrznego pierścienia wielokąta. Oto co mam do tej pory: import osgeo.ogr import glob path = "/home/woo/maps/" out = path + 'output.txt' file = open(out,'w') for filename in glob.glob(path + "*.shp"): ds = osgeo.ogr.Open(filename) layer1 = …

18 python  open-source-gis  ogr