Запитання з тегом «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 pair of …

20 python  gdal  gdalwarp  numpy 

Після використання gdalwarpдля проектування та вирівнювання до сітки (через -tap) декілька растрових я помітив, що вихідний растр був значно більшим, ніж оригінальний растр. Досить ретельний веб-пошук виявив цю проблему Trac : Френк Вармердам пояснив причину: "При уважному огляді різниця у розглянутому файлі полягає в тому, що gdal_translate використовує інтерфейс TIFFWriteScanline …

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

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

15 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 інструменти панелі меню, що випадає, внизу мали редаговане поле виклику, куди можна було вручну вводити команди, недоступні у gui (натискаючи кнопку олівця). У 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 0' …

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 865.827 …

10 python  raster  gdal  coordinate-system  gdalwarp 

Я трохи намагаюся генерувати плитки для зображення з високою роздільною здатністю. Поточне зображення, яке ми маємо, - це дуже велике (+ 20 Гб) зображення, збережене як файл GeoTiff. Я хотів би генерувати плитки за допомогою утиліти командного рядка gdal2tiles, а потім відкрити та переглянути її у Cesium, використовуючи постачальник зображень …

10 gdal  gdalwarp  cesium  gdal2tiles 

Я намагаюся спроектувати Растр. У R є projectRaster()функція цього (нижче повністю відтворюваного прикладу): # 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", я раніше звик до архпінгу. У будь-якому випадку, я не можу працювати зі своїм простим сценарієм, будь-яка пропозиція буде дуже вдячна! import qgis.core, qgis,utils QgsApplication.setPrefixPath("C:/OSGeo4W64/apps/qgis", True) QgsApplication.initQgis() …

9 gdal  pyqgis  gdalwarp