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

12
gdalmerge создает слишком большой файл

Я пытаюсь объединить 60 файлов геотифов в один огромный файл, чтобы иметь возможность создавать плитки из полученного огромного файла геотифов с помощью gdal2tiles. Процесс слияния работал, но по какой-то причине он создает файл размером 78 ГБ вместо ожидаемых 19 ГБ. Не совсем уверен, что там пошло...

12
написать OGR с пространственным полигоном, упрощенным с помощью gSimplify

Я использую gSimplify(пакет rgeos) для упрощения геометрии шейп-файла. Функция работает хорошо, но теперь я не могу записать вывод в новый шейп-файл. Я попробовал несколько способов: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') я получил объект должен быть...

12
установка GEOS, PROJ, GDAL / OGR в python virtualenv на Mac OS X

Я пытаюсь настроить изолированный python virtualenv для работы над ГИС-проектами на моем Mac OS X. Похоже, мне потребуется установить GEOS, PROJ, GDAL / OGR с сайта kyngchaos здесь http://www.kyngchaos.com/ софт / фреймворки, а не просто_установить его в мой virtualenv. Это точно? Загрузки,...

12
Могу ли я сохранить исчезающий альфа-слой при наложении мозаики на изображения с помощью gdal?

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

12
Должен ли я всегда использовать BIGTIFF = YES по умолчанию?

Я иногда пишу короткие сценарии / рабочие инструкции для коллег, которые не очень опытны в GDAL / программировании в целом. Чтобы избежать каких-либо проблем со сценариями, при объединении / преобразовании геотегов с GDAL я стараюсь просто использовать опцию BIGTIFF = YES независимо от того, какие...

12
TRI, TPI или Roughness

Я работаю над прогнозом выбора места, где я не смогу выполнять какие-либо наблюдения на месте. Ключевым компонентом процесса выбора является прочность неровной поверхности. Я планирую использовать GDALdem для создания производных карт прочности. Я читал как Уилсон, и соавт. Бумага 2007 года и...

12
Как получить доступ к метаданным Shapefile, используя OGR?

Я пытаюсь использовать ogrinfo, чтобы получить некоторые подробности о шейп-файле, который я скачал. В настоящее время единственный способ, которым я знаю, как это сделать, это загрузить его в QGIS и вручную щелкнуть по нему, чтобы найти любую информацию о нем, например открыть таблицу атрибутов. Я...

12
GDAL полигонизировать в Python, создавая пустой полигон?

У меня возникли проблемы с использованием функции Polygonize в Python. Пример поваренной книги для этого можно найти здесь . Соответствующая часть моего кода: sourceRaster = gdal.Open('myraster.tif') band = sourceRaster.GetRasterBand(1) bandArray = band.ReadAsArray() outShapefile = "polygonized"...

12
Доступ к растровым файлам базы геоданных (.gdb) с помощью GDAL?

Я работаю над инструментом, который преобразует растровые слои в массивы для обработки с помощью NumPy, и в идеале я хотел бы иметь возможность работать с растрами, которые упакованы в .gdb, без их экспорта (особенно если для этого требуется использовать ArcGIS или ArcPy). ). Это возможно с...

12
Преобразование растровых данных в полигональные

Я пытаюсь преобразовать растровые данные в полигональные. Ранее я использовал shapefile.py для создания точечных данных по отдельным широтам / долготам на растре. Но это создаст огромные данные, и их запрос станет большой проблемой. Таким образом, было бы лучше, если бы я мог преобразовать...

12
Есть ли способ правильно заставить GDAL вызывать исключения в Python?

Я работаю с GDAL в Python и хотел, чтобы GDAL вызывал исключения, чтобы я мог сделать с ними что-то полезное. Например, если файл нельзя открыть с помощью gdal.Open (), я бы хотел выдать исключение вместо gdal, возвращающего None. Я просто хотел бы получить способ фиксировать ошибки, а не...

12
GDAL RasterizeLayer не записывает все полигоны в растр?

Я пытаюсь записать шейп-файл в растр, используя RasterizeLayer из GDAL. Я предварительно создаю растр области интереса из другого шейп-файла, учитывая определенный размер пикселя. Этот AOI затем служит базой для всех последующих растеризаций (одинаковое количество столбцов и строк, одинаковая...

12
GDAL и Python: как получить координаты для всех ячеек, имеющих определенное значение?

У меня есть двоичная сетка Arc / Info - в частности, растр накопления потока ArcGIS - и я хотел бы идентифицировать все ячейки, имеющие определенное значение (или диапазон значений). В конечном счете, мне нужен шейп-файл точек, представляющих эти ячейки. Я могу использовать QGIS, чтобы открыть...

12
Python-эквивалент gdalbuildvrt

Есть ли способ выполнить ту же задачу, что и утилита gdalbuildvrt, используя привязки GDAL Python? До сих пор я не нашел никакого другого способа сделать это, кроме как создать vrt одного набора данных и вручную отредактировать xml. Я хотел бы создать VRT из нескольких растров (по сути, выполняя...

12
Gdal Dataset.ReadAsArray () аварийно завершает работу Python

Я использую Python 2.6.5 (32bit) с Numpy 1.3 и Gdal 1.9.1, установленными в Windows 7 64bit. Я пытаюсь прочитать набор растровых данных Imagine (.img) объемом 800 МБ в массив Numpy, чтобы выполнить некоторую растровую алгебру, но как только я запускаю следующий код, происходит сбой Python.exe. from...

12
PyQGIS сохраняет растр как визуализированное изображение, затем использует на нем инструменты GDAL

Я попытался использовать код из QGis Save Raster в качестве отрендеренного изображения . iэто растровый слой: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), pipelayer.height()) piperenderer = pipelayer.renderer() pipeprovider = pipelayer.dataProvider() crs...

12
Преобразование стека временных рядов растра GTiff в один NetCDF

Переезд из списка рассылки gdal-dev: В понедельник, 2 сентября 2013 года, в 19:09 Дэвид Шин написал: Привет список, я пытаюсь упаковать временные ряды растров GTiff с идентичной проекцией / экстентом / разрешением как один файл NetCDF для распространения. Я провел последний час, консультируясь с...

12
ArcGIS Raster Mosaic Dataset и GDAL VRT, когда использовать?

ArcGIS Desktop имеет набор растровых данных мозаики , QGIS и друзья имеют виртуальные растровые таблицы . Просмотр документации обнаруживает большое сходство. Они оба создают или включают альтернативное представление данных, которое находится где-то еще, без необходимости фактически создавать новое...

12
Невозможно суммировать растры из-за узловых пикселей

У меня есть набор растров (.vrt) с ежедневными данными влажности почвы. Я хочу суммировать пиксели всех растров, чтобы иметь меру по месяцам. Однако весь мир не покрывается каждый день, что приводит к значениям узлов в местах, где существуют меры для других дней. Что я хочу сделать, это суммировать...