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
2k views

How to get GDAL to create statistics for GTiff in Python

I regularly create my own GeoTIFF rasters with GDAL in Python, e.g.: from osgeo import gdal from numpy import random data = random.uniform(0, 10, (300, 200)) driver = gdal.GetDriverByName('GTiff') ds ...
7
votes
4answers
2k views

Are there GUI resources for gdal and ogr?

For those who are new to the programming side of GIS, I feel that a lot of great tools and resources are inacessible. I am interested in GUI based wizards for handling geoprocessing with libraries ...
14
votes
11answers
10k views

Python GDAL package missing header file when installing via pip

I asked this question on Stack Overflow before I realised this site existed. As this site it more appropriate, I'll ask it here too. I'm trying to install gdal from pip pip install gdal inside a ...
1
vote
1answer
226 views

gdal_retile - Out of range band requested

I'm using gdal_retile to create ImagePyramid from a set of GeoTIFF images. I'm using following command: gdal_retile.py -v -r bilinear -levels 8 -ps 2048 2048 -co "TILED=YES" -co "COMPRESS=JPEG" ...
5
votes
2answers
978 views

GDAL - Perform Simple Least Cost Path Analysis

I am investigating methods to perform a simple least cost path analysis with gdal. By simple, I mean using the slope of a dem as the only cost factor. I would prefer to do using the python or .net ...
2
votes
1answer
2k views

Convert area from square meters to degrees globally

I have a global raster file with 3600x1800 pixels. I don't have any information about the projection used. The raster contains one band with a quantity with the unit "per square meter". How can I ...
0
votes
3answers
450 views

raster data retrieval from postgis using c++

I have raster data in a postgis server (I am able to see it in qgis). When I retrieve it and form gtiff image out of it using gdal, the generated image is not correct (there is a coloring problem). I ...
3
votes
1answer
214 views

how to square the pixels of a GeoTiff

I have a GeoTiff that has been supplied to me with non square pixels (I suspect due to reprojection to lat lon at some point). However the next step in my processing chain requires a Arc ASCII grid as ...
5
votes
2answers
568 views

Tiling high resolution image using GDAL tools

I have a .tif image that I georeferenced using QGIS utility. It was originally in UTM with NAD83 as datum. I need to make tiles out of it and get a kml file to overlay onto google earth. Here's how I ...
1
vote
1answer
117 views

Can I avoid runtime raster interpolation by creating an image pyramid to match the zoom scales?

For image quality and performance reasons, it would be best to perform interpolation of my large TIFF image of the world (10GB) beforehand to all the zoom scales I need to display at so that MapGuide ...
3
votes
2answers
384 views

Eliminate small polygons from shapefile

I downloaded GTOPO30 *.DEM files, and extracted a contour line at 2000m altitude using gdal_contour. In the resulting shapefile, I have many small polygons in areas where the altitude is close to that ...
1
vote
1answer
471 views

GDAL: Raster Intersection

I have two raster layers. (A) is state-sized, 30m. (B) is county-sized, 3m. I'd like to get a new raster (C) which is the part of (A) which coincides with (B), preferably scaled to the same ...
3
votes
1answer
1k views

How to use gdal_contour in QGIS Python console?

Is there a way to use the gdal_contour utility with-in the python console of qgis ? Thank you
2
votes
1answer
907 views

Filtering a raster by pixel value ? (Qgis)

I have downloaded Globcover data in GeoTIFF format. I want to create a new grid which only has values where the Globcover data (the only band in the TIF file) has one certain value (==190). Can I do ...
10
votes
1answer
1k views

Understanding the Map Tile creation properties

I am trying to create a tile base map. I have an image that I want to create the tile from and I am using GDAL. I managed to create the tiles but when I view the result it seems that the tiled map ...
3
votes
1answer
824 views

how do I get the “outer mantle” of a polygon shape file? (without joining all the polygons)

I am interested in finding all the vertices of all polygons which are not on the inside the most outer polygoin that can be formed in a shape file: Here is how I get all the vertices: def ...
7
votes
2answers
3k views

How to georeference using GDAL tools?

I have a .tiff image that I am trying to overlay in google earth. For that I require to georeference it first. The image projection is UTM with NAD83 as datum. I know the lat/long of the four corners ...
7
votes
2answers
8k views

ECW for QGIS 1.8 in Ubuntu 12.04

I have followed various online help guides on how to make ecw for QGIS 1.8 work in Ubuntu 12.04. I had it all working in Ubuntu 10.04 but since migration to 12.04 no such luck. Is there a guide that ...
10
votes
1answer
541 views

OGR/GDAL threading results in low core utilization

I'm trying to process some raster data using ogr/gdal and I can't seem to get full utilization of all the cores on my machine. When I only run the process on a single core, I get 100% utilization of ...
2
votes
1answer
507 views

How to add GDAL Programs folder to ~/.bash_profile

I have recently installed on Mac OS X Lion 10.7.4 the GDAL Framework 1.9 with Python bindings from [http://www.kyngchaos.com/software:frameworks][1]. I need to add the following Programs folder to my ...
1
vote
1answer
1k views

TMS format specification

I want to convert large map file to tms format. how directory structure, zoom levels , tile size, xml file creation have to be dealt with. Please point to the relevant docs. And is there any gdal c++ ...
4
votes
1answer
864 views

GDAL: How to build an image pyramid starting from the base

Normally I would create an Image Pyramid using gdal_retile.py on one large image or a set of large images. I want to do the opposite. I have all the zoomed in tiles on the base level and I want to ...
5
votes
1answer
392 views

Mosaicing indexed colour rasters with independent colour tables

I have a large number of single band rasters that are colour indexed TIF files. I want to mosaic them using GDAL (or similar non-proprietary software with a Python API I can call as part of a bigger ...
3
votes
1answer
837 views

Serving 4 band GeoTIFF in GeoServer 2.1.4

I have an uncompressed GeoTIFF file (figure 1) with the following band information: Image Structure Metadata: INTERLEAVE=PIXEL Band 1 Block=5000x1 Type=Byte, ColorInterp=Red Band 2 Block=5000x1 ...
3
votes
1answer
317 views

gdalinfo gives no georeferencing information for a PostGIS raster

I have a PostGIS raster dataset which I am trying to access via GDAL. I've had various problems with making this work, so thought I'd try running the gdalinfo command to check that GDAL can actually ...
3
votes
1answer
368 views

GDAL minidriver and Tilecache

I'm trying to use GDAL minidriver to read tiles folders generated from Tilecache seeding. I have set this in my GDAL source <ServerUrl>http://myserver.com/test/${z}/${x}/${y}</ServerUrl>, ...
1
vote
1answer
378 views

Rasterize difference between QGIS and FWTools

I am wanting to use FWTools to rasterize a large number of Shapefile polygons. I have used the Rasterize tool in QIS and it works great. But when I use the generated code in FWTools the process does ...
3
votes
3answers
3k views

How to generate a shapefile containing an UTM grid (graticule)?

I'm looking for a shapefile containing an 1000m or 2000m UTM grid ("graticule"). All I was about to find on the web is a 10000m UTM grid (here: http://en.giswiki.net/wiki/UTM) and I am aware that a ...
5
votes
3answers
2k views

How to generate colorful images using GDAL?

I am using gdal_translate -of JPEG -scale -co worldfile=yes ndvi.tiff output.jpg command to convert tiff images to jpg and resize it. But it generates only black and white images. How can I generate ...
13
votes
3answers
7k views

ogr2ogr converting to shapefile - how to select feature types?

Shapefiles can only have one geometry type, but many other formats (kml, geojson) can have multiple types. When converting to shapefiles, is it possible to to tell ogr2ogr to create multiple ...
3
votes
1answer
410 views

How does latitude of origin work in Mercator

I'm trying to make a coordinate transformation using GDAL/osr python bindings. I have a longitude an a latitude and I want to transform them into the mercator coordinates. So I define the mercator ...
7
votes
1answer
809 views

Can QGIS display GeoRSS in real time?

QGIS has the ability to display GeoRSS by following these steps: Click on "Add Vector Layer..." icon Make sure the Data Source is "File" and enter the GeoRSS URL with "/viscurl/" prefix on the ...
5
votes
1answer
1k views

Reproject script in Python with GDAL

I'm having a lot of trouble with GDAL. Besides just the at times lacking documentation, there seems to be little support in Python. Anyways this is just a batch reproject script that goes thusly: ...
1
vote
1answer
421 views

how to use gdal command inside php exec() function

I need to develop a webportal which will display the images generated by a GIS software package. I my development I need to run gdal command inside php exec() function. but I don't know how to do ...
0
votes
3answers
931 views

How to run GDAL command in php page? [duplicate]

I want to develop a webportal that will show the QGIS generated images on dynamic website. For this reason I need to run gdal_translate command inside php page. but I don't know is it possible or not. ...
3
votes
2answers
3k views

how to convert qgis generated tiff images into jpg/jpeg using gdal command line utilities?

I want to develop a web portal that will show the images dynamically produced by qgis . but qgis provided tiff images are too big in size and it is not possible to show it on website using image ...
5
votes
1answer
619 views

How to create tiles with points from MongoDB?

I have a (potentially) large data set stored in MongoDB. Each record is just lat, lon and some float-point value. Each pair of lat and lon is unique, but the whole set is not covering the Earth, there ...
2
votes
2answers
216 views

Issue converting DEM in ASCII with LO17 (apparent) projection

I'm a newbie where GIS is concerned, so please don't hesitate to point me to a manual. We've got DEM files without any projection information. One is called 000426_LO17_200M_BD16_25.ORT and contains ...
2
votes
1answer
275 views

only 1 tile visible in GE after gdal_translate -of KMLSUPEROVERLAY

I have a 3-band ASTER scene in any of several formats (PCIDISK, GeoTIFF, Imagine), already reprojected to the geographic coordinate system. I can view and manipulate the scene successfully in PCI ...
4
votes
3answers
2k views

'NoneType' object has no attribute

I am new to python geospatial programming. i ran the following script and got the corresponding error message >>> import osgeo >>> import osgeo.ogr >>> shapefile = ...
8
votes
8answers
7k views

Splitting shapefile into separate shapefiles for each feature using QGIS/GDAL/SAGA?

I have a polyline shapefile with about 20 kinds of roads. I would like to know how can I create automatically a vector file for each one. I am aware of how to do it using ArcGIS for Desktop - see ...
2
votes
3answers
1k views

Batch script converting and resampling orthophotos

I would like to convert more than 300 orthophotos produced at scale 1:8.000 with a pixel size of 0.8m X 0.8m from Ecw format to Tiff format and after resampling pixel values. My aim is to batch ...
14
votes
4answers
2k views

The best language to program with GDAL

I'm going to start to create programs to handle geographical data from a PostGIS server, and going to use GDAL. According to your experience, with the best language to program with it? I can program ...
6
votes
2answers
1k views

Big Satellite Image Processing

Im trying to run Mort Canty's http://mcanty.homepage.t-online.de/ Python iMAD implementation on bitemporal RapidEye Multispectral images. Which basically calculates the canonical correlation for the ...
3
votes
1answer
209 views

Problem with precision when convert from WKT to shapefile

I have a problem with converting from WKT to shapefile in pygdal. In wkt I have fields with precision (for example real 7.2). When I create a shapefile I have field numeric with no precision (numeric ...
2
votes
2answers
4k views

GDAL: Convert specific RGB (GeoTIFF) color to transparent

I have an RGB GeoTIFF file and I want to convert a specific RGB colour to transparent. I would like to use GDAL for doing this, but external Python scripts are also an option (though I haven't learned ...
8
votes
4answers
3k views

Open source software to merge large shape files

Some background: I have 2 polygon shape files with approximately 500,000 records each. The two files are identical in terms of field structure. Using ArcGIS 10.0 I can'nt merge these two files using ...
3
votes
1answer
679 views

Cell calculations; QGIS,GRASS,GDAL,SAGA

I recently asked a question about trying to calculate the slope over distance (manually) using QGIS, but then I came up with another way of solving the problem, but I would have to get some help for ...
0
votes
1answer
977 views

Trying to Install on MAC OS Lion and getting Python error message

Brand new to this and struggling to install QGIS on MAC OS X with Lion, I am not a techy - but noticed in the install read me that /usr/bin needs to be added to the GDAL.frameworks path. Failing to ...
3
votes
1answer
390 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 ...