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
318 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
114 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
401 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
140 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
148 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
871 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
323 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
691 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
230 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
215 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
560 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
225 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
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
235 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
318 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
464 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
659 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
45 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
34 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
35 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

Qgis gdalogr:contour is not generating output

I am trying to use gdal:contour script in Qgis console and then script. Here is code which I tried write to console. raster = QgsRasterLayer('C:/Users/Bulva.user/Desktop/test_raster.asc', "DEM") ...
1
vote
0answers
18 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
21 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
29 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 ...
1
vote
0answers
62 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
40 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
89 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
34 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
68 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 ...
1
vote
0answers
31 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 ...
1
vote
0answers
35 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
78 views

mosaic rasters with gdal_merge.py

I want to mosaic 2 rasters "rasa" and "rasb" using gdal_merge.py . While running the code, it produces an error " No input files selected". Can somebody help me to find the error in the code? ...
1
vote
0answers
63 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
43 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
37 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
17 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
49 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 ...
1
vote
0answers
51 views

Storing non-orthorectified images in GeoTIFF

I'm reading GeoTIFF and GDAL documentations and I'm trying to understand what type of transformations are supported. GeoTIFF documentation says you can use ModelTransformationTag to define a ...
1
vote
0answers
18 views

Local Maxima (or Pour Point) Pixel Value extract in Raster

I am doing watershed & pourpoint using python GDAL. I have calculated Fill, Flow-Direction & Flow-Accumulation using Python & GDAL. The next is to get Pour Point in certain window ...
1
vote
0answers
53 views

error using gdal_translate with -projwin option

I have to crop geotiff image, I'm using gdal assembly and I'm trying to use the gdal_translate method If I crop my image using the -srcwin option all is fine gdal_translate -srcwin 0 0 400 400 ...
1
vote
0answers
29 views

Flipped (Y-axis) Raster from point cloud

I have been working on a project where I need to create rasters from irregularly spaced point clouds. I am resampling my point clouds to a regular 0.25 cm grid using a nearest neighbor method ...
1
vote
0answers
42 views

How to correctly rasterize a polygon which crosses 180E longitude?

I want to rasterize a polygon wich crosses 180E lon with gdal_rasterize. The polygon is in the WGS84 projection. I used the following command: gdal_rasterize -burn 1 -of GTiff -ot Byte -ts 1500 ...
1
vote
0answers
26 views

How delete feature from layer using Gdal binding java

I am trying delete one feature from layer using GDAL in java https://trac.osgeo.org/gdal/wiki/GdalOgrInJava But nothing happens. ogr.RegisterAll(); Driver driver = ogr.GetDriverByName("ESRI ...
1
vote
0answers
20 views

gdal_retile — force full tiles?

I have a set of large geotifs (DEM + orthoimagery) that I need to generate image pyramids from. So far gdal_retile from osgeo4w is working very well for my use case with one exception: I need all my ...
1
vote
0answers
51 views

undefined reference to `TIFFDefaultStripSize@LIBTIFF_4.0'

I am trying for the first time to use gdal in Ubuntu 15.10 with gcc. I call gdal.h in my program. At first I got the error "undefined reference to `GDALOpen'", so I check online and I found that I ...
1
vote
0answers
29 views

Creating tiling scheme to get multiple extents before gridding elevation data

Looking to create 1 degree by 1 degree tiles with pixel overlap to get extents before I start gridding. So basically I would like to have a shapefile of the tiling scheme with extents before having ...
1
vote
0answers
40 views

Generating continuous aspect surfaces from DEM with QGIS

I'm trying to extract aspect surfaces from different DEM sources at different resolutions in QGIS. I've had success working with SRTM-derived grids (3-arcseconds) and am trying to perform the same ...