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

learn more… | top users | synonyms

0
votes
0answers
26 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 ...
2
votes
1answer
34 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
12 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
10 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 ...
2
votes
2answers
31 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
58 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 ...
0
votes
0answers
21 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
78 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 ...
0
votes
0answers
33 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
39 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 ...
2
votes
1answer
78 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 ...
2
votes
1answer
48 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 ...
0
votes
1answer
46 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
49 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
26 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
0answers
49 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
34 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 ...
1
vote
1answer
118 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-...
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
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 ...
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
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 -...
0
votes
1answer
25 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 ...
4
votes
2answers
110 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
0answers
11 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", ...
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
0answers
41 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
32 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
86 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 ...
4
votes
1answer
107 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
79 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
44 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=...
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
46 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
47 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
44 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 ...
0
votes
0answers
32 views

intersect geotiff with bbox

In a project using gdal I need to crop a geotiff file with bbox in pixel or coordinates epsg:4326. This step is ok. The next step is to manage cases where the image size is "lower" than the bbox size ...
2
votes
1answer
43 views

Using gdalwarp to remove one row of pixels or add one row of pixels

I am trying to use GDAL_calc but my two rasters are off by one pixel in the y dimension. One is 256005 the other is 256004 so GDAL_calc fails with a multiplication. Is it possible in gdalwarp to ...
1
vote
1answer
115 views

GRIB to NetCDF with GDAL, specifying dimensions

I have a GRIB dataset representing temperature over North America, in Lambert Conformal Conic (with specified central meridian, latitude of origin, etc.). It has what should be understood as a ...
4
votes
3answers
39 views

gdalwarp to geographic coordinate system

I'm getting the folllowing Error when trying to re-project a raster from NAD_1983_Albers to GCS_North_American_1983. I used the proj4 defenition from going to Properties in QGIS of the layer whose ...
1
vote
1answer
64 views

GDAL: gdalwarp wrong transformation to “South_Pole_Lambert_Azimuthal_Equal_Area” EPSG:102020

I am trying to reproject a raster that is defined in EPSG:4326 to EPSG:102020. Command used gdalwarp -ot Float32 -s_srs EPSG:4326 -t_srs EPSG:102020 -r near -of GTiff -co COMPRESS=DEFLATE -co ...
0
votes
1answer
61 views

Reprojecting and retiling a large set of tiles

I am trying to understand the process to re-project and re-tile a set of raster (GeoTiff) but can't wrap my head around how to use the different GDAL tools to achieve my goal. More concretely, I ...
7
votes
1answer
238 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 ...
3
votes
1answer
107 views

Raster interpolation with gdalwarp, and replicating with python API

My problem is similar to this: Replicating result of gdalwarp using gdal Python bindings. I want to resample a raster using bilinear interpolation, which I've done successfully with gdalwarp by adding ...
1
vote
0answers
40 views

Retaining resolution of MODIS data while reprojecting using gdalwarp

Can anyone help me to use gdal(gdalwarp specifically) for reprojection? I am using it to reproject MODIS data : I originally have an hdf file - MCD43A4.A2005001.h08v05.005.2007351192333.hdf which has ...
11
votes
1answer
324 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(), ...
3
votes
1answer
84 views

Crop raster from vector problem

So I have a shapefile of a vector for Philadelphia's city lines and a raster image of Philadelphia county that I got from the lovely people from USDA:NRCS. I converted the raster from a sid file to ...