Questions tagged [gdalwarp]
A Geospatial Data Abstraction Library (GDAL) image mosaicing, reprojection, and warping utility.
741
questions
0
votes
0
answers
54
views
Cannot write GeoTIFF larger than 4GB with gdalwarp despite using BIGTIFF createopt
I am running a Bash script to create a series of NetCDF grid files from MB-System that contain gridded swath multibeam bathymetry data. Each grid is unprojected Geographic Lat-Long, about 15-20MB in ...
0
votes
0
answers
22
views
How to use gdal exceptions in .NET?
I'm converting a .vrt to a .tif with Gdal.wrapper_GDALWarpDestName(). I know the .vrt has been calculated to a size that causes an integer overflow error during warping because the command line ...
0
votes
0
answers
30
views
Warp result of GDAL translate jpg to geotiff using osgeo.gdal
I am attempting to warp a geotiff created from gdal.Translate where the original file was a jpg. I first obtain the ground control points from a co-located text file. The results of gdal.Translate ...
1
vote
0
answers
54
views
Warnings on converting NetCDF using gdalwarp
I have a curvilinear NetCDF file, which currently has a single time step (for a mean across time). The variable I am looking at is DIC at a single depth.
I am trying to convert this NetCDF file to ...
0
votes
0
answers
39
views
GDAL warp removing half of the map
I am using gdal warp to transform a geotiff between espg:4326 to epsg:3857, I have successfully done this previously on similar images but with the current geotiff it doesn't seem to be working. I am ...
1
vote
1
answer
71
views
Converting map from 0 to 360 to -180 to 180
this is second part of How to convert VICAR data to TIFF.
I was able to get the VICAR data into TIFF form by using VICAR2isis to make a .cub then gdal to make a tiff. (Yes i eventually found out I ...
1
vote
0
answers
56
views
Why use gdalwarp add -wo NUM_THREADS=ALL_CPUS is slow than don't use this option?
I try to test then reproject run time by use QGIS's toolbox-->Warp(reproject) on Windows 10, I find use gdalwarp add option -wo NUM_THREADS=ALL_CPUS is slow than don't use this option. this is the ...
3
votes
2
answers
188
views
Convert GeoPDF to 8-bit GeoTIFF using gdalwarp
I am using gdalwarp to convert GeoPDFs of processed maps to GeoTIFFs (for later stitching into a larger GeoTIFF) using the following command:
gdalwarp -t_srs EPSG:28356 -r cubic -cutline "...
0
votes
1
answer
62
views
Passing geometry dataset type to gdal.Warp() for a clip operation - order of operations and file type confusion
Trying to wrap my head around the correct order of operations and file type(s) required by gdal.Warp().
I have a valid geometry object created from the shapely object.unary_union.convex_hull method. ...
1
vote
0
answers
24
views
Bad values with gdalwarp and "aliased" coordinate transformation
I've strange values using gdallocationinfo on vrt with coordinate system transformation.
Here is my use case, maybe you can help to understand what I did bad or misunderstood
1) I've got a bunch of ...
1
vote
2
answers
131
views
QGIS 3.28 - cannot find GDAL commands / MacOS 12.5
I am having trouble linking QGIS to GDAL since a few versions (just updated to 3.28 hopping it would solve the problem but it didn't). For instance, a Wrap function returns this error:
GDAL command ...
1
vote
2
answers
89
views
Reprojecting each tif located in Folder Directory using GDAL in Python [closed]
Python beginner here. I am using python to reproject a folder of rasters (.tifs). Each raster has a projection that I want to change to WGS 1984 (espg:4326). I want to use an open-source GIS library ...
1
vote
1
answer
58
views
Avoid creating alpha band when using gdalwarp
I am using gdalwarp to reproject a greyscale GeoTiff tile from one projection to another so I can use it in a VRT. This is required as gdalbuildvrt won't accept files in multiple projections.
When I ...
0
votes
0
answers
22
views
gdalwarp() does not work to reproject OMI/Aura TOMS-Like Ozone data
I'm trying to convert the OMI/Aura TOMS-Like Ozone data provided in HDF5 to the GeoTIFF format using gdal_translate() in the R environment. gdal_translate() provided this unknown Extention.
so I want ...
4
votes
1
answer
89
views
gdalwarp sum resampling algorithm double counting at some specific resolutions
GDALWarp utility is double-counting my pixels in the sum algorithm for resampling in some cases. I don't believe it's something I'm obviously doing wrong.
I have a file, at 30 arcsecond resolution for ...
0
votes
0
answers
40
views
QGIS clip raster by mask layer: gdalwarp failed to start
I run into an error when I try to clip the land cover raster layer from ESA named ESACCI-LC-L4-LCCS-Map-300m-P1Y-2015-v2.0.7 for the Netherlands. I do this by saving the Netherlands polygon from the ...
0
votes
0
answers
144
views
GDAL Python - Warping Raster with Inaccurate GPS data (GCPs)
I'm currently working to georeference a raster (GTiff, output by Metashape) using inaccurate GPS points (serving as GCPs). I have some code that works well, in that it places the raster the in the ...
1
vote
1
answer
28
views
FAILURE : No target filename specified
path = '/mnt/Data/30_Stages_Encours/2023/ZonesHumides_Sarah/donnee_dezippee/sentinel2/2020/'
dir = os.listdir(path)
for i in dir :
# i = SENTINEL2X_20200415-000000-000_L3A_T31UES_C_V2-1
if i....
0
votes
0
answers
30
views
gdal.Warp reading 2d multichannel array incorrectly - loading bands (3) as the y axis
We have a python script that downloads NOAA GOES satellite NetCDFs, and processes the image datasets via Pillow\OpenCV. Then I convert it back to a GDAL dataset in order to translate and warp. With a ...
1
vote
1
answer
135
views
gdal.WarpOptions not accepting srcBands argument?
Defining as follows:
warpoptions = gdal.WarpOptions(format = 'GTiff', outputType = gdal.GDT_Float32, srcBands = [1], dstBands = [1], dstSRS = 'EPSG:4326', dstNodata = 0, resampleAlg = 'cubicspline', ...
4
votes
1
answer
227
views
How to get Amersfoort 28992 in QGIS?
I am currently analysing landscape change in the Netherlands, for which I am using maps which have Amersfoort 28992 as their standard CRS. However the most recent map uses Amersfoort 4289 instead.
For ...
0
votes
0
answers
55
views
Contours of a DSM render slowly
I am comparing two digital surface maps (dsms generated from drone photogrammetry) taken at two times. The contour lines render quite quickly for one (in seconds) and quite slowly for the other (1–2 ...
2
votes
1
answer
57
views
How to avoid gdalwarp converting to alpha band
I am running the following code in order to resample a single-band TIFF file, whose data are going from 0 to 100 with 255 being the nodata value
"gdalwarp -t_srs "+srs+" -tr "+pxsz+...
2
votes
0
answers
114
views
How to solve antimeridian split using GDAL
I'm fairly new to GIS and GDAL. But I need to solve the problem with anti-meridian. I have satellite images from Meteor-M2 that can span over -180 / 180 longitude.
I created .vrt file with correct GCP ...
0
votes
1
answer
47
views
Can't use Polygon to Clip Raster
I created a polygon to use to clip a DTM raster (DEMTiff), but the Clip by Mask tool gives an error:
GDAL command:
gdalwarp -overwrite -of GTiff -cutline temp.shp -cl temp -crop_to_cutline "D:\\...
1
vote
0
answers
35
views
Julia - Write KMZ file with transparency
I'm trying to save a .kmz file from a dataset representing a DEM of an area stored in a GeoTiff file.
My code is reading the DEM, transforming the CRS to UTM, and defining a subset based on ...
1
vote
0
answers
219
views
Run gdalwarp to rectify a raster (RPC embedded) with DEM but got error: Too many points failed to transform, unable to compute output bounds
I've been trying to rectify an unprojected raster (with RPC model) with digital surface model (DSM) by using gdalwarp command. Here is the command:
gdalwarp \
-co BIGTIFF=YES \
-co "TILED=YES&...
1
vote
1
answer
72
views
Warping ArcGIS 9.2 georeferenced TIFFs using GDAL warp
I have hundreds of very large TIFFs which were georeferenced 15+ years ago using ArcGIS 9.2. Viewing them in ArcMap 9.2 the images are properly located and sized/warped on the map. However, they ...
3
votes
2
answers
183
views
Command line/Python (Linux) GDAL Warp result much bigger than in QGIS (Windows) [duplicate]
I'm trying to run the following code to downsample a ~3GB GTiff while coregistering it to a reference extent:
def coreg_downsamp(source, destination, output_type="GTiff"):
#essentially ...
0
votes
0
answers
18
views
Example of explicit transform using gdalwarp
Could anyone give me an example for using the -ct option in gdalwarp? The instruction in GDAL docs is a little opaque to me. How about a urn:ogc:def:coordinateOperation:EPSG::XXXX?
0
votes
0
answers
134
views
Coordinate transformation with GDAL replicating Geotools Java package transform
We are developing an earth observation data processors utilizing GPU.
One of the operation trying to achieve is the coordinate transform. For some rasters the output CRS intended would be some local ...
2
votes
1
answer
178
views
Output filesize is too big when using gdal warp in Python
I'm trying to process a very large directory of smaller geotiff files into a single mosaic. The cumulative size of the entire set of files is < 300mb but there are about 6,000 individual files ...
0
votes
0
answers
71
views
Converting into correct format using GDAL
I am trying to use gdal to convert landcover .tiff files into .asc files in the AAIGrid format. I am using them for an application which seems to only accept formats where the xllcorner and yllcorner ...
0
votes
0
answers
111
views
Converting NetCDF to GeoTIFF WGS 84
I would like to convert the next NetCDF file to a GeoTIFF WGS 84.
NetCDF file gdalinfo output:
Driver: netCDF/Network Common Data Format
Files: /home/me/netcdf_file.nc
Size is 512, 512
Metadata:
...
0
votes
0
answers
148
views
Convert HDF4 MODIS to TIF
I need to convert an hdf4 file by MODIS, the MOD13A3 v006 (NDVI Index). In particular, I'm interested in the block "h18v04", which is north-west Europe. You can download the example file I'm ...
0
votes
0
answers
167
views
Shift after reprojection from EPSG:3035 to other EPSG CRS (25832 / 31467 / 4326)
I'm trying to reproject a raster layer (EPSG:3035) to EPSG:25832.
In QGIS I used gdalwarp with the following command:
gdalwarp -s_srs EPSG:3035 -t_srs EPSG:25832 -tr 30.0 30.0 -r near -of GTiff ...
0
votes
1
answer
69
views
Is it possible to warp a raster using -te in other projection?
I have some files in EPSG:4326 and some in "ESRI:102027".
I use warping with -te argument in EPSG:4326:
gdalwarp -overwrite -of GTiff -t_srs srs -te xmin ymin xmax ymax -te_srs "EPSG:...
1
vote
0
answers
144
views
How to convert elevation data from WGS84 ITRF2008 to NAD83 CGG2013a
I need to convert elevation data captured in WGS84 ITRF2008 to NAD83 CGG2013a.
I came across with GPS-H from NRCAN and downloaded the Geoid Model file (CGG2013ai08.byn) from here.
Although I use a ...
0
votes
0
answers
46
views
How to extend raster edges when rescaled using the warp tool in QGIS
I am rescaling a raster from a resolution of 100m to 2000m using the warp tool (method cubic spline) in QGIS, it looses a lot of information around the edges which I would like to retain. Is there a a ...
0
votes
0
answers
301
views
Vertical transformation of DEM from EPSG:3035 to EPSG:4326
I am trying to convert a DEM file from EPSG:3035 (ETRS89-extended / LAEA Europe) to EPSG:4326 (WGS 84) with gdalwarp. EPSG:3035 uses the European Terrestrial Reference System 1989 ensemble datum where ...
1
vote
1
answer
127
views
Warp Natural Earth raster without losing resolution
I'm trying to make a world map using Pseudo-Mercator projection (3857) and Natural Earth's 1:10m raster in QGIS. However, Natural Earth apparently doesn't like Mercator, so it cuts the extent of the ...
0
votes
1
answer
209
views
gdal.Warp is not creating a TIFF file (very big raster)
My raster size is more than 6 Gb. It was created from SRTM tiles, so I wasn't even able to save it as TIFF, hence I used vrt format.
Now I want to translate it to .map, so I use:
warpopt = gdal....
1
vote
1
answer
74
views
Trying to gdalwarp a TIFF results in mostly empty image with data at ends
I'm trying to use gdalwarp to project an FAA source chart in EPSG:9802 (Lambert Conic Conformal) to EPSG:3395 (World Mercator). The issue is that part of the data is east of longitude -180, and part ...
2
votes
0
answers
190
views
Speeding up WarpedVRT in rasterio (change CRS + resample)
Given a source GeoTIFF of some unknown resolution (better than the target), what is the quickest way to re-project it in-memory to the desired GSD? Masking in the input data should be preserved.
...
2
votes
1
answer
67
views
Calculating common area between shapefiles using GDAL
I'm looking for a way to get the intersection area between my shapefiles.
Basically, I got three shapefiles from my .tif files like this:
gdaltindex a.shp a.tif
gdaltindex b.shp b.tif
gdaltindex c.shp ...
0
votes
0
answers
210
views
Using Gdal warp on shapefiles with GDAL Java Bindings
I am trying to write a method in Java that transforms from one CRS to another using GDAL's Java Bindings. Specifically from 3857 to 4326 (and a few others will be common to my use cases).
Here is my ...
2
votes
1
answer
339
views
Summing two rasters of different extents in QGIS
In QGIS 3.24 I've got two rasters, one that is 1079x1079 and another that is 83x84, but they both have a very specific pixel size, 0.002083333333333333304,-0.002083333333333333304.
I'd like to add ...
1
vote
1
answer
1k
views
gdal.warp in Python script
I've a .tif file in EPSG:3426 and i'd like to convert it in UTM, zone 33, scale: meters.
In gdal I use:
gdalwarp -t_srs "+proj=utm +zone=33 +datum=WGS84 +units=m" image_input.tif image_utm....
2
votes
1
answer
452
views
Convert image from EPSG:4326 to EPSG:3857
I'm new to GDAL and map projections, I need to convert equirectangular (EPSG:4326?) projection to "pseudo Mercator" (EPSG:3857?) projection used by 90% web-based maps (which is simple this 0/...
2
votes
0
answers
174
views
Convert EPSG and vertical datum of a DEM file
I have DEM files that are in EPSG:4326 - WGS 84 and use latitude and longitude for coordinates. The elevation units are in meter but use the WGS84 ellipsoid with WGS84 elevation. I want to convert ...