Próbuję uchwycić wszystkie nie nakładające się wielokąty wskazane poniżej za pomocą Shapely (podane wielokąty A, B i C). Co więcej, mam nadzieję, że zrobię to bez iteracji, testowania przecięcia itp. Przyjęta odpowiedź na to pytanie wyraża metodę PostGIS, ale wydaje się, że „związek” oznacza różne...
Korzystam z danych GIS mojego hrabstwa do aktualizacji dróg w OpenStreetMap. Mam kopię, którą opublikowali w zeszłym roku, i jedną, która właśnie została opublikowana, i chciałbym znaleźć wszystkie LineStrings, które miały zmienione atrybuty lub geometrię. Pomoże mi to upewnić się, że wszystkie...
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. #...
Szukam innego, bardziej eleganckiego rozwiązania problemu statystyki przestrzennej. Surowe dane składają się ze współrzędnej xy dla każdego pojedynczego drzewa (tj. Przekonwertowanego na punktowy plik .shp). Chociaż nie zostało użyte w tym przykładzie, każde drzewo ma również odpowiedni wielokąt...
Krótko: chcę zmienić rozdzielczość rastra i wygładzić szarą rampę, jak pokazano na poniższych zdjęciach. Preferowane jest użycie GDAL, PIL lub Numpy. Opis: Za pomocą High Performance Geostatistic Library dzielę punkty na rastry o rozdzielczości wyjściowej 20 metrów . Nie chcę zmieniać...
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? --...
Czytałem ostatnio o kartogramach w tle i natknąłem się na interesującą alternatywę - mapy wartości według alfa . Koncepcja jest opisana bardziej szczegółowo tutaj, a pełny papier jest dostępny [nie otwarcie] [tutaj] 2 . Istnieje kilka szczegółów dotyczących implementacji na blogu indiemaps :...
Mam 2 pytania na temat analizy zestawu danych GPS. 1) Wyodrębnianie trajektorii Mam ogromną bazę danych zarejestrowanych współrzędnych GPS formularza (latitude, longitude, date-time). Zgodnie z wartościami daty i godziny kolejnych rekordów próbuję wyodrębnić wszystkie trajektorie / ścieżki, po...
Korzystam z poniższego kodu, aby znaleźć kraj (a czasem stan) dla milionów punktów GPS. Kod zajmuje obecnie około jednej sekundy na punkt, co jest niezwykle powolne. Plik shapefile ma 6 MB. Czytałem, że geopandy wykorzystują rtree do połączeń przestrzennych, co czyni je niesamowicie wydajnymi, ale...
Mam wielokąty ze współrzędnych w (kształtny python), które wyglądają tak POLYGON ((24.8085317 46.8512821, 24.7986952 46.8574619, 24.8088238 46.8664741, 24.8155239 46.8576335, 24.8085317 46.8512821)) Chciałbym obliczyć powierzchnię tego wielokąta w km ^ 2. Jaki byłby najlepszy sposób, aby to...
Mam duży plik rastrowy (129600 na 64800 pikseli) z globalnymi częściami wód (1 bit wartości 0 i 1) i próbuję wydobyć linie brzegowe oceanów i wód śródlądowych. Próbowałem z ArcGIS i QGIS przekonwertować z rastra na polilinię, ale to trwa wieki. Czy ktoś zna lepszy / szybszy sposób (Python lub R)...
Być może to trochę dziwne pytanie, ale pozwólcie, że wyjaśnię to krótko, zanim zacznę moje rzeczywiste pytania: Mikroskopia sił atomowych (AFM) to metoda, która w skrócie (i o mojej ograniczonej wiedzy) pozwala badaczom skanować obszary w mikro- i nanoskali. Działa poprzez „skanowanie” obszaru za...
Odwiedziłem i ponownie odwiedziłem stronę dotyczącą obsługi geometrii w książce kucharskiej PyQGIS: http://documentation.qgis.org/2.0/en/docs/pyqgis_developer_cookbook/geometry.html, ale wydaje się, że nie mogę dowiedzieć się, jak uzyskać wielokąt do rysowania z konsoli Python. Czy ktoś może...
Jak definiujemy połączenie SDE dla obszaru roboczego w skryptach Python?
Piszę proste narzędzie do przycinania partii wielopasmowych plików rastrowych geotiff do tego samego (mniejszego) obszaru. Korzystając z gdalwarp, mogę z łatwością przyciąć plik za pomocą pliku kształtowego przycinającego jeden wielokąt: gdalwarp -cutline clipper.shp -crop_to_cutline input.tif...
Zaczynam studia, które, mam nadzieję, doprowadzą do długiej kariery jako specjalista ds. Teledetekcji. Obecnie pracuję z ArcGIS dla niektórych aplikacji i uczę się ENVI dla innych. Zdałem sobie sprawę, że koniecznie muszę nauczyć się języka programowania i utknąłem w obliczu wyboru między IDL a...
Mam skrypt geoprzetwarzania, który uruchamia następujące ... arcpy.AddMessage("TEST...") Tak, tylko jeden wiersz kodu. Pusty plik .mxd zaczyna się dla mnie przy 391 KB, kiedy uruchamiam skrypt i zapisuję wzrost rozmiaru pliku. 391 kb 431 kb 451 kb 471 kb and so on... Wydaje się, że skrypt...
Regularnie tworzę własne rastry GeoTIFF z GDAL w Pythonie, np .: from osgeo import gdal from numpy import random data = random.uniform(0, 10, (300, 200)) driver = gdal.GetDriverByName('GTiff') ds = driver.Create('MyRaster.tif', 200, 300) band = ds.GetRasterBand(1) band.WriteArray(data) ds = band =...
Mam aplikację django z ulotką na froncie i muszę obliczyć bbox obiektu geojson, aby móc przekazać współrzędne do map.fitBounds() Próbowałem także niektórych bibliotek, a także próbowałem getBoundsmetod określonych obiektów w Ulotce, takich jak FeatureGroup, ale narzekałem, że granice nie zostały...
Przeprowadziłem badania nad użyciem R do dataminowania Twittera, ale tak naprawdę nie znalazłem odpowiedzi ani przyzwoitego samouczka na moje pytanie. Interesuje mnie pobieranie tweetów z Twittera z określonym hashtagiem w określonym czasie i wykreślanie lokalizacji tych tweetów na mapie w QGIS...