An image mosaicing, reprojection and warping utility

learn more… | top users | synonyms

1
vote
1answer
135 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
17 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
61 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
96 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
37 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 ...
6
votes
2answers
313 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
1answer
54 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 ...
0
votes
0answers
40 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 ...
2
votes
1answer
48 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
118 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
564 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
30 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
11 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
288 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
43 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
29 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
35 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
32 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
96 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
81 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
15 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 ...
0
votes
0answers
24 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
26 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
40 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
20 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
63 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
80 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
83 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
63 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
146 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
33 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
22 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
84 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
48 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
105 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
78 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
27 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
25 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
396 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
256 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 ...
1
vote
1answer
90 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
284 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 ...
3
votes
2answers
1k 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 ...
1
vote
1answer
78 views

How can I access the bands of a grib file with gdalwarp?

How can I access the subdatasets of a grib file? I am working with NLDAS data and I want to process only the first band of the grib file. I know how to do it with HDF data but I couldn't handle it ...
0
votes
0answers
6 views

Gdalwarp in maptilerpro command line

-crop_to_cutline Trying to use -wo command in Map Tiler Pro inside a myarguments.txt source input text file tried it many ways with no success
0
votes
0answers
35 views

GOES 13 EAST How to Georeferencing

I have a many image for Georeferencing, Goes 13 EAST this is metadata: GOES 13 Imager frame 78 at UTC 14:45:22.485 day 009 of 2015 Vis pixels : 4009, 24840 lines : 2621, 13440 Vis Lon : ...
1
vote
0answers
24 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
25 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 ...