Pracuję nad obliczeniowym modelem liczebności dzikich zapylaczy w krajobrazie. Sam model jest kompletny, a teraz mam problem z etapem przetwarzania końcowego. Mam raster zapylacza GDAL, który wygląda mniej więcej tak (jaśniejsze kolory oznaczają większą wizytę zapylacza na piksel): I mam plik kształtu OGR punktów reprezentujących przykładowe lokalizacje w krajobrazie: Próbuję …
Mam funkcję wielokąta i chcę być w stanie generować w niej punkty. Potrzebuję tego do jednego zadania klasyfikacyjnego. Generowanie losowych punktów, dopóki jeden nie znajdzie się w wielokącie, nie zadziałałoby, ponieważ to naprawdę nieprzewidywalny czas, jaki zajmuje.
Szukam sugestii, jak ulepszyć mój kod Python. Zwykle wydajność nie ma dla mnie znaczenia, ale teraz pracuję z plikiem tekstowym lokalizacji w USA z ponad 1,5 miliona punktów. Przy danej konfiguracji uruchomienie operacji w jednym punkcie zajmuje około 5 sekund; Muszę obniżyć tę liczbę. Używam trzech różnych pakietów GIS Pythona, …
Mam skrypt w języku Python, który dodaje pole atrybutu do pliku Shapefile, jeśli nie istnieje. Jest to łatwe do zrobienia w ArcGIS (graficznie lub przez Python), ale szukam czegoś, co nie zależy od ArcGIS. Próbowałem tego bez powodzenia z OGR, ponieważ mój plik Shapefile zawiera funkcje . Patrzyłem na pyshp …
Interfejs API geobazy danych Esri może być teraz używany z GDAL / OGR - http://www.gdal.org/ogr/drv_filegdb.html, jednak wydaje się, że każda osoba musi się zalogować do Esri i osobno pobrać pliki API. Na przykład jest to wymienione w poście na blogu tutaj oraz w notatkach kompilacji GDAL . Umowa licencyjna przed …
Czytając warstwę z połączenia OGR PostGIS, mogę uzyskać parametr SpatialReference warstwy, ale czy można uzyskać wartość EPSG? Czy jest na to jakaś dokumentacja? Na przykład: lyr = conn.GetLayerByName(tbl) # Where conn is OGR PG connection srs = ly.GetSpatialRef() print srs Zwroty: PROJCS["OSGB 1936 / British National Grid", GEOGCS["OSGB 1936", DATUM["OSGB_1936", …
Chcę przypisać projekcję do pliku shapefile w mojej internetowej aplikacji C #, do tego używam powiązań Gdal / OGR / OSR C # i dodaj odwołania do osr_csharp.dll i ogr_csharp.dll. Ale napotykam pewne trudności podczas pisania kodu do przypisywania projekcji do pliku kształtu. Chcę wiedzieć, czy jest do tego dokumentacja …
Moim celem jest zmodyfikowanie istniejącego pliku kształtu przez scalenie niektórych wielokątów. Po zaimportowaniu pliku kształtu i użyciu polecenia UnionSpatialPolygons otrzymuję kontur wielokąta, który chcę. Jednak jest to teraz obiekt SpatialPolygons, a nie SpatialPolygonsDataFrame, więc nie jestem w stanie wyeksportować go do pliku kształtu za pomocą writeOGR. Jak mogę obejść ten …
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 …
Próbuję zmienić wartość pola pliku kształtu. Nie zmienia to jednak wartości. Co jest nie tak z moim kodem? import ogr driver = ogr.GetDriverByName('ESRI Shapefile') fn = 'dist.shp' dataSource = driver.Open(fn, 0) layer = dataSource.GetLayer() feature = layer.GetNextFeature() dist = 233 while feature: feature.SetField("dist", dist) layer.SetFeature(feature) feature = layer.GetNextFeature() dataSource.Destroy()
Obecnie używam następującego skryptu, aby dodać niektóre dane atrybutów z tabeli do wielu pojedynczych plików kształtów: library(rgdal) specieslist <- read.csv("SpeciesList1.txt", header=F) attdata <- read.table("TestAtt.csv", sep = ",", header=T) for (n in 1:dim(specieslist)[1]) { speciesname <- specieslist[n,1] shp <- readOGR("Mesoamerica_modified_polygons", speciesname) shp$ENGL_NAME<-attdata[n,2] writeOGR(shp, "PolygonsV2", speciesname, driver="ESRI Shapefile") } Na końcu pojawia …
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 = …
Piszę plik kształtu po pewnym przetworzeniu; jednak w ostatnim kroku muszę dokonać wyboru według atrybutów z pliku kształtu. Korzystam z polecenia w powłoce w celu użycia go w skrypcie Python, gdy działa. ogr2ogr -f "ESRI Shapefile" -select * where ID="1" outfile.shp infile.shp Otrzymuję komunikat o błędzie: FAILURE: Unable to open …
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.