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

1
vote
0answers
44 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
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 ...
1
vote
0answers
36 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
48 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
45 views

How do I set nodata with custom color for gdal_translate and gdal2tilespy?

I start with a bunch of 1-bit geotiff files that I had to convert to 8-bit RGB tiff files to use with GDAL. Then I tried to set borders to red (content is in black and white so I can't use those) and ...
1
vote
0answers
37 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 ...
1
vote
0answers
60 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
204 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 ...
1
vote
0answers
83 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 ...
1
vote
0answers
100 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 ...
1
vote
0answers
120 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
0answers
91 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
0answers
87 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 ...
1
vote
0answers
80 views

GDAL and Python: How to “pad the extent” or otherwise handle band.ReadRaster() for area falling slightly outside the image bounds?

In a nutshell, I'm trying to create a bare-metal tiling utility using GDAL and Python. As you can see in the graphic below, which shows the perimeter of a raster and the OSM-structured tile index, ...
1
vote
0answers
82 views

Efficient computation of path profile with Numpy and GDAL

I have numpy array with SRTM elevation data and I need calculate path profiles - vectors with elevation data from one point (transmitter) to every points of my map. Here's my current solution: def ...
1
vote
0answers
102 views

GDAL in Java fails on ReadRaster()

To summarize, trying to read a raster keeps giving me an UnsatisfiedLinkError that makes sense, but looking into the compiled GDAL jar, I don't know how to fix. gdal.AllRegister(); TheFile = ...
1
vote
0answers
106 views

How to add new raster store and publish for ECW fromat in GeoServer 2.4-RC2?

I am using GeoServer 2.4 RC2 and Tomcat7 on windows xp OS. I have extracted the gdal-data.zip (2.4.0) and geoserver-2.4-SNAPSHOT-gdal-plugin into Tomcat 7.0\webapps\geoserver\WEB-INF\lib. Now ...
1
vote
0answers
70 views

How to join attributes by location with command line

I would like to joint the attributes of 2 shapefiles by location but with a command line (linux). This feature is easy in QGIS but it is VERY slow. Is there a way of doing this with ogr or pyqgis ? ...
1
vote
0answers
61 views

Contours generated from merged TIFs huge - much bigger than individually

I'm generating contours from four 5x5 degree SRTM DEMs in southeastern Australia. I'm using gdalwarp to reproject the DEM, then gdal_contour to generate contours. First, I generated contours for each ...
1
vote
0answers
159 views

How to generate JPEG tiles using gdal2tiles.py

I am using gdal2tiles.py to generate tiles and packing it into MBTiles using mb-util. gdal2tiles by default generates PNG tiles. When I convert PNG tiles to JPEG large black area around tiles are left ...
1
vote
0answers
408 views

How do I change from lat/long to a metric coordinate system in GRASS?

Yes, it's yet another transformation question. I am working with ASTER DEM data, which is LatLong WGS84. I want to process the DEM map into an STL file for use in a CFD simulation. The script I am ...
1
vote
0answers
339 views

How to properly convert EPSG:3857 to EPSG:4326 tif?

I tried using the answer from "Transforming EPSG:3857 to EPSG:4326" to do this but gdalinfo tells me that this didn't really work properly. Driver: GTiff/GeoTIFF Files: geo-map.tif Size is 1024, 512 ...
1
vote
0answers
122 views

GeoPDF Support in GDAL

I installed OSGeo4W using a Windows installer... it contains GDAL 1.7.3. I can't get GDAL to recognize PDF as a file type using commands in MSYS shell. Does this version of GDAL not support PDFs? ...
1
vote
0answers
60 views

How can I maintain MapInfo GRC file attributes in QGIS/GDAL?

I am trying to manage a classified MapInfo grid file in QGIS - I'm wondering if anyone has any experience with maintaining the attributes of the classified grid? QGIS seems to assign them just a ...
1
vote
0answers
239 views

Topography: How to convert raster into vector layers?

(nb:This project is part of the #Wikipedia #wikimaps project) 1. Question: How to create vector topographic maps from raster image ? How to convert raster topographic image into closed vector ...
1
vote
0answers
164 views

QGis Save Raster as Rendered Image

In QGIS 1.9.0 Master, when you right click on a raster in the Layers Panel and select "Save As", you can select the output mode to be "Raw Data" or "Renedered Image". When selecting the Rendered ...
1
vote
0answers
187 views

gdal2tiles.py problem - creating mini tiles

I have a similar problem to the issue posted here (gdal2tiles problem). I have 80 GeoTIFFs (4000x4000) with EPSG:27700 projection embedded and TFW and PRJ files. These I have built into a virtual ...
1
vote
0answers
136 views

problem with WFS Data, GDAL, Python, ArcGIS, & Defining Projection

I'm pulling in data in gml2 format with a EPSG#4267 from a wfs. I create a datastream in python and then run the following from the python window in ArcMap: data = datastream.read() #write the gml ...
1
vote
0answers
119 views

QGIS 1.8.0 on Linux treats NetCDF 4 as HDF5

I've been generating some large (28800x13920) grids of BYTE data representing the NALCMS2005 North American Land Cover dataset (original 250m projected data resampled to a WGS84 grid at 15 arc sec). I ...
1
vote
0answers
115 views

How to create multi-image single band single file geotiff?

I have a geotiff file that is single band but contains multiple subdatasets. I use gdal with numpy to open a specific subdataset and mask some of its values. Example: $ gdalinfo sample.tif ... ...
1
vote
0answers
259 views

Completely new to GDAL: please help me with the best way to study it

This is the first time that I post here so pardon me for asking such a beginner question. I'm completely new to GDAL and I'm going to use it to work with satellite images, DEMs, etc. I have some ...
1
vote
0answers
75 views

gdalwarp not clipping neatline properly

I have a PDF file (10mb) which contains a NEATLINE polygon in metadata. $ gdalinfo I5601_GeoPDF.pdf Driver: PDF/Geospatial PDF Files: I5601_GeoPDF.pdf Size is 5906, 3543 Coordinate System is: ...
1
vote
0answers
124 views

How can I combine these two commands? (Get a clipped raster ASC-II file from PostGIS)?

My objective is to clip a polygon from PostGIS raster table, and then save this raster (draw polygon shape) as an ASC-II in a folder. What I did now are 2 steps. 1st Step: Insert the clipped raster ...
1
vote
0answers
30 views

GDAL Plugin with a different version of GDAL aready built into QGIS

My qgis windows version is 1.8.0 Lisboa and use the gdaltools in version1.9.2. I want to use the plugin NetCDFBrowser, but this plugin need the gdal 1.10. I install 1.10 stand-alone gdal find here ...
1
vote
0answers
71 views

Python GDAL/OGR: check lines geometries connectivity

I'm working with shapefiles representing roads, so basically it's stored as (poly)lines. I need to perform a road connectivity checking. In other words, I want to ensure that roads meant to ...
1
vote
0answers
109 views

setup Mapserver WMS from tile files

I'm totally new to Mapserver and I've been banging my head trying to figure out how to mount a WMS server that will seed images of a small geographical area to a OpenLayers HTML page. I already ...
1
vote
0answers
301 views

Using gdal to convert from UTM to equidistant lat/lon

I am getting rather confused with converting a UTM gridded dataset into a equidistant/rectangular latlon gridded dataset. Can anybody help me? I have a UTM set with the following projection: ...
1
vote
0answers
105 views

How can I call gdal_translate from Java code?

How can I execute gdal_translate in Java? For example, will something like this work? gdal_translate "PG:host=localhost dbname='xyzy' schema=public' table='xyz'" "d:\xyx\ab.tiff"
1
vote
0answers
139 views

using GDAL/OGR (faster)

Eventually, I was able to get the GDAL installed and compiled correctly using Visual C++ 2010. I even tested some of its methods on some geo-referenced raster files, and it looks like that all is ...
1
vote
0answers
85 views

Normalize with gdal_translate

I do not understand the command line when I try to work with gdal_translate in QGIS. I found this at the forum: Normalize raster values to 0-255 And tried this: *"If not, from within QGIS choose ...
1
vote
0answers
141 views

How can I fix badly specified geographic coordinate system with units in arc-seconds?

Check out this depth-to-rock soil data. It is linked from here. The spatial reference system is a latlon, but the units are arcseconds. I could not reason out how to override the SRS even after ...
1
vote
0answers
153 views

Question about geotiffs

My colleague and I are georeferencing a bunch of .tif images. After each image, we'll click "update georeferencing". The resultant .tif files are different between our computers and I'm trying to ...
1
vote
0answers
107 views

Are there limitations with gdal virtual raster option?

I am new at Postgis and I'm groping... I pretend to rasterize some weighty vector layers in Postgis, do some map algebra and reclassification with them and visualize intermediate as well as final ...
1
vote
0answers
83 views

Coarse heat map placement using gdal

I have generated a coarse heatmap of where data is available by 1 degree x 1 degree cells that I want to overlay on top of a world image. Using a python script I merged two images (heat map and a ...
1
vote
0answers
47 views

How do I link an existing install of GDAL to FileGDB in Linux?

The distribution I use (CentOS) only has GDAL 1.8 available through elgis. I've downloaded the FileGDB API from ESRI and I'd like to be able to link it after an install. How can I do this?
1
vote
0answers
145 views

How to use GDAL to burn a selected value from one raster into another?

I think this should be a relatively simple process, but I'm not quite sure how to approach it. I have a shapefile containing 1000's of polygons. I also have 19 rasters. All of the polygons fall ...
1
vote
0answers
94 views

OSGB36 referenced PNG to Google Earth

I am creating a PNG that is referenced to the OSGB36 datum (1 pixel per 50m easting/northing that is colored according to the result of some calculation code) I want to display this on Google Earth ...
1
vote
0answers
270 views

OGR-Driver GML axis order handling (GDAL, QGIS)

By default the OGR-Driver will swap the axis order to longitude, latitude for geographic coordinate systems. "Since OGR 1.8.0, the GML driver has coordinate system support. This is only reported ...
1
vote
0answers
67 views

Size of output geotiff file not as in the .vrt file

I have the following vrt file (input.vrt) with GCPs in EPSG:3857 <VRTDataset rasterXSize="5409" rasterYSize="5408"> <Metadata /> <Metadata domain="IMAGE_STRUCTURE"> <MDI ...
1
vote
0answers
95 views

VSIMalloc3() and UnifiedSrcDensity mask errors when running gdal_warp in OSGEO4w shell

This is related to gdal.py as part of bat script in windows 7 I am getting ERROR 2: VSIMalloc3(): Out of memory allocating 854522136 bytes. ERROR 2: Out of memory allocating 119813400 bytes ...