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

У меня есть растровое изображение геотифа, которое имеет систему координат с долготой от 0 до 360. Горизонтальный центр изображения составляет 180 долготы. Смотрите изображение ниже: Я хочу преобразовать его в EPSG: 4326 SRS с дальностью -180 180. И я хочу, чтобы центр изображения был на...

Я пытаюсь обрезать свой растровый файл в соответствии с геометрией шейп-файла. Я использую следующий код gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif Но это приводит к черному цвету за пределами геометрии шейп файла. Я привожу пример здесь. В первом изображении, которое я хочу обрезать. Второе...

Я пытаюсь перепроектировать / повторно сэмплировать с привязками Python GDAL, но получаю немного другие результаты по сравнению с результатами из утилиты командной строки gdalwarp. Смотрите обновление ниже для более короткого примера Этот скрипт иллюстрирует подход Python: from osgeo import osr,...

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

Я хотел объединить около 20 файлов в формате 100 МБ каждый и использовал gdalwarp (я знаю о параметре виртуального файла): gdalwarp $(list_of_tiffs) merged.tiff Однако эта команда заняла очень много времени, и примерно через 15 минут она обрабатывала 5-е изображение из набора. Я прекращаю это....

Я пишу простую утилиту для обрезки пакетов многоканальных растровых файлов геотифов в одну и ту же (меньшую) область. Используя gdalwarp, я легко могу обрезать файл, используя шейп-файл с одним полигоном: gdalwarp -cutline clipper.shp -crop_to_cutline input.tif output.tif Однако фактическая...

В QGIS 2.XX инструменты в выпадающем меню были внизу с редактируемым окном вызова, в котором команды, недоступные в графическом интерфейсе, можно было вводить вручную (нажатием кнопки карандаша). В 3.0 инструменты панели меню совпадают с инструментами GDAL, которые есть в наборе инструментов...

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

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

Я использую gdalwarp для манипулирования плитками SRTM рядом с датой (т.е. 180 °, иначе антимеридиан). Плитки SRTM имеют очень небольшое (1/2 пикселя) перекрытие с меридианом. Вы можете увидеть это с помощью gdalinfo: gdalinfo S16W180.hgt Driver: SRTMHGT/SRTMHGT File Format Files: S16W180.hgt Size...

Я пытаюсь переназначить геостационарный конформный Ламберт, используя gdalwarp. Мои входные данные находятся в netcdf и в географических координатах (градусах), и я хотел бы вывести переназначенные данные в netcdf. Я создал соответствующий файл vrt для входных данных netcdf. Gdalwarp выведет файл...

Я хочу привязать растр с помощью pythonи GDAL. Мой текущий подход состоит в том, чтобы вызвать gdal_translateи gdalwarpиспользовать os.systemи ужасный список наземных контрольных точек. Я действительно хотел бы сделать это изначально внутри python. Это текущий процесс, который я использую: import...

Я немного изо всех сил пытался сгенерировать плитки для изображения с высоким разрешением, которое у нас есть. В настоящее время у нас есть очень большое (+ 20 ГБ) изображение, сохраненное в виде файла GeoTiff. Я хотел бы сгенерировать плитки с помощью утилиты командной строки gdal2tiles, а затем...

Я использую Python и QGIS 2.0. Я пытаюсь обрезать растры в папке одним полигоном. Это первый раз, когда я использую (скажем, «PyQGIS»), я раньше привык к Arcpy. В любом случае, мой простой сценарий не работает, любое предложение будет высоко оценено! import qgis.core, qgis,utils...

Я пытаюсь спроектировать растр. В R есть projectRaster()функция to this (ниже полностью воспроизводимого примера): # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project to newproj <-...