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

learn more… | top users | synonyms

0
votes
0answers
135 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 ...
1
vote
1answer
1k views

jpg to geotiff transformation with gdal_translate and gcps fails

I am trying to create a geotiff from a jpeg by inputting gcps with gdal_translate: gdal_translate -of GTiff -a_srs EPSG:4326 -gcp [1616.0 0 -121.459869532 38.5822533549] -gcp [0 0 -121.460081357 38....
1
vote
1answer
136 views

Issues with getting the height from SRTM gdalwarp WGS84(EPSG:4326) to Mercator(EPSG:3857)

I am trying to use the SRTM .hgt file to obtain the height of a location. My issue is that I wish to reproject the SRTM file from WGS84 to Mercator but the height returned is not consistent after ...
1
vote
0answers
50 views

What is the kernel width of the gdalwarp resampling methods?

gdalwarp provides 7 resampling methods. How many neighborhood pixels do contribute to upsampling respectively downsampling? I will try to guess. Please correct, if values are wrong: Upsampling: ...
4
votes
0answers
37 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 ...
0
votes
1answer
125 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
1answer
108 views

polygonize giving duplicated polygons

I am trying to reproject some ascii files and convert them into shapefiles. I am new to python but looking at some code online and help files I was able to write the following code to run on Python ...
0
votes
0answers
320 views

Reprojecting from Lambert Conformal Conic into EPSG:3857?

I'm trying to reproject a plain png image file into EPSG:3857 for use in an OpenLayers application. I'm pulling the image from the NWS site at: http://rapidrefresh.noaa.gov/HRRR/for_web/hrrr_ncep_jet/...
7
votes
1answer
2k views

Replicating result of gdalwarp using gdal Python bindings

I am trying to re-project/resample with the GDAL python bindings, but am getting slightly different results compared to those from the command line utility gdalwarp. See update below for shorter ...
0
votes
1answer
209 views

Crop using gdalwarp aspect ratio

I'm using Landsat 8 imagery, and the Indicar Landsat Processing Tools to produce NDVI imagery from the Landsat 8 source. That works fine. I then have a bounding box expressed as lat/lon coordinates, ...
3
votes
1answer
336 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 C:\...
0
votes
2answers
89 views

gdalwarp: Syntax

In GDAL, I have seen the sample line: gdalwarp $(list_of_tiffs) merged.tiff I would like to merge images to create a mosaic, but can't use gdal_merge (need an exe file), e.g to merge C:\Temp\image1....
0
votes
1answer
48 views

Query shapefile for crop_to_cutline using gdalwarp

I have a shapefile with several polygons. I would like to clip a series of rasters but only use a single feature from the shapefile feature class. I am using gdal_warp. There is an option to query the ...
1
vote
1answer
36 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 - 'http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO1/data/ice_surface/grid_registered/georeferenced_tiff/...
1
vote
1answer
110 views

MODIS resolution missmatch between MRT and GDAL processing

I've been using MRT (Modis Reprojectin tool) to extract EVI bands into a geotiff with geographic coordinates (epsg4326). I did not specify an output resolution to MRT. The output from MRT had a ...
2
votes
2answers
232 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
242 views

convert GTiff to AAIGrid results in cell size switch

I don't have a clue why the conversion of a tif to an asc has that kind of behavior: original tif metadata Size is 466, 427 Coordinate System is: PROJCS["ETRS89_UTM_zone_32N", ... Origin = (-1181164....
0
votes
0answers
110 views

gdalwarp for warping a Space View Perspective or Orthographic tiff

i need to warp an image with a Space View Perspective or Orthographic project. What is the parameter for the gdalwarp -s_srs "+proj=???? +ellps=WGS84 +datum=WGS84" -t_srs EPSG:4258 in.tiff out.tiff ? ...
1
vote
0answers
38 views

GDAL Warp cuttline not working

Has anyone ever got the cutline to work in C++? I have it working off the command line no problem, but it doesn't seem to do anything when I populate that argument in C++ as below: *//Create a ...
0
votes
0answers
46 views

Error in warping the merged hdf4 files

I've tried to merge a set of hdf4 files(unknown projection) and it run well. But when I tried to warp the output of the merge hdf4 files an error like this appears: ERROR 1: Unable to compute a ...
1
vote
1answer
255 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 don'...
2
votes
0answers
133 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 ...
4
votes
2answers
2k 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: ...
0
votes
0answers
148 views

gdalwarp hdf problem

I am working with some hdf files. The following command works on my Mac. command = gdalwarp -of GTIFF -tr 1000 1000 -tps -s_srs "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs " -t_srs "+...
1
vote
0answers
339 views

gdalwarp with NETCDF files

I have been looking for it for a while but couldn't find what exactly I need. I have WRF outputs in NETCDF format and I want to convert them to Geotiffs with gdalwarp. I think my -s_srs should be ...
2
votes
1answer
440 views

Creating a mosaic with gdalwarp - How to deal with no data and overlap

I want to create a mosaic of satellite imagery in .tiff format. I want to create an actual mosaic, not a vrt, because afterwards I will cut the image in to tiles that make sense for my users. Gdal ...
1
vote
0answers
81 views

hillshade and slopeshade geotiffs fail to convert to new projection using gdalwarp

I'm trying to convert my three DEM products (colour, hillshade and slopeshade; the latter two generated using gdaldem) from EPSG:4326 -> EPSG:3857 using gdalwarp, via the command: gdalwarp -of GTiff -...
3
votes
1answer
141 views

gdalwarp 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 ...
3
votes
1answer
762 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 ...
3
votes
1answer
151 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 ...
1
vote
0answers
182 views

Reprojecting and Clipping in one script using Python Subprocess module and GDAL

I am trying to reproject a .tif file then clip it to the boundaries of a shapefile using GDAL Warp. I want to do this in a single script, and have successfully used the Python subprocess module to ...
2
votes
1answer
162 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 ...
1
vote
1answer
104 views

Transforming Robinson to Equirectangular

I need convert this http://bioval.jrc.ec.europa.eu/products/gam/images/large/access-map.png (is it robinson?) to Equirectangular (ESPG:32662) How can i get this using gdalwarp?
3
votes
1answer
266 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 "...
2
votes
1answer
83 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 ...
3
votes
1answer
2k 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
2answers
247 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 ...
2
votes
1answer
249 views

gdalwarp '-tps' equilvalent in C# wrapper?

I am attempting to do the following from code: gdal_translate.exe --config GDAL_DATA "D:\gdal\bin\gdal-data" -of GTiff -a_srs EPSG:4326 -gcp 0 0 -108.724530873811 40.1376029956585 ..<truncated>....
1
vote
0answers
107 views

How to run cmd prompt in the background and show progress bar

I am warping files using gdal. Since I have numerous files to warp how will I run the cmd in background and just have a progress bar to indicate the percentage in warping the files? Because when I ...
1
vote
2answers
436 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
118 views

Geographic -> UTM Raster Rectification

I'm trying to use an image overlay in KML / Google Earth as the basis for adding extra info and returning this to Google Earth. In my utility app I can specify a rectangle to give me the north-south-...
1
vote
0answers
64 views

GDALwarp Creates Darker target Image sometimes

I'm yoyo and new to this forum. I'm relatively new to GIS and GDAL. Please accept some typos now and then... English is not my native language. Because of some problems with GDAL 1.7.0b2 in FWTools 2....
0
votes
0answers
54 views

GDAL_Translate not compressing when NoData value is set

I am trying to compress a .TIF file using GDAL_Translate or GDAL_Warp. My Raster is 32mb to start with. When I use GDAL_Translate with the COMPRESS=PACKBITS option, it compresses to 1.9mb. However ...
1
vote
1answer
253 views

How can I clip multiple geotiff tiles with a shapefile without greatly increasing their size?

I'm using gdalwarp to clip a number of geotiffs (tiles of a larger map) to a polygon described by a shapefile (which is about the size of the overall map, i.e. much larger than individual tiles) using:...
2
votes
1answer
320 views

Python and gdalwarp: change resolution of raster

I am new to GDAL and have been having some hard time figuring out how to change a raster's resolution and how to align it to another existing raster. I have come to the point where I have a raster ...
1
vote
1answer
307 views

Using gdalwarp with multiple elevation models

I'm trying to use gdal_warp to slice and convert some SRTM data. The data I'm interested in are in S30E152.hgt and S30E153.hgt. What I'd like to do is: Convert these DEM files from EPSG:4326 to EPSG:...
3
votes
1answer
256 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 ...
1
vote
0answers
119 views

Mapping from one plane on the other plane despite of masking regions

have a set of data given here, where in the first and second columns there are the sky coordinates (ra,dec) e.g. right ascension and declination, respectively and in the third and forth, the ...
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
1answer
2k 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 ...