An image mosaicing, reprojection and warping utility

learn more… | top users | synonyms

2
votes
2answers
2k 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 ...
13
votes
2answers
2k views

How to reproject raster from 0 360 to -180 180 with cutting 180 meridian

I have a geotiff raster image that has a coordinate system with longtitudes from 0 to 360. The horizontal center of the image is 180 longtitude. See image below: I want to transform it to EPSG:4326 ...
2
votes
2answers
679 views

How to re-project the EASE (Equal Area Scalable Earth) grid with a ~25 km cylindrical projection to WGS84 0.25 degree?

I have nc files. from the metadata ,the projection is cylindrical and the resolution is 25 km: indrical NC_GLOBAL#ease_resolution=25 ...
4
votes
2answers
2k 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 ...
3
votes
1answer
1k views

Re-project raster image from Mercator to Equirectangular

I have a raster image of the world map in Mercator projection, square 1:1 format. I want transform the image into an Equirectangular projection accurately. I assumed from my research that GDAL would ...
2
votes
1answer
240 views

gdalwarp won't export AAIGrid

I'm trying to use gdalwarp to reproject an ascii grid file. According to http://www.gdal.org/formats_list.html this shouldn't be a problem and indeed I can happily use gdal_translate to create grd ...
2
votes
1answer
312 views

How to remove strange gdaldem hillshade artifacts?

I've downloaded 9 chunks of 1 arc-second NED in IMG format from the National Map. I'm trying to create a hillshade from them. I'm using: GDAL 1.10.0, released 2013/04/24 First, I combine them and ...
2
votes
2answers
438 views

Cropping and transforming geoPDF with GDALwarp

I got some great advice on how to transform a geoPDF to a geoTIF on this post, which used the available downloads from this site. I would like to take it a little further and not only convert, but ...
2
votes
2answers
233 views

Problem with reprojecting raster from MGI 6 to WGS

When I reproject a raster from Mgi 6 to WGS 84 using Warp in Quantum gis it makes mistake ( about 380m). Is there any way to avoid this?
3
votes
1answer
190 views

Reprojecting a raster image using gdalwarp

I am trying to reproject an ECW Orthophoto from 25830 to 23030 srid. I tried the following but I am getting an error: C:\Program Files (x86)\FWTools2.4.7>gdalwarp -s_srs 25830 -t_srs 23030 ...
3
votes
2answers
157 views

gdalwarp leaves horizontal artifacts regridding from EASE-Grid (laea) to Polarstero (stere)

I'm trying to understand how to use gdalwarp to warp an image. I believe I'm doing the basics correctly, but I might be missing some gdalwarp options? Any help would be greatly appreciated. The ...
7
votes
1answer
1k views

gdalwarp, cutline and postgis

I can not figure out how I directly can read a polygon into gdawarp's -cutline argument. My current approach is: ogr2ogr -f "ESRI Shapefile" mydata.shp PG:"user=myuser dbname=gisdb" -sql "select * ...
5
votes
2answers
1k views

Converting Ordnance Survey raster maps to WGS84 webmap tiles

I would like to generate Google/Bing-like map tiles from Ordnance Survey's 1:25000 and 1:50000 raster maps, which are provided as 4000x4000 UK national grid tiles (GeoTIFF). A sample tile can be ...
3
votes
0answers
55 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 ...
3
votes
1answer
605 views

Georeferencing a raster using GDAL and Python

I want to georeference a raster using python and GDAL. My current approach is to call gdal_translate and gdalwarp using os.system and an ugly list of ground control points. I'd really like a way to do ...
3
votes
1answer
170 views

Holes between imagery tiles projected using gdalwarp

I am having the following issue when reprojecting imagery tiles using gdalwarp. These are 4096x4096 imagery tiles, which I'm gdalwarp'ing from WGS84 to a State Plane projection. When reprojected, ...
3
votes
2answers
1k views

Help making geotiff and reprojecting to EPSG:4326

My goal is to download a gif file from this webpage: http://aviationweather.gov/data/iffdp/2104.gif Then crop it to the known latlon bounds of 75N,45S,35W,70E, convert to geotif, and then project to ...
0
votes
1answer
614 views

Use subprocess.call with gdalwrap

I am now want to use subprocess.call with gdalwrap to clip raster from shapfile (in qgis 1.8.0) I use raster→projections→wrap(reproject)to generate command_line below gdalwarp -t_srs "+proj=tmerc ...
3
votes
1answer
186 views

Reprojecting MODIS Ice Sea Temperature products with gdalwarp

I have problems projecting MODIS Ice Sea Temperature (IST) products to EPSG:3031 (Antartic Polar Stereographic) using gdalwarp. The product I'm working with is MOD29E1D, documented here, p33: The ...
3
votes
1answer
512 views

How to reproject a raster to Winkel Tripel projection?

I'm trying to project a raster map (in ASCII grid format) from longlat to the Winkel Tripel projection. But I'm running into errors, for instance gdalwarp won't accept the proj-string. $ gdalwarp ...
2
votes
2answers
216 views

Reprojecting Natural Earth raster into natearth projection with gdalwarp causing black squares?

I'm trying to reproject the Natural Earth 2 raster layer (http://www.naturalearthdata.com/downloads/10m-raster-data/10m-natural-earth-2/) into the natural earth projection ...
2
votes
2answers
779 views

How do I know the vertical/horizontal units of my DEM? I'm creating a slope using gdaldem

I'm using gdaldem to create a slope map. I can apply a scale. The default is 1. The notes say: Ratio of vertical units to horizontal. If horizontal unit of source DEM is degrees (e.g WGS84 ...
1
vote
1answer
73 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 ...
1
vote
1answer
32 views

gdalwarp `-te -5.8 41 10 51.5 -ts 1980 0` not working with mercator?

I was working with native lat/log projection on etopo1: #donwload: curl -L -C - ...
1
vote
2answers
100 views

GDAL - Warping raster to new domain/resolution/projection causes data loss

I am attempting to warp a GTiff from a lat/lon projection to a higher (~2x) resolution Polar Stereographic grid. In the output GTiff any cells north of 64 deg lat become NaN instead of their actual ...
1
vote
1answer
108 views

Warping images with unknown projection

I have tried to used gdalwarp to project set of images(from hdf4 files to tif) before merging them using gdal_merge. But then the merge output has spaces/blanks. They're not merge properly, and I ...
1
vote
2answers
314 views

gdalwarp batch process

I have downloaded numerous hdf files from an ftp site. Half of them were level 1 and the rest were level 2. Lvl1 format is Afile.A1_AC.sample.hdf and the lvl 2 format is: Afile.A1_AC.sample.hdf. How ...
1
vote
1answer
244 views

WKT for local mine grid

I have been trying to create a WKT entry to use in gdalwarp etc. I need an affine transform from coordinates in MGA94 Zone 54 to our local mine grid. All efforts have so far failed, including using ...
1
vote
0answers
186 views

gdalwarp not clipping neatline properly

I have a PDF file (10mb) which contains a NEATLINE polygon in metadata. $ gdalinfo I5601_GeoPDF.pdf Driver: PDF/Geospatial PDF Files: I5601_GeoPDF.pdf Size is 5906, 3543 Coordinate System is: ...
1
vote
1answer
198 views

re-project from EPSG:4326 to polar projection

I have global georeferenced DEM tiff file in EPSG:4326 and want to re-project in EPSG:3411. I used this command: gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3411 global.tif polar.tif which resulted in: ...
0
votes
1answer
49 views

GDAL - Specify span of elevation values for conversion to Terragen format?

I'm attempting to use the gdal tools to convert some geotiff files I have to the Terragen format. I was hoping to use either gdal_translate or gdalwarp in a similar fashion to this post. (though I ...
0
votes
0answers
538 views

Reprojection problem with gdalwarp

I have a GeoTiff file with WGS84 datum projection and I'm trying to reproject it using gdalwarp. Using this command gdalwarp -t_srs '+proj=tpers +azi=19 +tilt=30 +lon_0=-109.5 +lat_0=23 +h=1274199.4 ...