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
208 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 ...
2
votes
0answers
499 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
134 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 <- ...
2
votes
0answers
323 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 ...
2
votes
0answers
116 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
405 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 ...
2
votes
0answers
142 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 ? ...
2
votes
0answers
149 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
2k 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 ...
2
votes
0answers
887 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 ...
2
votes
0answers
325 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 ...
2
votes
0answers
1k 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 ...
2
votes
0answers
701 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
237 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"
2
votes
0answers
221 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
571 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 ...
2
votes
0answers
229 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 = int(gdal....
2
votes
0answers
181 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 ...
2
votes
0answers
242 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
320 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
476 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
665 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
25 views

Is there any difference between gdal_proximity.py and ArcGIS Euclidean Distance?

The gdal_proximity.py description says: The gdal_proximity.py script generates a raster proximity map indicating the distance from the center of each pixel to the center of the nearest pixel ...
1
vote
0answers
9 views

how to prepare tiffs to create a netcdf file

I´d like to create a netcdf file containing a mean and stddev for some variables for 12 months. At this moment each geotiff (270MB) is either a mean or stddev of a month. Now my problem is that I ...
1
vote
0answers
20 views

Setting default value in FileGDB table column during GDAL/ogr2ogr export from PostGIS?

I'm doing an export for one table (table name: 'schaden_polygon') from a PostGIS-database (db-name: 'mydata') to ESRI FileGDB via ogr2ogr commandline-tool. The PostGIS-table contains a column 'm_uuid' ...
1
vote
0answers
20 views

Conditional statement ENVI

I am trying to do conditional statement in ENVI to calculate multiband image: result = (b1)^2 - (b2)^2 if (b1 eq 0) or (b2 eq 0) result * 0 else result * (result gt 0.0) Is there a way to combine ...
1
vote
0answers
40 views

Multiple HDF5 files in single GDAL virtual raster

Does anyone know how to generate a single GDAL virtual raster linking multiple HDF5 files to each other? Each file represents different swaths of the globe with different lat/lon coordinates and data ...
1
vote
0answers
16 views

Why can't get correct RGB img of Landsat 8 using np.dstack

I use np.dstack to stack Band4, Band3, and Band2 of Landsat8 data to get ture color image,I use code from Milos Miljkovic, but the result looks incorrect. import numpy as np from matplotlib import ...
1
vote
0answers
49 views

Merging features by admin or sovereingt

I followed Mike Bostock's 'Let's Make a Map' tutorial to make a GeoJSON map of the EU starting from a 1:10M cultural admin0 shapefile from Natural Earth. The problem that I'm facing if that some of ...
1
vote
0answers
25 views

Create Diference TIFF image from two TIFFs

I am trying to use gdal_calc to calculate the differences between two almost identical images. The problem is that the images are tiffs with three bands (RGB) and the result is difficult to interpret. ...
1
vote
0answers
42 views

ogr2ogr not running on Ubuntu 15.10

I'm trying to convert an osm.pbf file to a spatialite db but I can't get ogr2ogr to run. I have qgis 2.14.2 and ubuntu 15.10. I've reinstalled gdal-bin and I'm getting the following message from the ...
1
vote
0answers
54 views

How to Convert from GMT Line files to Shapfile for Qgis?

I have several line text-files of (used in GMT) saved in the format (this is a small sample): -111.29 42.28 -111.29 42.31 -111.29 42.33 > -111.37 41.99 -111.38 42.02 I am trying to convert this ...
1
vote
0answers
42 views

Converting TIFF (DIMAP) to .xyz file format with GDAL

We are attempting to convert a .tif file into .xyz file. We have been using gdal_translate on a linux box and this is the line we've used: gdal_translate -strict -of XYZ ABC.tif ABC.xyz -mo ABC.dim ...
1
vote
0answers
48 views

Is there any way to connect features (entities) of dxf to database in python?

I'm using gdal, ezdxf, sqlite python packages and sqlite ODBC driver. They work fine, and in AutoCAD I can easily manually connect to my dbase.sqlite (read only) using _DBCONFIGURE and _DBCONNECT ...
1
vote
0answers
23 views

GDAL compilation to .NET (C#) use

I want compile GDAL library for use in C# program. I use Microsoft Visual Studio 2013. First, I compiled GDAL C++. As a result I have DLLs and libraries (No problems) . Next, I tried to compile with ...
1
vote
0answers
24 views

How to create tiles from c# code

I need to create tiles directly from c# code. The input file is in *.tif format. Is there a way to do this simply? I tried GDAL, but I found it complicated, I don't understand how to cut tiff into ...
1
vote
0answers
31 views

Python GDAL open dataset in PostGIS with quotations

Error trying to use gdal to open a dataset in PostGIS with the following code: ds = gdal.Open("PG: dbname=nyc host=localhost user=postgres password=mypassword port=5432 mode=2 schema=public column=...
1
vote
0answers
77 views

Slow performance clipping raster with GDAL

I have a few large LiDAR DEMs that I want to merge. Each DEM has some glitches around the borders which need to be removed in order to get a seamless final product. I manually traced around each DEM ...
1
vote
0answers
44 views

ogr2ogr - Shapefile database SQL query; missing .cpg and .xml

I am using gdal and ogr2ogr to perform a SQL query on a database of shapefiles with the code below. The geoprocessing works and the appropriate shapefile is extracted from the query, but when I print ...
1
vote
0answers
97 views

Transformation from local projection to WGS84

I'd like to reproject point that is in local coordinate system to the WGS84, using gdal. But exception "No PROJ.4 translation for source SRS, coordinate transformation initialization has failed" ...
1
vote
0answers
37 views

Efficiently iterating through raster calculations with Rasterio or GDAL

I've got a very slow and round-about way to currently do this, but I struggle with the general concept of using windows in raster calculations and how to apply an affine or geotransformation between ...
1
vote
0answers
73 views

How to make NoData transparent?

I have been using python to turn a netcdf into a geotiff. Once I have the geotiff, I am using gdalwarp to transfer it into the right projection. Then use gdal_translate to get the tif into a vrt file....
1
vote
0answers
39 views

Gdalwarp UTM reprojection causes image distortions

I'm trying to re-project Sentinel-2 images to a common UTM zone in order to create a mosaic. For this I'm using gdalwarp with the following syntax: C:/OSGeo4W64/OSGeo4W.bat; gdalwarp -t_srs \"+proj=...
1
vote
0answers
37 views

gdalwarp won't run with “cubic” resampling method

I am having an issue with running the gdalwarp command from command prompt. Basically I have a raster WGS84 file, which I would like to reproject to spatial reference which is using the Azimuthal ...
1
vote
0answers
69 views

Troubling using MODIS package in R (mac). Error with gdal

I'm new to using MODIS data, and have been trying to set up a a way to use the HDF files produced by MODIS on R. I had trouble with gdal (could only download the package through system dialog), and ...
1
vote
0answers
47 views

How do I calculate the percentage of progress using raster data using Python?

How do I calculate the progress amount of two or more raster data being processed? I have seen examples of vector data but it becomes difficult to calculate the progress when no vector data are ...
1
vote
0answers
47 views

ogr2ogr (or other) simplification of polygons with shared boundaries

I am trying to simplify a collection of 100,000 or so polygons containing nearly 4 million vertices. The data is stored in a 2GB shapefile so there is no explicit representation of the topology ...
1
vote
0answers
21 views

How can I read the TIFFTAG_DATETIME using GDAL in Python

I am currently using GDAL from Python to read and transform coordinates and that works well. I would now like to be able to read some of the none geographic tags, for example, TIFFTAG_DATETIME using ...
1
vote
0answers
61 views

gdal_merge.py is not working with NODATA values

I am trying to merge a couple of GeoTiff files with gdal_merge.py. Example input: In the input data the nodata value is -32678. Hence I tried: gdal_merge.py -n -32678 -o DEM.tif -of GTiff na_DEM....