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
2answers
594 views

Image processing using Python, GDAL and Scikit-Image

I'm new here and I'm sorry if I do something wrong. I am struggling with a processing and hopefully I will be able to solve here. I work with Remote Sensing applied to Forestry, especially working ...
4
votes
1answer
265 views

PostGIS is not seeing any GDAL drivers

I've tried setting POSTGIS_GDAL_ENABLED_DRIVERS=ENABLE_ALL as an environment variable. As well as running SET postgis.gdal_enabled_drivers = 'ENABLE_ALL'; from psql. -- (With a pg_ctl restart after ...
2
votes
1answer
119 views

Implement r.cross / ArcGIS Combine with GDAL

I am using GDAL to process a lot of rasters. I then need to combine all the arrays to create a single output raster that represents all the unique combinations of values. I.E the same thing as r....
4
votes
1answer
387 views

Python GDAL: ogr.open or driver.open?

I’ve been using GDAL with python to manipulate shapefiles. When opening a source file I first set the driver to “ESRI Shapefile” before calling the driver’s "Open" method with the required shapefile. ...
0
votes
1answer
204 views

GDAL2Tiles crash

GDAL2Tiles seems to crash after generating tiles and no google viewers are generated, but only openlayers. This is the code I run: os.system("gdal_translate -of VRT -a_srs EPSG:4326" + ...
2
votes
3answers
2k views

Read HDF4 data using python

I am now working with VIIRS/NPP Active Fires by using python gdal. But I cannot read the data inside the files. Filename = "NPP_AVAF_L2.A2012019.0600.P1_03110.2014057125956.hdf" Here is some ...
2
votes
1answer
214 views

use gdal python bindings from the QGIS folder in a custom script

I have QGIS installed (on Windows) and can import this line from the QGIS Python console: from osgeo import gdal When I try the import in a python file it doesn't work Traceback (most recent call ...
1
vote
0answers
297 views

ECW support for QGIS 2.8 (Wien) and Ubuntu 14.10 (Utopic Unicorn) [duplicate]

Getting QGIS 2.8 (Ubuntu 14.10) ready for ECW files is really a pain. I've no terminal experience and it took me several hours to install the ECW JPEG 2000 SDK 3.3. What I did so far: #download ECW ...
1
vote
2answers
628 views

Shapefile to Mysql error using ogr2ogr

I'm following this tutorial. I installed the “GDAL Complete” Framework from kyngchaos.com and I tested the installation on Mac. I have GDAL 1.11.2, released 2015/02/10. When I am trying to run this ...
2
votes
1answer
385 views

Can gdal2tiles create tiles at zoom levels > 10?

I am trying to create tiles from tif images using the gdal2tiles command (at the OSGeo4W command line). I've been asked if it would be possible to create tiles at zoom level 12. Using the -z switch ...
0
votes
0answers
89 views

GDAL conversion to Terragen outputs file full of zeros?

I recently learned some gdal command line syntax in my attempt to convert some geotiff files to the Terragen format. During my first attempt I noticed a message about nodata values: gdalwarp -of ...
4
votes
0answers
37 views

Cropping EUDEM with GDAL creates zig-zag terrain

I have downloaded EUDEM (20 Gb) from http://www.eea.europa.eu/data-and-maps/data/eu-dem#tab-original-data When I try to crop it by GDAL using: gdalwarp -dstalpha -of EHdr -crop_to_cutline -cutline ...
0
votes
0answers
65 views

Optimizing shapefile lookup time?

So, I have a ~1GB shapefile which is running circles around a ~5MB shapefile when it comes to feature look-up time (using python, with GDAL): def check(name, lat, lon, attributeFilter:'SQL Query'=...
0
votes
1answer
125 views

GDAL - Specify span of elevation values for conversion to Terragen format?

I'm attempting to use the gdal tools to convert some geotiff files I have to the Terragen format. I was hoping to use either gdal_translate or gdalwarp in a similar fashion to this post. (though I ...
0
votes
0answers
594 views

Interpolation using GDAL

I am looking for some hints here. There is a shape file with point data of rainfall sampled at a regular grid, and there is a shape file with point data of soil attributes sampled at a regular grip ...
2
votes
3answers
987 views

GetRasterBand() method for gdal in Python

I was going through a tutorial book called Python Geospatial Development. On the chapter on using working with geospatial data in python there is an example of a script meant to handle and analysis ...
1
vote
0answers
369 views

gdal.jar and .dll placement with gdal.Allregister();

I built gdal and java bindings for gdal. My environment variables I added (Windows 7 64 bit): CLASSPATH = C:\Program Files\Java\JDK1.8.0_40\lib; C:\gdal-1.11.1\swig\java\gdal.jar the dlls (gdal111....
1
vote
1answer
273 views

Create web map tiles from large images

I have 70Gb of raster data from drone photography that I need to serve via Openlayers. The data is currently provided as georeferenced jpegs each is 175Mgb in size. What is the best strategy to cut ...
0
votes
1answer
108 views

polygonize giving duplicated polygons

I am trying to reproject some ascii files and convert them into shapefiles. I am new to python but looking at some code online and help files I was able to write the following code to run on Python ...
2
votes
0answers
94 views

3d model to point cloud

Excepting FME, is there a way to convert 3d model to point cloud?
2
votes
1answer
601 views

dissolve polygons of GeoJSON with gdal/ogr

I have a GeoJSON with only polygons in it. I have to dissolve (merge) the polygons with the same "DN"-value. I found out that you can do it with shapefiles like this: ogr2ogr outputfile.shp inputfile....
1
vote
1answer
110 views

extracting the outline of raster imagery using GDAL

Salams I want some help regarding developing a method using GDAL to extract and represent the footprint of any raster file (or group of raster data) in a vector format... any ideas please ??
0
votes
0answers
129 views

GDAL, Georeferencing

I have been trying to use GDAL library to add geoinformation to an image. From the documentation from GDAL web page I could find that I could use GDAlSetGeotransform(), for that I need a six ...
2
votes
1answer
2k views

Install gdal in python on Windows7

I'm having problems installing the gdal library for python. I use a Windows 7 and have WinPython 2.7.9.2 which runs Python 2.7.9. I'd like to use gdal + Python for raster manipulation. I've read ...
0
votes
0answers
320 views

Reprojecting from Lambert Conformal Conic into EPSG:3857?

I'm trying to reproject a plain png image file into EPSG:3857 for use in an OpenLayers application. I'm pulling the image from the NWS site at: http://rapidrefresh.noaa.gov/HRRR/for_web/hrrr_ncep_jet/...
0
votes
1answer
245 views

Ogr2ogr trying to import and convert shapefile to EPSG:3857

I'm attempting to import a shapefile (along with .prj file) to MySQL using Ogr2ogr and using the t_srs option set to EPSG:3857 (the default for Leaflet and Google Maps) but I'm getting odd results, eg:...
0
votes
0answers
205 views

Gdal2tiles issue when using -z for specifying zoom lvls in raster profile

Everything works fine when not using -z. But using -z and only wanting 0-1, i run into an issue that is best illustrated as below. C:\python\WinPython-64bit-2.7.6.4-gdal\python-2.7.6.amd64>python ...
2
votes
1answer
261 views

Mixing POINT and MULTIPOINT in PostGIS table

I'm trying to load a bunch of GeoJSON files into one PostGIS table. I get this kind of error: Warning 1: Geometry to be inserted is of type Point, whereas the layer geometry type is Multi Point. ...
2
votes
1answer
270 views

getting strange “grid” artifacts when creating hillshades using gdaldem

I'm trying to create hillshades from Czech elevation data (Fundamental Base of Geographic Data of the Czech Republic (ZABAGED®) - altimetry - grid 10x10 m) - a demo files is available here: http://...
3
votes
1answer
282 views

Why does projection WKT differ between osr.SpatialReference and gdalinfo

I am writing a simple geotiff using the GDAL python bindings, but the projection WKT reported by gdalinfo differs from the source. The script from osgeo import gdal, osr import numpy # Create ...
7
votes
1answer
2k 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 ...
1
vote
0answers
122 views

How to import ERDAS ECW in QGIS 2.6 using GDAL 1.11? [duplicate]

I am using Lubuntu 14.10, QGIS 2.6, and Ubuntu GIS(unstable). QGIS was used to use GDAL 1.10, and can import ERDAS ecw format. But after GDAL upgraded to 1.11, QGIS couldn't import ecw images. How to ...
0
votes
0answers
124 views

apply gdal to calculate raster that is already in qgis raster calculator

I am using the following code in qgis to calculate some raster.It is not giving the correct result as it is in arcpy.I want to apply same in gdal.Can anyone help me how to use the gdal api.Please give ...
1
vote
1answer
543 views

GDAL GetGeoTransform returning lat/lng instead of affine co-ordinates

I'm using the GDAL libraries via the c# wrapper and have hit a problem when using geotiff's where the co-ordinates are from a state plane zone. Previously I was having trouble getting the correct ...
3
votes
0answers
245 views

qgis 2.8.1 - processing gdal 1.11.2 error (mac, 10.10.2)

i've just recently installed qgis 2.8.1 on my mac. including all relevant packages from kyngchaos - i.e. GDAL Complete 1.11 framework package, Matplotlib Python module, numpy - and updated python to 2....
1
vote
0answers
470 views

create image out of geojson (gdal_rasterize?)

Short version: I need to know how I can convert a .geojson-file to an according .png-file (or any other image format) long version: I have a .geojson which was created by gdal_polygonize from a .png ....
0
votes
1answer
209 views

Crop using gdalwarp aspect ratio

I'm using Landsat 8 imagery, and the Indicar Landsat Processing Tools to produce NDVI imagery from the Landsat 8 source. That works fine. I then have a bounding box expressed as lat/lon coordinates, ...
1
vote
1answer
133 views

gdal RasterizeLayer gives different extent to arcpy PolygonToRaster

I have an existing arcpy script which uses PolygonToRaster_conversion to convert a polygon shapefile to raster (tiff). I am now exploring GDAL/OGR and have used the cookbook approach to carry out the ...
3
votes
1answer
336 views

Reprojecting a raster image using gdalwarp

I am trying to reproject an ECW Orthophoto from 25830 to 23030 srid. I tried the following but I am getting an error: C:\Program Files (x86)\FWTools2.4.7>gdalwarp -s_srs 25830 -t_srs 23030 C:\...
1
vote
1answer
180 views

Getting incorrect coordinates for ortho photo using GDAL libraries

I am very new to gis, and while most of the stuff I have is working well, this thing has me stumped. I'm using GDAL libraries and the c# wrapper to access various gis functions. I'm processing ortho ...
4
votes
1answer
585 views

Calculating NDVI with Rasterio

I'm using Rasterio to read GeoTIFF files from Landsat 8 and calculate NDVI into a new GeoTIFF file. My code looks like this: import numpy import rasterio import subprocess with rasterio.drivers(...
0
votes
2answers
89 views

gdalwarp: Syntax

In GDAL, I have seen the sample line: gdalwarp $(list_of_tiffs) merged.tiff I would like to merge images to create a mosaic, but can't use gdal_merge (need an exe file), e.g to merge C:\Temp\image1....
1
vote
3answers
679 views

GDAL build from source throwing OpenSSL errors on Ubuntu

I'm trying to get GDAL installed/up to date on my Ubuntu machine, and I can't seem to get it to work properly. I've been following this process: https://milkator.wordpress.com/2014/05/06/set-up-gdal-...
1
vote
1answer
512 views

convert raster to xyz without no data values

I want to convert a geotif to an xyz-file by gdal_translate. However, when I use the following statement: gdal_translate infile.tif outfile.xyz -of xyz the nodata values are always within the xyz-...
0
votes
0answers
102 views

ERROR 6: The PNG driver does not support update access to existing datasets

I was trying to overlay a shapefile to some raster images from a directory using gdal_rasterize. But then encountered this error ERROR 6: The PNG driver does not support update access to existing ...
0
votes
1answer
358 views

Coordinate Transformation in Gdal/Ogr (EPSG:4326 to EPSG:3068)

I think gdal/ogr(/proj4) doesn't recognize epsg:3068 right? Here is what I tried: public static void transformSRS(){ // creating EPSG:4326 SpatialReference sourceSRS = new ...
0
votes
1answer
193 views

Can't import libraries to rsaga

I get this error when I'm trying to import a raster to RSAGA: > rsaga.import.gdal(paste(rast_name, ".tif", sep=""),"RASTER", env= my_env) Module 'GDAL: Import Raster' not found in SAGA library '...
1
vote
1answer
186 views

ArcMap: Unknown Spatial Reference fix with GDAL/Python

I have a .tif file that I created with Python and when I open it in ArcMap it says it doesn't have spatial reference information. Is there a way to set the spatial reference with GDAL if I have a .prj ...
0
votes
1answer
129 views

connect gdal-java with tomcat

I am trying to run: SpatialReference poSourceSRS = new SpatialReference(); poSourceSRS.ImportFromEPSG(4326); String unit = poSourceSRS.GetAttrValue("UNIT"); System.out.println("unit: " + unit); on ...
0
votes
1answer
85 views

Cut a hole into a geotiff

I have two hillshading geotifs that I created using gdaldem hillshade. One covers a huge area with an OK resolution and the other has a much higher resolution for a smaller part within the huge area. ...