All Questions
528
questions
0
votes
0
answers
42
views
GDAL Warp extremely slow on Python API
I am trying to export one VRT to TIFF using GDAL Warp. My VRT has some scale, offset, reprojections, etc, but nothing too crazy.
I have observed that trying to perform the warp from the Python API is ...
1
vote
0
answers
31
views
GDAL - From jp2 image with rpc file to Orthophoto geotiff and keep the image rotation
We are trying to convert a jp2 image with rpc to orthoPhoto geo tiff. We do the following process. The problem is that in the ortho.tiff file the image don't preserve it's rotation angle from the jp2 ...
- 11
0
votes
0
answers
26
views
Add transparency to raster bands where did not exist before
My original radar image is a tiff image with one single band.
Metadata:AREA_OR_POINT=Area
Image Structure Metadata:
COMPRESSION=DEFLATE
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( ...
- 1
0
votes
1
answer
39
views
Stripes in map are missing after gdalwarp and gdalbuildvrt
I have a problem with some GeoTIFF:s that are in EPSG:3006 and I need them in EPSG:3857. The problem is that after gdalwarp, and gdalbuildvrt I get transparent stripes in the vrt map
This is one of ...
- 175
1
vote
0
answers
43
views
Access Violation following rasterdataset.Close() after GdalWarp, GDAL C# bindings
I have the following C# code to mosaic multiple rasters using GDAL C# bindings and GdalWarp method:
try
{
Gdal.AllRegister();
// Open all the raster files in ...
0
votes
1
answer
40
views
Transforming ortho (.tiff) file from one epsg to other epsg
I am using GDAL to transfrom ortho from one EPSG to other EPSG.
In my case I have used EPSG:32610 which is in (m) and transformed into EPSG:2227 which is in (us-ft). I am able to transform it using ...
0
votes
0
answers
55
views
Clip a NetCDF file using gdal_translate or gdalwarp
I'm trying to clip this .nc file using either gdal_translate or gdalwarp. However I don't think I'm getting the correct output. I'm able to clip the file in Python using either of the following ...
- 141
1
vote
0
answers
55
views
Creating Orthomosaic with GDAL
Platform: Windows 11
1. Aim trying to achieve: I have some very high quality multispectral data (drone multispectral imagery), from which I am trying to create an orthomosaic. My multispectral data ...
- 11
0
votes
1
answer
46
views
Convert RGB(A) GeoPDF to Grayscale GeoTIFF with GDAL
I'm using the GDAL command line tools gdal_warp and gdal_translate to convert RGB(A) GeoPDF images into GeoTIFF images. Unfortunately I found no config option that changes the RGB(A) color space to ...
- 5,857
1
vote
1
answer
32
views
Gdalwarp assigning zero values outside input raster extents
I'm attempting to combine several GeoTIFF files into a single mosaic raster using gdalwarp:
import os
os.environ['AWS_NO_SIGN_REQUEST'] = 'YES'
file1 = "https://dea-public-data-dev.s3-ap-...
0
votes
0
answers
21
views
Coordinate projection transformation Mars image from Ocentric Equirectangular to Transverse Mercator in Python
The task is to re-project the Mars image from Ocentric Equirectangular projection to Transverse Mercator.
I'm trying to do this in gdal in python.
In terms of Transverse Mercator projection, there is ...
0
votes
0
answers
32
views
Transformation of a point with GDAL
I know the place where the image was captured and I've made pixel to lat-lon correspondences (GCPs).
After running the following GDAL commands a .tif file is created.
gdal_translate -of GTiff -gcp {...
0
votes
0
answers
24
views
Concatenating several very tall JPG-in-TIFF images using GDAL
Preface: I am fairly new to GIS, and as such - I haven't wrapped my head around the terms. I'm trying to vertically concatenate several tall JPEG-in-TIFF files into a single file. I've attempted using ...
- 101
0
votes
0
answers
31
views
Translating (cropping) a GeoTIFF file to a specific extent but ending up with a different final extent
I'm trying to crop a larger GeoTIFF file to a smaller extent.
The darker GeoTIFF file in the middle has the extent I want to reproject to. I've used the following command to do the translation:
...
- 1
0
votes
0
answers
20
views
Produce ISEA "Icosahedral Snyder Equal Area" GEOTIFF/raster with gdalwarp
Is an icosahedral projection impossible to produce with PROJ/GDAL ?
Produce ISEA "Icosahedral Snyder Equal Area" GEOTIFF/raster with gdalwarp.
I am trying to use gdal tools to convert Plate-...
- 35
0
votes
0
answers
49
views
Reproject pseudo Plate-Carree to "Icosahedral Snyder Equal Area" using GDALWARP
Reproject pseudo Plate-Carree to "Icosahedral Snyder Equal Area" with gdalwarp.
I am trying to reproject a plate-carree to ISEA "Icosahedral Snyder Equal Area".
I suspect ISEA ...
- 35
1
vote
2
answers
79
views
Use gdal_translate on an orthographic PNG image of the Earth centered at Lat=+50N Lon=+5E to proper orthographic GeoTIFF format
How do I gdal_translate an orthographic PNG image of the Earth centered at Lat=+50N Lon=+5E to proper GeoTIFFf format ?
I cannot work out how to specify the limits of a full disc. (diameter of Earth ...
- 35
1
vote
0
answers
52
views
Cropping VRT with Shapefile to create new VRT with no dependecy from input
I'm new to work with gdal and I'm trying to crop a VRT-File with a shapefile to create a new VRT-File that is independent from the input.
I use gdalbuildvrt [...] to create the base VRT-File from ...
- 11
0
votes
1
answer
79
views
GDAL: GFS weather model reprojection issues
I'm trying to reproject a grib file using gdal.
The crib file has the following,relevant, info:
GEOGCRS["Coordinate System imported from GRIB file",
DATUM["unnamed",
...
0
votes
0
answers
110
views
Error trying to translate raster image from EPSG:3857 to EPSG:3068 with GDAL with subwindow
I'm new with this GIS stuff and starting work with GDAL.
I'm using this TMS dataset and trying to project to a specific subwindow like this:
<GDAL_WMS>
<Service name="TMS">
...
0
votes
0
answers
44
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 ...
- 535
0
votes
0
answers
55
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 ...
- 111
1
vote
0
answers
101
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 ...
- 11
0
votes
0
answers
67
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
159
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 ...
- 63
0
votes
1
answer
184
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. ...
- 109
1
vote
0
answers
35
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 ...
- 61
1
vote
2
answers
416
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 ...
- 160
1
vote
2
answers
263
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 ...
4
votes
1
answer
187
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 ...
- 323
1
vote
1
answer
35
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....
- 33
1
vote
1
answer
326
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', ...
- 13
0
votes
0
answers
92
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 ...
- 83
2
votes
1
answer
92
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+...
- 69
2
votes
0
answers
205
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 ...
- 121
0
votes
1
answer
55
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
1
vote
0
answers
42
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 ...
- 11
1
vote
0
answers
395
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&...
- 258
0
votes
0
answers
162
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 ...
- 103
0
votes
0
answers
88
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 ...
- 1
1
vote
0
answers
221
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 ...
- 1,745
0
votes
0
answers
452
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 ...
- 329
0
votes
1
answer
249
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
103
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 ...
- 113
2
votes
0
answers
278
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.
...
- 151
2
votes
1
answer
91
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 ...
- 31
0
votes
0
answers
268
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 ...
- 297
1
vote
1
answer
2k
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....
- 31
2
votes
1
answer
615
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/...
- 141
1
vote
0
answers
58
views
gdalwarp raster reprojection is generating pixelated edges
I have a raster image, which is satellite data, that I would like to reproject using gdal on the command line.
The files I am using can be found here:
the unprojected image (.png)
.wld file
.msinc ...
- 93