Вопросы с тегом «gdal»

После использования gdalwarpдля проецирования и выравнивания по сетке (через -tap) нескольких растров я заметил, что выходные растры были значительно больше, чем исходные растры. Достаточно тщательный поиск в Интернете выявил эту проблему Trac : Фрэнк Вармердам объяснил причину: «При тщательном рассмотрении разница в рассматриваемом файле заключается в том, что gdal_translate использует …

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

Итак, GDAL недавно добавил новую функцию, которая позволяет произвольно читать файлы S3 Bucket. Я хочу обрезать изображения GDAL из нескольких фрагментов изображения, не загружая весь файл. Я только видел очень скудную документацию о том, как настроить и получить доступ к корзине S3 через GDAL, и я немного запутался, как начать? …

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

Хорошо, немного глупый вопрос, но: Я слышал, как люди произносят библиотеку GDAL (Geospatial Data Abstraction Library) как goooo-dal. Моя естественная склонность состояла в том, чтобы произнести это «гее-дал», как в письме «г'-дал». Любое объяснение того, как Goo-Dal стал популярным? Был ли вовлечен Google?

18 gdal 

В Python, используя GDAL, я извлек проекцию растра в виде строки WKT следующим образом: 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]]' Используя строку WKT, я могу создать экземпляр SpatialReference следующим образом: src = osr.SpatialReference() src.ImportFromWkt(wkt) Это легко-peasy. Я могу достаточно легко извлечь несколько параметров проекции из srcзоны UTM и т. …

18 python  gdal  coordinate-system  python-2.7 

У меня есть растр (на самом деле USGS DEM), и мне нужно разделить его на более мелкие фрагменты, как показано на рисунке ниже. Это было достигнуто в ArcGIS 10.0 с использованием инструмента Split Raster. Я хотел бы метод FOSS, чтобы сделать это. Я посмотрел на GDAL, думая, что он обязательно …

18 raster  gdal  grass  dem  open-source-gis 

Этот пост является сообществом Wiki . Отредактируйте существующие ответы, чтобы улучшить этот пост. В настоящее время он не принимает новые ответы. Я занимаюсь разработкой некоторого кода на Python и хочу использовать привязки Python GDAL / OGR. Какой минимум мне нужно установить, чтобы использовать привязки Python GDAL и OGR? Какой самый …

18 python  gdal  ogr 

Я пытаюсь изменить значение поля шейп-файла. Однако это не меняет значение. Что не так с моим кодом? 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 

У меня есть шейп-файл с полигонами. И у меня есть глобальный растровый файл. Я хочу наложить полигоны шейп-файла на растровую сетку и рассчитать среднее значение растра для каждого полигона. Как я могу сделать это с помощью GDAL, записав результаты в шейп-файл?

18 python  raster  shapefile  gdal  spatial-statistics 

В настоящее время я использую следующий сценарий для добавления некоторых данных атрибутов из таблицы во множество отдельных шейп-файлов: 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") } В …

18 shapefile  gdal  r  ogr  attribute-table 

Я репроектирую растры в Python, используя GDAL. Мне нужно спроецировать несколько tiff из географических координат WGS 84 в WGS 1984 Web Mercator (Вспомогательная сфера), чтобы потом использовать их в Openlayers вместе с OpenStreetMap и, возможно, картами Google. Я использую Python 2.7.5 и GDAL 1.10.1 отсюда , и преобразовываю координаты, используя …

17 python  gdal  coordinate-system  epsg  web-mercator 

Я часто вижу в коде Python GDAL, что люди закрывают наборы данных в конце своего скрипта. Почему имеет смысл закрывать набор данных в Python GDAL? Есть ли последствия, если я не сделаю это? import gdal # open dataset ds = gdal.Open('test.tif') # close dataset ds = None

17 python  gdal 

Я готовлю данные для включения в узел геосети, где карты отображаются через геосервер. Проблема в том, что ячейки, у которых нет данных, отображаются в геотеке черными, и мне бы хотелось, чтобы цвет был белым или прозрачным. Я использую ArcView и не имею доступа ни к Spatial Analyst, ни к 3D-аналитикам. …

17 arcgis-desktop  geoserver  gdal  geotiff-tiff  geonetwork 

Я пытаюсь обработать некоторые биоклиматические растровые файлы, например, которые можно скачать с http://www.worldclim.org/current (набор биоклимов). Кажется, что для них установлены значения узлов в -3.4e+38соответствии с QGIS (если посмотреть на вывод gdalinfo, это так -3.39999999999999996e+38). Кажется, что инструменты gdal не могут справиться с этим значением нодаты, и qgis, похоже, тоже не …

17 qgis  raster  gdal 

Кто-нибудь знает, как проще всего установить GDAL 2.0 на MacOSX? Обычные подозреваемые бинарный дистрибутив kingchaos и homebrew, похоже, застряли на GDAL 1.11 Мне нужна поддержка свойств объектов TopoJSON, которые, кажется, были включены в 2.0 EDIT 2017/01/07: бинарный дистрибутив kingchaos теперь предоставляет дистрибутив GDAL 2.1

16 gdal  installation  ogr  mac  homebrew 

В настоящее время я работаю над тем, чтобы импортировать климатические данные CANGRID (предоставленные как Surfer Grid ascii, файлы .grd) в ArcGIS. Сетка размером 95 строк на 125 столбцов. Метаданные обеспечивают широту / долготу происхождения (нижний левый угол), размер ячейки (50 км), а также проекцию нот в виде полярной стереографии с …

16 arcgis-desktop  python  gdal  coordinate-system