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

0
votes
1answer
213 views

GDAL - projection from pixels to geo coordinates

I am using GDAL library in C++. On my dataset, I have obtained geo transform with GetGeoTransform to GT array (GCP is NULL). Than I use this values for this equation (Source): Xgeo = GT(0) + ...
3
votes
2answers
495 views

GDAL TMS (WMS) mini driver

I've a problem with GDALs TMS mini driver. I've created an xml file, which i can open in qgis. I am using local tiles on my computer which i prerendered with mapnik, the map covers entire country of ...
1
vote
2answers
856 views

GDAL/proj4 equirectangular projection coordinate system units assumed by OGRSpatialReference and GDALDataset::GetGeoTransform()?

How do I configure my OGRSpatialReference object to specify an equirectangular projection system http://en.wikipedia.org/wiki/Equirectangular_projection where the affine geo-transform matrix (provided ...
3
votes
2answers
726 views

combine raster with another raster mask file - gdal

I have 2 rasters in GeoTiff format. One has all of the land and water and the other is a mask file that shows either 0 for water or 1 for land. I basically want my output to just have the land values ...
0
votes
2answers
1k views

OS X Mavericks: accessing GDAL from Terminal issue

I successfully installed QGIS on my iMac with OS X 10.9.1 Mavericks using the installers provided by Kyngchaos. However, I'm having troubles accessing GDAL / ogr from terminal. It's giving me the ...
1
vote
1answer
580 views

What is the best way to ortho-rectificate image from the plane?

I am very new in GIS and I am trying to ortho-rectificate my images. But after putting it on the map I see that results are not precise. We have images taken from the airplane, it's longitude, ...
3
votes
2answers
312 views

Working with GDB in QGIS 2.0 - plain language instructions please!

I'm a QGIS user running v2.0 (Dufour) on Mac OSx10.8.5. I have some data in .gdb format and need to access this and work with it in QGIS. I'm a relatively savvy computer user, but not at all a ...
1
vote
1answer
148 views

Polygons not appearing when creating shapefile based on criteria using GDAL and Python (area calculated properly)

I am looping over the features of a shapefile containing multiple polygons. If certain criteria are met, I want to save the current feature in a new shapefile. I am doing this using GDAL 1.10.0 with ...
3
votes
1answer
354 views

Controlling output linewidth linestyle using gdal_rasterize shapefile shp to Geotiff

I have successfully combined a shapefile (.shp) and a GeoTiff using gdal_rasterize. The shapefile displays as one pixel wide in the GeoTiff, but I need it to be wider and more visible, say, 5 pixels ...
1
vote
1answer
415 views

using GDAL/python to stack georeferenced images of different sizes

I am in the process of porting a code I wrote in IDL (interactive data language) to python but am running into a bit of a problem that I am hoping someone can help me with. The code goes like this: ...
2
votes
2answers
509 views

GDAL_TRANSLATE Assumes Picture is Orientated North?

I used gdal_translate to change a jpeg to a geotiff. Attempting to use gdalwarp to create a mosaic I found the images skewed and this lead me to the observation in the title of the question. Below ...
1
vote
0answers
129 views

Strange black pixel on edge when generating with gdaldem -hillshade

I'm getting a strange result when I'm generating an hillshade with the gdaldem hillshade command. Here is my command (see the compute_edges parameter): ...
5
votes
1answer
709 views

Why close a dataset in GDAL Python?

I see quite often in Python GDAL code that people close datasets at the end of their script. Why does it makes sense to close a dataset in Python GDAL? Are there any consequences if I don't do it? ...
1
vote
1answer
268 views

GDAL tiling with a lot of MrSID files

I'm trying to generate tiles for a region with a surface of 85k km2. The data source are MrSID files (2734 files). I'm using gdal to achieve this task. It works fine with a small subset of MrSID, but ...
1
vote
1answer
911 views

Convert multiple GeoTIFF to one KMZ

I am using GDAL to convert GeoTIFF to KMZ. I am using the following command: gdal_translate.exe -of KMLSUPEROVERLAY *.tif out.kmz If I have one TIF file in the folder, it works ok. If I have more ...
1
vote
0answers
74 views

Error while retrieving a WCS coverage with GDAL

I'm currently using the OSGeo-Live 7.0 (Link) image to write a PyWPS process. Part of this process is to retrieve a WCS layer with GDAL (from a GeoServer instance). Since GDAL supports WCS out of the ...
0
votes
1answer
664 views

Using ECW with mapserver

I know that mapserver drop the support for ecw due to licensing issues. But if i understand well, as mapserver rely on GDAL, if GDAL support ECW, mapserver can serve ECW data. So I've compiled ...
4
votes
1answer
149 views

Suggestions for vector format for contour map

I've been trying to create a contour map of a DEM and I've been running into some technical problems. I want to avoid artifacts, so I'm using a DEM of the entire area. This means the size is rather ...
3
votes
2answers
278 views

in python, srs.ExportToProj4() returns an empty string, why?

I have a shapefile with the following .prj file: ...
2
votes
2answers
220 views

How to convert Point coordinates into CSV format?

The GDAL documentation describes how to convert a ShapeFile into CSV data. The layer create options (-lco parameter) allows to prepend the coordinates of Point features to the actual attributes ...
1
vote
2answers
158 views

Difference between OGRGeometry and Geometry in C++ and python

According to the documentation, there is a python version of OGRGeometry's C++ class (http://www.gdal.org/ogr/classOGRGeometry.html) which is http://gdal.org/python/osgeo.ogr.Geometry-class.html ...
6
votes
1answer
262 views

ogr2ogr csv conversion with latlon in one column

I have a csv with "lat\,lon" in a single column. I am trying to convert it with ogr2ogr to another format, and cannot figure out how to format the VRT file to accept this format. I cannot do ...
1
vote
2answers
636 views

OSX 10.9 + QGIS 2.0.1 - Missing osgeo module

Running the following: QGIS 2.0.1 (via Kyngchaos) OSX 10.9 GDAL/OGR 1.10.1 (via Homebrew). Recently, I've become unable to use the OpenLayers plugin, which is complaining that I'm missing the ...
1
vote
1answer
125 views

gdal_rasterize using a “string” as src_datasource

I would like to use a string as src_datasource for gdal_rasterize program. Is it possible to use something like a WKT definition to translate a single point to a raster? Somethink like ...
1
vote
0answers
100 views

Why does GDAL Translate not create any output in batch mode?

Has anyone tried using the Batch Mode function in GDAL translate - qgis2? I can process files individually ok, in my case MODISA HDF files. But the batch mode does not seem to work. Here is my set ...
1
vote
0answers
694 views

Importing Census TIGER Shape Files to SQL Server 2012?

I've been doing my research on how to import shape files into Microsoft SQL Server 2012. and am still unable to figure out how to do it after hours of playing around with it. I'll document what I've ...
2
votes
1answer
634 views

Cutting image using gdalwarp to extent of data

I am currently generating a raster from point data and applying inverse distance algorithm. This image is then cut with : gdalwarp -cutline <<src shape file>> -crop_to_cutline -srcnodata ...
3
votes
2answers
1k views

Help making geotiff and reprojecting to EPSG:4326

My goal is to download a gif file from this webpage: http://aviationweather.gov/data/iffdp/2104.gif Then crop it to the known latlon bounds of 75N,45S,35W,70E, convert to geotif, and then project to ...
0
votes
1answer
428 views

How to run Layer Stack in Batch Mode with GDAL?

My goal is process fifteen Landsat-8 scenes with three bands. I need to rescale the pixel depth from 16 bit to 8 bit, stack three bands and mosaic all them. This job must be making via open source ...
2
votes
1answer
287 views

How to remove strange gdaldem hillshade artifacts?

I've downloaded 9 chunks of 1 arc-second NED in IMG format from the National Map. I'm trying to create a hillshade from them. I'm using: GDAL 1.10.0, released 2013/04/24 First, I combine them and ...
1
vote
0answers
159 views

Uninstall ERDAS ECW/JP2 SDK Version 5.0 (Ubuntu)

I am using Ubuntu Linux as main OS when it comes to geoprocessing but I am not yet an expert using it. I recently needed QGIS to have ECW support, so I followed the instructions of Luke in the ...
2
votes
1answer
311 views

gdalinfo image.tif>stats.txt - returns empty text file (Windows 7 64 bit)

I am trying to output the gdalinfo image statistics results into a text file with the following command: gdalinfo image.tif>stats.txt The stats.txt file gets created, however it is empty. It ...
6
votes
2answers
646 views

Convert points to lines with Python GDAL

I have a dictionary with points. The points are derived from a raster (see image below). The points are not sorted and are not in order. The dictionary looks for example like this: pointDict = {0: ...
2
votes
5answers
984 views

How to export raster from PostGIS enabled DB as GTiff with gdal_translate

When I want to export a GTiff from my PostGIS DB with this code: gdal_translate -of GTiff PG":host='localhost' port:'5432' dbname='landslide' user='postgres' password='pass' schema='public' ...
5
votes
1answer
758 views

GDAL polygonize lines

I have a raster file containing roads (left image). I polygonize the raster file with GDAL (see the script below). In the end I like to have vector lines. However GDAL only gives me something like ...
1
vote
1answer
405 views

GDAL polygonize only certain raster values

I want to polygonize a raster file with GDAL. I only want to rasterize cells with a value of 0. Is it possible to specify that within the polygonize function or do I have to modify my raster before ...
4
votes
2answers
1k views

Issue Trying to create Zonal Statistics using Gdal and Python.

I am trying to create zonal statistics using Python and Gdal. I have a polygon shapefile and a raster file, and in order to do so, I am using a piece of code I found in StackExchange. The raster and ...
5
votes
1answer
558 views

Generate hillshading with GDAL, and respect nodata

I'm still trying to nut out this problem of merging overlapping geotiffs. From the original dataset (here, a 20m resolution digital elevation model of the state of Victoria), I have reprojected it ...
-1
votes
1answer
296 views

How to Geocode use GeoPy API combined with GDAL Python bindings?

My question is related to the 'How to Geocode 300,000 addresses on the fly?' topic that has been asked by Mehul Nahar. I'm interested in running the processes explained by capooti in response to mehul ...
0
votes
1answer
526 views

ogr2ogr DXF to PostgreSQL some LINESTRING coordinates mirrored to the negative longitude

I have been using ogr2ogr GDAL tool to migrate geospatial data from DXF file to PotgresQL database without much problems. But now there is a special problem I can not fix. I have a DXF file with ...
0
votes
1answer
236 views

Why do Shapefiles and GeoJSON behave differently in GDAL Python?

I have a test.shp. I loop over the file, buffer the features, and write the buffers to a new file. I tested writing the buffers to a Shapefile and to a GeoJSON. If I try to use the new buffer layer ...
1
vote
2answers
401 views

Calculate logarithm of a raster?

I'm trying to calculate the logarithm of each cell in a raster. Is there a way to do this with F/OSS GIS? gdal_calc.py only has basic operators and the QGIS raster calculator feature doesn't seem to ...
1
vote
0answers
210 views

Trying to build GDAL's ECW plugin in 64 bit Linux

I'm trying to build the GDAL's ECW plugin in Linux. I'm using the latest version of the SDK from Intergraph's site (5.0). It fails because it tries to compile against the 32 bit version of the ...
2
votes
0answers
353 views

Merge overlapping DEMs with irregular boundaries

I have a 10m DEM covering part of a state, a 20m DEM covering the whole state, and a 90m DEM covering the whole country. How can I merge them together, getting the best quality available in each area? ...
1
vote
1answer
123 views

Load Raster Values into Shapefile [duplicate]

I have a shapefile containing many points. Some of these points fall within an area for which I have raster data. I would like to load the values of the raster which underlie each point into the ...
1
vote
0answers
428 views

Can GDAL Retile run multi-threaded?

gdal_retile is currently using only one CPU of the available eight and it is taking a while (48 hours already) to retile 5TB of data.
1
vote
1answer
196 views

How to get WGS84 lat lon bounds from a VRT file

I have vrt files generated using gdalbuildvrt that reference GeoTiffs referenced using the NZTM coordinate system. Is there an easy way to determine the bounds as WGS84 latitude longitude ...
1
vote
0answers
235 views

GeoDjango with MongoDB

Has anyone had experience successfully deploying GeoDjango with MongoDB? Is it possible? If so, how is it as simple as as installing the non-relational fork of Django and setting MongoDB as the ...
2
votes
1answer
358 views

gdal_translate produces tif with different colors as original

I am using SDE as my raster data source with the ultimate goal to tile this dataset: I am running the following gdal_translate command to generate the initial VRT: gdal_translate -projwin 390605 ...
1
vote
3answers
541 views

KML to Shapefile using GDAL (LIBKML) - MultiGeometry data not transforming

I am trying to convert a KML file to Shapefile using GDAL (ogr2ogr). Regular KML files convert fine but KML files containing MultiGeometry features (eg: Geometry containing both Point and Polygon) do ...