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

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

31 coordinate-system  gdal  geotiff-tiff  proj  gdalwarp 

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

26 gdal  gdalwarp 

Я пытаюсь перепроектировать / повторно сэмплировать с привязками Python GDAL, но получаю немного другие результаты по сравнению с результатами из утилиты командной строки gdalwarp. Смотрите обновление ниже для более короткого примера Этот скрипт иллюстрирует подход Python: from osgeo import osr, gdal import numpy def reproject_point(point, srs, target_srs): ''' Reproject a …

20 python  gdal  gdalwarp  numpy 

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

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

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

14 gdal  gdalwarp  gdal-merge 

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

14 python  raster  gdal  clip  gdalwarp 

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

13 gdalwarp  qgis-3 

Я попытался использовать код из QGis Save Raster в качестве отрендеренного изображения . iэто растровый слой: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), pipelayer.height()) piperenderer = pipelayer.renderer() pipeprovider = pipelayer.dataProvider() crs = pipelayer.crs().toWkt() pipe = QgsRasterPipe() pipe.set(pipeprovider.clone()) pipe.set(piperenderer.clone()) pipedFile = os.path.join(tempfile.gettempdir(), safeLayerName + '_pipe.tif') print pipedFile file_writer …

12 qgis  gdal  pyqgis  gdalwarp  gdal-translate 

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

12 gdal  geotiff-tiff  gdalwarp  gdal-translate  gdalbuildvrt 

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

11 srtm  gdalwarp 

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

11 gdal  coordinate-system  proj  gdalwarp  netcdf 

Я хочу привязать растр с помощью pythonи GDAL. Мой текущий подход состоит в том, чтобы вызвать gdal_translateи gdalwarpиспользовать os.systemи ужасный список наземных контрольных точек. Я действительно хотел бы сделать это изначально внутри python. Это текущий процесс, который я использую: import os os.system('gdal_translate -of GTiff -gcp 1251.92 414.538 -7.9164e+06 5.21094e+06 -gcp …

10 python  raster  gdal  coordinate-system  gdalwarp 

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

10 gdal  gdalwarp  cesium  gdal2tiles 

Я пытаюсь спроектировать растр. В 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 <- "+init=epsg:4714" # using raster package to reproject pr1 <- projectRaster(r, crs = CRS(newproj), method …

9 raster  gdal  r  coordinate-system  gdalwarp 

Я использую Python и QGIS 2.0. Я пытаюсь обрезать растры в папке одним полигоном. Это первый раз, когда я использую (скажем, «PyQGIS»), я раньше привык к Arcpy. В любом случае, мой простой сценарий не работает, любое предложение будет высоко оценено! import qgis.core, qgis,utils QgsApplication.setPrefixPath("C:/OSGeo4W64/apps/qgis", True) QgsApplication.initQgis() CLIP= "C:/Users/unim/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/study_area_foscagno.shp" INPUT_FOLDER="C:/Users/unimi/Documents/Umberto/Universita/PhD/Guglielmin/Permafrost/Alta_Valtellina/Landsat_ita/LE71930282000259EDC00" OUTPUT= …

9 gdal  pyqgis  gdalwarp