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
1answer
231 views

Multiple versions of GDAL (32bit and 64bit) on same machine

I currently have GDAL installed on my Windows 7 64bit workstation. It is the 32bit version. I chose the 32 bit version because I wanted to use the GDAL Python bindings with a 32bit installation of ...
4
votes
1answer
609 views

Python equivalent of gdalbuildvrt

Is there a way to perform the same task as the gdalbuildvrt utility using GDAL Python bindings? So far I have not found any way to do this other than creating a vrt of a single dataset and manually ...
1
vote
0answers
333 views

OGR-Driver GML axis order handling (GDAL, QGIS)

By default the OGR-Driver will swap the axis order to longitude, latitude for geographic coordinate systems. "Since OGR 1.8.0, the GML driver has coordinate system support. This is only reported ...
2
votes
1answer
189 views

Tiling Rasters based on an existing Polygon Shapefile (Vector Grid)

I have a tile index representing the tile scheme used to tile Lidar data. I want to tile imagery rasters using the same tiling scheme. I have not found anyway using open source tools to tile a ...
3
votes
1answer
140 views

Two versions of gdal_merge and two diffrent effects

Is there any possibility that two diffrent versions of gdal (1.7.3 and I believe 1.9.0) can give me two diffrent results with given: two RGB files.Each twith 3 layers - input1.png and input2.png ...
5
votes
1answer
809 views

ESRI GDB Support in GDAL / ogr2ogr on Mac OS X 10.8

How do you incorporate support for ESRI's fgdb format into GDAL / ogr2ogr on Mac OS X 10.8?
6
votes
1answer
276 views

How to detect a GDAL/OGR 'Warning'?

When I run a batch script using a GDAL/OGR program, I can detect if an error occured by checking ERRORLEVEL for the value 1. When using GDAL/OGR in python or .NET, I get an exception when an error ...
3
votes
2answers
250 views

How to call gdal_retile from Python code?

I am trying to call gdal_retile from my Python code like this: import gdal_retile and in main: gdal_retile.main("-v -r bilinear -levels 4 -ps 2048 2048 -co \"tiled=YES\" -targetDir pyramid ...
1
vote
0answers
83 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 ...
7
votes
3answers
2k views

GDAL and Python: How to get coordinates for all cells having a specific value?

I've got an Arc/Info Binary Grid---specifically, an ArcGIS flow accumulation raster---and I'd like to identify all cells having a specific value (or in a range of values). Ultimately, I'd like a ...
3
votes
1answer
286 views

How to run OSgo4W shell script on Windows?

I have a lot of DEMs that I want to use to create hillshades, slopes and color relief from and would like to automate the processes using a shell script within OSgeo4W shell. Here is the code I used ...
16
votes
3answers
4k views

How to call gdal_translate from Python code?

Is it possible to use some gdal API to call gdal_translate from Python code? I do not mean simply executing the gdal_translate.py from the file system, but rather call it somehow in code so I do not ...
10
votes
2answers
1k views

What to do with -3.4e+38 nodata values?

I am trying to process some bioclimatic raster files, such as can be downloaded from http://www.worldclim.org/current (bioclim set). They seem to have nodata values set to -3.4e+38 according to QGIS ...
2
votes
0answers
173 views

projection discrepancy with gdal (1.9) and proj4

I boiled the problem down to that python example: import sys try: from osgeo import gdal from osgeo import osr except: import gdal import osr if __name__ == '__main__': version_num = ...
3
votes
1answer
400 views

ogr2ogr / GDAL remove features in FileGDB while in use by ArcGIS Server map service

I'm using ogr2ogr to pull features out of Oracle and put them into a File GeoDatabase. The FGDB is being used by an active ArcGIS Server map service. This is all done for performance / reducing DB ...
0
votes
0answers
48 views

How to select a subwindow from a source image in EPSG 3857?

I need to select a subwindow from an image projected to epsg 3857. If I use, for example: gdal_translate -projwin ${UL_x_3857} ${UL_y_3857} ${LR_x_3857} ${LR_y_3857} source_3857.tif ...
3
votes
3answers
668 views

DXF to georeferenced SHP

I know how to set projection and georeference bounds of scanned TIFF file, by using gdal, but I don't know how to do the same for vector file. For example I georeferenced TIFF file by using this ...
4
votes
3answers
525 views

Where can I find a GDAL distribution with DWG support?

I want to use the GDAL utility ogrinfo with dwg files. Where can I find an GDAL distribution with a support for DWG files (read only)? DWG support is moved out of the normal GDAL/OGR source ...
2
votes
2answers
167 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
1answer
774 views

“Please ensure that the GDAL_DATA environment variable is set and try again” on Kubuntu 12.10

I'm having a problem loading the 'convert os gml to ogr' plugin in QGIS 1.8 running on (k)ubuntu 12.10. Basically, clicking the plugin prompts the error: "Please ensure that the GDAL_DATA ...
6
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
109 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
138 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
250 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
674 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
171 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 ...
3
votes
1answer
255 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
246 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. ...
1
vote
2answers
473 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
757 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
275 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
207 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
218 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
275 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 ...
5
votes
2answers
443 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
248 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
591 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
460 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
402 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
603 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
578 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
179 views

gdal_retile question

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
417 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
194 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
421 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
109 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
711 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
388 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
333 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 ...