Stack Exchange Network

Stack Exchange network consists of 174 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [gdal]

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 ...

88
votes
13answers
94k views

Installing GDAL with Python on windows?

Can anyone explain how to install GDAL/OGR with Python on Windows? I have Windows Vista and I have tried following the information on the website and it does not seem to bind the exe files for me. ...
55
votes
20answers
60k 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 ...
45
votes
7answers
16k views

Should GDAL be set to produce GeoTIFF files with compression? Which algorithm should be used?

I have a folder of GIS data that consists mainly of GeoTIFF files. The whole set weighs in at about 1.2 GB. I noticed that if I pack the contents into a tarball, it smashes down to about 82 MB. I ...
40
votes
3answers
13k views

What is Lanczos resampling useful for in a spatial context?

GDAL includes a resampling method beyond the normal mix of nearest neighbor, bilinear, cubic and splines: "Lanczos windowed sinc resampling". I understand that its a convolution filter, but unlike ...
39
votes
5answers
27k 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.exe from the file system, but rather call it somehow in code so I do not ...
39
votes
1answer
19k views

Getting pixel value of GDAL raster under OGR point without NumPy?

I am working on a computational model of the abundance of wild pollinators across a landscape. The model itself is complete, and I am now struggling with a post-processing step. I have my GDAL ...
35
votes
3answers
32k views

How can I get a shapefile from a postgis query?

My data is in a PostGIS database. I would like to generate a shapefile from a query. How can I do that?
32
votes
2answers
24k views

How do I check which version of GDAL I have and how many?

I have installed in my computer the following software: PosgresSQL 9.3.2.3 PostGIS 2.1.1.1 QGIS 2.0.1 Dufour (this one installed through the OSGeo4w package) QGIS 1.7.0 Wroclaw As far as I know all ...
31
votes
3answers
3k views

Creating transparent hillshade?

I get elegant transparent hillshades via a combination gdal and convert. Compare to grey-based hillshade commonly used, such transparent hillshades are very cool because they can be placed between the ...
29
votes
6answers
3k views

Learning resources for GDAL?

I quite often get stuck trying to solve particular problems, or getting syntaxes right. Apart from the details provided in the links on GDAL Utilities, does anyone know of any 3rd party resources/...
27
votes
4answers
33k views

Installing GDAL on Linux Ubuntu Server?

apt-get install python-gdal installs version 1.6.0 and this is already installed. I would like to install the 1.9.0 version. Do I have to de-install 1.6.0 first and then compile manually or is there ...
27
votes
4answers
27k views

Writing numpy array to raster file

I'm new to GIS. I have some code that converts infrared images of Mars into thermal inertia maps, which are then stored as 2D numpy arrays. I've been saving these maps as hdf5 files but I'd really ...
27
votes
1answer
513 views

How to calculate distance to a feature with gdal_proximity?

I'm using gdal_proximity to find the distance to the nearest major river across the USA (lower 48 states). I've projected the NHD+ network flowlines to Conus Albers (epsg:5070), selected rivers with ...
26
votes
4answers
30k views

How to fully load a raster into a numpy array?

I have been trying to check my filters on DEM raster for pattern recognition and it is always resulting in missing last rows and columns(like..20). I have tried with PIL library, image load. Then with ...
25
votes
4answers
47k views

Splitting shapefile into separate files 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 ...
25
votes
4answers
6k 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 ...
24
votes
3answers
24k views

Clipping raster with vector layer using GDAL

I have installed GDAL using the Osgeo installer. How can I clip a raster layer with a vector layer programmatically? Is there any GDAL API available which can help me with this? I am using Python.
23
votes
1answer
3k views

Can File Geodatabase API be Redistributed?

Esri's file geodatabase API can now be used with GDAL/OGR - http://www.gdal.org/ogr/drv_filegdb.html however each individual seems to have to login to Esri and download the API files individually. For ...
22
votes
3answers
19k views

How to get raster corner coordinates using Python GDAL bindings?

Is there a way to get the corner coordinates (in degrees lat/long) from a raster file using gdal's Python bindings? A few searches online have convinced me that there is not, so I have developed a ...
22
votes
3answers
19k views

Selecting feature types when using ogr2ogr to convert to shapefile?

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 ...
22
votes
3answers
3k views

How are GeoTIFF pyramids/overviews standardised?

While pyramids/overviews are not part of the GeoTIFF standard, many tools support creating them. For example vips/nip2, Orfeo Toolbox (otb) and ossim which all promise some support for creating those. ...
22
votes
7answers
11k views

How to cut a 'tiled image' from a raster dataset using GDAL?

I am working on an app that needs to create tiles (jpeg) from a given raster dataset. My initial inclination was to use gdal2tiles.py, but its performance seems like it could be improved. That leads ...
21
votes
3answers
21k views

gdalwarp cutline along with shapefile

I am trying to clip my Raster file according to a shapefile geometry. I am using the following code gdalwarp -cutline INPUT.shp INPUT.tif OUTPUT.tif But it is resulting a black color outside the ...
20
votes
4answers
5k 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 ...
19
votes
4answers
16k views

Using ogr2ogr to convert 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 ...
19
votes
1answer
4k views

How to rename field names in a shapefile from the commandline?

I need to rename the field names (attributes) in a shapefile. Is there an easy way to do this from the command line?
18
votes
8answers
64k views

Installing GDAL and OGR for Python?

I am developing some code in Python, and I want to use the GDAL/OGR python bindings. What is the minimum I need to install in order to use the GDAL and OGR python bindings? What is the easiest way to ...
18
votes
2answers
7k views

Smoothing/interpolating raster in Python using GDAL?

I am developing in Python and using GDAL from OSGEO to manipulate and interact with rasters and shapefiles. I want to take a shapefile that has point features and interpolate it into a surface raster....
17
votes
3answers
27k views

How can I get the proj4 string or EPSG code from a shapefile .prj file? [duplicate]

I have to work with a shapefile defined in an unusual projection, which my GIS software does not know about. How can I get the proj4 definition or the EPSG code for the projection?
17
votes
1answer
1k views

Why do people pronounce GDAL (Geospatial Data Abstraction Library) like goo-dal?

Ok, a bit of a silly question, but: I've heard people pronouncing the library GDAL (Geospatial Data Abstraction Library) like goooo-dal. My natural inclination was to pronounce it geeee-dal, as in ...
17
votes
2answers
12k views

API documentation for Gdal/Ogr with C#

I want to assign projection to shapefile in my web C# application, for this i am using Gdal/OGR/OSR C# bindings and add osr_csharp.dll and ogr_csharp.dll references. But i am facing some difficulty ...
17
votes
2answers
7k views

How to access feature classes in file geodatabases with Python and GDAL?

I am trying to access a vector dataset in a ESRI File Geodatabase using Python + GDAL. I have successfully compiled GDAL with the file geodatabase API. The FileGDB driver is working correctly since ...
17
votes
3answers
14k views

How do I change geotiff “no data” color to white/transparent rather than black?

I am preparing data for inclusion into a geonetwork node, where the maps are displayed via geoserver. The problem is that cells that have"no data" are displaying as black in the geotiff and I am ...
17
votes
1answer
9k views

Replicating result of gdalwarp using gdal Python bindings

I am trying to re-project/resample with the GDAL python bindings, but am getting slightly different results compared to those from the command line utility gdalwarp. See update below for shorter ...
17
votes
2answers
4k views

Reprojecting WGS 1984 Web Mercator (EPSG:3857) in Python with GDAL

I am reprojecting rasters in python using GDAL. I need to project several tiffs from geographic WGS 84 coordinates to WGS 1984 Web Mercator (Auxiliary Sphere), in order to use them later in Openlayers ...
16
votes
5answers
5k views

Splitting raster into smaller chunks using GDAL?

I have a raster (USGS DEM actually) and I need to split it up into smaller chunks like the image below shows. That was accomplished in ArcGIS 10.0 using the Split Raster tool. I would like a FOSS ...
16
votes
3answers
15k views

Calculate lat lon bounds for individual tile generated from gdal2tiles

I have multiple tile sources generated using gdal2tiles I'd like to present on the same map. So when presenting a tile I need to determine which source to serve it up from, checking against the bounds....
16
votes
3answers
3k views

Generating contours from DEM using GDAL?

We have just received a large set of DEMs at work and I would like to generate contours from them. The DEMs have a resolution of 1m and a size of 1kmx1km. Output from gdalinfo: Driver: AAIGrid/Arc/...
16
votes
3answers
12k views

How to georeference a web mercator tile correctly using gdal?

As an example I will take the following tile http://a.tile.openstreetmap.org/3/4/2.png and save it as "4_2.png". The WGS84 coordinates of this tile can be calculated or read there by clicking the ...
16
votes
3answers
1k views

Geotransformation for polar stereographic?

I'm currently working to import CANGRID climate data (provided as Surfer Grid ascii, ".grd" files) into ArcGIS. The grid is 95 rows by 125 columns in size.Metadata provides lat/lon of origin (lower ...
15
votes
1answer
19k views

GDAL/Python: How do I get coordinate system name from SpatialReference?

In Python, using GDAL, I've extracted a raster's projection as a WKT string as follows: wkt = dataset.GetProjection() # wkt is 'PROJCS["GDA_1994_Transverse_Mercator",GEOGCS["GDA_1994",DATUM["GDA_1994"...
15
votes
1answer
5k views

How to change the field value of a Shapefile using GDAL/OGR?

I am trying to change the field value of a shapefile. However it doesn't change the value. What's wrong with my code? import ogr driver = ogr.GetDriverByName('ESRI Shapefile') fn = 'dist.shp' ...
15
votes
2answers
4k views

Adding band to existing GeoTiff using GDAL?

This seems like it should be straightforward but I can't for the life of me figure it out. Given an existing GeoTiff containing a single band, how can you add another band to that GeoTiff? Is the ...
15
votes
2answers
5k 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 (...
15
votes
4answers
7k views

Smoothing/reinterpolating raster with GDAL?

Short: I want to change the resolution of a raster and smooth the gray ramp like shown in the images bellow. The preference is to use GDAL, PIL or Numpy. Description: I'm kriging points into rasters ...
14
votes
8answers
6k views

how to overlay shapefile and raster?

I have a shapefile with polygons. And I have a global raster file. I want to overlay the shapefile's polygons onto the raster grid and calculate the mean raster value for each polygon. How can I do ...
14
votes
1answer
3k views

File size inflation normal with gdalwarp?

After using gdalwarp to project and align-to-grid (via -tap) a number of rasters I noticed that the output rasters were significantly larger than the original rasters. A fairly thorough web search ...
14
votes
1answer
8k views

Georeferencing using GDAL?

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 ...
14
votes
5answers
18k views

Installing GDAL on MacOSX?

Does anybody know what's the easiest way to install GDAL 2.0 on a MacOSX? The usual suspects kingchaos binary distribution and homebrew seem to be stuck on GDAL 1.11 I need support for TopoJSON ...
14
votes
1answer
628 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 ...