GDAL (Geospatial Data Abstraction Library) is an open source translator library & set of command-line utilities for processing raster spatial data formats. As a library, it presents a single abstract data model to the calling application for all supported formats. It also comes with a variety of ...

learn more… | top users | synonyms

7
votes
3answers
2k views

How to extract slope from topographic raster map?

I am working on a project for a machine learning class wherein I extract features from a topographic raster map and then classify the traversability of the terrain. One of the features I need to ...
1
vote
0answers
118 views

VSIMalloc3() and UnifiedSrcDensity mask errors when running gdal_warp in OSGEO4w shell

This is related to gdal.py as part of bat script in windows 7 I am getting ERROR 2: VSIMalloc3(): Out of memory allocating 854522136 bytes. ERROR 2: Out of memory allocating 119813400 bytes ...
4
votes
1answer
150 views

How to echo gdal process to screen and pipe to text file in Windows?

gdal_warp and similar give a nice progress summary when run and I want this output to show in the osgeo shell and also pipe to a text file. The issue is when I add the pipe the output doesn't show on ...
2
votes
1answer
272 views

gdal.py as part of bat script in windows 7

I need to run the following bat file from a MS4W shell which has gdal_retile.py as the final option. I have added C:\ms4w\python\gdal to my pythonpath but within a pyscripter shell I can not import ...
3
votes
2answers
763 views

set geotiff 0 values to nodata

Is there a way to set all the 0 value pixels on a raster to nodata using gdal or qgis?I world like to use fillnodata to fill some holes, but holes are set to 0,so I need to delete them first
2
votes
0answers
176 views

Localized Linear Histogram Match using python

I would like to use the local linear histogram match for filling the missing pixels of a Landsat 7 image with another older image (SLC Gap-Filled Products). Do you have any idea how can I do it by ...
4
votes
1answer
287 views

How can I select which interpolation technique will be used in a gdal command? (gdal_fillnodata.py)

I work on Landsat 7 slc-off images and I would like to use the gdal_fillnodata.py command to fill-in the missing pixels. The problem that I have is that the gdal_fillnodata.py uses the "inverse ...
2
votes
2answers
255 views

Why all WFS and WPS are CGI based?

Even though we have a lot of web technologies now, why are all WFS and WPS still cgi based? PHP allows tp read from postgres or other db's and HTML5Canvas or SVG can render the maps on the web. ...
2
votes
2answers
537 views

error message when clipping raster with vector layer

When I try to crop a raster using QGIS Clipper tool, I get the following error message: ERROR 1: Attempt to create 0x0 dataset is illegal,sizes must be larger than zero.. Why? The code QGIS generated ...
3
votes
1answer
857 views

GeoTiff convert to PNG with transparency

I need to convert a GeoTIFF image to PNG format with "NoData" pixels in the original image being transparent in the output. Is there any tool particularly in ArcGIS ArcToolbox or GDAL for this ...
5
votes
2answers
285 views

How to modify raster pixel inside a polygon?

I would like to set to novalue all the pixels of my raster inside a polygon. I'm trying the gdalwarp command line with -cutline option. But this command extract the values of my raster inside the ...
3
votes
1answer
217 views

Bug in gdalwarp's simple nearest neighbor? Or more likely I need some help

I'm trying to understand how to use the GDAL utilities a little better. I think I've defined something incorrectly or I've found a bug in gdalwarp. My bet is that I've defined something incorrectly, ...
3
votes
1answer
248 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 ...
11
votes
1answer
288 views

How to use GDAL with PyPy?

I am trying to optimize my python code using pypy. While I am able to import most libraries and run succesfully, I am not able to import GDAL libraries. The same libraries when imported in Python 2.7 ...
6
votes
2answers
470 views

Read an ArcGIS Feature dataset using GDAL

I am new to GDAL in Python and would like to use it to create a list of populated layers in a ArcGIS mdb Feature datasets. I have the code below that provides a list of all the populated layers but ...
2
votes
0answers
258 views

Image Pyramid + Geoserver + Google Earth slow with errors

I have hundreds of GB of somewhat small tif files. Here is gdalinfo of one of them: F:\gis\geotiffs>gdalinfo finaloutputSSAM211Jun181306hfwaterfall0.tif Driver: GTiff/GeoTIFF Files: ...
3
votes
2answers
635 views

gdal doesn´t support my *xyz file

I have just exported a tif DEM to XYZ format with SAGA GIS. When I check this with gdalinfo it states that the format is unsupported? Seeing as I converted this with saga gis which used gdal I find ...
3
votes
1answer
491 views

Is it possible to produce a decent looking shaded-relief map at 50K-scale-or-larger using SRTM3 DEM?

I'm trying to make a decent-looking shaded-relief map of a particular area in 50K-scale-or-larger using open-source tools. The best (and free) option I have for DEM data for that area right now is the ...
1
vote
1answer
421 views

How to extract a WCS layer with gdal_translate?

There is a certain WCS layer that I'd like to save it to my local computer for analysis. According to GDAL WCS driver draft: "Accessing a WCS server is accomplished by creating a local service ...
9
votes
2answers
648 views

Python, GDAL and building raster attribute tables

I have an integer raster for which I would like to build a raster attribute table using Python and GDAL. I can create a GDAL raster attribute table in Python as follows: >>> rat = ...
4
votes
2answers
607 views

How to programmatically create exports from PostGIS database into different file formats (e.g. shapefile, dxf, dgn,…)?

Is there a kind of wrapper/library written in python which would enable me to programatically create exports from postgis database into different file formats (most common being, shapefile, dxf, ...
11
votes
2answers
2k views

How to split shapefile per feature in Python using GDAL?

is it possible to split a shapefile per feature in python? (best would be a solution where i can temporarily save the resulting vector objects to the memory instead to disk). The reason: I want to ...
2
votes
1answer
188 views

gdal_retile question [closed]

I have been trying to get gdal_retile working with the --optfile command. At first it did not work because in my files.txt I have '-' characters in the filename. I tried adding quotes around the ...
3
votes
1answer
440 views

Get the geotransform in Mercator given a pixel size away from the origin

I have a radar image with the following information: The projection is Mercator (Mercator_1SP in GDAL) The projection origin is 0º lat 0º lon At a given geographic coordinates (e.g. lon=10º ...
2
votes
2answers
209 views

Driver Error : DTM to Raster using GDAL_GRID

I am looking to interpolate a DTM (eventually over 4,000 DTMs) that is in a personal geodatabase format. I would like to do this in as few step as possible, eventually writing a batch script to do ...
5
votes
1answer
449 views

Slow FileGDB Performance in GDAL

I am attempting to write many ASPRS LAS point files into an ESRI File Geodatabase using my compilation of GDAL 1.9.2. The FileGDB driver for GDAL/OGR appears to be incredibly slow when writing large ...
1
vote
1answer
115 views

keeping the coordinate system of raster files in the resulting raster file after operation with numpy

I'm trying to read as array two raster files using ReadAsArray in GDAL/python, perform some map algebra operation on it using numpy and then write the resulting raster file in GTiff format using ...
1
vote
1answer
741 views

DLL load error with Gdal polygonize in QGIS 1.8 (windows 7) [duplicate]

Possible Duplicate: Qgis polygonize problem - DLL load fail [I had already posted this question, but was closed and redirected to a supposedly duplicate question. I'm re-posting since I ...
5
votes
1answer
410 views

Where can I find GDAL bindings for Python 3.3?

I have been using GDAL 1.9.x on Python 2.7 and took the bindings from here. Is there a place to get GDAL bindings for Python 3.3 or is GDAL not officially supported in Python 3.3 yet?
1
vote
2answers
367 views

How to produce a mean raster file based on other raster files? [duplicate]

Possible Duplicate: Tools to flatten the jp2 to a single band (average the bands)? I have a number of raster files from a same area, that coordinates of corresponding pixels perfectly match ...
3
votes
1answer
374 views

Problem editing raster and saving with projection and transformation intact

I am trying to make a custom tool to Arcmap using Python. It is supposed to read in a raster DEM, a point shape file and a distance, and output another raster. The output raster should be ...
2
votes
2answers
357 views

--optfile not working with gdal_retile?

I have the following gdal_retile command I am trying to run, but it errors on the --optfile: gdal_retile.py -v -r bilinear -levels 8 -ps 512 512 -co "tiled=YES" -targetDir pyramid --optfile files.txt ...
1
vote
0answers
287 views

ArcGIS 10.1 crashes when a button in an addin invokes GDAL DLL through SWIG CSharp wrapper

I am developing an addin with a button that tries to do coordinate transformation using GDAL OSR. When I create a new CoordinateTransformation object ArcMap crashes and close. If a debug against a ...
3
votes
2answers
536 views

Why does Ogr2Ogr Postgres import fail?

I have a lot of S57 files that need to be imported to PostgresDB. The process takes over 2 hours. In order to improve performance, I have tried using PG_USE_COPY to use copy instead of DB inserts. I ...
6
votes
2answers
617 views

How to add GDAL to existing Python installation?

I would like to download and install GDAL for Python version 2.6 [MCS v.1500 32 bit(Intel)] on win32 that supports reading FileGDB's with Python. But I am not sure how to go about this. I found one ...
1
vote
0answers
80 views

resampling errors

I have to reproject my images, but it contains some resampling errors. Any one has any idea on significance of resampling error. Are they neglectable? Is there any good documentation to refer to? I'll ...
3
votes
1answer
704 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 ...
4
votes
4answers
3k views

How do I use ogr2ogr to convert a GML to shapefile in Python

I am trying to convert a GML to an ESRI shapefile using ogr2ogr utility in a python script. I have successfully installed the GDAL/OGR package via osgeo but I am now struggling to find/understand ...
4
votes
1answer
329 views

How to use the OGR/GDAL configure options in QGIS python environment

I’m looking for a solution how to use the OGR/GDAL configure options (--config) in the qgis python environment. http://trac.osgeo.org/gdal/wiki/ConfigOptions Setting special configure options is ...
3
votes
2answers
179 views

How can I create a portable script for coordinate conversions to run from USB?

I would like to be able to convert Lat long values in wgs84 to the swedish RT90, the values are in CSV files under headings X and Y. I would like to have this portable, an eg runned in a script in ...
5
votes
1answer
491 views

Error Installing libgdal1 Package

I'm trying to install the libgdal1 package on Ubuntu 12.04, but the apt is returing this error: dpkg: error processing /var/cache/apt/archives/libgdal1_1.9.1-2~precise4_amd64.deb (--unpack): ...
1
vote
1answer
555 views

How can I add all the 4 corners as GCP in an image that I want to reproject?

I am working on a Modis - oceancolor image and I want to reproject it to wgs84. I have the upper left/right and lower left/right values from the metadata of the image and I want to add them all in the ...
2
votes
1answer
375 views

Flattening JP2 files to single band (i.e. averaging bands)?

I have created thousands of png tiles using ArcGIS and have merged groups of them into large jp2 files for serving via GeoServer. I need to flatten these to a single band so that I can apply a colour ...
4
votes
1answer
1k views

GEOS and Shapely installation on Windows

When trying to install Shapely on my Windows 64bit computer, I cannot get the GEOS library to work. So far, I have run the OSGeo4W installer from which I installed GDAL (I believe the geos library ...
3
votes
3answers
1k views

Should I use gdal_retile.py or gdal2tiles.py?

I have many GeoTif files that I want to display in google earth with Geoserver. I am wondering if I want to use gdal_retile.py or gdal2tiles.py to create the image pyramids for them to be usable in ...
12
votes
2answers
946 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
364 views

Latest version of ogr2ogr fails with Layer Creation

I am able to successfully run ogr2ogr to add S57 data to PostgresDB. The version I use is GDAL 1.7.0b2, FWTools 2.4.7 and release date 2010/01/19 I upgraded my version to GDAL 1.9.1 and release date ...
3
votes
1answer
534 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" ...
5
votes
1answer
426 views

How to read Geotiff tags and pixel per pixel values using GDAL C#

I'm new here and I'm new at programing with GDAL C#. My question is how I can get the value of SRTM GEOTIFF pixel to compare with it's neighbors, and how can I get the values of Geotiff tags.
4
votes
3answers
1k views

gdal_merge to merge hundreds of png files

I need to merge about 760 png files into one and am looking at using gdal for this. This is related to the following question Getting a constant background color for transparency ...