An image mosaicing, reprojection and warping utility

learn more… | top users | synonyms

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
17 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
63 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
21 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: ...
3
votes
1answer
84 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
59 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 ...
6
votes
3answers
329 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 ...
0
votes
0answers
19 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
55 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
42 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 ...
1
vote
1answer
153 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
19 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
1answer
64 views

Issues at 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 ...
0
votes
2answers
99 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 ...
3
votes
0answers
44 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 ...
1
vote
1answer
57 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 ...
2
votes
1answer
93 views

Properly Converting from UTM to Lat/Long with gdalwarp

I have a file in GeoTiff format in UTM coordinates. I convert the file to LatLong, then back to UTM, and the coordinates aren't what I started with. I have the printout of the process below. ...
6
votes
2answers
134 views

How to convert the geographical projection of an extremely large file (GDAL)

My Task I need to convert 9 x 14Gb GeoTiff files to WGS84 geo projection. The rough size in pixels is about 200k x 100k. The data-type is byte and it is not a subject to interpolation, it is some ...
2
votes
2answers
625 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 ...
1
vote
1answer
35 views

Transform from Lambert Azimuthal “old” to ETRS89_LAEA

I'm currently struggling to convert a small .asc dataset with no defined projection. The projection that has to be first applied is Lambert Azimuthal "old" projection: In PROJ.4 +proj=laea +lat_0=48 ...
0
votes
0answers
14 views

Rasterizing a SSURGO shape file with GDAL — Choosing right number of Pixels

I am trying to rasterize a SSURGO shape file into UTM coordinates with 5m resolution. My first step gives me a GeoTiff with 3000x3000 pixels in lat/lon. gdal_rasterize -a MUKEY -ts 3000 3000 -l ...
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 ...
2
votes
1answer
298 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 ...
1
vote
1answer
54 views

Tilemill and GeoTiff not working together

Loading in a geotiff with projection of EPSG:4326 isn't fully rendering: The original: The original in Tilemill: When using gdalwarp to change the projection from EPSG:4326 to EPSG:3857, ...
0
votes
1answer
33 views

gdal_translate or gdalwarp to define projection

Which of these is the best approach to define (aka. set) a projection for a tiff? 1) gdal_translate -a_srs EPSG:32182 in.tif out.tif 2) gdalwarp -s_srs EPSG:32182 -t_srs EPSG:32182 in.tif out.tif ...
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 - ...
0
votes
2answers
58 views

gdalwarp changing value in nodata pixels

I'm using gdalwarp to convert MODIS NDVI mosaics from Sinusoidal to epsg:4326. To save space I'm using VRT files until I reach the final product. Here is my workflow. Mosaic NDVI Modis tiles using ...
0
votes
0answers
45 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
132 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 ...
1
vote
1answer
87 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 ...
0
votes
0answers
16 views

Relating to GeoLocation tool of Geospatial PDF Document

I created a PDF file with OGC_BP as GEO_ENCODING option like shown in the below command. The below tif file’s Datum and Projection details are WGS84 /Geographic (lat/lon). ./gdal_translate.exe -of ...
1
vote
0answers
29 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: ...
0
votes
0answers
28 views

gdalwarp does not keep the original color-ramp

I am trying to clip a PNG file with the following command in the OSGeo4W Shell: gdalwarp -cutline coastline.shp -dstnodata "143 186 219" original.png clipped.tif This is how the original png looks ...
0
votes
1answer
45 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 ...
3
votes
0answers
21 views

Cropping eudem with GDAL creates zig-zag terain

I have downloaded EUDEM (20Gb) 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
65 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 ...
3
votes
1answer
91 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
0answers
93 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: ...
0
votes
1answer
73 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
164 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 ...
0
votes
2answers
40 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 ...
0
votes
1answer
23 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
2answers
92 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
52 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
1answer
116 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
90 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 = ...
0
votes
0answers
33 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 ? ...
0
votes
0answers
26 views

Downsample raster average map with control on the average algorithm

I need to downsample raster average map with control on the average algorithm to 1/4 from original size
2
votes
2answers
416 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 ...
3
votes
1answer
281 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 ...