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
40 views

Using ogr2ogr for clipping vector files

I am looking for alternative (more efficient) methods for clipping large feature class datasets. When using the 'Clip' tool in ArcGIS, often times the tool will crash after many hours of processing. ...
2
votes
0answers
36 views

How to use GDAL RasterizeLayer (python binding) -where flag

I am trying to use GDAL RasterizeLayer python binding, but I can not find a way to use the -where flag. The 'pure' GDAL command I am trying to run is: gdal_rasterize -b 1 -burn 0.00001 -where 'DN=...
2
votes
0answers
63 views

Is there a good example of using GDAL_WMS to add an authenticated TMS layer to QGIS?

I've been trying to add a TMS layer in QGIS using a GDAL_WMS xml configuration file. I added UserPwd to the xml and this seems to provide the authentication required by GeoServer but I'm struggling ...
2
votes
0answers
27 views

failed gdal script from QGIS Georeferencer using Thin Plate Spline (TPS)

I'm trying to georeference a photo using QGIS Georeferencer plugin (this is an oblique photograph. I know I should be using a proper georectification routine, but my current work now is just to simply ...
2
votes
0answers
182 views

gdal/python: tif not recognised a supported file format

I'm trying to open a tif using gdal, read it as a array, do some math, and write it back as a tif. However, when I try to gdal.Open(tif) in python, I get the error: RuntimeError: /att/gpfsfs/...
2
votes
0answers
61 views

Why azimuth result is wrong for this pyqgis code?

In this code I want to obtain azimuth from a single point to each pixel center, but when run it, in raster created azimuth result for each pixel is wrong, what is correct code? from osgeo import gdal,...
2
votes
0answers
66 views

Why this pyqgis code cause crash dump in QGIS?

In ArcGIS there is gp.EucDirection to calculate Euclidean allocation, distance, and direction around a point, Is there equivalent function in pyqgis or gdal? I wrote a script with help of gene and ...
2
votes
0answers
28 views

No Exception Handling in GDAL/OGR java bindings: ogr.UseExceptions() has no effect

When a GDAL/OGR 2.0.1 swig/java operation encounters an error, the message is written to stderr but no exception is thrown within the application -allowing it to continue, which is highly undesirable....
2
votes
0answers
107 views

display of raster in QGIS - cells are not rotated

when viewing rotated rasters in QGIS - the overall raster is rotated (as designed), though each grid cell is not. (the nodata value doesn't appear to pick up as the default in QGIS either - -2.9999e+...
2
votes
0answers
189 views

NetCDF to Shapefile with Python and GDAL?

I need to save my outputs, which are netcdf files or panda python tables, to shapefiles. They cannot be rasters because the software I am loading them into accepts shapefiles only. I have googled ...
2
votes
0answers
61 views

Export File Geodatabase Raster to TIFF or IMG

Is there a license-free technology that can export rasters from an ESRI File Geodatabase to a file-based raster (e.g. GeoTIFF or IMG file)? ESRI has the FileGDB API for "reading" the contents of a ...
2
votes
0answers
151 views

How to import netcdf in qgis

I have been trying to impot netcdf files following this guide online http://www.ggiuliani.ch/download/netcdf_qgis_GG.pdf into qgis 2.10 with no success. ALthough the file is imported, the dimensions ...
2
votes
0answers
70 views

ogr2ogr ST_Difference operation on OGRVRTDataSource

I'm trying to subtract from one shapefile the portions of another which overlap using the following: land.shp: subunits.shp lakes.shp COOKIE=subunits CUTTER=lakes; \ echo '<OGRVRTDataSource>...
2
votes
0answers
290 views

Comparing the use of GDAL to build raster pyramids or overviews versus Arcmap 10.1

I have a large set of imagery that I created a mosaic from in arcmap 10.1. I want to use GDAL to create the overviews because arcmap has issues building the pyramids and never seems to complete. I ...
2
votes
0answers
102 views

Output formats not showing in heatmap plugin (QGis 2.10.1 in ubuntu 14.04)

I just upgraded to QGIS 2.10 (using the debian repos) and the heatmap plugin doesn't show any formats in the dropdown menu. I'm able to open every raster format I have at hand and, of course, the ...
2
votes
0answers
120 views

gdal_merge shifts rasters after processing, layer stack is slightly off

I have experienced that gdal shifts rasters. The outputs of gdal_warp are often skewed, have different pixel sizes and so forth. A known error which has also been discussed here. But now, I had to ...
2
votes
0answers
75 views

Multiple layer intersection

I am trying to find a way with qgis/python to get the intersection area of more than 2 layers. All the tools in Qgis take just the 2 input layers. ='( For example in the image below, I have 4 layers ...
2
votes
0answers
58 views

Is GDAL/OGR (specifically CoordinateTransformation) thread-safe?

I'm using the Java bindings if it matters. I have a heavily multi-threaded program and this is one of the bottlenecks that I have in a synchronization block. If I could remove the synchronization, it ...
2
votes
0answers
302 views

Loading Raster with SharpMap (GDAL) doesnt work

i have problems loading rasterfiles (.tiff) with SharpMap I tried initializing a SharpMap GDAL Layer but nothing happens. It doesnt even throw and error when debugging after i set a path for the file ...
2
votes
0answers
118 views

BlueMarble is appearing strange

I am fairly new to WorldWind and want to make WorldWind Server 1.5 server blue Marble Imagery to worldwind 2.0. I am using WMSLayerManager applet to test my server. Blue marble appears strange on ...
2
votes
0answers
248 views

Zonal computations: computing raster pixel counts within polygons using GDAL versus rasterstats

I'm computing some zonal statistics, and was comparing the results obtained using the GDAL cookbook example code here with the rasterstats zonal_stats function. I want to compute the mean value per ...
2
votes
0answers
127 views

GDAL equivalent to ArcGIS Union?

I have a single shapefile full of overlapping polygons. The ArcGIS Union tool will generate individual, non-overlapping polygons with that geometry. Is there a way to do this with GDAL/OGR or another ...
2
votes
0answers
92 views

Reverse parameters of a projection map

I'm user of a soaring simulator named "Condor Soaring". This simulator is at the center of a community of enthusiastic people who build sceneries (landscapes) which represents quite accurately real ...
2
votes
0answers
102 views

How to re-grid a raster using area-overlap weighted averaging?

I have two rasters r and r2 with different resolutions. I want to re-sample (re-grid) r2 to be the same as r using area-overlap weighted averaging:If an r grid cell is covered by more than 1 ...
2
votes
0answers
49 views

GDAL merge; Why add 0.1?

In the implementation of gdal_merge: https://svn.osgeo.org/gdal/trunk/gdal/swig/python/scripts/gdal_merge.py a value of 0.1 is added to the conversion from geo-referenced coordinates to pixel ...
2
votes
0answers
93 views

3d model to point cloud

Excepting FME, is there a way to convert 3d model to point cloud?
2
votes
0answers
2k views

How do I install the GDAL bindings for C#/.NET 4.5.1?

I'm trying to use the C# gdal bindings but I get an error in my C# code which says "Unable to load DLL 'gdal_wrap': The specified module could not be found. (Exception from HRESULT: 0x8007007E)". I'...
2
votes
0answers
183 views

Using OGR2OGR with LIBKML (Examples please)

I am trying to figure out how to use GDAL (ogr2ogr) to export spatial data from SQLServer to KML. The standard KML driver works, but it is not very robust. LIBKML has so much promise (layers, styling, ...
2
votes
0answers
68 views

How to change values in a raster file at particular positions?

I have a collection of coordinates which comes from a QGSGeometry object. The shape is a polygon. example: qgsGeometry [[(75.2172,17.6566), (75.8499,17.9613), (75.5218,17.2113), (75.2172,17.4574), (...
2
votes
0answers
86 views

Help with Python/GDAL documentation

The instructions from the url http://gdal.org/gdal_tutorial.html give the following example: dataset = gdal.Open( filename, GA_ReadOnly) with the following description: It's interpretation is ...
2
votes
0answers
138 views

QGIS polygonize GDAL DLL load failed

I'm using QGIS 2.6.1 and i've been having problems using the polygonize function. I keep getting a "DLL load failed" import error. I've searched older threads which unfortunately weren't able to help....
2
votes
0answers
148 views

ogr2ogr fileGDB translating points

I am converting data from spatialite to fileGDB (and reverse) for a project using QGIS and ArcGIS. When I convert using ogr2ogr (or QGIS) from spatialite to fileGDB, the points are translated ...
2
votes
0answers
113 views

QGIS Problem with extract and warp modules in model builder

I have problems using the gdal modules extract and warp in the model builder. They run correctly when started without the modelbuilder, but when I integrate them in the modelbuilder they won't work. ...
2
votes
0answers
65 views

What is the issue with geopdf format?

I'd like to add the geopdf format as output of a library that I am working on. After have looking around, I have not been able to find a (c, c++, python) library that meet my license requirements (...
2
votes
0answers
97 views

How to Use GetGeomtype Gdal Python

I am just wondering what I am doing wrong when trying to use the GetGeomType. I need to focus on separating layers from a pix file. Two categories are the bitmaps(which are vectors) and then the ...
2
votes
0answers
163 views

Ignore empty tiles in gdal2tiles.py when generating tiles

I am in the process of generating a large set of tiles from a .vrt file composed of many images. The image is following a road section, being very long and small strip of images and due the shape of ...
2
votes
0answers
140 views

Any solution to mosaic Modis MOD13Q1 using batch file?

I have a Modis MOD13Q1 product for Pakistan, four tiles cover my complete area and I have data consisting of 14 years, after every 16-day interval. I want to mosaic them using a batch file which ...
2
votes
0answers
206 views

Improve gdaldem hillshade output (lighten up flat parts)

I have a hill-shading that I made using gdaldem hillshade which is already pretty good. I used the option -alt 60 to get nice contrast, but the problem is, that with lowering the light-source, also ...
2
votes
0answers
55 views

Error in polygonising raster data with QGIS on Mac

I'm using QGis Chugiak on a Mac and I'm trying to converting a raster layer in a vector data by: Raster/Conversion/Polygonize. So, I select the layer needed from the drop-down menu, I give a name to ...
2
votes
0answers
241 views

Memory leak with GDALOpen?

I'm very new in using GDAL, and trying the tutorial given here. But with a piece of code with only GDALRegister_GTiff(); GDALDataset *poDataset = (GDALDataset *) GDALOpen(argv[1], GA_ReadOnly ); ...
2
votes
0answers
1k 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 ...
2
votes
0answers
616 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 ...
2
votes
0answers
262 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 ...
2
votes
0answers
459 views

Why does python reads only the no data value from GeoTiff file?

I had a number of tiff files with a country's map. They had different rows, columns or different projections etc. I needed to make them with same projection, spatial resolution and same number of rows ...
2
votes
0answers
351 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 ...
2
votes
0answers
114 views

How to preserve GCP's after gdalwarp

I'm trying to figure out how I can preserve my GCP's after a gdalwarp. I need the georeferenced GCP's data, so I can calculate the Root Mean Squared Error. Is there a way to preserve this data, and ...
2
votes
0answers
263 views

Convert a classification image to aggregated polygon in shapefile using gdal/python

My goal is convert a 2D numpy array to shapefile. The matrix results from a classification, so the numbers have a meaning. I solved my problem, but the solution is somehow involuted and complicated. ...
2
votes
0answers
93 views

Append ArcGIS feature layer data using GDAL/OGR API

ArcGIS servers are configured with a maximum records count that they're allowed to return following a REST query. This maximum is usually 1000 records, which means that GDAL/OGR clients should ...
2
votes
0answers
472 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 ...
2
votes
0answers
57 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, ...