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

2
votes
0answers
141 views

How to fix “Raster layer: Provider is not valid” error when loading PostGIS rasters?

I compiled and installed the current QGIS master as usual and I imported some raster data into my postgres 9.3 database with postgis 2.0.2 extension. Whenever I try to import the raster data into qgis ...
0
votes
1answer
110 views

POSTGIS raster/Java with GDAL bindings - Load byte array in Gdal Dataset Object

I would like to load a raster as java.awt.image.BufferedImage (or equivalent class that would allow to perform some display) out of my PostGIS database. My Raster is inDb (loaded with raster2pgsql) ...
0
votes
1answer
40 views

How to resolve a 'Received NULL pointer' ArgumentOutOfRangeException for OSGeo.OGR.Geometry.TransformTo?

I am attempting to project a geometry using the GDAL c# bindings and I'm getting a ArgumentOutOfRangeException ('Received NULL pointer') exception on the Transform methods of Geometry. I'm also ...
0
votes
0answers
74 views

MODIS products reprojection problem

I am working with Modis Data (MOD02, MOD11, MOD03, MOD35). In particular I am looking for a solution to reproject the subdatasets below, into albers equal. Because I am using NLCD dataset which uses ...
2
votes
1answer
123 views

How to preserve original resolution of scanned 1-bit PDF raster with gdal_translate?

I am using gdal_translate to convert non-spatial PDFs to non-spatial TIFFs (intended for subsequent manual georeferencing). This is my gdal_translate command: gdal_translate -of GTiff -sds -co ...
0
votes
1answer
78 views

gdalwarp increase resolution

I currently have a 87Mb Compressed GeoTiff with a resolution of 50m. I need to interpolate this to 2.5m I have used gdalwarp to achieve this using a cubicspline, but it took 12hrs and resulted in a ...
1
vote
2answers
112 views

Unsuccessfully trying to iteratively run gdal_polygonize in a batch file [closed]

Thanks in advance for your help with this. I've been reading, experimenting, and banging my head against the wall for a couple days on this, and could really use some assistance. This is my first ...
0
votes
0answers
140 views

GDAL HDF not recognized as supported file

I am having the same issues as in this thread; however, I'm new to programming and need more advice on how to get gdal to open hdf files. I downloaded hdf from MODIS, which I believe are hdf 4. I have ...
1
vote
1answer
124 views

GDAL - How to get gemetry type of a MapInfo TAB file?

I am using the OGR MapInfo File driver to read a MapInfo TAB file. However, the function GetGeomType() returns 0 which means 'wkbUnknown` even if the features in the TAB are only multi-polygons. Is ...
0
votes
0answers
54 views

need help creating layer with geom type wkbGeometryCollection with gdal python

I am having trouble creating a layer with geom type as geometry collection. Thus I am unable to output to a shapefile. I have attached code below. dstLayer is null. This does not happen if geom type ...
2
votes
2answers
233 views

Error exporting data from Oracle to shp using ogr2ogr

I have an oracle table which I would like to query/make a selection. I would then like to export the selection to a shp file. I'm using gdal and ogr2ogr. Following the guidelines here and here I ...
0
votes
1answer
72 views

gdalbuildvrt error when mosaicking SRTM data with different NoData values

I'm attempting to mosaic together 9 SRTM tiles downloaded from CGIAR to produce a DEM of Borneo. I'm using: gdalbuildvrt Borneo.vrt BorneoSRTM/*.tif where BorneoSRTM contains the 9 tif files to be ...
0
votes
0answers
44 views

Import GRD and count days with temperature below zero, using python

Minimum temperature of a year is given in GRD file. I have to read that file and then calculate the frost day (less than zero degrees celcius) from it in python.
3
votes
3answers
89 views

QGIS Mosaicing with null values in images

I have a set of landuse landcover data which needs to be mosaiced. Whenever i try using QGIS Raster>miscellenous>merge tool. giving no datavalue as '0', it returns either (A) ' Traceback (most ...
0
votes
0answers
83 views

aggregate counts on a grid and determine grid adjacency

So here is the problem I am working on. I have a bunch of lat, long coordinates for event occurrences. I want to create a grid and then count how many points fall into each grid pixel. Once I have ...
3
votes
2answers
182 views

How to add Georeferenced PDF as layer to QGIS 2.0?

I have some georeferenced pdf's that I would like to open in Qgis. I have read that GDAL version 1.8 and higher will support pdf, and I have GDAL 1.10 installed with Qgis 2.0 on windows 7, But when I ...
0
votes
1answer
148 views

create polygon adjacency matrix using python

I am trying to create an adjacency matrix from a set of polygons. That is, I have a bunch of polygons and I want to identify which polygons have a common edge or "touch" each other. Once I find this ...
1
vote
0answers
87 views

Problem with polygon overlay in raster stats using gdal

I am using the raster_stats module in python to calculate the mean carbon value for different sets of polygons. Here is the code that I am using: africa_sat = ...
0
votes
0answers
82 views

Convert Shp to KML error org2org Error4

I am attempting to convert Shp file to Kml using command line org2org application. I have used this in the past to convert more simple shp files and was sucessful. However I get an error message when ...
3
votes
1answer
150 views

Are the functions in OGRPolygon accessible through the python bindings?

I'm currently trying to use the OGR python bindings to process some shapefiles. I want to obtain the number of interior rings in a polygon feature. Unfortunately I'm unable to use the documented ...
0
votes
0answers
58 views

Problems rendering tilestache served mbtiles in QGIS

I came across Geoffery Hing's walk through on using Tilestache to serve up mbtiles to QGIS. I followed the configuration and I was able to view tiles in the browser without any problems. However, when ...
0
votes
2answers
382 views

using Python to calculate NDVI, output as a Geotiff

I need a script that can calculate NDVI from two separate input Geotiff files, then output the results as a Geotiff. I wrote following script, most of which is borrowed from various on-line sources ...
1
vote
1answer
85 views

3D functions with SFCGAL in PostGIS

I have this setup: "POSTGIS="2.1.2 r12389" GEOS="3.3.8-CAPI-1.7.8" PROJ="Rel. 4.8.0, 6 March 2012" GDAL="GDAL 1.9.2, released 2012/10/08" LIBXML="2.7.8" LIBJSON="UNKNOWN" TOPOLOGY RASTER" on ...
3
votes
3answers
429 views

Python - Find the largest extent in a list of extents (tuple format)

I am trying to find a way (simplest is the best) to retrieve the largest extent in a list of tuple extents in python. The list is looking like that: extents = [(-180.0, -90.0, 180.0, 83.624), ...
1
vote
2answers
141 views

PostGIS to GeoJSON with GDAL: how to get lat/longs?

I'm trying to convert line features from PostGIS to GeoJSON files, as follows: from psycopg2 import * from subprocess import call conn = connect(dbname='gis') cur=conn.cursor() cur.execute("SELECT ...
3
votes
1answer
315 views

Re-project raster image from Mercator to Equirectangular

I have a raster image of the world map in Mercator projection, square 1:1 format. I want transform the image into an Equirectangular projection accurately. I assumed from my research that GDAL would ...
2
votes
1answer
100 views

How can I extract only band 4 from a geotiff multi bands to create geotiff or AAIgrid with this band using Python?

I have a geotiff file with many bands... I will like to create a .asc file or tiff file by converting only one band from the geotiff file. I don't find the solution by using python... With ...
1
vote
0answers
37 views

How do I split a DEM at a certain elevation?

I have a DEM that includes bathymetry and terrain elevation and I would like to split it at the coastline. Ideally, I would end up with one tif that's land and one tif that's bathymetry. I have ...
0
votes
2answers
126 views

Get srid from a shapefile using ogr

I'm writing a script to import shp to postgis, using python, osgeo and psycopg2 ( see my previous question: Import shp to Postgis using Python and ogr ) I'd know if there is a way to get by osgeo ...
0
votes
0answers
102 views

Make a raster with pixel scale 100 by 100 meter with gdal_rasterize

I'm asking here, because, I've got great results last time. Here it is the next of my adventure with gdal_rasterize, I want to make a raster picture where each pixel correspond to a square of 100 by ...
2
votes
3answers
196 views

Merge 30000 polygon shapefiiles in a logical manner in ArcGIS model Builder

I have 30000 shapefiles (a1.shp to a30000.shp). They are all polygons. Each polygons has one and only one polygon feature in it. I need to Merge the polygons into 1000 shapefiles. My output would be ...
1
vote
1answer
36 views

How can I specify the ranges in the gdaldem's color configuration file?

I want to render an IMG file according to a specific renderer, and export to a jpeg File. I have asked a question about this, and got a satisfactory answer. The issue now is that I need to write a ...
1
vote
0answers
31 views

using gdal2tiles 2 tile a png longlat. keep getting wrong coordinate box in google kmls

I have tried lots of variations, but can't seem to get gdal2tiles to do what it used to do. Prior version? images are 60N to 60S and 180W to 180E 9000 pix x 3000 lines. gdal2tiles keeps outputting ...
1
vote
1answer
108 views

How to call ogr2ogr from command line with Python 2.7 from ArcGIS 10.2 for Desktop?

I am trying to access ogr2ogr in the command line, and this is the message I receive: 'ogr2ogr' is not recognized as an internal or external command, operable program or batch file. This is ...
2
votes
2answers
277 views

Conversion of point grid to a raster with open source tools (Gdal, QGIS …)

I'm looking to do something which seem simple to me, but I haven't found a solution since 2 days I looked for. I've got a Postgis DB with lot of points (WGS84) spaced of 100 meters which cover a ...
0
votes
0answers
51 views

How to rasterize a point layer? Cell value should be 1 if feature attribute is 1

I’m a fairly experienced GIS user, but I’m new to QGIS. I’m using version 2.2, and I’m trying to convert a point data layer into a raster data layer using the Rasterize tool. The point data layer has ...
1
vote
0answers
90 views

Decompress a LZW-compressed GeoTIFF

I have a LZW-compressed GeoTIFF file and would like to split it into several smaller regional files and change the format to png. I tried to do the spliting with gdal_translate -srcwin. I only get a ...
0
votes
1answer
25 views

how to calculate the average temperature of each band of a .grd file?

I am having a GRD file in which the temperature at different longitude and latitude are given in a 32x35 array for one day(band),i am having 365 such bands and want to calculate the average of all ...
3
votes
2answers
136 views

Memory and Processing Efficient Approach to Read a Timeseries in Python

I am trying to store a timeseries of high resolution raster images containing NDVI values into a dictionary for calculations. I am curious as to which approach may the most efficient when loading ...
1
vote
1answer
191 views

Split Large GeoTIFF to smaller regions with python and gdal

I have quite a large GeoTIFF-file (compressed 50MB but unpacked 6GB) and want to split it into several smaller regional files. Therefor I am using python. As I am a new gdal user I do not really have ...
0
votes
1answer
169 views

gdal polygonize in python error, blank polygon created

I am having trouble using the Polygonize function in python. The cookbook example for this can be found here. The relevant portion of my code is: sourceRaster = gdal.Open('myraster.tif') band = ...
0
votes
1answer
83 views

How can I reproject a layer in pyqgis?

I have a problem. I'm currently trying to get the extent of a vector layer in latitude and longitude coordinates in python. The layer is currently in EPSG:3857, so .fullExtent() is not returning what ...
0
votes
1answer
137 views

Error with Gdal 1.10 install on Ubuntu 13.10

I've installed GDAL 1.9 then upgrade to 1.10 on Ubuntu 13.10 and then now it's error (gdalwarp can work but gdal_translate is fail) gdal_translate ERROR 1: gdal_translate was compiled against GDAL ...
0
votes
1answer
133 views

Gdal creating raster rotated wrong way

I have been trying to read in a stack of images for a timeseries and created an output raster that contains the max pixel values of the timeseries. The arrays were created using numpy img_stack = ...
1
vote
1answer
241 views

QGIS DEM (Terrain models) Mode Slope creates wrong data [duplicate]

I have a DEM raster file downloaded from the USGS site. The file can be downloaded here. I am trying to create a slope raster in percent with the DEM tool in QGIS. The GDAL command is the following: ...
2
votes
3answers
124 views

How to apply a specific Rendering to an Image and export to JPEG?

We have a specific process which generate an interpolated single band 32-bit Float Raster in Erdas Imagine (.IMG) format. This process generates an image every hour. We wish to apply a specific ...
1
vote
0answers
72 views

GDAL C/C++ API and filtered pyramid tile generation

I have been tasked with adding functionality to an existing C++ developed toolset whereby a source dataset (typically a single GeoTIFF) needs to be split into tiles at various resolutions ('zoom' ...
4
votes
2answers
149 views

Converting shapefile to raster: Character attributes

I've got a shapefile that I would like to convert into a raster, but the trick is that the attribute I want a raster version of is character data. I'm familiar with geoprocessing using GDAL (in the ...
0
votes
1answer
144 views

GDal problem “no suitable image found” when I run “ogr2ogr” command

I am trying to run "ogr2ogr" command macbook$ sudo ogr2ogr -f "PostgreSQL" PG:"dbname=test user=postgres" /Users/macbook/planning.gdb But I am getting this error: dyld: Library not loaded: ...
0
votes
0answers
18 views

GDal problem: “Unable to open datasource” when I run “ogr2ogr” command [duplicate]

I am trying to load gdb file to postgis, I created a database "test", then run a command using "ogr2ogr". sudo ogr2ogr -f "PostgreSQL" PG:"dbname=test user=postgres" planning.gdb FAILURE: Unable to ...