An image mosaicing, reprojection and warping utility

learn more… | top users | synonyms

12
votes
2answers
945 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 ...
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 * ...
6
votes
1answer
477 views

How do I make gdalwarp set target extents to -180 -90 180 90?

I'm using gdalwarp with the -te option and it doesn't appear to be working like I expected. I'm running a command like: gdalwarp -of GTiff -te -180 -90 180 90 -tr 0.1 0.1 -t_srs EPSG:4326 input.tif ...
5
votes
3answers
2k views

gdalwarp cutline along with shapefile

I am trying to clip my Raster file according to a shape file geometry. I am using the following code gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif But it is resulting a black color outside the ...
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 ...
5
votes
1answer
164 views

How does QGIS reproject data?

Projection systems are always troubling and I cam across and issue yesterday. I created a buffer and did a select by location using that buffer, however it was producing the wrong results. It would ...
5
votes
1answer
244 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 ...
4
votes
1answer
2k views

What's the difference between gdalwarp and gdal_merge for mosaicing

I wanted to merge around 20 dem tiffs each around 100MB, and I used gdalwarp (I know about virtual file option): gdalwarp $(list_of_tiffs) merged.tiff However this command took very long time, and ...
4
votes
1answer
1k views

How do I solve a “No PROJ.4 translation” error with gdal?

I've installed gdal using the brew utility on OS X so the proj and geos dependencies were also installed. I'm now trying to warp a GeoTIFF so that I can import it into Google Earth. I am getting a ...
4
votes
1answer
210 views

GDALWARP and OGR2OGR / CS2CS not agreeing on transformation

I am trying to transform some data (both raster and vector) from WGS84 to a Croatian Coordinate system defined by: +proj=tmerc +lat_0=0 +lon_0=15 +k=0.9999 +x_0=5500000 +y_0=0 +ellps=bessel ...
3
votes
2answers
1k 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
3answers
1k views

Should I use gdal_retile.py or gdal2tiles.py?

I have many GeoTif files that I want to display in google earth with Geoserver. I am wondering if I want to use gdal_retile.py or gdal2tiles.py to create the image pyramids for them to be usable in ...
3
votes
2answers
605 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 ...
3
votes
2answers
37 views

Gdal: clipping a raster with another raster

I'm writing a simple utility to crop batches of multi-band geotiff raster files to the same (smaller) area. Using gdalwarp, I can easily crop a file using a single-polygon clipping shapefile: ...
3
votes
1answer
703 views

Coordinate system of a MODIS file to be introduced in GDAL for transformation

I am working on MODIS data using GDAL. When I extract coordinate system information of the files I got this: Coordinate System is: PROJCS["unnamed", GEOGCS["Unknown datum based upon the custom ...
3
votes
1answer
472 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 ...
3
votes
1answer
403 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 ...
3
votes
1answer
43 views

Gdalwarp with “out of globe” values in sinusoidal projection

I want to project a global image (source ERA INTERIM) from WGS 84 to sinusoidal, using gdal 1.11 and bilinear resampling. It work fine BUT I have values everywhere on the rectangular extent of the ...
3
votes
2answers
75 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 ...
3
votes
2answers
344 views

informations about “gdaltransform” utility of GDAL

I have a doubt about GDAL reprojection. I use "gdalwarp" utility to reproject an image in a specific projection system (for example WGS84 UTM49). "gdalinfo" on original image displays this corner ...
3
votes
1answer
534 views

How to convert a geotiff from WGS84 to 3857 with gdalwarp?

I'm trying to convert images with WGS84 projection to 3857 but warp is not giving me a output I would expect. This is my command: gdalwarp -s_srs "+proj=merc +datum=WGS84" -t_srs "EPSG:3857" ...
3
votes
1answer
76 views

Generating a GDAL script through the QGIS Georeferencer

For several months, I have been georeferencing files manually using the Georeferencer plugin within QGIS. I am now attempting to automate this process using GDAL. The Georeferencer provides a ...
3
votes
1answer
247 views

gdal_warp of multiple rasters keeping maxium intensity value

I have multiple rasters which overlap eachother geospatially but do not completely occupy the same space. I can combine them with gdalwarp like: gdalwarp -s_srs EPSG:4326 *.tif larger_output.tif ...
3
votes
1answer
263 views

Projecting GCP points into translated file

I have file A.png and file B.png. I translate file A.png so i could translate it into B.png (creating multilayer file). Is there any easy way to get new GCP points coordinates so i could combine those ...
3
votes
1answer
351 views

Adding metadata to a geographic gridded data-set and projecting it

I asked a question before that nobody answered...but I kind of figured it out myself. But I still remain with some things that I just don't seem to get right. Can anybody please find the time to help ...
3
votes
1answer
639 views

Why does gdalwarp return incorrect corner coordinates for my geotiff?

I have a png source image that's in Lambert Azimuthal Equal Area and I have the long/lat extent. So I ran gdal_translate to create tie points and convert to gtif. Then I ran gdalwarp to reproject to ...
3
votes
1answer
217 views

Bug in gdalwarp's simple nearest neighbor? Or more likely I need some help

I'm trying to understand how to use the GDAL utilities a little better. I think I've defined something incorrectly or I've found a bug in gdalwarp. My bet is that I've defined something incorrectly, ...
2
votes
3answers
132 views

Use GDAL command line to copy projections

I have two images, a georeferenced Tiff and a second Tiff without a reference. The second image is a copy of the first, but it lost its reference. Now, I want to copy the reference from the first Tiff ...
2
votes
1answer
72 views

Clip raster by shapefile in parts

Please help write .sh script command to clip raster by polygonal shapefile. Shapefile have name field with unique name of polygon feature. On the output must be rasters with name of polygons in ...
2
votes
2answers
369 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 ...
2
votes
1answer
1k 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 ...
2
votes
2answers
126 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
1answer
97 views

File size inflation normal with gdalwarp?

After using gdalwarp to project and align-to-grid (via -tap) a number of rasters I noticed that the output rasters were significantly larger than the original rasters. A fairly thorough web search ...
2
votes
1answer
306 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
1answer
91 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 ...
2
votes
1answer
81 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, ...
2
votes
2answers
68 views

Cropping with GDALWarp in QGIS

I followed this post, the first answer, to clip a USGS Quad in order to remove the borders so that I can tile them. Basically it utilized the lat/long at the corners in a .txt file to import a ...
2
votes
1answer
49 views

does Pseudo Mercator Projection causes problem while clipping rasters

I tryed to clip a Raster file using an ESRI shapefile using gdalwarp command, but i always get an empty output file or a wrong clip with wrong location. this the command bellow: gdalwarp -dstnodata ...
2
votes
1answer
460 views

Cutting image using gdalwarp to extent of data

I am currently generating a raster from point data and applying inverse distance algorithm. This image is then cut with : gdalwarp -cutline <<src shape file>> -crop_to_cutline -srcnodata ...
2
votes
2answers
172 views

Coordinate system of a raster file to be introduced for transformation

I am working on raster data using GDAL. When I extract coordinate system information of the files I got this: PROJCS["Geostationary_Satellite", GEOGCS["GCS_Unknown datum based upon the custom ...
2
votes
1answer
102 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 ...
2
votes
1answer
185 views

How to correctly resize raster GIS images to a given px width?

I work on ETOPO1, a 466MB file, and want to move to SRTM (~20GB Earth). My workflow involve several time-consuming processing steps so I strongly want to work on a minimum raster. It usually want to ...
2
votes
1answer
42 views

Invert gdalwrap cutline

I have a big raster image with medium details, but for some regions I do have far more details as a seperate raster image. Now I'd like to display both of them in my map that I render with Tilemill. ...
2
votes
1answer
402 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 ,here is one file https://echange-fichiers.inra.fr/get?k=LUfqiyO1Dl7az5bihnS. from the metadata ,the projection is cylindrical and the resolution is 25 km: ...
2
votes
1answer
384 views

Technique using GDAL to overlay / fit one PNG map to directly overlay on another?

I have one .PNG basemap with known North edge of Latitude = 52.48, a South edge of Latitude = 21.94, an east edge of Longitude = -61.87 and a west edge of Longitude = -129.37. And I have another PNG ...
2
votes
1answer
300 views

Crop collar from chart

I need advice on how to complete the following task: I have aeronautical charts covering parts of the US (example: http://aeronav.faa.gov/content/aeronav/sectional_files/Albuquerque_89.zip - Each ...
2
votes
0answers
25 views

Problem when converting to LCC projection

I am very much so a beginner.. but here goes. I have two maps that should line up with each other, but their projections are different. I've tried to use the gdalwarp utility to match them up, but ...
2
votes
1answer
44 views

GDAL Warp to clip a .tif: Works fine standalone, won't work in Python

I am trying to use the subprocess module to call gdal warp and clip a .tif using a shapefile. When I try the command in GDAL it works fine, but when I attempt to call it into a Python script it will ...
2
votes
1answer
76 views

Re-project from Polar projection to EPSG:4326

I simply want to project a polar stereographic raster (EPSG:3411) to Mercator (EPSG:3857) or WGS84 (EPSG:4326). The data I is from NOAA. $gdalwarp -s_srs EPSG:3411 -t_srs EPSG:3857 northpsg.20141027 ...
2
votes
1answer
35 views

warping and outputting coordinates of GeoTiff data: advice and error checking

So I have a massive TIF file (projected in Lambert azimuthal equal-area), and I'm trying to get lat/lon coordinates for each point of the grid. My first step is to warp the projection to the WGS84 ...