Stack Exchange Network

Stack Exchange network consists of 174 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

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

1
vote
1answer
48 views

Error using QGIS 3.2.2 Warp (Reproject) tool with custom coordinate system

I'm using the QGIS Warp (Reproject) tool to reproject a raster from WGS84 to a Lambert Conformal Conic custom projection in QGIS 3.2.2. The call to the GDAL/OGR console is the following: gdalwarp -...
0
votes
1answer
19 views

gdalinfo and other applications not running through Anaconda Prompt

I am using GDAL 2.0.3(installed through wheel) in Anaconda 4.2.0 with Python 3.5 installed. I am able to import GDAL in my code. But whenever I am trying to run gdalinfo and other applications(such ...
0
votes
1answer
34 views

GDAL - Unable to convert to Web Mercator using gdalwarp

I am making geoTIFs of weather data (for a Mapbox overlay) and I am converting them to web mercator for use on the map. The images were generated from the same source file and both images have the ...
1
vote
1answer
44 views

Efficient downsampling with GDAL?

I have a 350go GeoTIFF file with DEM data for the entire world. (a mix of SRTM, EUDEM, etc. ...) This TIFF file is tiled with a block size of 256*256. I would like to downsample this file to a lower ...
1
vote
0answers
17 views

reprojection and write grib files

I'm working on GRIB file, In input, I have the file with data from 0 to 360°. As gdalwarp does not write GRIB file (only gdal_translate), I use gdal warp to projct data from [0,360] & [-90;90] to ...
1
vote
1answer
27 views

Why doesn't QGIS georeferencer support gdal scripting for projective transformation?

I am georeferencing drone images with QGIS, and the transformation method that gives the best results is the projective one. As I have many images, I would like to write a gdal script to automate the ...
2
votes
1answer
22 views

GDAL processing of GRIB2 raster prints hundreds of unsuppressable warnings

I am processing some GRIB2 raster files from NCEP (weather model stuff). An example file would be: http://nomads.ncep.noaa.gov/pub/data/nccf/com/hiresw/prod/href.20180906/ensprod/href.t00z.conus.mean....
0
votes
0answers
57 views

GDAL interpolation ignores edges with nodata vals

I have a grib file with weather data in 0.25-degree resolution and I want to upsample its higher resolution. I used gdalwarp: gdalwarp -t_srs EPSG:3857 -of GTiff -wo SOURCE_EXTRA=1000 -wm 2000 --...
1
vote
1answer
21 views

gdalwarp: differentiate between empty-in-source and outside-of-cutline?

I'm using gdalwarp to clip a polygon out of a map layer using --crop_to_cutline. The original map layer has a lot of empty pixels (well, transparent / alpha=0, to be precise). As it stands, my output ...
1
vote
1answer
42 views

Very large VRT file generated by gdal after cutline

In Python, I generate a virtual raster using gdal warp with a cutline as per below gdal.Warp(warpedFilePath, ds, format = 'VRT', cutlineDSName = csvFilePath, srcNodata = 0, dstAlpha = True, ...
1
vote
0answers
11 views

Creating polar-stereographic projection of HDF data subset?

I have an HDF4 file. I would like to create a geotif file containing a polar-stereographic projection of one of the HDF4 data sets. I've tried using gdalwarp. This works for some projection types, ...
2
votes
1answer
52 views

GDAL reprojection 4326 to 3857 creates gap

I have an EPSG:4326 tif covering the world that I need to reproject to EPSG:3857 in order to use it with Mapnik to create a map. When displaying the original tif and leaving Mapnik doing the ...
0
votes
0answers
19 views

gdal.Warp workingType option default

I have been using gdal.Warp() in Python 3.6 for some time now and have just noticed the workingType option (see here for gdal.WarpOptions()). From its definition: workingType --- working type (...
0
votes
0answers
20 views

remove white outline artefacts after gdalwarp/gdal2tiles

I am struggling with the similar problem as described here but with some differences: every raster I project twice using GDALWarp - into their initial projection (raster + .tfw) and then to EPSG:3857 (...
0
votes
1answer
14 views

GDALReprojectImage (C++) quiet

I am using the GDALReprojectImage method from gdalwarper.h. Is there a way to suppress the output? The function has GDALWarpOptions * psOptions as an input but there in no quiet option. I have seen ...
0
votes
1answer
33 views

Extracted GeoJSON from TIFF

If there are a way to do call("gdal_contour -i 10 -f GeoJSON".split(' ') + [file1, file2]) as Python call method of some library? Not via command line call and intermediate file storing?
1
vote
1answer
31 views

gdal_translate and gdalwarp do not close when converting an ecw file to GeoTIFF

I'm trying to convert a portion of a fairly large ECW file to a GeoTIFF using either gdal_translate or gdalwarp. The problem I'm having is that the GeoTIFF is generated but gdal_translate and ...
1
vote
0answers
21 views

gdal: WMTS is working with 256x256 tiles, but not with @2x (512x512) tiles

I'm working on using imagery from a web map tile service which I've cached locally as a gdal dataset. I have both 256x256 tiles and 512x512 tiles; the latter are the "retina" versions of the former; i....
0
votes
0answers
9 views

rasterio.warp.reporject: works incorrectly on big regions?

I want to reproject a big region of territory to ortho projection. For example, I take a raster of whole Earth in longlat projection and reproject it to ortho (make a view of whole Earth in ortho ...
0
votes
1answer
25 views

Using gdalwarp to merge bands into a mosaic created by gdal_merge?

I have created a mosaic with gdal_merge.py using the following command: gdal_merge.py -createonly -separate -co PHOTOMETRIC=RGB -o merged.tif LE07_L1TP_046030_20180810_20180810_01_RT_b3.tif ...
2
votes
0answers
58 views

gdal_warp with huge tif: integer overflow

I retrieved the SRTM 1" data and merged it in one huge tif. Now I want to warp it. It's working fine for low resolution (eg 5000, 1000, 500), but it doesn't for 90m resolution. I'm using this command: ...
0
votes
2answers
35 views

gdalwarp not clipping to index grid

I am using the 4 corners of a grid to clip a raster to it with gdalwarp with the following code: pkdumpogr -i index.shp -n left1 -n bottom1 -n right1 -n top1| while read Block_ID left1 bottom1 right1 ...
0
votes
1answer
17 views

gdalwarp for clipping images using multiple polygon in a bash shell for loop

C:\OSGeo4W64\bin\gdalwarp.exe -dstnodata 0 -q -cutline D:\hillsborough_Images_work\split_Polygon\120-140k\Hillborough_FINAL_140000_region_IGD_120570003d29ed0dc.shp -crop_to_cutline -dstalpha -tr 0....
0
votes
0answers
19 views

gdal: apply orthorectification to ogr dataset

I'm trying to process Digital Globe L2A image with gdal utilities (e. g. this one: https://www.digitalglobe.com/resources/product-samples/washington-dc-usa) I am using SRTM digital elevation model to ...
0
votes
0answers
18 views

What transformation is used for EPSG 28409 in GDAL?

I used gdalwarp to convert a raster file in EPSG:4326 to one EPSG:28409 I believe that using this method doesn't apply a full transformation as there is no option to select between the numerous ...
0
votes
0answers
25 views

Output arrays from Rasterio masks of Sentinel-1 prodcuts from the same shapefile are not the same shape

I am masking Sentinel-1 data using rasterio.mask.mask. I had first reprojected all Sentinel-1 products I have downloaded using command line gdalwarp as Sentinel-1 data is georeferenced but no ...
0
votes
0answers
47 views

Raster color after crop image by mask with rasterio

I have stacked 4 bands of Sentinel which have 10 meter of spatial resution. I I did it with command above: command = 'gdal_translate SENTINEL2_L1C:{scene_file_name}.SAFE/' \ 'MTD_MSIL1C....
1
vote
0answers
47 views

Problems with S3 georeferencing using Python-GDAL [closed]

I have to work with many Sentinel 2 and 3 images, with two different algorithms. I have to treat both images separately but only with Python and GDAL. Concerning Sentinel-2 images, I download them on ...
0
votes
0answers
12 views

GDAL WARP where geotif goes from west to east hemisphere [duplicate]

I'm using gdalwarp to reproject an elevation geotif downloaded from http://srtm.csi.cgiar.org/SRT-ZIP/SRTM_V41/SRTM_Data_GeoTiff/srtm_70_22.zip Using the following command gdalwarp -s_srs EPSG:4326 ...
0
votes
1answer
42 views

GDAL-Georefrencing already Georefrenced single band Raster

I am trying to georeference already georeferenced single band raster data using GDAL with four GCP's. Input image is having Coordinate system GCS- WGS1984 (EPSG:4326). Following Command is run on ...
0
votes
1answer
55 views

FAILURE: No target filename specified

Each time I run gdalwarp -te 212051.9319179999874905 214051.9319179999874905 183962.2946890000021085 Quada.tif Clipped.tif I get the following error: FAILURE: No target filename specified. Any ...
0
votes
0answers
17 views

Geo-tagging each pixel within an image using GDAL or Python and creating a jpg

I have a large multi band image, in ENVI format that was taken using an aerial platform. Three files are generated out of that, ENVI file, hdr file and a nav file. The system provides lat/long co-...
-1
votes
1answer
102 views

Gdal Translate Creating Non-Standard File from netCDF

I am trying to convert netCDF file from NASA Ocean Color website such as https://oceandata.sci.gsfc.nasa.gov/cgi/getfile/A2018001180000.L2_LAC_OC.nc now the file converts fine using the answer to ...
0
votes
0answers
15 views

Workpath for Tiled GeoTIFF Warp and Worldfile

I am attempting to identify a workpath (or script) for warping TIFFs (and TFWs) and saving them out. Have tried diff methods "batch mode" and pointed it to the relevant folder etc. But I just get an ...
0
votes
0answers
15 views

Does GDALWarp not support writing MBTiles? [duplicate]

I'm trying to create MBTiles from GeoTIFF datasets (UTM33 projection). I have tried the following command gdalwarp -of MBTiles -t_srs EPSG:3857 -co "QUALITY=80" tiffs\*.tif dataset.mbtiles The ...
0
votes
0answers
81 views

difference between gdalwarp in gdal and in R

I have got some issues using reprojection tool in R, which I already mentioned it HERE. Also read some topics like this but still have some problems understanding the different outputs from R ...
2
votes
1answer
52 views

Use gdal to set multiband raster values to no data

I need to set a 3 band multiband raster (GeoTIff) where band 1 = 0 AND band 2 = 0 AND band 3 = 0 then this is no data (0,0,0 is No Data). Not OR in each band but AND across the 3 bands. I think it is ...
2
votes
1answer
63 views

Transforming Lambert Conformal Conic to EPSG:3857

I have a map file in a tiff format that is not geo-referenced, but I know the following projection information: +proj=lcc +lon_0=10.856 +lat_0=51.322 +lat_1=42 +lat_2=57 +datum=WGS84. I would like to ...
2
votes
1answer
68 views

Gdal merging every two images

I have two sets of MODIS imagery, changed and re-projected to GEOTIFF: 1.includes the south part of my region (100 images) 2.includes the north part of my region (100 images) I want to merge them into ...
3
votes
0answers
46 views

Match Projection & Extent Using GDAL

I am looking to match both the projection and extent of two separate .tif files. The file I want all other .tif files introduced into raster calculations is as follows: The other raster image covers ...
2
votes
0answers
53 views

Reprojecting EPSG:5514 to EPSG:3857 for use in Leaflet

I have one JPG (170 MB) and JGW file of a small area projected in EPSG:5514 (Image provider told me.) JGW file: 0.1250000 0.0000000 0.0000000 ...
2
votes
2answers
53 views

How to estimate the width/height size of output image in a target CRS without gdalwarp?

I have a hard problem to find the width/height of a 2D GeoTiff file which is reprojected from a source CRS (e.g: EPSG:4326) to target CRS (e.g: EPSG:3857) without running gdalwarp. Because, the file ...
0
votes
0answers
13 views

When using gdalwarp it reorintates the map to grid north

I have a Geopdf map geo-referenced with GDA94 datum and UTM projection. The map is orientated to true north and so the grid lines are at an angle (away from the central meridian). When using gdalwarp ...
0
votes
0answers
17 views

QGIS/GDAL: Clip Raster to Vector omits cells

I am trying to clip a raster (blue) by a vector shape (green) in QGIS (which seems to use gdalwarp as a backend). The result is the black raster in the image below and as you can see there are some ...
0
votes
0answers
41 views

What transformer should I use for `gdal.Warp`?

The documentation for gdal.Warp (http://gdal.org/python/osgeo.gdal-module.html#Warp) lists the following alluring options: tps --- whether to use Thin Plate Spline GCP transformer rpc --- whether ...
1
vote
1answer
68 views

Re projecting Tiff file from the Mollweide Projection to Mercator Projection using Gdal 2.2

I'm currently looking at re-projecting a tiff file I have of population density. The file is using The Mollweide Projection. I have come to understand that Gdal would allow me to do this. I have ...
0
votes
0answers
39 views

Reproject NetCDF with 'lat' and 'lon' 2D-arrays

I have a netcdf with 'lat' and 'lon' 2D-arrays and would like to reproject into a GeoTiff. link to netcdf file: https://1drv.ms/u/s!AoFTayz0dxdxj5wPQkjj5Hv6hXC6yQ $gdalinfo goes15.2018.060.183018....
2
votes
0answers
52 views

Batch convert geotiff to jpg but ensure nodata is white

I've been using gdal_translate to convert geotiffs of rotated paper maps to jpg+wld, but nodata sections would always turn up black (i.e., where the original data has been rotated). So based on other ...
1
vote
0answers
27 views

gdalwarp -r max and -r average giving the same (incorrect) result

I am resampling a raster from ~500m to ~1km in using two different methods: max and average. However, both results are the exact same and not what I'd expect for either. gdalwarp -ot Byte -tr 926....
3
votes
0answers
64 views

Does Gdal read temporary files?

I am trying to change geographic reference system images from landsat 8. I created one method to do it in my class. I used the temporary directory as strategy to convert files and clean all row ...