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

4
votes
2answers
627 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
189 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
462 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
232 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
466 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
125 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
773 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
448 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
395 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
385 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
390 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
295 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
598 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
675 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
2answers
802 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 ...
6
votes
4answers
4k 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
357 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
187 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
516 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): ...
2
votes
1answer
673 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
404 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
2answers
2k 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
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
390 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
561 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" ...
6
votes
1answer
520 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 ...
5
votes
2answers
1k views

Shape-file encoding problem in QGIS 1.9.0 (built with GDAL 1.9.2)

I have googled a lot but failed to find the answer. My configuration is: Ubuntu 11.10 Quantum GIS - 1.9.0-Master 'Master' (caae7b7) GDAL 1.9.2, released 2012/10/08 Shape-files refuse to open in ...
5
votes
1answer
415 views

GDAL.RasterizeLayer not ported to the C# bindings?

I am trying to use the OSGeo.GDAL.Gdal.RasterizeLayer(…) from GDAL 1.9.2 and it requires some SWIGTYPE parameters: OSGeo.GDAL.Gdal.RasterizeLayer( OSGeo.GDAL.Dataset dataset, int bands, ...
13
votes
3answers
4k views

Writing numpy array to raster file

I'm new to GIS, and I'm lost with something I feel ought to be relatively simple. I have some code that converts infrared images of Mars into thermal inertia maps, which are then stored as 2D numpy ...
2
votes
2answers
337 views

Batch process of GDALWARP in FWtools

I am using FWtools for pre processing of time-series MODIS data, I want to mosaic series of two scenes. But I don't know the batch code for gdalwarp. I am using windows. Kindly help me out.
12
votes
2answers
1k views

How to equalise colour, brightness and contrast across several scanned paper maps

I am a vectors guy mostly, but for a current project I need to work with a bunch of old scanned paper maps (ww2 bomb damage maps for London, if your interested!) We scanned and georeferenced the ...
5
votes
1answer
1k views

What is the best practice for reading data from a raster band with Java and GDAL?

After using GDAL with Python for a year or so now, I'm getting back into Java and trying to figure out how to use GDAL with it. I am wondering if anyone has had experience with reading data from a ...
1
vote
2answers
574 views

Proper gdaltrasform syntax to transform a list of coordinates?

this is probably a very simple question, but I'm very new to GIS so I just need one last nudge. I've got a sense of the Linux utility gdaltransform and I am feeding it coordinate pairs in ...
2
votes
1answer
398 views

How do I convert EPSG:26986 coordinate pairs to EPSG:4326 from the MassGIS web service, the Linux command-line, or PHP?

I'm a GIS newbie, normally a PHP web applications developer. I am trying to map all boundaries of the MA state representative districts. So I just need to convert polygon coordinates from the ...
1
vote
1answer
443 views

Converting arcpy based python script to OGR for semi-automated Metadata creation

I have created a script to automatically generate base ANZLIC standard metadata using arcpy based coding. Can anyone help me convert it to OGR/GDAL python? This is so that it can be run in QGIS and ...
1
vote
1answer
224 views

How to choose suitable pyramid levels and tile size?

I have a big tiff file and i want to create tiles with gdal_retile. However I am still concerned about the pyramids levels and the tiles' size. how can I choose suitable pyramids levels (for example ...
1
vote
1answer
407 views

Why do glob.glob and gdal_translate not work for me?

I have a little problem with my python script. I wrote it to convert MODIS data from .hdf into .tif. It will convert all files in a folder using a loop. I do not know which command I must enter to ...
1
vote
0answers
74 views

Create 6 band IMG file with GDAL out of two JPG files, one with WorldFile coordinates

I have two JPG files 1 ortomap: http://imageshack.us/photo/my-images/703/ortofoto1.jpg/ with WorldFile: 0.683772574816354 0.000000000000000 0.000000000000000 -0.683739485024190 ...
1
vote
1answer
349 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
399 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 ...
5
votes
1answer
154 views

How to Add FileGDB driver to PyGDAL installation

I installed ArcGIS 10.0 and it installed Python 2.6.5 with it. Later, I installed pyGDAL using the windows installer within that Python installation. I cannot find a specific place where GDAL is ...
4
votes
1answer
223 views

Update PostGIS after updating GDAL?

I am using PostGIS 2.0 and I want to use the raster functionality through mapserver. It seems as if the latest official GDAL release 1.9.1 has an unstable PostGIS raster driver, as pointed out in a ...
0
votes
1answer
210 views

PostGIS raster support

I have the following requirements to manage raster data 1) Convert raster-maps from BSB format to PNG or GeoTIFF 2) Cut the borders of the map 3) Merge individual maps to form a contiguous coverage ...
4
votes
1answer
581 views

gdal_merge produces incomplete mosaic

I am using GDAL 1.91 on a Windows 7 machine. I have been mosaicing large numbers of rasters together and most of the time this process works exactly as expected. However a very few of my mosaics ...
5
votes
2answers
1k views

How to use gdal utilities to subset from a raster?

How to you use gdal to subset a smaller raster from a bigger one given the bounding coordinates? EDIT: Can you subset a region from a vrt ?
1
vote
1answer
171 views

gdal_rasterize -at function still won't capture covered pixels properly

I'm currently using the gdal_rasterize tool in bash. My code looks like this: gdal_rasterize -l $(basename $file .shp) -burn 1 -a_nodata 0 -te -180 -90 180 90 -tr 0.1 0.01 -at -ot Byte $file ...