Pracuję nad DEM w obszarze, w którym nie ma dużej ulgi. Chciałbym „wypalić” sieć rzeki do DEM, aby móc dokładnie obliczyć akumulację i długość przepływu. Korzystam z ArcGIS Desktop
Pracuję nad DEM w obszarze, w którym nie ma dużej ulgi. Chciałbym „wypalić” sieć rzeki do DEM, aby móc dokładnie obliczyć akumulację i długość przepływu. Korzystam z ArcGIS Desktop
czy możliwe jest podzielenie pliku shapefile na funkcję w Pythonie? (najlepiej byłoby rozwiązaniem, w którym mogę tymczasowo zapisać wynikowe obiekty wektorowe w pamięci zamiast na dysku). Powód: Chcę użyć funkcji gdal rasterizeLayer z kilkoma różnymi podzbiorami pliku kształtu. Ta funkcja wymaga...
Mam plik kształtu linii reprezentujący sieć dróg. Chcę zrasteryzować te dane, a wartości wynikowe w rastrze pokazują całkowitą długość linii, które mieszczą się w komórce rastra. Dane są w projekcji British National Grid, więc jednostkami będą metry. Idealnie chciałbym wykonać tę operację przy...
Próbuję wypalić plik kształtu do rastra przy użyciu RasterizeLayer GDAL. Wstępnie tworzę obszar zainteresowania rastrem z innego pliku kształtu, biorąc pod uwagę określony rozmiar piksela. Ten AOI służy następnie jako podstawa dla wszystkich kolejnych rasteryzacji (ta sama liczba kolumn i wierszy,...
Chcę utworzyć raster o rozdzielczości 25 metrów × 25 metrów, w którym każda komórka zawiera odległość do najbliższej linii brzegowej, obliczoną od środka komórki. Aby to zrobić, wszystko, co mam, to kształtna linia brzegowa Nowej Zelandii . Próbowałem wykonać samouczek Dominica Roye'a dotyczący...
W moim obszarze roboczym jest 100 SpatialPointsDataFrames. Załadowałem je w ten sposób: filenames <- list.files(path="", pattern="XYhectareTravelTimes_ez+.*shp") for(i in filenames){ filepath <- file.path("/",i) assign(i, readShapePoints(filepath)) Nazywa się je
Przekształcam wektor na raster w R. Jednak proces ten był zbyt długi. Czy istnieje możliwość włączenia skryptu w proces wielowątkowy lub przetwarzanie GPU, aby zrobić to szybciej? Mój skrypt do zrasteryzowanego wektora. r.raster = raster() extent(r.raster) = extent(setor) #definindo o extent do...