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

learn more… | top users | synonyms

1
vote
1answer
223 views

Transform an image with projection EPSG:4326 to EPSG:3857

I have an image with known boundary and projection EPSG:4326 and want to overlay it to openstreetmap using Leaflet.js. I have tried using function imageoverlay, but the image is not aligned with the ...
1
vote
0answers
36 views

getting started designing satellite projection raster

How do I get started designing a map with a satellite projection like this:http://bl.ocks.org/joyrexus/6619251. So far using d3 I've created a map of North Africa w/ topojson and tiles from Mapbox but ...
3
votes
1answer
147 views

How to ‘invert’ or ‘reverse’ a vector cutline when gdalwarp'ing a raster?

I have a raster of North America and I'd like to remove the larger lakes from it, replacing them with the same no-data value as the original raster. So I downloaded Natural Earth's 10 meter lakes ...
1
vote
1answer
13 views

Create single map image from two images with different color tables

Is it possible to merge the color tables of two images when using gdal to create a single image? Or is there another way to merge said images not using a color table?
3
votes
0answers
84 views

Downsampling for beginners using QGIS

I apologise from the outset. I need a full run through of working gdalwarp to downsample a raster to match another i.e. How to resample GeoTIFF images to the same resolution?. Using QGIS 2.12. Using '...
2
votes
1answer
68 views

gdalwarp leaves big hole at 180° longitude when warping from 4326 to custom AEqD

I'm trying to reproject the Natural Earth I 1:50m raster into an Azimuthal Equidistant projection centered on my hometown. As you may know, the source raster covers from -180° to 180° latitude, and -...
1
vote
0answers
96 views

gdalwarp UTM Zone 1 to WGS84 - fails to transform

I have a global dataset consisting of tiles in UTM projection that I want to mosaic and therefore need to reproject to a common projection. Can I solve my reprojection error with options to gdalwarp ...
1
vote
2answers
181 views

gdalwarp - merge two images, but only the no-data values

I am preparing tiff elevation files. For one particular country I have a 10 meter elevation model, but the data is present only for the specific country - so when the image is sliced into smaller ...
0
votes
0answers
21 views

how to change whole raster folder?

I used to gdalwarp to reproject raster tiff file, but I did it one by one How can I manage to do all raster in my folder (about 70 file) ? Is it possible? I am using ubuntu-64bit
0
votes
0answers
80 views

How to transform any GTiff to Lambert Azimuthal Equal Area with GDAL?

I am trying to transform a GTiff EPSG:4326 +proj=longlat +datum=WGS84 +no_defs and a GTiff EPSG:31982 +proj=utm +zone=22 +south +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs both to EPSG:...
2
votes
2answers
150 views

projection conversion from local to UTM

I have some images at a local coordinate system with the following specifications: TIFF header file: -gdalinfo im.tif Driver: GTiff/GeoTIFF Coordinate System is: PROJCS["Eastern Ridge Project Grid ...
1
vote
1answer
38 views

resample etopo.grd with gdal_warp gives many null values

I'm having a problem with resampling ETOPO NetCDF GRD to a bigger pixelsize for using it with Generic Mapping Tools. I'm using gdal_warp for the resampling (just for testing with pixelsize of 1°): ...
3
votes
1answer
104 views

Mass reprojection of raster files in a loop [Linux]

I have a folder with numerous .asc raster files that each accompanied by a corresponding .prj file. The file basename is the same for both files. The projection however can vary between all of them (...
1
vote
0answers
98 views

GDAL - replacing rgb value (0,0,0) with (nodata,nodata,nodata)

I have a set of rasters in which I want to replace any rgb value of (0,0,0) with (nodata,nodata,nodata) I've tried using gdal_translate -a_nodata "0,0,0" in.tif out.tif But this seems to convert ...
1
vote
0answers
26 views

Transform stereographic projection using GDAL

How can I read/wrap a JPG image without any embedded projection/coordinate system to an autoreferenced GEOTIFF with Mercator or similar coordinate system using GDAL? The image is 500x500 px and I ...
1
vote
2answers
42 views

gdalwarp change spatial resolution

I just realized that I have Geo-TIFF's in the samme coordinatesystem but different resolution. DOP20 (20cm=1px) and DOP10 (10cm=1px) mix. Now wanted to use gdalwarp to transform the DOP20 TIFF into ...
2
votes
2answers
242 views

How to re-project a netcdf file?

I want to re-project and resample a netcf file but I am getting a usage problem: gdalwarp -t_srs EPSG:3410 -ts 1383 586 NETCFD:"C:\data\product.nc":gsfc output.tiff
1
vote
0answers
35 views

Changes in MOD09A1 quality mask as processed with GDAL

I'm analyzing MOD09A1 spectral data using GDAL/Python. The quality mask is a uint32 dataset and want to pull out the value of the first bit to mask the spectral data (1 == poor data; 0 == good data). ...
1
vote
1answer
61 views

gdal_warp some geotiffs change SRS

Hello I#m quite new to GDAL. I have some Geotiffs that I want reproject. I tried that with one file and it works great, but because the TIFF has to be slightly rotated for the conversion I have some ...
1
vote
0answers
109 views

how to reproject an image programatically in c# using gdal

I have to reproject an image programmatically in c# using gdal library. For the moment i'm a little bit lost in gdal and i don't really now which methods i have to use. I think i will have to ...
2
votes
1answer
93 views

Reproject “unknown” coordinate system to LCC 2SP across lon 180

Starting with a grib file with limited projection information, I am trying to import it as a GeoTiff in Geoserver 2.7.x. The original gdalinfo is: Size is 218, 81 Coordinate System is: GEOGCS["...
0
votes
1answer
140 views

Use lamber 93 imageOverlay with leaflet on top of OSM tiles layer

I'm trying to use a lambert 93 image overlay on top of an Open Street Map title layer. var map = L.map('map'); L.tileLayer('http://{s}.tile.osm.org/{z}/{x}/{y}.png', { }).addTo(map); var coords = { ...
0
votes
0answers
93 views

Define band specific metadata of a GeoTIFF file using GDAL (gdalwarp/gdal_translate)

Does anybody know how to change or set the "Description" option/tag (band specific metadata) of a GeoTIFF file using GDAL? To specify what I mean, this is an example of gdalinfo return from a GeoTIFF ...
1
vote
1answer
42 views

How to re-project only one band from a hdf?

I have this file which contains: gdalinfo product.hdf Coordinate System is `' Metadata: calendar=standard center=gsfc Subdatasets: Corner Coordinates: Upper Left ( 0.0, 0.0) Lower Left ( ...
1
vote
1answer
150 views

Clipping a coloured raster with polygon shp in gdalwarp changes original values, why?

I need to clip my raster based on a polygonal shapefile, basically leave raster values inside the polygon shape (inside red perimeter in the figure below) as shown here: Therefore, I use: gdalwarp ...
1
vote
0answers
66 views

Why does gdalwarp's cubicspline method expand into nodata areas?

The image below shows part of a peninsula on the Australian island of Tasmania. The black pixels are 9" climate data with GDA94 (EPSG:4283) SRS, projected by QGIS on the fly to Australian Albers (EPSG:...
2
votes
1answer
214 views

Mosaic Imagery using GDAL Edited

I have a couple of large sets of imagery I need to mosaic using GDAL. The sets are currently in tiff format with 18 tiles totaling about 20gb in size. I have attempted to mosaic using both the ...
0
votes
1answer
54 views

Raster extents incorrect for vrt when changing projection with -a_srs flag

I want to generate a .vrt in epsg:3857 from a series of raster tiles in epsg:4326, and am trying to do so in one step using gdalbuildvrt with the -a_srs flag: gdalbuildvrt -a_srs epsg:3857 ortho_3857....
3
votes
1answer
54 views

Wedge of missing data reprojecting to stereographic

I am trying to reproject an equirectangular projection image into stereographic around the north pole of mars using gdalwarp and I am not sure what is happening in my output image. The input image ...
1
vote
2answers
102 views

Use gdalwarp to crop area with simple coordinate system

I'm new to GIS and trying to cut image with gdalwarp I need image projected like google map and I need to input coordinates of area of interest in latidute longitute system, like: from 45.3 10.2 to ...
0
votes
2answers
45 views

How do I run gdalwarp on a number of csv of coordinates?

I have a CSV of Ground Control Points in NAD83 UTM 12N and I need to convert them to WSG84. Is it possible to do this with gdalwarp or something similar?
0
votes
0answers
73 views

GDALwarp: Re-projecting from EPSG<4326> to EPSG<3857>; problem: “double pixels”

I would like to change the projection from a png file that shows river pixels that exceed a certain threshold. Therefore I am following this workflow: (note: it is part of a larger R script) 1) ...
3
votes
0answers
155 views

How do I specify an outsize for gdal translations to an ENVI file?

I'm trying to create an elevation map using the output from gdal_translate. When I run gdal I can convert to both PNG and ENVI and combine the results to make something like the following image. ...
1
vote
1answer
111 views

Using Gdalwarp to Change Resolution of File in Degrees to Metres

I am using the command: gdalwarp -t_srs WGS84 -te 144 -36 150 -32 -crop_to_cutline AustTopo250_SI55.tif I need to change the pixel resolution to 50m x 50m but the option is given as -tr xres yres: ...
1
vote
1answer
86 views

Strange effect in gdalwarp output

I am using gdalwarp in python to extract and re-project reflectance bands from MODIS level 1B (MOD021km). An example of the command that I am using is the following: gdalwarp -of GTIFF -tps -t_srs ...
4
votes
1answer
334 views

Difference between gdalwarp and projectRaster

I am trying to project a Raster. In R there is the projectRaster() function to to this (below a fully reproducibly example) : # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40,...
1
vote
1answer
241 views

Unable to warp HDF5 files

I have numerous HDF5 files(.nc) which I need to batch process warp using one of the gdal utilities gdalwarp. When I tried to warp the files an error occurred: INPUT: gdalwarp -geoloc -te 109.975 3....
0
votes
0answers
46 views

Add hillshade to naip imagery

I've been searching for quite a while now. Is there a way to merge a hillshade image with an aerial image with gdal? I'm just trying to create a hillshaded aerial image.
1
vote
2answers
77 views

Black image using gdalwarp

Im a geomatics student; Currently, I try to geo-referencing a tif image with the QGIS GUI, and then I load and save the Gdal script. With the QGIS GUI, I managed to get a decent picture, then I ...
3
votes
1answer
35 views

gdalwarp reprojection causes missing data

I have tried to reproject and resample an SRTM DEM file from WGS84 to a cylindrical equal area projection. The problem is there is missing data as seen in the image below. The command I use is: ...
5
votes
0answers
204 views

can I preserve a fading alpha layer when mosaicing images with gdal?

I am mosaicing some images with gdal and would like to improve the final result by using a fading / gradual alpha layer towards the edge of each image to remove the sharp edges in the middle of the ...
1
vote
1answer
142 views

How to change the projection and resolution of a file using gdal?

I have a file with these information: lat_min -60.0 lat_max 80.0 lon_min -180.0 lon_max 179.9 lon_sampling 0.1 degree lat_sampling 0.1 degree samples = 3600 lines = 1401 I want to re-project the ...
2
votes
1answer
1k views

Properly Converting from UTM to Lat/Long with gdalwarp

I have a file in GeoTiff format in UTM coordinates. I convert the file to LatLong, then back to UTM, and the coordinates aren't what I started with. I have the printout of the process below. ...
0
votes
0answers
128 views

Create GeoTIFF from JPG and JFW files with gdal

I have a JPG image and the corresponding geospatial information on a JFW wolrd file in a ETRS89 geodesical system and UTM coordinates (EPSG:25830) 5.0000000000 0.0000000000 0.0000000000 -5.0000000000 ...
6
votes
2answers
242 views

How to convert the geographical projection of an extremely large file (GDAL)

My Task I need to convert 9 x 14Gb GeoTiff files to WGS84 geo projection. The rough size in pixels is about 200k x 100k. The data-type is byte and it is not a subject to interpolation, it is some ...
1
vote
1answer
112 views

Transform from Lambert Azimuthal “old” to ETRS89_LAEA

I'm currently struggling to convert a small .asc dataset with no defined projection. The projection that has to be first applied is Lambert Azimuthal "old" projection: In PROJ.4 +proj=laea +lat_0=48 +...
1
vote
1answer
106 views

Tilemill and GeoTiff not working together

Loading in a geotiff with projection of EPSG:4326 isn't fully rendering: The original: The original in Tilemill: When using gdalwarp to change the projection from EPSG:4326 to EPSG:3857, ...
0
votes
1answer
135 views

gdal_translate or gdalwarp to define projection

Which of these is the best approach to define (aka. set) a projection for a tiff? 1) gdal_translate -a_srs EPSG:32182 in.tif out.tif 2) gdalwarp -s_srs EPSG:32182 -t_srs EPSG:32182 in.tif out.tif ...
1
vote
2answers
516 views

gdalwarp changing value in nodata pixels

I'm using gdalwarp to convert MODIS NDVI mosaics from Sinusoidal to epsg:4326. To save space I'm using VRT files until I reach the final product. Here is my workflow. Mosaic NDVI Modis tiles using ...
0
votes
0answers
134 views

Cutting a shape from a VRT file with gdal warp

I have a couple of overlapping tif files which I want first to merge and then to clip a shape file from the result tif. The problem is that the tifs are pretty big and when the tiffs are overlapping ...