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

6
votes
0answers
290 views

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

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 ...
6
votes
0answers
984 views

How to create a 4096x4096 png from a 1440x720 data set using GDAL?

I'm having a bear of a time trying to do something that seems fairly straightforward. I'm trying to take a netcdf file from ftp://eclipse.ncdc.noaa.gov/pub/OI-daily-v2/NetCDF/2011/AVHRR-AMSR/ and turn ...
4
votes
0answers
394 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 ...
4
votes
0answers
308 views

How to read Geotiff tags and pixel per pixel values using GDAL C#

I'm new here and I'm new at programing with GDAL C#. My question is how I can get the value of SRTM GEOTIFF pixel to compare with it's neighbors, and how can I get the values of Geotiff tags.
4
votes
0answers
263 views

How to store Large Raster Images of Equipotential Surface Generation Using Monte Carlo Simulation?

I have a large (1GB+) raster that I am using as the base for the generation of an equipotential surface. This is to assess the accuracy of DEMs derived from another source. In situ observation is ...
3
votes
0answers
175 views

Reprojecting WGS 1984 Web Mercator (EPSG:3857) in Python with GDAL

I am reprojecting rasters in python using GDAL. I need to project several tiffs from geographic WGS 84 coordinates to WGS 1984 Web Mercator (Auxiliary Sphere), in order to use them later in Openlayers ...
3
votes
0answers
229 views

Fortifying re-projected shapefile for ggplot2 in R

I have a shapefile of US county boundaries which I got from the NHGIS. I can load the file using the rgdal package and see information about it: library(rgdal) shp_1860 <- readOGR("shp", ...
3
votes
0answers
77 views

Can GDAL Perform Vertical Transformations?

Can GDAL perform vertical transformations? If the answer is yes, please indicate which amongst these listed pages can explain how exactly to accomplish vertical transformations? If there is some ...
3
votes
0answers
78 views

Why are my raster mean and standard deviation -1?

I have a set of about 25 rasters whose descriptive statistics I have been examining. The data in the rasters were calculated in Python and the rasters were created using Python-GDAL. In all but one ...
3
votes
0answers
236 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 ...
3
votes
0answers
176 views

How to create raster attribute table to a raster data using gdal python?

I classified a remote sensing data and saved it as a raster using Gdal python. I want to create a raster attribute table for this classified result.. my classified raster has five classes with ...
3
votes
0answers
181 views

Build Raster Attribute Table with GDAL Utilities

Is there a way to build a raster attribute table using one of the GDAL utilities? I know you can view it using gdalinfo and I have successfully created one using Python, but it would be so much ...
3
votes
0answers
175 views

Bug in gdalwarp's simple nearest neighbor? Or more likely I need some help

I'm trying to understand how to use the GDAL utilities a little better. I think I've defined something incorrectly or I've found a bug in gdalwarp. My bet is that I've defined something incorrectly, ...
3
votes
0answers
507 views

How to create mask bands in GDAL and Python?

I have small program in Python that uses gdal and numpy to perform a logarithmic scale of some data into the RGB 0-255 range. Now this input data has a nodata value set to mask the valid values, and ...
3
votes
0answers
321 views

How to compile GDAL with HDF4? (for dummies)

I need GDAL with both Shapefile and HDF4-Support, if possible in 64 bit. Since the normal GDAL-versions from http://www.gisinternals.com/sdk/ have Shapefile, but no HDF4 and the version from FWTools ...
3
votes
0answers
313 views

How to save table from PostgreSQL to MapInfo File?

I have two tables: 1 with geometry in WGS 2 with geometry in Local System I tried this: ogr2ogr -f "MapInfo File" output2.tab PG:"host=localhost user=postgres password=123456 dbname=postgis" -sql ...
2
votes
0answers
92 views

Per-pixel calculations on a raster stack using GDAL

In R, it's relatively trivial to perform per pixel calculations based on a raster stack (e.g., get std.dev for each pixel on a 12 layer GeoTIFF). Unfortunately, the speed is less than desirable when ...
2
votes
0answers
20 views

GeospatialPDF javascript access

Using javascript in GeospatialPDF, is it possible to zoom in to an extent, say when the document is opened the first time ? To what extent can I use javascript to access the geospatial data? pan, ...
2
votes
0answers
57 views

Zonal Statistics using Gdal and Python

I recently came across a loop_zonal_stats function by @ustroetz found here: Issue Trying to create Zonal Statistics using Gdal and Python.. I am using a (single) polygon shapefile and a raster file to ...
2
votes
0answers
217 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 ...
2
votes
0answers
66 views

C# GDAL DEM Speed

My 64-bit, Win 7 C# application using GDAL creates a number of constant-bearing paths with 1000+ points along each path. At each point, I want to look up the 4 closest elevations for interpolation. ...
2
votes
0answers
98 views

Why is GDALPolygonize so much slower than ArcGIS Raster to Polygon?

I am attempting to polygonize a raster using GDALPolygonize() in a Python script. The script began polygonizing yesterday at 5pm and is still polygonizing now at 9:30 am. I have no clue how far ...
2
votes
0answers
40 views

Tabulate Area using GDAL and Numpy

How could I perform a Tabulate Area using GDAL and Numpy? I am trying to summarize how many pixels of each class in a Landcover raster are within each polygon of a Census Blocks shapefile. I've ...
2
votes
0answers
108 views

How to convert a raster (array) to a NetworkX graph and do a * search?

I want to use a slope-raster for a a_star path analysis in NetworkX. I am using Python for this project. My test raster looks like this First I read in the raster with GDAL input_raster = ...
2
votes
0answers
121 views

Buffering a Raster with GDAL and Numpy

I am using GDAL to convert a raster dataset into a numpy.array and then using numpy to buffer the data in the raster. For example, I need to know all the regions within 100ft of a school point ...
2
votes
0answers
85 views

Displaying DEM raster in QGIS

Using Raster>Analysis>DEM to display hillshade produces the message "gdaldem.exe had stopped working" so I asked the community. I'm new to QGIS, and I now realise after a little experimentation that I ...
2
votes
0answers
176 views

GDAL: Geotransform - rotated rasters (and ArcGis)

I am working on a GDAL driver for a proprietary grid-format. In some cases, the grids are rotated in order to get smaller file sizes. In my GDAL geotransform array, I do the following (showing a ...
2
votes
0answers
315 views

How do I solve Gdalwarp error 'too many points failed to transform' for remapping Geostationary to Lambert conformal?

I am trying to remap from Geostationary to Lambert conformal using gdalwarp. My input data is in netcdf, and is in geographic coordinates (degrees) and I would like to output the remapped data to ...
2
votes
0answers
239 views

ogr.GetDriverByName('FileGDB') returns None

I installed: gdal-110-1600-x64-core.msi GDAL-1.10.0.win-amd64-py2.7.msi gdal-110-1600-x64-filegdb.msi from www.gisinternals.com I can open my gdb on the command line: C:\data>ogrinfo ...
2
votes
0answers
296 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 ...
2
votes
0answers
39 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 ...
2
votes
0answers
87 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 ...
2
votes
0answers
99 views

gdal_grid failing to set nodata

Using gdal_grid for spatial interpolation I have run across a problem: the nodata argument is ignored. The only solution I have found for this is to then use gdal_translate but it is rather cumbersome ...
2
votes
0answers
665 views

gdal ImportError in python on Windows

Trying to import gdal for Python 2.7.3 on Windows XP: >>> import gdal Traceback (most recent call last): File "<stdin>", line 1, in <module> File ...
2
votes
0answers
137 views

projection discrepancy with gdal (1.9) and proj4

I boiled the problem down to that python example: import sys try: from osgeo import gdal from osgeo import osr except: import gdal import osr if __name__ == '__main__': version_num = ...
2
votes
0answers
152 views

Localized Linear Histogram Match using python

I would like to use the local linear histogram match for filling the missing pixels of a Landsat 7 image with another older image (SLC Gap-Filled Products). Do you have any idea how can I do it by ...
2
votes
0answers
222 views

Image Pyramid + Geoserver + Google Earth slow with errors

I have hundreds of GB of somewhat small tif files. Here is gdalinfo of one of them: F:\gis\geotiffs>gdalinfo finaloutputSSAM211Jun181306hfwaterfall0.tif Driver: GTiff/GeoTIFF Files: ...
2
votes
0answers
204 views

How do I use scipy.ndimage.gaussian_filter1d on an arbitrary axis?

I have a GDAL raster that looks like this: And I would really like to blur this raster along an arbitrary transect. I see that scipy.ndimage has a gaussian filter that allows me to blur the raster ...
2
votes
0answers
54 views

Tiled image transforming to new profile

i have a set of tile which was tiled profile of mercator and proper to tms scheme mechanism. unfortunately i dont have any tif or georeference table or anything for this. all was tiled with mercator ...
2
votes
0answers
417 views

GDAL tiling DTED (dt2) dataset as gdal2tiles for 3D Maps

i am searching some webgl technologies to create 3d map with digital elevation data. for web tech. i have found ReadyMap SDK to create 3d map from TMS... in ReadyMap for elevation data it uses sth ...
1
vote
0answers
21 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 = ...
1
vote
0answers
27 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 ...
1
vote
0answers
33 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 ...
1
vote
0answers
31 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' ...
1
vote
0answers
42 views

Ohio State Plane to Lat/Lon

I've looked into a few other posts/libraries as far as converting between State Plane (Easting/Northing) to get a Lat/Lon value, however I've run into some trouble. I saw here ...
1
vote
0answers
39 views

How to use GDAL on Heroku?

I have a geospatial python app that I want to use on Heroku. This is my requirements.txt file: Flask==0.10.1 Jinja2==2.7.2 MarkupSafe==0.19 Werkzeug==0.9.4 gunicorn==18.0 itsdangerous==0.23 ...
1
vote
0answers
20 views

Transformation of OGR style data

I have been using GDAL/OGR in a project and it is going quite well. I have been re-projecting data from various source formats into a single CRS (using the C# bindings), but one issue I have come ...
1
vote
0answers
62 views

Broken PostGIS support GDAL

I'm having some problems with Postgres/PostGIS after rebuilding my GDAL library:( I wanted to add MrSID support in GDAL so rebuilt it with the MrSID SDK. Everything compiled ok and I am now able to ...
1
vote
0answers
73 views

Workflow with raster in Qgis 2.2

I'm trying to build a workflow with the Qgis 2.2 processing toolbox. I want to clip a raster layer with a vector layer using the GDAL tool "clip with a layer mask". But it doesn't work. The error ...
1
vote
0answers
19 views

Problem defining the tiling starting point

We are tiling a image datasets. However, we want to specifically define custom windows. To do this, we are using -te option during warping. The warped output is then tiled. The problem though is that ...