Questions tagged [gdalwarp]

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

Filter by
Sorted by
Tagged with
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 ...

1
2 3 4 5
12