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
1answer
87 views

getting strange “grid” artifacts when creating hillshades using gdaldem

I'm trying to create hillshades from Czech elevation data (Fundamental Base of Geographic Data of the Czech Republic (ZABAGED®) - altimetry - grid 10x10 m) - a demo files is available here: ...
3
votes
1answer
68 views

Why does projection WKT differ between osr.SpatialReference and gdalinfo

I am writing a simple geotiff using the GDAL python bindings, but the projection WKT reported by gdalinfo differs from the source. The script from osgeo import gdal, osr import numpy # Create ...
3
votes
1answer
80 views

Replicating result of gdalwarp using gdal Python bindings

I am trying to re-project/resample with the GDAL python bindings, but am getting slightly different results compared to those from the command line utility gdalwarp. See update below for shorter ...
0
votes
0answers
90 views

How to import ERDAS ECW in QGIS 2.6 using GDAL 1.11?

I am using Lubuntu 14.10, QGIS 2.6, and Ubuntu GIS(unstable). QGIS was used to use GDAL 1.10, and can import ERDAS ecw format. But after GDAL upgraded to 1.11, QGIS couldn't import ecw images. How to ...
0
votes
0answers
29 views

Subtract parts of polygon that overlap with ocean

I have some seismic data in format of shapefiles/geojson. The data has no relevance in the oceans so I am looking to strip away those parts of the polygon. What I am currently doing is using the ...
0
votes
0answers
59 views

apply gdal to calculate raster that is already in qgis raster calculator

I am using the following code in qgis to calculate some raster.It is not giving the correct result as it is in arcpy.I want to apply same in gdal.Can anyone help me how to use the gdal api.Please give ...
1
vote
1answer
109 views

GDAL GetGeoTransform returning lat/lng instead of affine co-ordinates

I'm using the GDAL libraries via the c# wrapper and have hit a problem when using geotiff's where the co-ordinates are from a state plane zone. Previously I was having trouble getting the correct ...
2
votes
0answers
136 views

qgis 2.8.1 - processing gdal 1.11.2 error (mac, 10.10.2)

i've just recently installed qgis 2.8.1 on my mac. including all relevant packages from kyngchaos - i.e. GDAL Complete 1.11 framework package, Matplotlib Python module, numpy - and updated python to ...
1
vote
0answers
74 views

create image out of geojson (gdal_rasterize?)

Short version: I need to know how I can convert a .geojson-file to an according .png-file (or any other image format) long version: I have a .geojson which was created by gdal_polygonize from a .png ...
0
votes
0answers
31 views

Using GDAL with a line layer and points layer, find next point along line layer

With the points layer and line layer, I'm trying to create a graph for path finding. I've got a line layer that I've dissolved, thinking that might make this easier. The points in the points layer lie ...
0
votes
0answers
104 views

How to read Sentinel-1 tiff images from ESA scihub with Python (and GDAL)?

I can read georeferenced tif images of Sentinel-1 downloaded from PolarView with Python / GDAL but I cannot read the SAFE format from ESA. I need it because PolarView (very nice web interface) keeps ...
0
votes
1answer
63 views

Crop using gdalwarp aspect ratio

I'm using Landsat 8 imagery, and the Indicar Landsat Processing Tools to produce NDVI imagery from the Landsat 8 source. That works fine. I then have a bounding box expressed as lat/lon coordinates, ...
1
vote
1answer
47 views

gdal RasterizeLayer gives different extent to arcpy PolygonToRaster

I have an existing arcpy script which uses PolygonToRaster_conversion to convert a polygon shapefile to raster (tiff). I am now exploring GDAL/OGR and have used the cookbook approach to carry out the ...
3
votes
1answer
146 views

Reprojecting a raster image using gdalwarp

I am trying to reproject an ECW Orthophoto from 25830 to 23030 srid. I tried the following but I am getting an error: C:\Program Files (x86)\FWTools2.4.7>gdalwarp -s_srs 25830 -t_srs 23030 ...
1
vote
1answer
80 views

Getting incorrect coordinates for ortho photo using GDAL libraries

I am very new to gis, and while most of the stuff I have is working well, this thing has me stumped. I'm using GDAL libraries and the c# wrapper to access various gis functions. I'm processing ortho ...
1
vote
1answer
192 views

Calculating NDVI with Rasterio

I'm using Rasterio to read GeoTIFF files from Landsat 8 and calculate NDVI into a new GeoTIFF file. My code looks like this: import numpy import rasterio import subprocess with ...
0
votes
2answers
33 views

gdalwarp: Syntax

In GDAL, I have seen the sample line: gdalwarp $(list_of_tiffs) merged.tiff I would like to merge images to create a mosaic, but can't use gdal_merge (need an exe file), e.g to merge ...
1
vote
2answers
138 views

GDAL build from source throwing OpenSSL errors on Ubuntu

I'm trying to get GDAL installed/up to date on my Ubuntu machine, and I can't seem to get it to work properly. I've been following this process: ...
0
votes
1answer
101 views

convert raster to xyz without no data values

I want to convert a geotif to an xyz-file by gdal_translate. However, when I use the following statement: gdal_translate infile.tif outfile.xyz -of xyz the nodata values are always within the ...
0
votes
0answers
33 views

ERROR 6: The PNG driver does not support update access to existing datasets

I was trying to overlay a shapefile to some raster images from a directory using gdal_rasterize. But then encountered this error ERROR 6: The PNG driver does not support update access to existing ...
0
votes
1answer
81 views

Coordinate Transformation in Gdal/Ogr (EPSG:4326 to EPSG:3068)

I think gdal/ogr(/proj4) doesn't recognize epsg:3068 right? Here is what I tried: public static void transformSRS(){ // creating EPSG:4326 SpatialReference sourceSRS = new ...
0
votes
1answer
84 views

Can't import libraries to rsaga

I get this error when I'm trying to import a raster to RSAGA: > rsaga.import.gdal(paste(rast_name, ".tif", sep=""),"RASTER", env= my_env) Module 'GDAL: Import Raster' not found in SAGA library ...
1
vote
1answer
69 views

ArcMap: Unknown Spatial Reference fix with GDAL/Python

I have a .tif file that I created with Python and when I open it in ArcMap it says it doesn't have spatial reference information. Is there a way to set the spatial reference with GDAL if I have a .prj ...
0
votes
1answer
39 views

connect gdal-java with tomcat

I am trying to run: SpatialReference poSourceSRS = new SpatialReference(); poSourceSRS.ImportFromEPSG(4326); String unit = poSourceSRS.GetAttrValue("UNIT"); System.out.println("unit: " + unit); on ...
0
votes
1answer
40 views

Cut a hole into a geotiff

I have two hillshading geotifs that I created using gdaldem hillshade. One covers a huge area with an OK resolution and the other has a much higher resolution for a smaller part within the huge area. ...
1
vote
0answers
39 views

QGIS Polygonize (GDAL) problem - Null value in attribute

Is it possible to use Polygonize in QGIS (GDAL) so that a line shapefile is converted to a polygon shapefile with an attribute(s) in the line file carried into the polygon file. I am trying to build ...
2
votes
0answers
135 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)". ...
0
votes
1answer
48 views

change projection in an ascii (x,y,z) file

I am having an ASCII file (lat,long,depth) in the Mars projection. I need this file to be converted into WGS-84 projection. How can I do the conversion for the ASCII files.
2
votes
2answers
89 views

Install GDAL frameworks and QGIS — redundant?

Just wondering if GDAL tools are actually already within QGIS -- was it redundant for me to install them on my OS X Yosemite independently? Recently I used GDAL tools to explode .kml data for use ...
0
votes
0answers
45 views

writeOGR alters MultiPolygon holes

I'having trouble using writeOGR to write GeoJSON files, when there are MultiPolygon objects. Here is an example: I'm creating a SpatialPolygonsDataFrame object from GeoJSON, with a single feature made ...
1
vote
0answers
54 views

python-mapscript on OSX Yosemite

I am keen to try out the mapserver export plugin from QGIS 2.6 This needs python-mapscript which comes from the mapserver package. This can be downloaded from ...
0
votes
3answers
94 views

Extract raster values gdal looping rasters

I have this code, most of them getting from this forum. I need to loop over a folder and sub-folders, or just to run a list (putting previously all the raster inside, I don't mind one way or the ...
1
vote
2answers
72 views

How do I access GeoTransform array from gdal on the command line?

I'm using the gdal library via OSGeo for windows, and I was just wondering how to access the geotransform array for a .ecw file.
0
votes
1answer
90 views

How to convert SPOT NDVI hdf file to tiff file using gdal_translate command

I have SPOT NDVI images. I want to convert NDVI hdf to tiff. I have used the following code to get the information of hdf file: gdalinfo 20060501_NDVI.HDF I get the following output Driver: ...
0
votes
1answer
153 views

ogr2ogr unable to open KML datasource

I downloaded some Wikimapia data using the wikimapia.api. Everything is A-OK since I've opened this in Google Earth and everything is in place. However, I need to turn it into a shapefile, or get the ...
2
votes
1answer
50 views

Georeferencing uneven/irregularly gridded rasters

I'm trying to take a set of NetCDFs and turn them into georeferenced datasets using GDAL. However, the only Geospatial data is a set of 2-D lat/lon arrays. The projection type is not lat/lon. THe ...
0
votes
1answer
26 views

gdaldem hillshade stop working with new projection?

Parallel to gdalwarp `-te -5.8 41 10 51.5 -ts 1980 0` not working with mercator? , my hillshade works with lat/long input : gdaldem hillshade input.tmp.tif shadedrelief.tmp.tif -s 111120 -z 5 -az ...
0
votes
0answers
60 views

Create tiles for on Google Maps from PNG images

I have some PNG images in 2048, 4096, 8192 and 16384 pixel size to generate tiles for zoom level 3,4,5 and 6. My data are plotted in Mercator projection and centered over a particular area of clat ...
1
vote
1answer
32 views

gdalwarp `-te -5.8 41 10 51.5 -ts 1980 0` not working with mercator?

I was working with native lat/log projection on etopo1: #donwload: curl -L -C - ...
2
votes
1answer
131 views

Why processing.runalg(“gdalogr:rasterize”…) is not creating out files when running in the QGIS 2.6.0 console?

I want to use the QGIS python scripting functions to turn to rasterize a point vector using many different dimensions. Right now my code looks like this: import processing RasterPixelSize = ...
1
vote
0answers
38 views

How can I improve image quality when using GDAL to generate Google Earth superoverlays?

I'm using the free version of Google Earth, so in order to create superoverlays of my large TIF images I have to use additional software. Google suggests GDAL. Using the instructions Google ...
0
votes
0answers
87 views

gdal ReadAsArray returning NoneType

I'm trying to read in several rasters as array using gdal in a python script, but I keep getting NoneType returned and can't figure out why. The process is this: 1. Use gdalbuildvrt ...
1
vote
2answers
71 views

Using GDAL to get vertices of a multipolygon

I am successfully using GDAL to get the vertices of a polygon from a ShapeFile using .NET code. This link showed me how to do that: ...
1
vote
0answers
61 views

netCDF GDAL C# - time:unit?

I'm trying discover how to read the time:unit variable information information from a netCDF file and have run out of ideas. Without this info i cannot auto-magically process the files by date (which ...
0
votes
1answer
67 views

How to use gdal_grid with VRT?

Hi folk's it seems to be that I'm to stupid to use gdal_grid. I' have file called dgm-10-epsg-32633.csv with equidistant samples 10x10m in this form. X,Y,Z 386700,6068670.000,.000 ...
2
votes
1answer
91 views

Difference between osgeo gdal/ogr python module and GDAL C++ library

I am developing application in python and I use the osgeo gdal/ogr python bindings for gdal, but its not clear to me what exactly is the difference between the gdal/ogr in python and the gdal library ...
1
vote
1answer
48 views

MODIS resolution missmatch between MRT and GDAL processing

I've been using MRT (Modis Reprojectin tool) to extract EVI bands into a geotiff with geographic coordinates (epsg4326). I did not specify an output resolution to MRT. The output from MRT had a ...
1
vote
1answer
319 views

convert vector (point shape) to raster (tif) using python gdal lib in qgis

I'm trying to convert a vectorlayer (point- shapefile) to a raster-dataset (tiff) using python gdal-library in qgis. On Graphical UserInterface I simple choose the menue ...
0
votes
1answer
44 views

Write single lon/lat coordinates in TIFF with GDAL

Is there a way to write the georeferentiation value by value in a tiff file using GDAL?
1
vote
1answer
36 views

Installing GDAL with HDF4 driver on Debian system

I am working on a Debian 7 Wheezy system using Anaconda and working from Spyder. I've had a hard time installing a version of GDAL that reads HDF4 drivers. Could someone walk me through the ...