Questions tagged [gdalwarp]

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

Filter by
Sorted by
Tagged with
0
votes
0answers
30 views

Constructing boundary polygon around raster image when mask layer is ill-defined using GDAL with Python?

I have a raster image and I want to construct a new vector layer that has the boundary of the raster, using python. Below is a sample of the image. The problem is that the nodata value is set to "...
0
votes
0answers
15 views

How to reproject GeoTIFF from WGS84 to UTM with gdalwarp?(win 10) [on hold]

I used code like below: gdalwarp -t_srs "+ proj = utm + zone = 1 + datum = WGS84" geo.tif output_utm.tif and I caught a problem : Translating source or target SRS failed: + proj = utm + zone = 1 + ...
1
vote
0answers
29 views

Gdalwarp orthorectification WorldView-3 not using RPC projection properly

I have a high-resolution WorldView-3 image that I am trying to orthorectify using gdalwarp -rpc -to "RPC_DEM=dem.tif" input.tif output.tif However the output is not aligned with Bing or Google maps ...
2
votes
0answers
26 views

GDAL Warp Memory Datasource as Cutline

I am attempting to use a memory vector to crop a raster using GDAL warp but I get the following error: ERROR 1: Cannot open <osgeo.ogr.DataSource; proxy of <Swig Object of type '...
3
votes
0answers
45 views

“Clip Raster by Mask Layer” produces zero values (and not Nodata) outside mask boundaries

This is my first question here. I have a big raster of the world's nighttime lights and want to clip it down to the Adm0 of a country: In order to do so, in QGIS 2.18 what I do is using the Clipper ...
0
votes
1answer
24 views

Error using gdalwarp

I have this NetCDF file of the Greenland ice sheet velocities but the coordinate system is NSIDC Sea Ice Polar Stereographic North (y: -655750 to -3355750 and x: -639750 to 855750 meters). I wanted to ...
0
votes
0answers
14 views

Gdal_calc changing the values incorrectly

I am trying to re-project a file to a uniform grid EPSG:4326 and for that I am using the command: gdalwarp -t_srs EPSG:4326 -srcnodata -32767 HDF5:"LC08_L1TP_029037_20190611_20190611_01_RT_L2GEN_2.nc"...
1
vote
0answers
30 views

Error when reprojecting raster (warp) QGIS [duplicate]

When trying to reproject a raster from degrees to meters in warp (QGIS) I receive the following error: Processing algorithm… Algorithm 'Warp (reproject)' starting… Input parameters: { 'DATA_TYPE' : ...
1
vote
1answer
37 views

Values change after re-projecting the NetCDF files to TIFF file

I am trying to re-project a NetCDF file to a uniform grid EPSG: 4326 and convert it to a TIFF file. However, after the GDAL re-projection the values in the file change, i.e. the values are not ...
0
votes
1answer
63 views

Transforming a image to a shape defined by image corners in earth coordinates - Python

I have a source satellite image. It is a rectangular shape x by y. I also have the stored lon/lat coordinates of the corners and the center point in the xml meta file. Is there a way with the gdal or ...
1
vote
1answer
45 views

How to convert a TIFF file to HGT file?

In my project I needed to convert TIFF file to SRTM data and came up with following warning/error. gdal_translate -of SRTMHGT Terrain.tif Terrain.hgt Output: Warning 1: The source projection ...
2
votes
1answer
54 views

Averaging overlapping rasters with gdal

I am working on creating a class which will merge several georeferenced rasters into one using different strategies, essentially taking average, max, min where the images are overlapping. So far I've ...
2
votes
0answers
36 views

Is there a difference in using Haversine distance vs re-projecting coordinate frame?

I have a GeoTIFF raster file in the EPSG 4326 (WGS 84) geographic coordinate system. I want to crop a 20 km x 20 km square around a particular (lat, long) coordinate, and save that as a new raster. ...
1
vote
1answer
20 views

Trying to make calculation with gdal_calc.py?

I have a tif image that I'm trying to run gdal_calc.py. In the source file the values are between 1,118-9,704 and I want to divide all the values by 10,000 to get results between 0-1. This is the ...
0
votes
0answers
11 views

Reproject Sinusoidal Big Raster to EPSG 4396

Today reproject MODIS Sinudoidal to EPSG 4396, in my older GDAL and Qgis 2.14 easy. But today no understand why not and sending this bugs!! david@dummy13:/mosaiced_tiles$ gdalwarp -s_srs "+proj=sinu +...
1
vote
1answer
50 views

Transforming NetCDF file in longlat projection into TIF in EASE 2.0 using gdal_translate and gdalwarp

I am very new in this. My goal is to convert 1 subband from .nc files (which I have several thousand if it matters) stored in longlat projection into .tif files in EASE 2.0 projection. I do that in 2 ...
2
votes
0answers
70 views

Saving osgeo.gdal.Dataset as GeoTIFF or JPEG with world files using GDAL from Python?

I am trying to take a plain JPEG with no spatial references, apply ground control points (GCPs) to it, and output a georeferenced GeoTIFF or JPEG with world files. I found some code that seemed like ...
0
votes
0answers
12 views

gdal_warp cannot specify target extents

This is related to this previous post. As was trying to implement the said answer and inserting it to my code, but the code doesn't seem to work. #importing all the necessary modules import os, ...
0
votes
1answer
52 views

How do I crop a TIFF down to a specific region?

I have two TIFFs that have the same pixel size. One of them is geolocated inside the other. I would like to crop the larger one down to the exact same extent as the smaller one. Is this possible with ...
1
vote
1answer
110 views

Using GDALwarp for reprojecting netCDF file?

I am trying to re-project a netCDF file using the gdalwarp command. However, the gdalwarp is giving me a black output. Below is my code: import gdal import osr import netCDF4 from osgeo.gdalconst ...
0
votes
0answers
30 views

gdalwarp : I lose some pixels when I crop

My problem: When I crop with gdal gdalwarp -te minX minY maxX maxY input_raster.tif output_raster.tif I expect an output shape of (maxX - minX, maxY - minY) = (10980,10980) BUT : I get an output ...
1
vote
0answers
63 views

Why gdal does not find gcps coordinates on SENTINEL images?

I use gdalwarp with python to overlay polygons on ENVISAT radar images. Everything goes well for this type of images but when I try with SENTINEL images, gdal does not find the georeferencing. It ...
1
vote
1answer
33 views

gdalbuiltvrt output is different than the source images, when using dstalpha?

I have warped some images from various state plane coordinate to EPSG 3857, but the output of gdalbuildvrt is messed up (see image). I've tried dozens of different combinations of gdalwarped and ...
0
votes
1answer
17 views

CGIAR DEM + gdalwarp — projection to Cartesian is not happening

I have downloaded a 30-by-30 tile from http://srtm.csi.cgiar.org/srtmdata/ Then gdalinfo on the obtained single *.tif-file would yield Driver: GTiff/GeoTIFF Files: cut_n30e000.tif ...
1
vote
1answer
32 views

UTM Projection of DEM-data with gdalwarp

I have downloaded a certain tile from https://www.eorc.jaxa.jp/ALOS/en/aw3d30/data/index.htm, and after the requisite logistics (unzipping, putting all *.AVE_DSM.tif files into a single folder), I ...
0
votes
1answer
54 views

Gdal.warp broken in version 2.3.3?

I have two versions of Python gdal running in different environments. One is slightly older (2.1.1) and gdal.Warp seems to work just fine there, the newer version doesn't. Has anyone else had this ...
2
votes
1answer
49 views

GDALwarp HDF5 Image With Unknown Projection

I have a NetCDF file that contains a series of HDF5 bands. I am looking to project one of these bands contained within this file. I've attempted to use the following command to achieve this: gdalwarp ...
0
votes
0answers
21 views

GDALWARP ERROR 6: Only single band dataset supported for non Byte datatype

I'm trying to cut an image using GDALWARP and I have the error: ERROR 6: Only single band dataset supported for non Byte datatype I think my image is of a single band and data type Float32. The result ...
1
vote
1answer
51 views

How do I transform an EPSG:3857 world map to an EPSG:4326 ECE world map?

This follows up a question posted over seven years ago. The results look very similar to the original inquiry without an answer or update I could test successfully. I made a custom stylized EPSG:3857 ...
0
votes
0answers
26 views

Specifying dimension for variable in NetCDF input file using GDAL?

I have a data which is in .nc format. I want to extract a variable which is in XxYxZ dimensions for example 4032x496x3 inside the netcdf file. But the problem is when using gdal_translate it is ...
0
votes
0answers
14 views

How can I convert MISR(multi-angle image spectro-radiometer) level 2 AOD (MIL2ASAE) swath dataset to grid data

I am able to convert as batch the 'Aerosol_Optical_Depth' data set inside MISR(multi-angle image spectro-radiometer) level 2 (MIL2ASAE) swath data to grid data,which is in .nc format, and in 2 ...
0
votes
0answers
95 views

GDALWarp and Target Aligned Pixels Not Effective

I have two sets of TIFFs of different resolutions (100m and 1000m), different alignment, and varying data types with corresponding no data values. I need to downsample the 100m TIFFs to match the ...
0
votes
0answers
34 views

Matching Mollweide extent?

I'm trying to match two rasters to same size, so i can do cell to cell analysis on them. The first cell That I'm trying to match is in Mollweide projection, with 250 m resolution. Here is the gdalinfo ...
1
vote
0answers
31 views

Split raster to many irregular polygons [closed]

I'm splitting a national soils raster file into regional (county) rasters, and organizing the soils classifications into a database keyed by region, eventually plot- or farm-level detail. At the ...
2
votes
1answer
36 views

gdalwarp - correct for a slight spatial misregistration?

Some Sentinel-2B band images have slight spatial misalignments as I explained on the step forum. https://forum.step.esa.int/t/sentinel-2-tci-10m-level-2-images-slight-spatial-misregistrations/14235 ...
2
votes
1answer
58 views

sampling error using gdalwarp on a subsetted VRT?

I am currently trying to prepare an SRTM subset including height conversion from EGM96 to WGS84 using gdalwarp. I encountered some strange shifts in part of the image. Here is a reproducible example: ...
0
votes
1answer
271 views

GDAL and reprojecting a NetCDF file

I have a .nc file that you can download here. This is satellite imagery. I want to re-project this and save it to a new .nc file, without saving to raster. Just NetCDF to NetCDF only. Is this ...
-1
votes
1answer
71 views

Python gdal.Warp doesn't create a file

I would like my code to write a tiff file from GRIB while changing the resolution. So I take file1 (GRIB fmt) and translate it to tiff. The file appears in the folder. Then, I make a warp, but no ...
1
vote
0answers
100 views

Gdal raster translate and wrap

I am trying wrap my raster with GCP control points using gdal_translate and gdalwrap tools. I've read at gdal mailing list about wrap with affine transformation, so I thought it is easy to done using ...
1
vote
0answers
52 views

Project a subset of GeoTIFF image from UTM 32651 to EPSG:4326 by bounding box in EPSG:4326?

I have a GeoTIFF file which is in UTM 32651 and I need to develop a program, using GDAL library version 1.x (e.g., not using gdalwarp/gdal_translate binary tools) which takes a subset bounding box in ...
1
vote
0answers
26 views

gdalwarp looses pixel information

I have an ungeoreferenced .tif of which I know the spatial information. However, when I run the following command it looses all unique values. All pixels have the value 32, which is the upperleft ...
0
votes
0answers
65 views

[Resolved) Resampling ALOS images with Python

I'm trying to resample an ALOS image of 12.5 m resolution to 1 m resolution with cubic interpolation using GDAL warp. Here is the current code used in pyzo import os from os import path import ...
0
votes
0answers
46 views

No transformation when applying gdal.warp

I am trying to convert a GeoTiff resolution and projection by applying gdal.warp. Unfortunately, there is no change on the output. > transformation = gdal.Warp('/vsimem/transfotest.tif'.encode('...
0
votes
0answers
91 views

gdal2tiles (or some steps before) produces erroneous tiles from a mosaic TIFF

I have for some days now tried to create a image pyramid from a set of around 1400 TIFF files with no luck. I can't get my head around what is the cause of the problem. I will try to explain as good ...
1
vote
1answer
57 views

Resampling two raster files

I have two rasters, one for road network and another one for the DEM of the country. I have to overlay the road network raster on top of DEM raster so as to get the elevation of roads. But my rasters ...
2
votes
0answers
63 views

Is there a GDAL function similar to OpenCV's WarpPerspective?

I't trying to mosaic a GeoTiff raster against another. Initially, they were both georeferenced in WGS-84 CRS. Using OpenCV, I calculated a homography transformation (a 3x3 projection matrix), that ...
2
votes
0answers
59 views

gdal_translate does not crop the GeoTIFF file [closed]

When I use gdal_translate to crop the GeoTiff file, the output file is saving the same as the input file. Nothing changed. The size needs to decrease at the boundaries I specify. What can I do? from ...
2
votes
1answer
69 views

How to set correct geo information on pleiades images with RPC model?

Now I have pleiades raster (img.tif) and a RPC model (rpc.XML). In the rpc.XML, there are the following lines which can be used to reset geo information. <RFM_Validity> <...
1
vote
1answer
80 views

Using GDAL to strip nodata in border regions only?

I am trying to create a script that given a tiff file with islands of data (which can be disjoint) and black pixels of nodata, can generate polygons (or new tiffs) that only wrap around the islands of ...
0
votes
0answers
35 views

Gdalwarp producing odd results around 0 degree longitude

I am warping ECMWF grib data (full globe, 0.25 degree resolution) to match the SEVIRI geostationary satellite projection. I use this command: gdalwarp -overwrite -ot Float32 -t_srs "+proj=geos +a=...