An image mosaicing, reprojection and warping utility

learn more… | top users | synonyms

4
votes
2answers
2k 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 ...
1
vote
1answer
144 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
140 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
31 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
536 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
89 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
1k 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
171 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
389 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
865 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
232 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
155 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
169 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
0answers
462 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
493 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 ...
1
vote
1answer
317 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
88 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 ...
2
votes
2answers
215 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
244 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
108 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
1k 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
275 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 ...
4
votes
1answer
3k 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
391 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
92 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
289 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
187 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 ...
1
vote
2answers
527 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
318 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 ...
12
votes
2answers
1k 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
581 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" ...
3
votes
1answer
456 views

How to reproject a raster to Winkel Tripel projection?

I'm trying to project a raster map (in ASCII grid format) from longlat to the Winkel Tripel projection. But I'm running into errors, for instance gdalwarp won't accept the proj-string. $ gdalwarp ...
1
vote
1answer
372 views

Georeference Raster File with WorldFile into an aerial photo using GDAL [duplicate]

Possible Duplicate: How can I transform a georeferenced JPG to a georeferenced TIF? UPDATE:Corrected link I must note that I am new to GIS so pleas be understanding. PROBLEM: I have a ...
2
votes
1answer
413 views

Technique using GDAL to overlay / fit one PNG map to directly overlay on another?

I have one .PNG basemap with known North edge of Latitude = 52.48, a South edge of Latitude = 21.94, an east edge of Longitude = -61.87 and a west edge of Longitude = -129.37. And I have another PNG ...
1
vote
2answers
207 views

How to apply v.transform on a set of vector layers (GRASS 6.4)?

I have 64 vector layers, all of them must be warped using the v.transform module; the pointsfile is the same for all vectors. How can I apply the transformation recursively to all layers? I'm not ...
1
vote
1answer
142 views

Cannot get gdalwarp to work

I am having some trouble getting gdalwarp to work. I want to clip my input raster in.img to the boundary of a polygon stored in a shapefile polygon.shp and save it out as out.img. The data type of ...
0
votes
0answers
27 views

Convert Lambert png images to Mercator projection [duplicate]

Possible Duplicate: How should I use GDAL to perform a transformation? Is there any way of converting the output .png images of a WRF system which are in Lambert projection to Mercator? I ...
0
votes
1answer
822 views

gdalwarp results in incorrect extent

I have a GeoTIFF that's 172800x86400 on the MODIS "250 meter" Sinusoidal grid (more like 231.7, but who's counting?). (You can find its metadata at the bottom of this post.) This bounty of resolution ...
1
vote
1answer
135 views

Open street map overlay white background

I have been working through trying to get the NOAA ridge radar image to act as an overlay. Looking at other entries on the site I have worked with the image and have it projecting correctly but the ...
1
vote
0answers
80 views

gdalwarp with -srcwin

I am trying to cut a window (in pixels) from a geotiff file and output in a different size. gdal_translate allows that but gives no option to set resampling, while gdalwarp allows resampling while it ...
1
vote
1answer
491 views

Gdal tiling error for EPSG: 3857

i have lots of data which was coordinated with gdal some programatic way at below: gdal_translate -of GTiff -a_ullr %s %s %s %s -a_srs EPSG:3857 %s/%s %s/%s.tif then for tiling them all i am ...
6
votes
1answer
502 views

How do I make gdalwarp set target extents to -180 -90 180 90?

I'm using gdalwarp with the -te option and it doesn't appear to be working like I expected. I'm running a command like: gdalwarp -of GTiff -te -180 -90 180 90 -tr 0.1 0.1 -t_srs EPSG:4326 input.tif ...
5
votes
2answers
1k views

Converting Ordnance Survey raster maps to WGS84 webmap tiles

I would like to generate Google/Bing-like map tiles from Ordnance Survey's 1:25000 and 1:50000 raster maps, which are provided as 4000x4000 UK national grid tiles (GeoTIFF). A sample tile can be ...
1
vote
1answer
306 views

gdalwarp - strange result

I have an input jpg image (sat thumbnail) and I know a number of gcps for this image. I want to reproject the jpg to have a north up image. This is what I do: ./gdal_translate -of GTiff -a_srs ...
3
votes
1answer
692 views

Why does gdalwarp return incorrect corner coordinates for my geotiff?

I have a png source image that's in Lambert Azimuthal Equal Area and I have the long/lat extent. So I ran gdal_translate to create tie points and convert to gtif. Then I ran gdalwarp to reproject to ...
7
votes
1answer
1k views

gdalwarp, cutline and postgis

I can not figure out how I directly can read a polygon into gdawarp's -cutline argument. My current approach is: ogr2ogr -f "ESRI Shapefile" mydata.shp PG:"user=myuser dbname=gisdb" -sql "select * ...