Pytania otagowane jako gdal

Po użyciu gdalwarpdo wyświetlania i wyrównywania do rastra (przez -tap) pewnej liczby rastrów zauważyłem, że wyjściowe rastry były znacznie większe niż oryginalne rastry. Dość dokładne wyszukiwanie w sieci ujawniło ten problem z Trac : Frank Warmerdam wyjaśnił przyczynę: „Po dokładnym sprawdzeniu różnica w tym pliku jest taka, że ​​gdal_translate używa …

19 gdal  coordinate-system  geotiff-tiff  gdalwarp  file-size 

GDAL niedawno dodał nową funkcję, która pozwala na losowy odczyt plików segmentu S3. Chcę przyciąć obrazy GDAL z wielu kafelków obrazu bez konieczności pobierania całego pliku. Widziałem tylko bardzo rzadką dokumentację, jak skonfigurować i uzyskiwać dostęp do segmentu S3 przez GDAL i jestem trochę zdezorientowany, jak zacząć? Czy ktoś byłby …

19 raster  gdal  gdal-merge  amazon-ec2  amazon-s3 

Ok, trochę głupie pytanie, ale: Słyszałem, jak ludzie wymawiają bibliotekę GDAL (Geospatial Data Abstraction Library) jak goooo-dal. Moją naturalną skłonnością było wymawiać to geeee-dal, jak w literze „g”-dal. Jakieś wyjaśnienie, jak popularność goo-dal stała się popularna? Czy w grę wchodziła goo-gle?

18 gdal 

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 raster (właściwie USGS DEM) i muszę go podzielić na mniejsze części, jak pokazano na poniższym obrazku. Dokonano tego w ArcGIS 10.0 za pomocą narzędzia Split Raster. Chciałbym, aby to zrobić metoda FOSS. Patrzyłem na GDAL, myśląc, że na pewno by to zrobił (jakoś z gdal_translate), ale nic nie mogę …

18 raster  gdal  grass  dem  open-source-gis 

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 

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()

18 gdal  ogr  fields-attributes 

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 

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 …

18 shapefile  gdal  r  ogr  attribute-table 

W kodzie Python GDAL dość często widzę, że ludzie zamykają zestawy danych na końcu skryptu. Dlaczego warto zamykać zestaw danych w Pythonie GDAL? Czy są jakieś konsekwencje, jeśli tego nie zrobię? import gdal # open dataset ds = gdal.Open('test.tif') # close dataset ds = None

17 python  gdal 

Ponownie projektuję rastry w Pythonie za pomocą GDAL. Muszę rzutować kilka tiffów ze współrzędnych geograficznych WGS 84 na WGS 1984 Web Mercator (Sfera Pomocnicza), aby użyć ich później w Openlayers razem z OpenStreetMap i być może mapami Google. Korzystam z Pythona 2.7.5 i GDAL 1.10.1 stąd i przekształcam współrzędne za …

17 python  gdal  coordinate-system  epsg  web-mercator 

Przygotowuję dane do włączenia do węzła geonetwork, w którym mapy są wyświetlane za pośrednictwem geoserver. Problem polega na tym, że komórki, które nie mają danych, są wyświetlane jako geotiff w kolorze czarnym, a ja chciałbym, aby kolor był ustawiony na biały lub przezroczysty. Korzystam z ArcView i nie mam dostępu …

17 arcgis-desktop  geoserver  gdal  geotiff-tiff  geonetwork 

Próbuję przetworzyć niektóre bioklimatyczne pliki rastrowe, takie jak te, które można pobrać ze strony http://www.worldclim.org/current (zestaw bioklimatyczny). Wydaje się, że mają ustawione wartości nodata -3.4e+38zgodnie z QGIS (patrząc na wynik gdalinfo, to -3.39999999999999996e+38). Wydaje się, że narzędzia gdal nie są w stanie poradzić sobie z tą wartością nodata, a qgis …

17 qgis  raster  gdal 

Czy ktoś wie, jak najłatwiej zainstalować GDAL 2.0 na MacOSX? Zwykła podejrzana dystrybucja binarna kingchaos i homebrew wydaje się utknąć na GDAL 1.11 Potrzebuję wsparcia dla właściwości funkcji TopoJSON, które wydają się być zawarte w 2.0 EDYCJA 2017/01/07: Dystrybucja binarna kingchaos zapewnia teraz dystrybucję GDAL 2.1

16 gdal  installation  ogr  mac  homebrew 

Obecnie pracuję nad importem danych klimatycznych CANGRID (dostarczanych jako Surfer Grid ascii, pliki „.grd”) do ArcGIS. Siatka ma 95 rzędów na 125 kolumn. Metadane zapewniają lat / lon pochodzenia (lewy dolny róg), rozmiar komórki (50 km), a także odwzorowanie notatek jako stereograficzne biegunowe z południkiem środkowym (110 stopni W) i …

16 arcgis-desktop  python  gdal  coordinate-system