Skip to main content
15 votes

What's the difference between gdalwarp and gdal_merge for mosaicing

I just happened to come across this question and a potential answer when looking for something else. gdal_merge.py uses nearest neighbor resampling. If you want control over the resampling used, ...
dave_h's user avatar
  • 176
14 votes
Accepted

Why does reprojection with gdalwarp change pixel size?

gdalwarp is doing the right thing: preserving total resolution of your image by changing the pixel-size. WGS 84 / Pseudo-Mercator projection is heavily distorted when moving away from the equator. ...
pLumo's user avatar
  • 6,519
12 votes

Reprojecting a raster to match another raster using GDAL

http://www.gdal.org/gdalwarp.html says -t_srs srs_def: target spatial reference set. The coordinate systems that can be passed are anything supported by the OGRSpatialReference....
bugmenot123's user avatar
  • 10.9k
10 votes
Accepted

How to create transparent tif images using GDAL?

Use gdaldem to add color to your tif file. gdaldem color-relief input.tif style.txt -alpha output.tif style.txt looks like this, each line having up to five values. First one is the raster's pixel ...
S. Thiyaku's user avatar
10 votes
Accepted

GDAL Warp to in-memory raster

The best solution for you would be the VSIMEM filesystem which lets you save outputs of gdal utilities into a filesystem stored in memory. gdal.Warp('/vsimem/reprojected.tif', ds, srcSRS=in_proj, ...
onakua's user avatar
  • 690
10 votes
Accepted

Averaging overlapping rasters with GDAL

The following approach worked pretty well. First I build virtual raster. gdalbuildvrt raster.vrt -srcnodata 0 -input_file_list paths.txt paths.txt is file with following content: a.tif b.tif Then ...
Владислав Мокроусов's user avatar
9 votes
Accepted

How do I change the pixel size in GDAL without changing the extent?

Here's a quick script I put together to correct the cell size of a folder of GeoTiff rasters: import os, sys from osgeo import gdal from osgeo import gdalconst WorkingFolder = sys.argv[1] # first ...
Michael Stimson's user avatar
9 votes
Accepted

How to replicate gdalwarp -cutline in Python to cut out tif image based on the extents.kml file

Using the GDAL Python bindings (GDAL >= 2.1), it should be: from osgeo import gdal input_raster = "path/to/rgb.tif" output_raster = "path/to/rgb_output_cut.tif" input_kml = "path/to/extent.kml" ds =...
Antonio Falciano's user avatar
8 votes
Accepted

gdalwarp can't find libgrass_dgl or libgrass_vector files

I had the same issue. Might be a bug in gdal configuration or something. To fix it, do this: sudo ln -s /usr/lib/grass72/lib/libgrass_* /usr/lib/
Aleksey Vitebskiy's user avatar
8 votes
Accepted

Artifacts in gdalwarp output (from polar stereographic to azimuthal equidistant)

I would suggest trying to add the -et (error threshold) option with a lower thresholds than the default. If you try a lower -et threshold the horizontal artifact should disappear Or try changing ...
whyzar's user avatar
  • 12.1k
8 votes

How to fix positive North-South Pixel size in GeoTIFF with GDAL

I had the same problem. I used gdalwarp to correct the NS resolution of xyz files and converted them to tif (in one step). Here is my solution in bash: for file in *.xyz; # adjust NS resolution ...
sehan2's user avatar
  • 457
8 votes

Crop raster in memory with python GDAL bindings

With GDAL 2.2.0+ use the VSIMEM filesystem (http://www.gdal.org/gdal_virtual_file_systems.html). It allows you to treat blocks of memory as a file within the virtual file system. This version uses ...
onakua's user avatar
  • 690
8 votes
Accepted

gdalwarp cutline without using a shapefile

It should be possible by using the -csql parameter http://www.gdal.org/gdalwarp.html and SQLite dialect http://www.gdal.org/ogr_sql_sqlite.html. Gdalwarp does not have an option for selecting the SQL ...
user30184's user avatar
  • 67.3k
8 votes
Accepted

Gdal warp cutline overhang

The solutiuon was to add the option: -wo CUTLINE_ALL_TOUCHED=TRUE
Francis's user avatar
  • 957
7 votes

What's the difference between gdalwarp and gdal_merge for mosaicing

gdal_merge.py loads all files into memory before processing them. therefore it is not able to process large files if your memory is small. see here
joaoal's user avatar
  • 891
7 votes

Convert MODIS HDF file in sinusoidal projection into GeoTIFF using Python GDAL

You can either leave the -s_srs empty or use +proj=sinu +R=6371007.181 +nadgrids=@null +wktext gdalwarp -of GTIFF -s_srs '+proj=sinu +R=6371007.181 +nadgrids=@null +wktext' -r cubic -t_srs '+proj=...
Loïc Dutrieux's user avatar
7 votes
Accepted

How to use gdal2tiles on a custom tiff image recieved from providers to generate tiles

I think all you need to do is reproject it using: gdalwarp -co TILED=YES -co COMPRESS=DEFLATE -t_srs EPSG:3857 image.tif newImage.tif and then tile it: gdal2tiles.py newImage.tif If your file is ...
Ian Turton's user avatar
  • 82.2k
7 votes
Accepted

How to make Processing algorithms use default parameter values in QGIS 2.14?

It turns out you can actually use default parameter values calling Processing algorithms from PyQGIS using a different syntax. I didn't find it in the docs, but that's what we have GIS.SE for :D. ...
Germán Carrillo's user avatar
7 votes
Accepted

Gdalwarp dstnodata problems

Since GDAL 2.1 it's possible to use GDAL/OGR utilities as library functions. So this is an alternative way to accomplish this task: from osgeo import gdal input_raster = r"c:\blah\someinputtiff.tif" ...
Antonio Falciano's user avatar
7 votes
Accepted

File size inflation using GDAL

Your existing data is compressed. You're getting an increase in file-size because the default is uncompressed and you aren't specifying a compression option (-co compress=LZW). However, it's a ...
user2856's user avatar
  • 67.1k
7 votes

What can I do with gdalwarp that I can't do with gdal_translate, and vice-versa?

In some ways gdalwarp is a generalization and improvement on gdal_translate. gdalwarp can reproject rasters, this covers much of what gdal_translate can do (resize, crop, resample, convert format etc)....
mdsumner's user avatar
  • 8,206
7 votes
Accepted

How to get Amersfoort 28992 in QGIS?

You hit the already classical problem: do not change the layers's CRS, instead use Menu Raster > Projections > Warp (reproject) (in case of raster layers). Alternatively, you could have the same ...
Babel's user avatar
  • 73.1k
6 votes
Accepted

Convert MODIS HDF file in sinusoidal projection into GeoTIFF using Python GDAL

You might get unexpected results because your dataset crosses the 180° meridian. As a consequence, the tile is squeezed around the globe when reprojected to WGS84. To avoid that, you have to cut the ...
AndreJ's user avatar
  • 76.7k
6 votes

Gdal: clipping a raster with another raster

The solution in Python directly, without making shape: import gdal from gdalconst import GA_ReadOnly data = gdal.Open('img_mask.tif', GA_ReadOnly) geoTransform = data.GetGeoTransform() minx = ...
XavierCLL's user avatar
  • 380
6 votes
Accepted

Understanding gdalwarp reprojection

Just display your raster in a GIS software in UTM 32N: You see that the raster is rotated in the new CRS, and the former Lower Left point is not the most southward point anymore. Since gdalwarp ...
AndreJ's user avatar
  • 76.7k
6 votes
Accepted

how can 2 geotiff rasters be aligned and forced to have same resolution in python using gdal?

Select which image you want to use as reference and check the pixel size of that with gdalinfo. Then warp the second image to have the same pixel size with gdalwarp http://www.gdal.org/gdalwarp.html ...
user30184's user avatar
  • 67.3k
6 votes

How to edit the GDAL/OGR console call in QGIS 3.0?

This can be done using Python console. First copy (or write) a GDAL command and edit in a text viewer. For instance: gdaldem hillshade C:/MyInputFile.tif C:/MyOutputFile.tif -of GTiff -b 1 -z 1.0 -s ...
Zoran's user avatar
  • 302
6 votes

Use gdal to set multiband raster values to no data

The syntax is documented on the manual page https://www.gdal.org/gdalwarp.html as If more than one value is supplied all values should be quoted to keep them together as a single operating system ...
user30184's user avatar
  • 67.3k
6 votes

Clipping GeoTIFF with shapefile?

In QGIS, you can use the iterator of the clip tool to iterate over the polygon shapefile, and it will clip the GeoTiff raster based on each polygon inside the polygon shapefile. For example, the ...
ahmadhanb's user avatar
  • 41.1k

Only top scored, non community-wiki answers of a minimum length are eligible