NumPy to moduł Pythona, który dodaje obsługę dużych, wielowymiarowych tablic i macierzy, a także dużą bibliotekę funkcji matematycznych wysokiego poziomu do działania na tych tablicach.
Jestem nowy w GIS. Mam kod, który konwertuje obrazy Marsa w podczerwieni na mapy bezwładności cieplnej, które są następnie zapisywane jako tablice dwuwymiarowe. Zapisuję te mapy jako pliki hdf5, ale naprawdę chciałbym je zapisać jako obrazy rastrowe, aby móc je przetwarzać w QGIS. Przeszukałem wiele wyszukiwań, aby dowiedzieć się, jak …
Próbowałem sprawdzić filtry w rastrze DEM pod kątem rozpoznawania wzorców i zawsze powoduje to brak ostatnich wierszy i kolumn (np ...20) . Próbowałem z biblioteką PIL, ładowanie obrazu. Następnie z numpy. Dane wyjściowe są takie same. Myślałem, że coś jest nie tak z moimi pętlami, podczas sprawdzania wartości w tablicy …
Próbuję ponownie rzutować / ponownie próbkować za pomocą powiązań Pythona GDAL, ale otrzymuję nieco inne wyniki w porównaniu do tych z narzędzia wiersza poleceń gdalwarp. Zobacz krótszy przykład aktualizacji Ten skrypt ilustruje podejście do Pythona: from osgeo import osr, gdal import numpy def reproject_point(point, srs, target_srs): ''' Reproject a pair …
Chciałbym nauczyć się korzystać z tablic NumPy w celu optymalizacji geoprzetwarzania. Duża część mojej pracy dotyczy „dużych zbiorów danych”, w których geoprzetwarzanie często zajmuje kilka dni, aby zrealizować określone zadania. Nie trzeba dodawać, że jestem bardzo zainteresowany optymalizacją tych procedur. ArcGIS 10.1 ma wiele funkcji NumPy, do których można uzyskać …
Zacząłem uczyć się, jak manipulować danymi LAS w pythonie i chciałem zobaczyć, jak inni obsługują pliki LAS. Chciałbym przeczytać punkty (używam tablicy numpy) i odfiltrować klasy 1 i 2 (niesklasyfikowane i uziemione) do osobnej tablicy. Mam następujący kod, ale nie mogę odfiltrować punktów. # Import modules from liblas import file …
Aby uzyskać wartość pojedynczego piksela w punkcie rastra za pomocą rasterio, istnieje przykład: https://github.com/mapbox/rasterio/pull/275 Czy jednak istnieje bezpośredni interfejs API w rasterio (a nie cli), którego można użyć do wyodrębnienia wartości w jednym punkcie rastra? -- EDYTOWAĆ with rasterio.drivers(): # Read raster bands directly to Numpy arrays. # with rasterio.open('C:\\Users\\rit\\38ERP.tif') …
Korzystam z ArcGIS 10.1 i chcę utworzyć nowy raster na podstawie dwóch istniejących rastrów. RasterToNumPyArray ma dobry przykład, który chcę, aby dostosować. import arcpy import numpy myArray = arcpy.RasterToNumPyArray('C:/data/inRaster') myArraySum = myArray.sum(1) myArraySum.shape = (myArray.shape[0],1) myArrayPerc = (myArray * 1.0)/ myArraySum newRaster = arcpy.NumPyArrayToRaster(myArrayPerc) newRaster.save("C:/output/fgdb.gdb/PercentRaster") Problem polega na tym, że …
W dokumentacji gdal_calc podano kalkulator rastrowy wiersza poleceń ze składnią numpy . Później jest kilka przykładów, w których w jednym z nich: gdal_calc.py -A input.tif --outfile = result.tif --calc = "A * (A> 0)" --NoDataValue = 0 - oznacza ustawione wartości od zera i poniżej na zero Niestety nie ma …
Mam dane rastrowe w tabeli postgres, którą chcę wprowadzić do Pythona jako tablicę numpy. Korzystam z psycopg2, aby połączyć się z bazą danych. Jestem w stanie pobrać dane, ale są one zwracane jako ciąg (prawdopodobnie serializowany plik binarny). Czy ktoś wie, jak pobrać ten ciąg i przekonwertować na tablicę numpy? …
Używam Pythona 2.6.5 (32-bitowego) z Numpy 1.3 i Gdal 1.9.1 zainstalowanym na Windows 7 64-bitowym. Próbuję odczytać zestaw danych rastrowych Imagine (.img) o pojemności 800 MB do tablicy Numpy, aby wykonać algebrę rastrową, ale gdy tylko uruchomię następujący kod, Python.exe ulega awarii. from osgeo import gdal g = gdal.Open(r'path\to\dataset', gdal.GA_Readonly) …
Próbuję nauczyć się obsługi przetwarzania obrazu za pomocą teledetekcji za pomocą powiązań i numpy w Pythonie GDAL. Jako pierwszą próbę czytam plik geotiff Landsat8, wykonuję prostą manipulację i zapisuję wynik w nowym pliku. Poniższy kod wydaje się działać dobrze, z tym wyjątkiem, że oryginalny raster jest zrzucany do pliku wyjściowego, …
Próbuję wyodrębnić naturalne wartości przerwania z rastra za pomocą Pythona. Zamierzony proces idzie: Arcpy obiekt rastrowy do NumPy Array (przy użyciu RasterToNumPyArray ) do naturalnych wartości przerwania (przy użyciu funkcji PySAL Natural Breaks ). import arcpy, pysal from pysal.esda.mapclassify import Natural_Breaks as nb # code to create greenIndex arcpy Raster …
Jak mogę uzyskać rzutowane współrzędne, a także rzeczywiste wartości pikseli na tych współrzędnych z pliku GeoTiff, a następnie zapisać je w tablicy liczbowej? Mam plik arsenci020l.tif, a jego współrzędne są w metrach. Poniżej znajduje się skrócona wersja gdalinfo, którą na nim uruchomiłem. ~$ gdalinfo arsenci020l.tif Driver: GTiff/GeoTIFF Files: arsenci020l.tif arsenci020l.tfw …
czy ktoś z powodzeniem nakłania pandy do pracy z instalacją QGIS w Pythonie w systemie Windows? Użyłem instalatora OSGeo4W do zainstalowania QGIS2.2, który jest dostarczany z Python 2.7 i Numpy 1.5.1. Każda rozsądnie nowa wersja pand wymaga Numpy 1.6.1 lub nowszej. To wydaje się moim głównym problemem. Próbowałem: - Budowanie …
Używamy plików cookie i innych technologii śledzenia w celu poprawy komfortu przeglądania naszej witryny, aby wyświetlać spersonalizowane treści i ukierunkowane reklamy, analizować ruch w naszej witrynie, i zrozumieć, skąd pochodzą nasi goście.
Kontynuując, wyrażasz zgodę na korzystanie z plików cookie i innych technologii śledzenia oraz potwierdzasz, że masz co najmniej 16 lat lub zgodę rodzica lub opiekuna.