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
108 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
179 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
174 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
257 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
279 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
58 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
492 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
38 views

Unable to call GDAL functionalities from Python

I am trying to call GDAL functionalities from the following script I made in the past. import os, fnmatch from subprocess import call call("dir", shell=True) # define the variables inFolder= ...
1
vote
0answers
34 views

Error running polygonize from QGIS!

I'm trying to vectorize a raster image and when I'm running polygonize within QGIS 2.6.1. I'm getting the following message: Traceback (most recent call last): File ...
1
vote
0answers
31 views

In qgis using python, how can I get multiple pixels from a raster based on a vector point file

I just started working with qgis, and for the life of me I cannot figure out how to solve this problem: I have several points of interest, these are stored in a shape file (just simple points). I ...
1
vote
0answers
18 views

Python Band.SetNoDataValue(), sets all bands at once

Can someone tell me what I am doing wrong? I want to set bands 1 and three to 255, and band 2 to 0. The code below only should only set band 2, but gdalinfo shows that all bands get set to 0. Ubuntu ...
1
vote
0answers
28 views

“fill in” .vrt with missing locations in a mosaic

I have ~2000 NASA SRTM .hgt files. I wanted to join them together in a single mosaic, and to save disk space I used GDAL to create a .VRT rather than a very large GeoTiff. This is read by Mapnik. ...
1
vote
0answers
29 views

different data range of raster file in mapserver and qgis

a NetCDF file with only one band with a data range from [0,1.1]. IDV & QGIS identifies it correctly, however gdalinfo finds the data range [1,26505]. Or could it be possible that gdalinfo is right ...
1
vote
0answers
43 views

Extract polygons of PNG with GdalOgrInJava

I have PNG pictures with polygons of different colours on them. What I plan to do is to extract these Polygons using the polygonize function (see here) Then I will try to encode the polygons in ...
1
vote
0answers
40 views

How to add CADRG RPF A.TOC Raster Data Source in Geoserver 2.6

I'm currently trying to add CADRG as a datasource in Geoserver 2.6. I've installed the necessary GDAL plugins (1.9.2) for Geoserver 2.6. I placed my RPF folder in C:\Program Files\Apache Software ...
1
vote
0answers
22 views

Python 2.6 Module for GDAL 1.9.2, ImportError

I work on an internal network(IE I cant get the latest 1.11 stuff) and have GDAL 1.9.2 setup on a Red Hat Linux Enterprise Server 6. My linux box has python 2.6.6 however the python binding I received ...
1
vote
0answers
40 views

Reprojecting and Clipping in one script using Python Subprocess module and GDAL

I am trying to reproject a .tif file then clip it to the boundaries of a shapefile using GDAL Warp. I want to do this in a single script, and have successfully used the Python subprocess module to ...
1
vote
0answers
22 views

SQLite Driver available with R's rgdal:writeOGR?

Does anybody know how I can get OGR's SQLite Driver to work with R's rgdal package? When listing the drivers in R with the ogrDrivers() command the SQLITE driver is missing..
1
vote
0answers
63 views

Using gdal ReadAsArray or Arcpy.RasterToNumPyArray on python 2.6 32bit without it crashing

I am trying to run ReadAsArray on an entire raster and on a single band, python will crash without an error message when trying to do either. Is there a way I can create a 3d array without having to ...
1
vote
0answers
37 views

How to run cmd prompt in the background and show progress bar

I am warping files using gdal. Since I have numerous files to warp how will I run the cmd in background and just have a progress bar to indicate the percentage in warping the files? Because when I ...
1
vote
0answers
68 views

NASA SRTM: precipitous contour lines in TileMill

I was trying to make a nice looking print map of Bali (Indonesia) in TileMill. After creating a terrain layer (using this walkthrough) I noticed that the island's coastline was precipitous. Is ...
1
vote
0answers
64 views

Batch GDAL_Calc multiple band calculation

I have a question concerning gdal_calc. I´m running windows 7 64b and am working from the OSGeo4W shell. I have multispectral images that need to be prepared for analysis but they are not saved in a ...
1
vote
0answers
46 views

Use Gdal2tiles to generate 512px tiles

I'm a newbie in gdal2tiles & python but I would like to modify gdal2tiles to create 512px tiles for retina display (screen). Can someone suggest how I can do this?
1
vote
0answers
29 views

Sending commands into FGDB

I have a large .gdb of all blocks in the United States, available from Census. I would like to create a point shape file of all census tract centroids in the US. I have QGIS with FileGDB enabled, as ...
1
vote
0answers
63 views

GDAL WMS in QGIS 2.4

I'm trying to add in QGIS a layer shown here : propluvia When I watch in my console I see this request corresponding to the layer I need. The server is not able to respond to a getCapabilities ...
1
vote
0answers
45 views

Converting MDB File With Unicode Characters

I'm using C# wrapper of GDAL and I'm having a problem converting MDB file that contains unicode characters. Lets take for example conversion to KML file, non-English characters are replaced by '?' ...
1
vote
0answers
39 views

Rendering a geotiff block with gdal and gdk-pixbuf - Some flaws in the result

I currently try to render a geotiff file with GTKmm\GDKmm (gdk-pixbuf). The file is a tiled geotiff with a block size of 256*256 and three rasterbands. The shortened gdalinfo output is: Driver: ...
1
vote
0answers
113 views

Correct way to extract pixel data using GDAL. Current method gives repeated but changing value.

I am trying to open a spectral file and move bands of data from GDAL to OpenCV's Mat data structures. As of now I am just trying to make sure that I can correctly access each band using GDAL and pull ...
1
vote
0answers
66 views

ogr2ogr “Error: pointer being freed was not allocated” on mac mavericks

I've successfully installed GDAL via homebrew on mac But while I'm using command ogr2ogr, keep getting error below $ ogr2ogr -f geoJSON tl_2010_20_tract10.json tl_2010_20_tract10.shp ...
1
vote
0answers
65 views

How to read .adf (TIN) map files?

I have a set of maps in few formats : TIN ├── [4.0K] info │   └── [2.6K] arc.dir ├── [ 708] log ├── [374K] M-33-45-D-d-4.asc └── [4.0K] m-33-45-d-d-4tin ├── [ 96] tdenv.adf ├── ...
1
vote
0answers
41 views

Clip Aerial Image Pair by Overlapping Area?

I have a set of aerial images which I am trying to clip by their overlapping areas, for use in a more involved program I am writing in Python with Python Imaging Library (PIL). Is there a command ...
1
vote
0answers
98 views

Rasterize shapefile using another raster's size and resultion

Rasterizing a vector (shp) I face the following problem: I want my output file (the rasterized vector) to have the same resolution and dimensions as an existing raster. If I understand the interface ...
1
vote
0answers
115 views

MySQL+GDal+MapServer+OpenLayers WebMercator problem

I have a problem/question what is the most efficiant way to go from WGS84 (epsg:4326) to webmercator. Outline of the structure right now: MySQL Database: with point Latitude/Lontiude Data [~200 ...
1
vote
0answers
72 views

How can you calculate the union/merge two polygons in ogr2ogr?

I want to calculate the union/merge two overlapping polygons using ogr2ogr (as I have 300+ to do). I want the output to be the same as you would get if you used the 'Merge shapefiles to one' command ...
1
vote
0answers
50 views

GDAL Color distortion around nodata value stitching

I am trying to extract JPEGs from SID images using GDAL. I ran into a an issue where black lines were present and figured out that it was a nodata issue, so I added -src/dst/vrt nodata arguments and ...
1
vote
0answers
110 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
39 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
38 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
0answers
89 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
111 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
121 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
49 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
121 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
158 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
29 views

ogr2ogr giving “ERROR 1: Failed connect to 192.168.4.18:8181; Operation timed out” error?

This is my command, which was working OK yesterday: ogr2ogr -f "GeoJSON" ../geojson/myfile.geojson ../data/gml/myfile.gml -s_srs EPSG:27700 -t_srs EPSG:4326 The output is this: ERROR 1: Failed ...
1
vote
0answers
89 views

How can I draw features on GeoTIFF in Python?

I have a set of data (coords + value) and want to turn it into a heatmap in GeoTIFF. Should I calculate pixels and set them as array, or there are shortcuts? This is what I do now: driver = ...
1
vote
0answers
31 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 ...
1
vote
0answers
78 views

Extra row and column in a raster after GDAL grid import

I'm exporting an ESRI grid to TIFF format and after that I'm reading this exported raster file into SAGA GIS. However when Saga reads this file in there appears to be an additional row and column in ...
1
vote
0answers
78 views

Write raster with its attribute table and symbology

I am classifying lots of rasters and I'd like to automate the legend creation and symbology displaying for easy manipulation. Here is the kind of raster I'm dealing with: library(raster) r <- ...
1
vote
0answers
65 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 ...