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

0
votes
0answers
19 views

How to fix GDAL_Merge problem?

I've got a lot of modis images (MOD15A2 productos -- LAI) separated by years. Inside each year folder there are the folders with the days i.e 049 folder for DOY 049. What I've done so far is to ...
1
vote
1answer
51 views

Can't write decimal data into shape (shp) attribute table with Python

I am trying to create a shapefile with Python GDAL/OGR library. I add to it a real number field, but when I write into it doesn't appear as it should but rounded. I will explain better with the code: ...
1
vote
0answers
12 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 ...
0
votes
1answer
33 views

Analyze height data using a digital elevation map

i ran the following code on python 2.7: import sys, struct from osgeo import gdal from osgeo import gdalconst minLat = -48 maxLat = -33 minLong = 165 maxLong = 179 dataset = ...
1
vote
1answer
18 views

Keep “no data” tiles as “no data” w/ proximity (raster distance) operation

Working in QGIS 2.8.1 Wien. I have a raster file with three classes: "No data" for tiles outside of the country of interest's bounds "0" for tiles that correspond to roads "1" for all other tiles. ...
1
vote
1answer
41 views

Convert to lat-lon coordinates GDAL

I have a netCDF file and would like to convert it to LAT-LON coordinates. 1) I do not know the projection in the file but from the website of the data provider, the data is on a polar-stereographic ...
0
votes
0answers
18 views

Project NetCDF file using GDAL or QGis

Can I project a NetCDF file with the help of GDAL or QGis? The file info is: gdalinfo c:\tomcat-6.0.36\webapps\ncWMS\layers\RAIN_MM_1hour_ac.nc Driver: netCDF/Network Common Data Format Files: ...
-1
votes
0answers
14 views

How to use longitude and latitude coordinates to make a mobile network coverage map using mapbox [duplicate]

I am new to mapping and currently I am trying to make a map showing coverage of a network provider here in my country. The map I am trying to achieve should look similar to that at ...
4
votes
3answers
18 views

gdalwarp to geographic coordinate system

I'm getting the folllowing Error when trying to re-project a raster from NAD_1983_Albers to GCS_North_American_1983. I used the proj4 defenition from going to Properties in QGIS of the layer whose ...
1
vote
0answers
33 views

GDAL: gdalwarp wrong transformation to “South_Pole_Lambert_Azimuthal_Equal_Area” EPSG:102020

I am trying to reproject a raster that is defined in EPSG:4326 to EPSG:102020. Command used gdalwarp -ot Float32 -s_srs EPSG:4326 -t_srs EPSG:102020 -r near -of GTiff -co COMPRESS=DEFLATE -co ...
2
votes
0answers
33 views

How to make a network /tower coverage map in in mapbox

I am new to mapping and currently I am trying to make a map showing coverage of a network provider here in my country. The map I am trying to achieve should look similiar to that at ...
2
votes
0answers
50 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 ...
2
votes
2answers
41 views

GDAL: *.nc rasters, opened as “netCDF/Network Common Data Format” under Linux and as “ HDF5/Hierarchical Data Format Release 5” under Windows

I have a set of rasters that are in the "*.nc" format and I noticed that gdalinfo (and subsequently other GIS packages based on GDAL like QGIS) is opening them with two different drivers. Linux: ...
1
vote
1answer
19 views

How to extract a 5 bands geotiff file into RGB file?

I have a five bands geotiff file for Red, Green, Blue, Red-edge and Infra-red, and want to use gdal_translate to extract bands Red, Green and Blue gdal_translate -b 1 -b 2 -b 3 -mask "none" ...
1
vote
1answer
60 views

ogr2ogr Shapefiles to Postgres SRS problem

I am trying to import to my Postgresql (9.1+129ubuntu1) + PostGIS (1.5.3-2) database a Shapefile with polygons with the ogr2ogr tool from Gdal-Bin package (1.7.3-6ubuntu3). I know the origin SRS of ...
0
votes
0answers
26 views

Error clipping extent of shapefile with boundary box with ogr2ogr from python [closed]

Trying to clip polygon shapefile with boundary box in ogr2ogr called from python gives error: "FAILURE: Unable to open datasource 'C:\UserData\save.shp' with the following drivers." then lists ...
1
vote
0answers
30 views

Why is GDAL/QGIS creating Geospatial PDF's with erroneous spatial data?

I have a GeoTiff that I wish to convert into a geospatial PDF for use in Avenza. I attempted to use gdal_translate, and also used the QGIS's Raster->Conversion->Translate which apparently just builds ...
2
votes
1answer
49 views

Tearing of polygons using ggmap and readOGR

I'm trying to incorporate a shape file of US states with ggmap and keep getting "tearing" of my polygons. I tried changing group=id to group=group which helped, but didn't seem to resolve the ...
0
votes
0answers
12 views

Copy specfic layer to new datasource using OGR Swig Java bindings

Is it possible to efficiently copy a specific layer from one datasource to a different datasource using GDAL/OGR Java bindings? I know it is possible to create a new layer in the target data source ...
2
votes
0answers
48 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 ...
0
votes
0answers
23 views

gdalwarp error cropping raster

I'm getting an error from gdalwarp using the following code: gdalwarp -cutline input.shp -crop_to_cutline -dstalpha input.tif out.tif Error: *** Error in 'gdalwarp': free(): invalid pointer: ...
2
votes
1answer
58 views

Is it possible to do parallel processing in GDAL and QGIS?

It seems when processing with large raster dataset in QGIS, It causes crash dump, for example in QGIS failing to complete clip or intersect and crash dump in qgis , Is it possible to use gdal parallel ...
3
votes
1answer
73 views

How can I use GDAL Merge with multiple files in one folder in python 2.7?

I have multiple files in one folder which are called, for example: 049.h11v11.hdf, 049.h11v12.hadf, 049.h11v13.hdf, 065.h11v11.hdf, 065.h12v12.hdf, and so on. The first 3 numbers are the julian day of ...
3
votes
1answer
29 views

Rasterize vector file using gdal.RasterizeLayer() without losing attribute information in python

Is it possible to rasterize a vector file with gdal.RasterizeLayer(), saving attribute as Raster value ? Like gdal_rasterize can do with : gdal_rasterize -a attribute -ts rasterXsize rasterYsize -l ...
2
votes
0answers
16 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 ...
1
vote
2answers
31 views

Error Installing GDAL python package on a MAC

I'm trying to install the [GDAL python package][1] on a Mac (OSX Yosemite) with python 2.7.9 and GDAL 1.11.3. When I try $ sudo easy_install GDAL I get a number of errors ending with: error: Setup ...
1
vote
0answers
24 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 ...
3
votes
2answers
38 views

Converting Python GDALDataset to C++ GDALDataset

I'm currently writing a C++ extension to python to speed up a raster data viewer which is running extremely slowly with our large datasets. I've managed to get a small speedup by writing a C++ ...
4
votes
0answers
33 views

How to unproject a UTM ESRI Shapefile over WGS84 extent boundaries

Background: I am working on a collaborative project in Python that uses GDAL/OGR to work with rasters and vectors. There is one function that attempts to project one ESRI Shapefile's coordinate ...
1
vote
2answers
116 views

GDAL Translate error

I am trying to translate a .tif raster to an xyz one in Python The code is this: import sys from osgeo import gdal gdal.AllRegister() driver = gdal.GetDriverByName('Local Disk') dataset = ...
1
vote
0answers
33 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 ...
1
vote
2answers
25 views

gdalinfo for Modis layer “NDVI”

I am trying the following command to get info on the NDVI MODIS layer stored in HDF file. I always receive an error. the command and error as below: gdalinfo ...
0
votes
1answer
28 views

GDAL utility in Bash for loop (Ubuntu) [closed]

In Windows GDAL terminal (in Windows, GDAL has a specific terminal/prompt) I usually do loops and it works like a charm. What I'm trying to do the same in a Shell terminal in Ubuntu with the ...
0
votes
1answer
39 views

Reprojecting and retiling a large set of tiles

I am trying to understand the process to re-project and re-tile a set of raster (GeoTiff) but can't wrap my head around how to use the different GDAL tools to achieve my goal. More concretely, I ...
1
vote
0answers
17 views

Access MEM data from GDAL?

I don't have any answers (yet) from from previous question on GDAL memory raster + OpenLayers, so i'll try to ask an easier question :) Anyone knows how to access GDAL memory with command lines when ...
3
votes
1answer
83 views

Convert KML to GeoJSON with ogr2ogr and remove KML special attributes

I'm trying to convert KML to GeoJSON (or some other formats) using ogr2ogr (look at this question). It's all working fine, except of one problem: I'm always receiving these metadata attributes with ...
1
vote
1answer
37 views

gdalwarp resampling problems

I'm trying to clip and resample a raster to match the dimensions and resolution of another raster. Here's my gdalwarp command: gdalwarp -t_srs EPSG:4326 -te 102.3375307079206 14.6874037407838 ...
-1
votes
1answer
55 views

Receiving GDAL Clipper Error from QGIS? [duplicate]

I realise there has been quite a lot of discussion on this topic and have already read through the suggestions on this forum of how to rectify this problem, but still nothing has worked. I am ...
7
votes
0answers
139 views

gdalwarp producing empty rasters when called from python script, but not via command line

I'm trying to aggregate some raster data (all .tif) up to a coarser resolution (from .05 degree to .25 degree) using gdalwarp in python, but the command is not working. Instead of getting an output ...
0
votes
0answers
25 views

GDAL memory raster + OpenLayers

I'm trying to improve my use of GDAL for creating and sending images to an OpenLayers map, but I have some memory issues. Here is the thing : 1. JS send a query to server with map parameters (bbox, ...
0
votes
0answers
6 views

mapnik.Gdal(file=my_latlon_raster_file_path) returns a datasource with an empty all_features() list

When I get my mapnik layer datasource like that for a raster file: layer = mapnik.Layer(my_raster_layer_name) layer.srs = my_srs layer.datasource = mapnik.Gdal(file=full_path_to_my_raster_file) and ...
1
vote
2answers
82 views

Export/Convert .tif + tfw into GeoTiff with tagged Metadata (coordinates)

What is the correct command line for using warp_management in Python (ArcGIS 10.1), to convert several georeferenced .tif's (or folders) with XML world files to GeoTIFF's with permanently written ...
1
vote
0answers
20 views

Changing Dataset Property

I 'm working on GDAL C++. But I have some trouble to change GDALDataset property. I have converted "jpg" to "bmp" an image. But I want to do read an image(jpg) to gdaldataset then change dataset ...
2
votes
1answer
38 views

Raster interpolation with gdalwarp, and replicating with python API

My problem is similar to this: Replicating result of gdalwarp using gdal Python bindings. I want to resample a raster using bilinear interpolation, which I've done successfully with gdalwarp by adding ...
2
votes
3answers
64 views

ASCII file with latitude, longitude and data to GeoTiff using Python

I have a .asc file with three columns of data: longitude, latitude and a data value. Each latitude and longitude point represents the location of a centriod of a 1/8 degree unprojected grid. I want ...
4
votes
0answers
41 views

OS terrain 50: .vrt from zip files

I tried to follow what was suggested in How to use OS Terrain 50 in QGIS, but I got the following problem: in spite of being generally able to generate .vrt files using gdalbuildvrt from unzipped ...
2
votes
1answer
29 views

Convert *.asc OSGB36 files to WGS84 GeoTIFF using gdal/ogr2ogr

I am trying to obtain elevation readings using gdallocationinfo (as per other datasets I have) from a set of .asc files containing elevation data (plus an associated .prj file), using WGS84, but I ...
1
vote
0answers
40 views

Write stacks of GeoTIFF image to NetCDF by chunk size

I converted GeoTIFF image stacks to NetCDF file with writing data by Index order (layer by layer) (this was a main solution). The problem is, it writes slowly (for 25 images - 1Gb takes 90 seconds). ...
1
vote
1answer
29 views

GRID-file: get max extent for all z-values above given value

I have a GRID-file with a quite large extent. For an automatic workflow I need the maximum extent/bounding-box for all z-values that are over a certain range. For example: Get the maximum extent for ...
1
vote
1answer
20 views

Over writing one large DEM over another in QGIS or GDAL

I want to over write/mosaic one large DEM (around 3gb file, Geotiff) over another. What is the best way to do this? Both cover roughly the same region, one is based on 2m res LiDAR, the other 1m res ...