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

learn more… | top users | synonyms

4
votes
0answers
36 views

Cropping EUDEM with GDAL creates zig-zag terrain

I have downloaded EUDEM (20 Gb) from http://www.eea.europa.eu/data-and-maps/data/eu-dem#tab-original-data When I try to crop it by GDAL using: gdalwarp -dstalpha -of EHdr -crop_to_cutline -cutline ...
1
vote
0answers
38 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 ...
3
votes
1answer
81 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 ...
3
votes
1answer
96 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 (...
6
votes
3answers
3k views

how to for loop a folder to batch clip rasters by polygon using python and QGIS

I am using python and QGIS 2.0. I am trying to clip rasters in a folder by one polygon feature. It's the first time for me using (let's say) "PyQGIS", I was used to arcpy before. Anyway, I don't get ...
1
vote
1answer
206 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 ...
3
votes
1answer
140 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
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 ...
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
75 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
64 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
94 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
157 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
73 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:...
0
votes
1answer
93 views

Reprojecting coordinates using gdalwarp and OpenLayers?

I have an interpolated raster image with EPSG:27700 (British National Grid), after I warp it using gdalwarp into EPSG:3857 (WGS84 Web Mercator (Auxiliary Sphere)): gdalwarp -s_srs EPSG:27700 -t_srs ...
2
votes
2answers
144 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
36 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°): ...
6
votes
3answers
466 views

How to stop gdalwarp creating world-spanning outputs near the dateline?

I'm using gdalwarp to manipulate SRTM tiles near the the dateline (i.e. 180°, aka the antimeridian). SRTM tiles have a very slight (1/2 pixel) overlap with the meridian. You can see this using ...
1
vote
0answers
86 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
25 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
40 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
232 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
0
votes
1answer
134 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 = { ...
1
vote
1answer
58 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
33 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
0answers
105 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
90 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["...
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
133 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
63 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
482 views

GDAL tiling with a lot of MrSID files

I'm trying to generate tiles for a region with a surface of 85k km2. The data source are MrSID files (2734 files). I'm using gdal to achieve this task. It works fine with a small subset of MrSID, but ...
2
votes
2answers
3k views

GDAL2Tiles: MapTiles from BSB/KAP are Switched

As an example I have a BSB file from NOAA which is a navigational map. The map has a SRS of WGS84 and a Mercator projection. I initially translate the file from BSB format to GTiff format using the ...
1
vote
2answers
100 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 ...
2
votes
1answer
199 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 ...
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 ...
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) ...
2
votes
2answers
1k views

Different NoData values when processing a GeoTIFF in QGIS and bringing it to ArcGIS

I'm trying to process a GeoTIFF in QGIS 1.8 using SEXTANTE and GRASS via the r.neighbors tool. I'd like to process the image, then be able to utlize it in ArcGIS. The problem I'm running into is ...
3
votes
0answers
149 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
85 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 ...
1
vote
1answer
103 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: ...
4
votes
1answer
299 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
226 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
45 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.
0
votes
0answers
126 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 ...
1
vote
1answer
404 views

How do I set nodata with custom color for gdal_translate and gdal2tilespy?

I start with a bunch of 1-bit geotiff files that I had to convert to 8-bit RGB tiff files to use with GDAL. Then I tried to set borders to red (content is in black and white so I can't use those) and ...
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: ...
0
votes
2answers
114 views

GDAL clips raster at wrong position

I try to clip raster data with GDAL (in fact I am using QGis, which produces me a command line for GDAL) The raster has been transfored to EPSG:3857, which is also the projection of the shape file ...
1
vote
1answer
132 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 ...