An image mosaicing, reprojection and warping utility

learn more… | top users | synonyms

1
vote
0answers
29 views

Problem defining the tiling starting point

We are tiling a image datasets. However, we want to specifically define custom windows. To do this, we are using -te option during warping. The warped output is then tiled. The problem though is that ...
1
vote
1answer
281 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 ...
1
vote
1answer
111 views

WKT for local mine grid

I have been trying to create a WKT entry to use in gdalwarp etc. I need an affine transform from coordinates in MGA94 Zone 54 to our local mine grid. All efforts have so far failed, including using ...
0
votes
3answers
296 views

Leaflet.js displaying incomplete base map

I am working on a tile server with a leaflet.js UI. I have successfully created the tiles for my base map and the tiles for the overlay. The result can be seen here (zoom in further to see overlay). ...
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 ...
1
vote
2answers
601 views

Cropping Geotiff with gdalwarp

I downloaded a geotiff from the recently published 25m-DEM of the European Union (see: http://courses.neteler.org/eu-dem-new-digital-surface-model-at-25m/). One tile has about 5GB so I would like to ...
1
vote
2answers
238 views

how to georeference my png using gdal

Hello I have created my own fantasy map, I have the map in a full size png and the map broken down into tiles ( with no geo data ) I installed MS4W and tried to use it to get the georeference data ...
1
vote
1answer
290 views

How to do geometric correction of remote sensing data in R?

My challenge is to replace most of the remote sensing commercial tools with R code (I'm teaching assistant at the University and I would like to move to open source R code to show the added values to ...
1
vote
1answer
169 views

GDAL tiling with a lot of MrSID files

I'm trying to generate tiles for a region with a surface of 85k km2. The data source are MrSID files (2734 files). I'm using gdal to achieve this task. It works fine with a small subset of MrSID, but ...
0
votes
0answers
44 views

How to increase the width and height of an image using gdalwarp?

I had a input tif file of size 936x736, pixel size <30.000,-30.000> and Block=936x10. When I increase the width and height to 1872x1572 using gdalwarp -ts 1872 1572 input.tif output.tif then ...
2
votes
2answers
281 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 ...
3
votes
2answers
494 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 ...
0
votes
0answers
127 views

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 ...
4
votes
1answer
178 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 ...
1
vote
1answer
131 views

How to get WGS84 lat lon bounds from a VRT file

I have vrt files generated using gdalbuildvrt that reference GeoTiffs referenced using the NZTM coordinate system. Is there an easy way to determine the bounds as WGS84 latitude longitude ...
2
votes
2answers
951 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 ...
0
votes
1answer
99 views

How to consider more than one no data value in utility gdalwarp of GDAL

I have to do a bilinear interpolation on an image with two no data values, using the utility gdalwarp of GDAL (http://www.gdal.org/gdalwarp.html), so these two values don't must be considered in ...
1
vote
0answers
29 views

How to reproduce the gdal utility programs in Java by GDAL Java bindings [duplicate]

Can anyone tell me how can I reproduce the gdal utility programs gdal_translate, gdalwarp and gdal2tiles in java by GDAL Java bindings? Thanks!
1
vote
1answer
119 views

How to duplicate this QGIS workflow with gdal/GRASS? (reproject and crop a raster)?

I have a number of rasters (rgb images) that I want to re-project and crop. I have established a workflow for doing this in QGIS, but I'd like to do it using gdal modules, because the output is being ...
2
votes
0answers
29 views

Workflow to assign and reproject VRT [duplicate]

I have the following workflow. generate vrt from imagelist file(with paths to tifs/grids) check if coordinate transformation is required, if so convert to target SRS Create tile using gdal_retile ...
1
vote
0answers
434 views

How to properly convert EPSG:3857 to EPSG:4326 tif?

I tried using the answer from "Transforming EPSG:3857 to EPSG:4326" to do this but gdalinfo tells me that this didn't really work properly. Driver: GTiff/GeoTIFF Files: geo-map.tif Size is 1024, 512 ...
0
votes
1answer
79 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 ...
2
votes
1answer
953 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 ...
5
votes
1answer
158 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 ...
3
votes
2answers
319 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 ...
2
votes
0answers
540 views

How do I solve Gdalwarp error 'too many points failed to transform' for remapping Geostationary to Lambert conformal?

I am trying to remap from Geostationary to Lambert conformal using gdalwarp. My input data is in netcdf, and is in geographic coordinates (degrees) and I would like to output the remapped data to ...
0
votes
1answer
147 views

gdalwarp cant find libproj.dlib

Thought I had finally got Python, gdal and QGIS to play nicely together but just tried to do a batch transform using a bash script to call gdalwarp and it claimed it couldn't read the libproj.dlib ...
1
vote
0answers
103 views

gdalwarp not clipping neatline properly

I have a PDF file (10mb) which contains a NEATLINE polygon in metadata. $ gdalinfo I5601_GeoPDF.pdf Driver: PDF/Geospatial PDF Files: I5601_GeoPDF.pdf Size is 5906, 3543 Coordinate System is: ...
1
vote
1answer
128 views

re-project from EPSG:4326 to polar projection

I have global georeferenced DEM tiff file in EPSG:4326 and want to re-project in EPSG:3411. I used this command: gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3411 global.tif polar.tif which resulted in: ...
0
votes
1answer
97 views

how to stretch GMT raster with gdal to work with google

How can I stretch a PNG raster coming from GMT (Generic Mapping Tools 4.5.6) in such a way that does align with google satellite or bing? such raster is in wgs84 and crosses two UTM zones (17 and 18 ...
0
votes
0answers
336 views

Reprojection problem with gdalwarp

I have a GeoTiff file with WGS84 datum projection and I'm trying to reproject it using gdalwarp. Using this command gdalwarp -t_srs '+proj=tpers +azi=19 +tilt=30 +lon_0=-109.5 +lat_0=23 +h=1274199.4 ...
0
votes
1answer
328 views

Use subprocess.call with gdalwrap

I am now want to use subprocess.call with gdalwrap to clip raster from shapfile (in qgis 1.8.0) I use raster→projections→wrap(reproject)to generate command_line below gdalwarp -t_srs "+proj=tmerc ...
0
votes
1answer
235 views

How to prevent Qgis clipper to clamp nodata values?

I'm trying to crop some old georeferenced cartography in order to have a continuous map (every map has its blank borders that overlaps). Usign raster->extract->clipper and setting no value to -32768 ...
2
votes
0answers
74 views

Need help tuning gdalwarp to match source cell size and alignment

I was hoping gdalwarp could merge and crop MODIS tiff images such that the original cell size and alignment is maintained. I tried the command below, but cell size is off slightly and cell alignment ...
1
vote
2answers
184 views

Problem with reprojecting raster from MGI 6 to WGS

When I reproject a raster from Mgi 6 to WGS 84 using Warp in Quantum gis it makes mistake ( about 380m). Is there any way to avoid this?
1
vote
0answers
199 views

How can I fix badly specified geographic coordinate system with units in arc-seconds?

Check out this depth-to-rock soil data. It is linked from here. The spatial reference system is a latlon, but the units are arcseconds. I could not reason out how to override the SRS even after ...
1
vote
0answers
94 views

Coarse heat map placement using gdal

I have generated a coarse heatmap of where data is available by 1 degree x 1 degree cells that I want to overlay on top of a world image. Using a python script I merged two images (heat map and a ...
4
votes
1answer
908 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 ...
3
votes
1answer
257 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
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 ...
3
votes
1answer
338 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 ...
1
vote
0answers
82 views

Size of output geotiff file not as in the .vrt file

I have the following vrt file (input.vrt) with GCPs in EPSG:3857 <VRTDataset rasterXSize="5409" rasterYSize="5408"> <Metadata /> <Metadata domain="IMAGE_STRUCTURE"> <MDI ...
3
votes
1answer
217 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 ...
2
votes
2answers
166 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 ...
4
votes
0answers
228 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
2answers
443 views

How to surround a clipped raster with transparency? (gdalwarp)

I am trying to perform the following task, but it does not seem to work: I have got a raster file (raster.tif) containing an aerial photo, which I want to clip using a Shape file (clip.shp) that ...
3
votes
1answer
634 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 ...
2
votes
1answer
283 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 ...
11
votes
2answers
843 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 ...
3
votes
1answer
495 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" ...