Questions tagged [gdalwarp]
A Geospatial Data Abstraction Library (GDAL) image mosaicing, reprojection, and warping utility.
592
questions
0
votes
0answers
38 views
Clipping raster using multiple shapefiles using gdal in Python
I am currently using batch process in QGIS based on clip raster by extent to clip geotiff raster layers based on different shapefiles. I am trying to clip many raster layers from the same region with ...
0
votes
0answers
15 views
Select first valid datasource in VRT
When reading data from a VRT, gdal appears to interpolate values from all available rasters. Is there a way to configure the VRT in such a way that it selects the first valid raster in the dataset, or ...
0
votes
1answer
43 views
Converting JP2 to GeoTIFF QGIS 3.16 Win 64
I have hundreds of USGS NAIP imagery GeoTIFFs that are in JP2 format, EPSG:3857.
I am using QGIS 3.16 to convert them to GeoTIFFs ("Rendered Images") as EPSG:4326, but I can only process one ...
0
votes
0answers
24 views
Going from degrees to meters in projection conversion using GDAL?
I have raster data in a Lambert-Conformal conic projection created from a NetCDF file by running sub-processes in a Python app that runs the following (Linux) command line argument:
gdal_translate -ot ...
0
votes
0answers
13 views
How can I use gdal warp on a .sid to create an .tif of the same resolution?
I’m trying to convert a .sid to .tif in OSGeo and was left this script by a former co-worker. I’m not very savvy in scripting, so I’m unsure of my translation.
First, I set the directory. The ...
0
votes
0answers
13 views
Merging mbtiles databases
I have a large set of hires GeoTif images.
I am trying to build a .mbtiles database of these images. However, using gdalbuildvrt and osgearth_conv makes me hit a wall. I have also tried using gdalwarp ...
2
votes
0answers
74 views
SRTM HGT to GeoTIFF conversion using GDAL
I have (for example) SRTM-1 file n48e083.hgt with resolution 3601x3601.
What is right way to convert it to GeoTIFF?
gdal_translate -srcwin 0 0 3600 3600 -a_ullr 83 49 84 48 n48e083.hgt n48e083.tif
or
...
0
votes
0answers
15 views
Runtime error GDALWarp command in C++ programme
I am trying to use GDALWarp in VS2019: A snapshot of the code is
#include <gdal_utils.h>
#include <gda_priv.h>
#include <cpl_string.h>
int main() { GDALDataSetH* Dataset=GDALOpen();...
2
votes
1answer
61 views
Reprojecting raster changes pixel scale (makes it darker)
I have an HH polarised SAR image as follows:
This is how it looks in QGIS in the default EPSG:4326 projection.
Now I would like to reproject this raster to ESRI:102020 (South Pole Lambert Azimuthal ...
0
votes
0answers
22 views
Converting a UTM .bil file to a lat/long TIFF
I have many 3 arcsec (90m) SRTM DEM files in BIL format in UTM projection that I need to convert into 1 arcsec (30m) lat/long TIFFs. I have tried both gdal_translate and gdalwarp to no success.
First ...
4
votes
0answers
68 views
Reproject WGS84 raster file to topocentric LTP-ENU using GDAL with a custom proj string
I have some raster files which are given in WGS84. As an example, here are some gdalinfo for one of those:
Driver: GTiff/GeoTIFF
Files: dem_r.tif
Size is 2000, 911
Coordinate System is:
GEOGCRS["...
2
votes
0answers
23 views
Calculating inverse polynomial transforms for pixel sampling when map georeferencing
In map georeferencing, one typically chooses a set of ground control points that represent pixel locations in the map image along with their real world geographic coordinates, which are then used to ...
1
vote
0answers
83 views
How can I align the pixels with using GDAL and Python?
I am trying to use Python/GDAL to align different GeoTIFFs that I warped (gdal.Warp) to the same SRS. My problem is one of the images is rotated, the pixels do not overlap perfectly, which is ...
6
votes
1answer
75 views
Difference between QGIS export and gdalwarp
Does QGIS export/save as use a different projection system/process than gdalwarp? I am working with MASIE ice coverage data available here. My goal is to automate the process of reprojecting the daily ...
2
votes
1answer
28 views
Using GDALwarp with not-to-scale maps
(This question was asked on another forum 3 months ago but received no response. I'd really welcome some help getting over the problem)
My input map is a JPG scan of a not-to-scale map. I would like ...
0
votes
0answers
16 views
Resizing a Raster (Row and column values) using GDAL in Spyder
I'm trying to extract the LST values from MODIS daily LST product for 3 months. So I'm Using Timeseries analysis to extract multi raster values to a point shapefile. The main problem here is that the ...
0
votes
0answers
27 views
How to use GDAL to write lat,lon,z from NetCDF to GeoTIFF?
I have the following NetCDF file:
Format:
netcdf4_classic
Global Attributes:
history = 'FERRET V7.43 (optimized) 9-Jun-20FERRET V7.43 (optimized) 19-Jun-20FERRET V7.43 (...
3
votes
0answers
49 views
Converting coordinate system produces distorted image
I have a GeoTIFF file produced from a LAZ file. I would like to convert its coordinate system to the WGS84. However, after the conversion the resulting image is skewed (has a perspective distortion).
...
0
votes
0answers
8 views
Performing ellipsoid to spheroid conversion on a DTM using gdal
I have a DTM that is referenced to an ellipsoid (equatorial radius is not equal to the polar radius) and I want to convert the elevation values to the spheroid (equatorial radius = polar radius).
How ...
2
votes
0answers
28 views
Maximize CPU and memory usage when using gdalwarp
I'm using gdalwarp to clip and warp a large (but not huge) GeoTIFF (6.2 GB total size). The GDAL command I am using is:
gdalwarp --config GDAL_CACHEMAX 80% -wm 80% -multi -wo NUM_THREADS=ALL_CPUS -...
0
votes
0answers
24 views
gdalwarp “RPC_DEM” argument not working
I'm trying to orthorectify an image with the following commands:
gdalwarp -rpc -to RPC_DEM=dem.tif input.tif output.tif
gdalwarp -rpc input.tif ortho_output.tif
I'd expect the "ortho_output.tif&...
0
votes
0answers
12 views
How to change the resolution of NetCDF data and start at a specific point in GDAL
I have temperature data which is in NetCDF format and lat=1801 lon=3600 and I want to convert it to a NetCDF file with lat=360 and lon=720. In other words, I want to change the resolution of my data ...
0
votes
1answer
82 views
Reprojection with gdal.Warp doesn't work
I am trying to reproject a raster. I have tried gdal.Warp and it doesn't work.
RuntimeError: Cannot find coordinate operations from 'EPSG:6257' to 'EPSG:6686'.
Working in Python 3, Windows 64bit
...
0
votes
0answers
107 views
COG Driver not recognized in GDAL 3.2.0
I want to generate a Cloud Optimized Geotiff (COG) using Python (with Anaconda environment). I wanted to use GDAL 3.2.0 since it has a COG Built in Driver since 3.1. Initially I cannot update GDAL to ...
0
votes
1answer
75 views
Converting shapefile to raster by passing specific column/variable using GDAL in Python
Followed the steps from this link to create raster from shapefile using GDAL
Convert an OGR File to a Raster
but it is returning an empty raster.
Can someone tell how to define a particular column ...
2
votes
1answer
50 views
GeoTIFF is stretched sideway when translating CRS to EPSG:4326 from EPSG:32654
I want to convert a GeoTIFF's coordinate reference system to EPSG:4326 from EPSG:32654. But the output image is stretched to the side.
I have tried two different methods (gdalwarp and rasterio), but ...
0
votes
1answer
40 views
gdalWarp : Assigning more than one nodata value
I am working with a single band raster file in pyQgis and want to assign multiple Nodata values while reprojecting. I tried doing this but it took only the first value 32761 as Nodata value:
gdal.Warp(...
1
vote
1answer
90 views
Issue creating tiles from GeoTIFF which crosses the 180th Meridian
UPDATE:
I've included more information and the scripts that were used to generate the output.
Example results can be found here at the links below and are all generated from the same original GeoTIFF. ...
1
vote
0answers
27 views
GDAL warp & GDAL contour inconsistent results
I am creating 2 output contours based on the same input DEM, one is clipped to a smaller grid feature (green line) while the other is clipped to a larger irregular polygon (blue line).
The problem is ...
5
votes
0answers
97 views
Georeference and warp a drone image based on sensor orientation and location
I have been flying some drone surveys over the ocean and need to properly project and georeference my images. I have all the information I think I need: lat, lon, altitude, yaw, pitch, and roll along ...
1
vote
1answer
30 views
Using gdal.Warp() on netCDF file. Simple case fails
I'm using the python API to GDAL to try to do a simple geographic warp of a netCDF4 satellite file to a small sub region. My python program fails to produce any output but does not throw any ...
0
votes
1answer
82 views
More efficient processing of many XYZ tiles into a merged response
I've built a mobile mapping device for my Search & Rescue group using a Raspberry Pi (3b+) and its primary purpose is to serve a cache of XYZ tiles via its web server, over its WiFi network. The ...
1
vote
1answer
93 views
gdalwarp does not use lanczos resampling on image with alpha channel
I want to warp Geotiff image of shaded relief (RGBA) in projection EPSG:8353 to the projection EPSG:3857 to use it in Mapnik web tiled map. I want to use lanczos resampling, but the result image looks ...
2
votes
1answer
38 views
Warping a Polar Stereographic image to UKNG projection
I have a North Polar Stereographic image over the UK. The top left and bottom right coordinates are: (-17.644 E, 59.683 N) (9.231 E, 45.25 N).
I want to reproject this to UKNG (aka BNG) projection.
...
0
votes
0answers
34 views
Remapping a raster map using GDAL in Python
I have a numpy array of data in Python that is a processed version (outliers set to zero, for example) of data from an Albers Conical United States (US) map:
PROJCS["Albers Conical Equal Area&...
0
votes
0answers
15 views
Why is GDALWARP making straight lines in hi-res ortho imagery look sawtooth? How to Fix?
I downloaded six hi-res (6") imagery tiles from USGS. Used gdal_translate to add State Plane Coordinates (SPC) projection to the tif. Used gdal_merge to combine the six GTiffs into a single GTiff....
0
votes
1answer
74 views
How to get the EGM96 and EGM2008 elevations of a list of WGS84 lat/lng coordinates from a file using GDAL?
I'm trying to do the following:
Given a list of WGS84 latitude / longitude pairs in a text file (one pair per line), I'd like to compute the EGM96 and EGM2008 elevations corresponding to those points ...
2
votes
1answer
115 views
Reprojecting raster (GeoTIFF) using PyQGIS in standalone script
I'm trying to reproject a GeoTIFF file within a PyQGIS standalone script.
Using QGIS 3.14 the steps would be:
resulting in the following process history and the output_raster.tiff being successfully ...
2
votes
0answers
28 views
How to use gdalwarp to reproject a raster without changing band labels?
Procedure
I used QGIS to reproject a raster data. QGIS does this via gdalwarp like so:
gdalwarp -overwrite -t_srs EPSG:4326 -of GTiff /path/to/input.tif /path/to/output.tif
Comparison of input and ...
0
votes
0answers
17 views
inconsistent behavior of identical gdalwarp operation on two different GTiff files (one works, the other has no effect)
I have 2 GTiff images in WGS84 EPSG:4326 CRS that I am trying reprojecting to EPSG:3857. The files are available from dropbox:
File 1: file_01.tif 3600 x 1800, 3 Band RGB Pixel interleaved
File 2: ...
0
votes
1answer
41 views
Georeference raster in WGS84 and project in North azimuthal in QGIS
I have a raster (JEPG) image from a full planet surface (in this case Jupiter).
I use QGIS and its plugin Georeferencer to add control points on the image. I use the four corners. Then I use Gdal to ...
0
votes
0answers
18 views
gdal.Warp ignores some pixels when reprojecting it
I have a tif in EPSG3035, pixel size 100x100m (very important for the study).
I use gdal.Warp to reproject the tif to EPSG4326 for getting lat/long of pixels along with their values. i am forcing the ...
0
votes
0answers
40 views
How to sum pixels of two overlapping rasters with different extents
I have two DEM rasters which overlap and have different extents. I would like to create a new raster which sums the pixels where there is an overlap.
I tried using gdal_calc first by doing A+B but ...
0
votes
0answers
28 views
Gdalwarp - Clip Raster by Mask - Missing Pixels (Cutline)
I am clipping two - partly overlapping - geotiffs with Gdalwarp by the same shapefile, wich I want to merge in a next step. This works well, but after merging the rasters, I noticed that some pixels ...
2
votes
0answers
62 views
gdal warp reproject creates artifact in raster
I’m trying to re-project a tif-file using gdal warp:
fn = "\dtm.tif"
out = "\dtm_25832.tif"
gdal.Warp(out, fn, dstSRS="EPSG:25832")
The output file is re-projected, but ...
2
votes
0answers
44 views
Converting vertical datums GeoTIFF to GeoTIFF
Anyone know where I've gone wrong here?
gdalwarp input: geotiff that is in NZVD2009 [EPSG:4440]
gdalwarp output: geotiff that is in NZGD2000 [EPSG:4167]
(gdal) C:\>gdalwarp -s_srs "+proj=...
0
votes
1answer
69 views
TypeError when using gdal.Warp in Jupyter Notebook
I try to use the gdal python binding to project a raster in World Mollweide.
import gdal
map_raster_proj = pm.getTmpDir() + Path(map_raster).stem + '_proj.tif'
input_ = gdal.Open(map_raster)
gdal....
1
vote
1answer
81 views
gdalwarp DEM with different resolution extents
I'm trying to convert DEM source (LIDAR is original source, in ASC format) to a custom internal format that has specific resolution and alignment values. I also want to generate different output ...
1
vote
0answers
35 views
Speed up gdalwarp
I'm cropping a large amount of geotiff images by use of geojson masks:
gdalwarp \
-of GTiff \
-co COMPRESS=JPEG \
-co TILED=YES \
-co BLOCKXSIZE=512 \
-co JPEG_QUALITY=85 \
-co NUM_THREADS=ALL_CPUS \
...
2
votes
1answer
148 views
Cropping raster file with GDAL warp C++ API
I need to crop a DEM raster file with a shapefile. This shapefile has many polygons. I need to calculate several variables after the crop operation.
The algorithm is written in C++ so I want to call ...
