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

1
vote
2answers
448 views

How to calculate Image boundary / footprint of satellite images using open source tools?

I'm referring to this question&discussion: How do I create a shapefile showing footprints of Rasters? I need to create the polygon outlines of several single layer raster images, not the ...
5
votes
1answer
711 views

Python - OGR: Transform coordinates from meter to decimal degrees

I have a point in a shapefile. The shapefile is in a WGS_1984_Web_Mercator Projected Coordinate System. I want to get the coordinates of that point in decimal degrees with GDAL/OGR in Python. I am ...
-5
votes
1answer
124 views

Problem installing GDAL and fiona module in arcmap10? [closed]

I am getting message like this Traceback (most recent call last): File "", line 1, in File "C:\Python26\ArcGIS10.0\lib\site-packages\gdal.py", line 2, in from osgeo.gdal import ...
0
votes
1answer
18 views

Is there a way to tell the supported input formats of libmap.dll without the associated .exe?

I have only compiled Mapscript_csharp and supporting dlls used in a .NET application - not the executables. Is there a way to tell what input / output formats are supported by these dlls?
0
votes
0answers
285 views

Reprojection problem with gdalwarp

I have a GeoTiff file with WGS84 datum projection and I'm trying to reproject it using gdalwarp. Using this command gdalwarp -t_srs '+proj=tpers +azi=19 +tilt=30 +lon_0=-109.5 +lat_0=23 +h=1274199.4 ...
0
votes
0answers
458 views

How do I add elevation data to a GeoJSON line from a GPX file?

I understand that elevation can be added as an optional third value in Coordinates (lat, long, elev) in a geojson file, but I can't find a way to actually do that from a gpx file. I've used ogr2ogr to ...
1
vote
1answer
29 views

GDAL: What is the best way to test if a CPLString is not empty?

If I define a CPLString this way: CPLString osString; And I want to check if osString is not empty, what would be the best way? My bet: if (osString.c_str() == NULL) But I find it a little ...
3
votes
2answers
268 views

informations about “gdaltransform” utility of GDAL

I have a doubt about GDAL reprojection. I use "gdalwarp" utility to reproject an image in a specific projection system (for example WGS84 UTM49). "gdalinfo" on original image displays this corner ...
1
vote
1answer
154 views

Python GDAL: find out spikes and pits in DEM

I'm working with on a DEM (Geotiff) in Python GDAL, and I'd like to find out the spikes or pits in my terrain. I have no idea how to proceed, does anyone have good algorithm to help me out? Thanks, eo ...
3
votes
0answers
308 views

How to set GDAL/cs2cs/proj to use NTv2 grid shift file when using EPSG codes?

I am using the GDAL API for Python for a script that takes only EPSG codes for inputs to specify source and target coordinate systems. While the cs2cs command seems to handle quite well ...
1
vote
0answers
87 views

Python GDAL/OGR: check lines geometries connectivity

I'm working with shapefiles representing roads, so basically it's stored as (poly)lines. I need to perform a road connectivity checking. In other words, I want to ensure that roads meant to ...
3
votes
1answer
115 views

Where can I find a python GDAL 1.8.0 package?

I have python 2.7.4 installed on my RHEL 5 machine along with gdal-devel-1.8.0 and gdal-1.8.0. I tried to install the GDAL 1.9.1 package available on python's page ...
0
votes
1answer
367 views

Installation of GDAL 1.10 in Windows 8 and Ubuntu 12.04 [duplicate]

Is there an easy way to install GDAL/OGR 1.10 in Windows 8 and Ubuntu 12.04.
1
vote
0answers
119 views

setup Mapserver WMS from tile files

I'm totally new to Mapserver and I've been banging my head trying to figure out how to mount a WMS server that will seed images of a small geographical area to a OpenLayers HTML page. I already ...
3
votes
2answers
595 views

circle detection on binary image

My image is gray scale image which i have converted to binary by smoothing and then edge detection method. Now i want to place the circle /ellipse on edge detected. I don't know how to proceed ...
0
votes
2answers
325 views

Why gdalwarp cannot make transparent black nodata?

I have a MODIS ariel image geotiff and I want to remove black nodata pixel (this is Aerosol Of Thickness) 1 band image. I've done like this C:>convert -transparent black ...
7
votes
1answer
1k 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 ...
1
vote
1answer
134 views

python gdal_polygonize error

I am trying to use gdal_polygonize to create a polygon from a binary raster. My code is: os.system('gdal_polygonize.py ' + filename + ' -f ESRI Shapefile ' + shapefile) where filename is a string ...
1
vote
1answer
120 views

OGR UnionCascaded gives strange lines (flat polygons)

I'm using OGR to perform a union using UnionCascaded. All my polygons have been created by GDAL polygonize which ensures that all polygons follow pixel edges and should therefore align nicely. But I ...
0
votes
1answer
255 views

CSV and GDAL 1.10 in QGIS

It seems the GDAL upgrade from 1.9 to 1.10 changed the way CSV files are handled so that I now have a non-functioning QGIS. App hangs trying to import csv data file when previously had no problems ...
1
vote
0answers
391 views

Using gdal to convert from UTM to equidistant lat/lon

I am getting rather confused with converting a UTM gridded dataset into a equidistant/rectangular latlon gridded dataset. Can anybody help me? I have a UTM set with the following projection: ...
0
votes
0answers
132 views

QGIS GDAL polygonize - how to set scale

I've got some old contour plots, paper copy only, that I'd like to be able to put into QGIS. I get as far as a successful polygonize into shp and it looks fine, but the scale is wrong. How might I ...
1
vote
0answers
118 views

How can I call gdal_translate from Java code?

How can I execute gdal_translate in Java? For example, will something like this work? gdal_translate "PG:host=localhost dbname='xyzy' schema=public' table='xyz'" "d:\xyx\ab.tiff"
1
vote
1answer
571 views

GDAL bindings problem with ArcMap 10.1 Python tool

I'm really getting mad with this. ArcGIS sometimes can be so just...irritating! Forgive me if this question is by any means duplicating this one or similar ones, but I really can't find my way out of ...
11
votes
1answer
256 views

Assign projection to a satellite image

Today I stumbled upon a picture taken from the ISS looking down Europe and it is really beautiful. So, I started thinking I could add more to the picture. Like world borders, cities, road network, ...
0
votes
1answer
217 views

Is it possible to use GDAL to combine (rather than mosaic) rasters

Is there a GDAL command which can combine multiple rasters so that a unique output value is assigned to each unique combination of input values, like the combine command in ArcGIS? gdal_merge seems ...
0
votes
1answer
91 views

how to stretch GMT raster with gdal to work with google

How can I stretch a PNG raster coming from GMT (Generic Mapping Tools 4.5.6) in such a way that does align with google satellite or bing? such raster is in wgs84 and crosses two UTM zones (17 and 18 ...
2
votes
3answers
191 views

How to test if geometry is inside another?

I need to tell whether geometries from a shapefile are inside other geometries or not. For example, I want to test if there are trees (represented as points in trees.shp) inside the urban areas ...
0
votes
1answer
157 views

Local variance image in python using gdal and a running window approach

I want a local variance image with a 3x3 of a geospatial raster image using python. My approach so far was to read in the raster band as an array, then using matrix notation to run a moving window and ...
3
votes
1answer
1k views

convert from HDF to GeoTiff

I am having some problems with the conversion from HDF to GeoTiff. I am following the steps described in this webpage. The problem is that I can't find what is the NDVI subdataset name (XXXX) in ...
12
votes
1answer
535 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?
1
vote
1answer
1k views

Python GDAL: Save array as raster with projection from other file

I have an array of data, and for each datapoint I know the latitude and longitude. I'd like to save it as a GTiff with the same projection as other rasters I have. This is what I've tried so far, but ...
2
votes
1answer
253 views

OGR - Weird GeomType -2147483645 on polygon shapefile

I am trying to import a shapefile (multipolygon) in an app with ogr but I got an error about invalid geom type. If I inspect the data with ogr like this: datasource = ogr.Open("c:\\temp\\data.shp") ...
3
votes
1answer
419 views

Export KML ExtendedData with OGR2OGR

I would like to export with ogr2ogr the data I have in a kml as ExtendedData <Placemark> <ExtendedData> <SchemaData schemaUrl="#15540_punti"> ...
0
votes
1answer
406 views

OGR OSM-driver unable to open datasource

I want to read an OSM file with OGR. I use Version GDAL 1.9.1, released 2012/05/15. To read a SQLITE file works fine: $ ogrinfo test.sqlite INFO: Open of `test.sqlite' using driver `SQLite' ...
1
vote
4answers
406 views

Get vector features inside a specific extent

The problem: I have a very large vector layer with many features and a much smaller raster layer in a defined region. I want to get only those vector features which are within the raster layers ...
2
votes
2answers
380 views

Why do I get different output from gdal2tiles vs. maptiler?

I am trying to generate map tiles using gdal2tiles for an area much larger than MapTiler can handle. However, my tile caches are not showing up on the maps. I checked gdal2tiles by running the same ...
7
votes
1answer
1k 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 ...
2
votes
0answers
95 views

efficiently slicing tiles with gdal C API

Im working a project using the gdal C api(in a block of objective-c) to slice TMS style png tiles. The code currently works, but I'm running into some trouble optimizing to run fast enough for use on ...
1
vote
0answers
155 views

using GDAL/OGR (faster)

Eventually, I was able to get the GDAL installed and compiled correctly using Visual C++ 2010. I even tested some of its methods on some geo-referenced raster files, and it looks like that all is ...
2
votes
1answer
110 views

gdal_retile.py tile naming scheme

Is there a way to define custom tile names during a gdal_retile.py process, other than the default _01_01. For instance, suppose I want to start naming my tiles at 100_100, is this possible?
0
votes
1answer
147 views

How to save image to disk with with gdal_translate?

I do: import os from gdal import * filename = 'c:/rast.jpg' outputTif = 'c:/rastOut.tif' os.system('gdal_translate filename outputTif') but it returns "1". gdal_translate doesn't save the image ...
1
vote
1answer
81 views

gdal_retile.py and gdal_translate Requirements

I basically need to use two GDAL utilities: gdal_retile.py and gdal_translate, on a Solaris 9 machine. I have tried compiling GDAL from sources, but this seems not to be working. So, since I have ...
0
votes
2answers
171 views

Handling MODIS NoData cells with gdalwarp

I'm using the gdalwarp tool to convert MODIS .hdf into GeoTiff files. I want to set the value of the nodata cells so that when I import the .tif image into Grass there is no black edge. I'm using the ...
1
vote
1answer
856 views

Python GDAL: Write new raster using projection from old

If I read in a raster image as an array, then make some changes to the values in the array, how do I then save the array as a raster with the same projection information as the original array? In ...
1
vote
2answers
725 views

Merging many rasters using Python

I have a list of many many rasters. I want to have a simple and easy to understand script that takes this and merges it into one big mosaic. I found a similar question but the accepted answer adivises ...
0
votes
1answer
225 views

gdaladdo positive and negative aspects?

I read and I read opinions and topics regarding this tool but I am still unclear if I should us this or not. In my project I have thousands off tiff files which need to be put on a map. The problem ...
0
votes
1answer
497 views

How to clip raster by multiple polygons in multi rasters?

Using QGIS, I want to clip a raster by multiple polygons contained in a single shapefile. I read that gdalwarp can clip my raster by multiple polygons BUT the output is a single merged clipped raster. ...
3
votes
1answer
105 views

Differences in raster extents in gdalinfo, SAGA & Quantum after resampling

I recently resampled a raster file using SAGA's Resampling module. The aim was to widen the extent of the new raster to two decimal places, where the original raster has three decimal places. The ...
1
vote
1answer
891 views

QGIS can't find osgeo after GDAL upgrade

Mac OS 10.6.8: I installed QGIS 1.8.0 with GDAL 1.8; then another app requires GDAL 1.9 and installed it; now QGIS can't find osgeo [python-gdal]. (I can't find any trace of GDAL 1.8 left in my ...