A Geospatial Data Abstraction Library (GDAL) image mosaicing, reprojection, and warping utility

learn more… | top users | synonyms

1
vote
1answer
22 views

Use gdalbuildvrt then gdal_translate changes Corner Coordinates

when I use: gdalwarp -dstalpha -multi -t_srs epsg:900913 b172.tif b282.tif b770.tif b771.tif b772.tif b832.tif w_25k.tif I can get one merged tif then go through gdal2tiles for tiling But if I use ...
2
votes
2answers
52 views

Clip raster with shapefile and set a value outside the mask

I want to clip a raster using a shape file with the command: gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif from this answer. However I want to set the pixels outside the mask to a particular ...
11
votes
1answer
335 views

PyQGIS save raster as rendered image, then use GDAL tools on it

I have tried to use the code from QGis Save Raster as Rendered Image. i is a raster layer: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), ...
1
vote
3answers
617 views

Using very large image formats with GDAL

I would like to use gdalwarp and gdal_translate with some very large raster images (313,040 x 313,040 px) and I'm wondering what the best format or method is to do this, if it's even possible. I have ...
0
votes
0answers
31 views

How to crop orthophoto using GDAL

I have orthophoto(GeoTiff) and 4 GPS coordinates that represent rectangle. How can I crop orthophoto by using GDAL? (I need to remove areas outside of rectangle) UPDATE: I tried gdal_translate ...
0
votes
1answer
42 views

GDAL Warp giving garbage data on edges when transforming to WGS 84

I want to change one tiff file to USGS DEM file. For that I am doing following steps: Transform the tiff file to WGS84 coordiantes using Gdal Wrap utility gdalwarp -t_srs EPSG:4326 -dstnodata -...
2
votes
1answer
35 views

Understanding how `gdalwarp` works under the hood

I try to understand how gdalwarp works. Here is a reference to Affine GeoTransform, but I can't understand the terminology pixel/line is it x,y coordinates of image? does pixel width and height ...
2
votes
1answer
20 views

Does GDALWARP have support for OpenCL speedup to orthorectification?

So I am attempting to orthorectify some images very large images. The process can take several hours so currently I am relying on using 8 threads to perform the speed up. However, the GPU does not ...
0
votes
1answer
13 views

How to project a large nonortho satelite image on a much smaller high res DEM?

I am currently attempting to orthorectify an image using GDAL. I have a very detailed DEM that I would like to rectify the image on top of. However, the DEM is smaller than the image itself; this is ...
7
votes
1answer
243 views

gdalwarp producing empty rasters when called from python script, but not via command line

I'm trying to aggregate some raster data (all .tif) up to a coarser resolution (from .05 degree to .25 degree) using gdalwarp in python, but the command is not working. Instead of getting an output ...
2
votes
2answers
42 views

QGIS clipper tool misaliagnment

I'm using QGIS 2.14 and GDAL 1.11.3 and I'm running into a weird error with the clipped outputs. I created a shapefile in QGIS and used it to clip my image. Here's a simplified version of the ...
2
votes
1answer
60 views

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

I need to change the pixel size of a GTiff in Gdal without changing the extent. When I use gdal_translate -tr xres yres src_dataset dst_dataset It changes the extent. I do not need a rescaling or ...
1
vote
1answer
253 views

How can I clip multiple geotiff tiles with a shapefile without greatly increasing their size?

I'm using gdalwarp to clip a number of geotiffs (tiles of a larger map) to a polygon described by a shapefile (which is about the size of the overall map, i.e. much larger than individual tiles) using:...
0
votes
0answers
24 views

gdalwarp not clipping VRT to correct extent

I'm trying to clip a projected .vrt raster to an extent expressed in GCS coordinates and it's not doing what I would expect. Here is my command: gdalwarp -te -133.08813624 56.69042536 -132.69405069 ...
3
votes
2answers
86 views

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

I have been struggling a bit to generate tiles for a high-resolution image that we have. The current image we have is a very large (+20GB) image, saved as a GeoTiff file. I would like to generate ...
2
votes
1answer
79 views

Trouble using “Clip Raster by Mask Layer” algorithm QGIS

I'm having issues using the "Clip Raster by Mask Layer" algorithm to clip a raster by a shapefile. When I run the algorithm, I get an error saying that the mask layer could not be loaded, even though ...
0
votes
0answers
34 views

using gdal to convert from NAD83 to WGS84

I'm trying to use gdal in OSGEO4W to convert from NAD83 to WGS84. I'm very new to this so if you can, be as thorough as possible with your answers. The command I'm running is gdalwarp --config ...
0
votes
0answers
40 views

Unable to generate tiles using gdal2tiles.py for a single band gray image

I am using and playing around with Cesium at the moment, and I am enjoying the functionality so far. I am trying to provide my own imagery to Cesium with the TMS imagery provider: imageryProvider: ...
0
votes
0answers
11 views

gdalwarp fails based on output filename? does it have to do with my using -t_srs?

I am so confused by what is going on with my gdal warp. From my terminal I can enter a command for it: $gdalwarp -t_srs '+proj=ortho +datum=WGS84' test_imagery/Maui/OBS_IMAG.nitf testout3.tif ...
1
vote
2answers
529 views

gdalwarp changing value in nodata pixels

I'm using gdalwarp to convert MODIS NDVI mosaics from Sinusoidal to epsg:4326. To save space I'm using VRT files until I reach the final product. Here is my workflow. Mosaic NDVI Modis tiles using ...
4
votes
2answers
118 views

Convert MODIS HDF file in sinusoidal projection into GeoTiff with

I'm looking for a way to convert MODIS LST data from HDF to Tiff format, and reproject from sinusoidal to WGS84 - all using python GDAL. An example file is here (I only care about band 0) : http://...
0
votes
1answer
47 views

Understanding GDAL ERROR 1: Too many points… failed to transform?

How can one trace the source of the following error when using gdalwarp? (I'm very new to GIS, GDAL and python!) Warning 1: Unable to compute source region for output window 0,0,5789,4748, ...
0
votes
0answers
53 views

Patch processing of multiple HDF5 files using Python and Gdal (gdalwarp)

I am trying to extend this process http://geoinformaticstutorial.blogspot.com/2014/06/reading-and-map-projecting-amsr-2-data.html wherein a set of raster data with separate geolocation arrays are ...
1
vote
1answer
27 views

gdalwarp reprojects only a portion of the raster

I´m trying to re-project a raster from epsg:4326 to epsg:3035 to work with European based datasets. I have been using this command: gdalwarp -srcnodata 9.95924e+36 -dstnodata 9.95924e+36 -s_srs epsg:...
1
vote
1answer
138 views

mask raster with shapefile using gdalwarp/gdal

I am struggling to use gdalwarp with cutline option. The following command keeps failing and I do not understand why : gdalwarp -cutline shapefile.shp file.tiff newfile.tiff Warning 1: Self-...
1
vote
0answers
53 views

How can I prevent gdal_translate from rounding coordinates when clipping?

I want to clip a raster (a VRT of several Ascii files) to a set of coordinates. Ideally I'd like to use gdal_translate to do this, but I only get an accurate result when using gdalwarp. gdal_translate ...
1
vote
0answers
35 views

Error executing algorithm Clip raster by mask layer 'NoneType' object has no attribute 'GetGeoTransform'

Processing modeler: GDAL "clip raster by mask layer" tool I have a raster layer that I wish to be clipped by a polygon layer (buffer zone). I intend to repeat to process many times in a more ...
0
votes
0answers
23 views

gdalwarp -t_srs EPSG:XXXX failed on data overlaping “Antimeridian”

I have program which process some image from satellite (crop and projection). For data overlaping "antimeridian" my process failed. This is the FootPrint of my image: Corner Coordinates: Upper Left ...
0
votes
0answers
11 views

Include partial grid cells with gdalwarp -cutline

I am using gdalwarp (v 1.10) to crop a large raster to the size required by my application using the command line utility with the following command: gdalwarp -cutline border.shp -crop_to_cutline ...
0
votes
0answers
30 views

Creating superoverlay from geotiff covered whole globe, shouldn't of?

I am new to this GIS world and found this blog post: https://developers.google.com/kml/articles/raster#step-1-download-and-install-gdal Following its instructions I used my own geotiff (ntif) image ...
5
votes
0answers
209 views

can I preserve a fading alpha layer when mosaicing images with gdal?

I am mosaicing some images with gdal and would like to improve the final result by using a fading / gradual alpha layer towards the edge of each image to remove the sharp edges in the middle of the ...
0
votes
1answer
29 views

NED (NAVD88) to WGS84 heights

I need to convert the 1/3 arc second DEM files from the NED into WGS84 ellipsoid heights. The documentation says that the NED is referenced to NAVD88. I've tried using gdalwarp as follows: gdalwarp ...
0
votes
0answers
12 views

Resampling using gdalwarp gives me no data

I have an image with the following information. Driver: GTiff/GeoTIFF Files:test_DEM.tiff Size is 91, 73 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", ...
4
votes
1answer
108 views

Georeferencing Himawari-8 in GDAL (or other)

Does anyone know the appropriate project information to georeference (and hence reproject) Himawari-8 AHI data in the IR? The images are 5500x5500 pixels. Currently I try these: ulx=-5570248.832537 ...
1
vote
0answers
46 views

Gdalwarp UTM reprojection causes image distortions

I'm trying to re-project Sentinel-2 images to a common UTM zone in order to create a mosaic. For this I'm using gdalwarp with the following syntax: C:/OSGeo4W64/OSGeo4W.bat; gdalwarp -t_srs \"+proj=...
6
votes
2answers
5k views

What's the difference between gdalwarp and gdal_merge for mosaicing

I wanted to merge around 20 dem tiffs each around 100MB, and I used gdalwarp (I know about virtual file option): gdalwarp $(list_of_tiffs) merged.tiff However this command took very long time, and ...
3
votes
0answers
37 views

How to blend edges of aerial imagery with gdalwarp

I need to fade/blend the edges of a series of aerial images with gdal. The closest thing that I've been able to accomplish is using the -cutline and -cblend switches of gdalwarp to get the edges to ...
0
votes
1answer
56 views

Raster extents incorrect for vrt when changing projection with -a_srs flag

I want to generate a .vrt in epsg:3857 from a series of raster tiles in epsg:4326, and am trying to do so in one step using gdalbuildvrt with the -a_srs flag: gdalbuildvrt -a_srs epsg:3857 ortho_3857....
0
votes
0answers
43 views

How to convert raster in UTM to lat/lon using GDAL (or other tool)

I have a .tiff file given with UTM coordinates that I'd like to convert to a .nc (NetCDF) file with latitude/longitude coordinates. I've tried with gdal_translate and gdal_warp but these did not work. ...
2
votes
1answer
33 views

Corners of GeoTiff's not lining up after doing gdalwarp EPSG:3006 -> EPSG:4326

I have some GeoTiff's that are in EPSG:3006 (SWEREF99 TM) atm, and I'm trying to convert them to EPSG:4326 (WGS84) because some other software needs it in this projection. But after doing the ...
1
vote
0answers
90 views

Slow performance clipping raster with GDAL

I have a few large LiDAR DEMs that I want to merge. Each DEM has some glitches around the borders which need to be removed in order to get a seamless final product. I manually traced around each DEM ...
1
vote
0answers
80 views

How to make NoData transparent?

I have been using python to turn a netcdf into a geotiff. Once I have the geotiff, I am using gdalwarp to transfer it into the right projection. Then use gdal_translate to get the tif into a vrt file....
1
vote
0answers
39 views

gdalwarp won't run with “cubic” resampling method

I am having an issue with running the gdalwarp command from command prompt. Basically I have a raster WGS84 file, which I would like to reproject to spatial reference which is using the Azimuthal ...
1
vote
1answer
40 views

gdalwarp output shifted 15 KM

This is similar to this question regarding shifts in raster output, but the difference I'm experiencing is fairly large-- the output raster is shifted 15 KM to the south. Here's an image: When ...
1
vote
1answer
20 views

Error using gdalwarp in EPSG:4326

I tried to change the projection of AMSR-2 in Python using: 'gdalwarp -geoloc -t_srs EPSG:4326 C:\GW1AM2_201301311114_050A_L1SGBTBR_1110110.vrt C:\test100.tif' It said there are something wrong ...
1
vote
1answer
54 views

gdalwarp causing shift in pixels

Can somebody help me with gdalwarp issue? I am having an input raster and would like to clip it with polygons. After running the below code, while overlaying output raster on input raster, i see there ...
0
votes
0answers
49 views

How to improve poor image quality of “gdalwarp” operation result?

I'm trying to improve the result of gdalwarp operation. The image width and height are 6732x4721 GCP's are: -gcp 0 0 38.005301407753 55.048725423557 -gcp 6732 0 38.184943751144 55....
0
votes
0answers
45 views

gdalwarp leaves black boundary of cropped geotiff with cutline

This is a new question based on a response on the thread white page border when converting geoPDF to geoTiff with gdal_translate I am using gdalwarp as described to do a cutline of geotiff images ...
1
vote
2answers
78 views

Black image using gdalwarp

Im a geomatics student; Currently, I try to geo-referencing a tif image with the QGIS GUI, and then I load and save the Gdal script. With the QGIS GUI, I managed to get a decent picture, then I ...
3
votes
1answer
141 views

gdalwarp to clip a .tif: Works fine standalone, won't work in Python

I am trying to use the subprocess module to call gdal warp and clip a .tif using a shapefile. When I try the command in GDAL it works fine, but when I attempt to call it into a Python script it will ...