Tagged Questions

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
0answers
616 views

Why is GDAL ECW compression much poorer than ECW compression from ERDAS

We use ECW compressor shipped with ERDAS Imagine to export Orthoimagery for our different departments. Because we are facing problems with the Imagine license server, we had to choose GDAL (built with ...
6
votes
0answers
378 views

What is the correct syntax for using gdal_rasterize with a PostGIS datasource [closed]

I have a multipolygon layer in PostGIS and I would like to use gdal_rasterize to generate a raster. I cannot find the correct syntax to specify PostGIS as the data source. $ gdal_rasterize -a ...
3
votes
1answer
345 views

gdal_merge.py with single band raster

I'm trying to mosaic the True Marble world imagery 250M geotiffs using gdal_merge.py though one of them doesn't have RGB bands like the rest and gdal_merge fails (same number of bands is listed as a ...
0
votes
1answer
424 views

How to export only one band from an image using GDAL?

I have a multi-band GeoTiff, and I want to extract only the first band and write a new image consisting of only that band. How do I do this using GDAL?
1
vote
1answer
368 views

Python GDAL/OGR: WGS84 & meters units

I'm dealing with shapefiles projected in: GEOGCS["GCS_WGS_1984", DATUM["WGS_1984", SPHEROID["WGS_84",6378137,298.257223563]], PRIMEM["Greenwich",0], ...
2
votes
0answers
413 views

Use GDAL/Python to add rasters (tiff) to an ESRI File Geodatabase

I'm using GDAL python to create a bunch of tiff files from a netcdf file. def modelout_to_geotif(netcdfile, outdir, VarName): """ Given the path to the netcdf output, produces geotiffs for ...
1
vote
1answer
740 views

Python GDAL: degrees to meters without reprojecting

I am trying to get senseful results in my Python script. I want to measure distance between features in a shapefile. I'm working with the SRS: GCS_WGS_1984, so my units are in degrees. But for the ...
2
votes
1answer
389 views

Technique for splitting raster image into tiles

I have written a standalone application which renders a globe similar to GoogleEarth/WorldWind which is currently able to display raster imagery from WMS servers. The native geodesy of this ...
1
vote
1answer
1k views

How to transform World data from EPSG 4326 to EPSG 3857

I have tryed different opensource tools to transform a world data shapefile from 4326 to 3857 but I always have an error. Im not sure, but is it because of Antartica? I have tryed with GDAL and QGIS ...
3
votes
0answers
57 views

How do I get the domain from an OGRSpatialReference?

I would like to know what the projected bounds are for a spatial reference using OGRSpatialReference. Is there some way to use the GetProjParm method to get this information? An example of what ...
0
votes
1answer
101 views

What's regular blocking arrangement ? (GDAL)

I am using GDAL tools to retrieve PostGiS raster. The code(command line) I am using is as follows : gdal_translate -of AAIGrid "PG:host=localhost dbname='database' user='postgres' password='admin' ...
1
vote
1answer
201 views

How to set a spatial filter with Python/GDAL?

Is there a way with GDAL (Python) to set a spatial filter (bounding-box) on a raster before I process it? I want to read it as an array. BandReadAsArray(band, xoff=0, yoff=0, win_xsize=None, ...
1
vote
0answers
34 views

GDAL Plugin with a different version of GDAL aready built into QGIS

My qgis windows version is 1.8.0 Lisboa and use the gdaltools in version1.9.2. I want to use the plugin NetCDFBrowser, but this plugin need the gdal 1.10. I install 1.10 stand-alone gdal find here ...
3
votes
2answers
872 views

Getting started with gdal/ogr

I found some batch code using ogr2ogr in order to convert and reproject vector layers. The difficulty for me is to know what to install and which gdal commandline should I use, and how to configure ...
2
votes
1answer
126 views

Using the C++ GDAL API how do I set the GDAL_DATA path?

If I want to change the directory of the GDAL_DATA path or point to a different version of a csv file for EPSG numbers how do I do that using the C++ API?
1
vote
2answers
536 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
941 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
142 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
361 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
590 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
31 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
334 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
172 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
430 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
97 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
124 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
432 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
136 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
750 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
407 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 ...
8
votes
1answer
2k 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
169 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
135 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
258 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
544 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
146 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
135 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"
2
votes
1answer
728 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
301 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
232 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
103 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
224 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
206 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 ...
4
votes
1answer
2k 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 ...
14
votes
1answer
705 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?
2
votes
1answer
2k 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
272 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
494 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
468 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' ...