An image mosaicing, reprojection and warping utility

learn more… | top users | synonyms

3
votes
0answers
32 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
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 ...
2
votes
1answer
45 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. ...
0
votes
0answers
39 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 ...
6
votes
2answers
117 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 ...
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 ...
1
vote
1answer
28 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 ...
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
28 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 ...
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
29 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
95 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
80 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
14 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 ...
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
39 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
62 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
78 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: ...
3
votes
1answer
79 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
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
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
1answer
47 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 ...
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
75 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
26 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
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
34 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 ...
1
vote
1answer
89 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 ...
2
votes
0answers
58 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 ...
3
votes
2answers
351 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
62 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
106 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 ...
1
vote
1answer
178 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
58 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 ...
2
votes
1answer
81 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
254 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
79 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
91 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
95 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
50 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
108 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
48 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 ...