An image mosaicing, reprojection and warping utility

learn more… | top users | synonyms

4
votes
3answers
18 views

gdalwarp to geographic coordinate system

I'm getting the folllowing Error when trying to re-project a raster from NAD_1983_Albers to GCS_North_American_1983. I used the proj4 defenition from going to Properties in QGIS of the layer whose ...
1
vote
1answer
37 views

GDAL: gdalwarp wrong transformation to “South_Pole_Lambert_Azimuthal_Equal_Area” EPSG:102020

I am trying to reproject a raster that is defined in EPSG:4326 to EPSG:102020. Command used gdalwarp -ot Float32 -s_srs EPSG:4326 -t_srs EPSG:102020 -r near -of GTiff -co COMPRESS=DEFLATE -co ...
0
votes
0answers
25 views

gdalwarp error cropping raster

I'm getting an error from gdalwarp using the following code: gdalwarp -cutline input.shp -crop_to_cutline -dstalpha input.tif out.tif Error: *** Error in 'gdalwarp': free(): invalid pointer: ...
0
votes
1answer
40 views

Reprojecting and retiling a large set of tiles

I am trying to understand the process to re-project and re-tile a set of raster (GeoTiff) but can't wrap my head around how to use the different GDAL tools to achieve my goal. More concretely, I ...
7
votes
0answers
142 views

gdalwarp producing empty rasters when called from python script, but not via command line

I'm trying to aggregate some raster data (all .tif) up to a coarser resolution (from .05 degree to .25 degree) using gdalwarp in python, but the command is not working. Instead of getting an output ...
2
votes
1answer
40 views

Raster interpolation with gdalwarp, and replicating with python API

My problem is similar to this: Replicating result of gdalwarp using gdal Python bindings. I want to resample a raster using bilinear interpolation, which I've done successfully with gdalwarp by adding ...
1
vote
0answers
27 views

Retaining resolution of MODIS data while reprojecting using gdalwarp

Can anyone help me to use gdal(gdalwarp specifically) for reprojection? I am using it to reproject MODIS data : I originally have an hdf file - MCD43A4.A2005001.h08v05.005.2007351192333.hdf which has ...
7
votes
0answers
173 views

PyQGIS save raster as rendered image, then use GDAL tools on it

I have tried to use the code from QGis Save Raster as Rendered Image. i is a raster layer: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), ...
2
votes
1answer
62 views

Crop raster from vector problem

So I have a shapefile of a vector for Philadelphia's city lines and a raster image of Philadelphia county that I got from the lovely people from USDA:NRCS. I converted the raster from a sid file to ...
1
vote
1answer
75 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
30 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
55 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
46 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
36 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
43 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
67 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
17 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
39 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 ...
2
votes
2answers
102 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
26 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
50 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
48 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
20 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
32 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
149 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
26 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
46 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
39 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
51 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: ...
0
votes
1answer
81 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
51 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
32 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
57 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
46 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 ...
2
votes
1answer
131 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
43 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 ...
3
votes
1answer
48 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
68 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
42 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
58 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
75 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
65 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
80 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 ...
3
votes
1answer
158 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, ...
1
vote
1answer
162 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 ...
0
votes
0answers
34 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
1answer
59 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
27 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: ...
4
votes
0answers
114 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 ...