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

3
votes
2answers
148 views

Gdal: clipping a raster with another raster

I'm writing a simple utility to crop batches of multi-band geotiff raster files to the same (smaller) area. Using gdalwarp, I can easily crop a file using a single-polygon clipping shapefile: ...
1
vote
0answers
68 views

Unable to call GDAL functionalities from Python

I am trying to call GDAL functionalities from the following script I made in the past. import os, fnmatch from subprocess import call call("dir", shell=True) # define the variables inFolder= ...
2
votes
0answers
71 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 ...
2
votes
1answer
121 views

Merge tif files with python

I'm also having problems with running gdal_merge.py within another python script and I cannot run it. I already checked a similar problem here. My example is as follows: First I set my workspace and ...
2
votes
1answer
220 views

Orfeo tool box (OTB) in Ubuntu 14.04 and SAGA incompatibility (gdal versions)

I need to test segmentation of a 12000 x 12000 pixels scene. I successfully compiled OTB in Ubuntu 14.04 and ran the process. But when I tried to start SAGA, I got a startup error. I found that gdal ...
1
vote
0answers
62 views

Error running polygonize from QGIS!

I'm trying to vectorize a raster image and when I'm running polygonize within QGIS 2.6.1. I'm getting the following message: Traceback (most recent call last): File ...
2
votes
0answers
112 views

Raster to polygon using R [closed]

I am trying to convert a raster object to a polygon object using R but I am running into some problems using my Windows 7 machine. I originally tried to use rasterToPolygons (package: raster) but it ...
1
vote
2answers
34 views

EPSG Equidistant cylindrical

I need to know what is the epsg code for equidistant cylindrical because i have to georeference an image Aquarius using gdal_translate
3
votes
1answer
100 views

artifacts in SRTM contour layers processed with GDAL tools and rendered with mapnik

Here is my workflow for importing individual SRTM tiles into postgis using the GDAL utils and rendering a png using mapnik: gdal_fillnodata.py -srcnodata -32768 gdal_contour -i 10 -a ele -snodata ...
0
votes
1answer
112 views

Read dted with gdal in c++

I need to find the elevation of lat and long point using a dted map. I read that gdal can read this kind of map. I didn't find a tutorial. I read about gdal here but it's not very "user-friendly". Can ...
2
votes
1answer
67 views

How to convert a RGB tiff into pseudocolor with GDAL etc?

In QGIS I have made a singleband pseudocolor image out of one of the bands from a RGB GeoTiff. I choose one band and then make a linear color interpolation (see attached image). I now want to do ...
0
votes
0answers
35 views

Problem with Leaflet Boundary Canvas plugin

I'm trying to identify the source of a problem that I'm having, I cannot show a tile layer(created with gdal2tiles) inside a polygon generated in Leaflet, using the plug-in ...
1
vote
1answer
83 views

In qgis using python, how can I get multiple pixels from a raster based on a vector point file

I just started working with qgis, and for the life of me I cannot figure out how to solve this problem: I have several points of interest, these are stored in a shape file (just simple points). I ...
0
votes
0answers
159 views

Description of parameters of GDAL SetGeoTransform

Can anyone help me with parameters for SetGeoTransform? I'm creating raster layers with GDAL, but I can't find description of 3rd and 5th parameter for SetGeoTransform. It should be definition of x ...
1
vote
0answers
70 views

Why can't gdal save CRS information to GeoJSON?

gdal fails to add CRS definitions to GeoJSON files in some cases. Try this, for example, from any linux terminal that has gdal installed: # download the data wget ...
1
vote
1answer
51 views

Slanted raster along the line? Wind erosion model and windbreak effectivity

I'm working on wind erosion modelling and windbreak potential. Imagine windbreak, which has some values of optical porosity (in %) and this values determinates the effectivity of windbreak in specific ...
2
votes
1answer
27 views

Does GDAL has any substitute function for the ExtractByMask in arcpy

I am now using arcpy: arcpy.sa.ExtractByMask(weight_raster,elev) both weight_raster and elev are raster data in different size. I want to use functions in GDAL to replace this arcpy methods. Does ...
0
votes
0answers
11 views

How to get linestringe from kml into shp or postgis [duplicate]

Im looking into extracting linestrings from a kml file into a shp file or postgis db. Each linestring would be a feature/row. How would I go about it usinge ogr2ogr or python? I'd like to get ...
1
vote
1answer
55 views

Gdal+Geos point in polygon test fails under MinGW64

I have to build a grid of points and remove the points that falls inside a polygon. Polygons are defined in a ESRI shapefile, points are in a list. This is how I make the test: OGRLayer *layer = ...
2
votes
2answers
164 views

merging raster aerials in QGIS issue

I follow the instructions here very closely and yet every time QGIS gets caught in a loop. It takes a few small rasters totally less than 100mb and creates HUGE gb files and will not stop until ...
1
vote
0answers
77 views

Old GDAL/OGR version with QGIS from ubuntugis in 14.10?

Installing QGIS 2.6 from ubuntugis installs an old GDAL/OGR (ogrinfo --version GDAL 1.10.1, released 2013/08/26) - it should be 1.11 for support of the Open FileGDB connector in QGIS 2.6 So I don't ...
2
votes
0answers
71 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 ...
1
vote
1answer
78 views

How to save changes in Feature-Geometry with GDAL/OGR and Python? (Segmentation fault)

I have a problem with my code: driver = ogr.GetDriverByName('ESRI Shapefile') dataSource = driver.Open('merge.shp', 0) if dataSource is None: print 'Could not open file' sys.exit(1) layer = ...
2
votes
0answers
72 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. ...
0
votes
1answer
27 views

d.rast -o failed after successfully imported.

I imported a geotiff raster using GUI grass, it was successful. (Thu Nov 20 15:55:35 2014) r.in.gdal input=C:\Users\Administrator\grassdata\kavaja\JPEG_8.tif output=JPEG_8 Projection of input dataset ...
0
votes
0answers
49 views

How to create a 8Bit unsigned geotiff with gdal?

I'm trying to create a 8 bit geotiff with this code. Is there a way to use unsigned 8 bit? I only found the GDT_UInt16? import numpy as np from osgeo import ogr, osr, gdal from osgeo.gdalconst ...
0
votes
1answer
59 views

Python buffer in a line using GDAL

I am trying to make a buffer from a line using gdal. I use the example in the gdal Cookbook import ogr, os def createBuffer(inputfn, outputBufferfn, bufferDist): inputds = ogr.Open(inputfn) ...
1
vote
1answer
83 views

How do I convert a TIF to EPSG:3857

I'm trying to convert this TIF to EPSG:3857, but when I run: gdalwarp -of GTiff -s_srs epsg:4326 -t_srs epsg:3857 -co "tfw=yes" original.tif converted.tif I just get a black TIF do you know what I ...
0
votes
5answers
66 views

Aligning Raster Over the Polygon

I have a raster image which have to be aligned over the polygon. Both the raster and polygon have the same coordinate system. The figure shows the raster image which has to be fitted over the RED ...
0
votes
0answers
28 views

Projection of COSMO-EU?

I have data points of COSMO-EU data stored in WGS84 coordinates. Now I try to transform them into coordinate system where the points becomming arranged to a regular grid. Does anybody knows the Proj4- ...
1
vote
1answer
99 views

Using GDAL's GetProjection information to make a coordinate conversion in pyproj

I'm opening a geoTIFF file in Python using gdal. Upon reading with the GetProjection() method, I find the following information ...
0
votes
0answers
64 views

GDAL virtual file system, how to create your own

I just found a blob post about, GDAL virtual file system, and I was wondering if I could create my own and how to do so? I want to create a virtual file system for reading directly from blob storage. ...
1
vote
2answers
127 views

Convert a grid file from one set of coordinate systems to another

I have a grid file in a certain coordinate system. I'd like to change it so I have the same grid in geographic lat/long based on WGS84. This seems like something that gdal can tackle quite easily ...
1
vote
0answers
79 views

Python Band.SetNoDataValue(), sets all bands at once

Can someone tell me what I am doing wrong? I want to set bands 1 and three to 255, and band 2 to 0. The code below only should only set band 2, but gdalinfo shows that all bands get set to 0. Ubuntu ...
2
votes
1answer
68 views

Determine spatial average of several tiff files

I am relatively new to GIS. I have a set of tiff files containing a raster each. I would like to average them, though their extents are slightly different. They each have their geospatial data ...
3
votes
2answers
196 views

Piping data to ogr2ogr

Is it possible to directly pipe data to ogr2ogr? I'm trying to curl a GeoJSON file and pipe it to ogr2ogr to covert it to another format. Without using a pipe, I was able to achieve this by first ...
1
vote
1answer
52 views

Mapping NetCDF CF Mercator parameters to GDAL

I have coordinates on a mercator grid projection from a NetCDF (CF) file. The metadata attributes in the file include : standard_parallel, longitude_of_projection_origin, earth_shape (spherical), ...
3
votes
2answers
237 views

How to create a raster tiff (dem / heightmap) from a mesh of xyz points?

I have to create a terrain heightmap(raster tif preferably) from mesh of xyz ungridded points(in geographical coordinates), I tried taking min and max of bbox of mesh and making width & height by ...
0
votes
1answer
41 views

Merging separate country features together into one region

I am trying to creating a map of the UK where Scotland and England are merged into one shape and not separate shapes for each country. As a minimum working example, I'm following Mike Bostock's Let's ...
1
vote
0answers
64 views

Fiona write crashes

I am trying to write a new shapefile with fiona but after creating the features, i keep getting this output after the program crashes: The program '[6668] python.exe' has exited with code -1073741819 ...
1
vote
1answer
87 views

Unsupervised raster classification using GDAL

Is there an equivalent in GDAL to the Arcpy ISO data unsupervised classification tool, or a series of methods using GDAL/python that can accomplish this?
0
votes
0answers
22 views

Change “World Size” with gdal2tiles.py

When I use the flag "-p 'geodetic'", the output is a world that is 512x256px at zoom level 0. In order for it to overlay with other datasets, I need it to be 640x320px at zoom level 0. I am looking ...
2
votes
1answer
65 views

Which algorithm does gdal_polygonize use?

I would like to know which algorithm (i.e. Sobel operator ) is used to extract the polygons in gdal_polygonizer. I searched here and here but could not find it. Any help would be appreciated.
2
votes
1answer
315 views

Convert GeoPDF with GDAL

I'm attempting to convert a USGS US Topo GeoPDF in Windows. I'm using GDAL 1.11.1. My initial attempts were with this Python script, which automates the process by feeding a command to the Windows ...
1
vote
0answers
40 views

“fill in” .vrt with missing locations in a mosaic

I have ~2000 NASA SRTM .hgt files. I wanted to join them together in a single mosaic, and to save disk space I used GDAL to create a .VRT rather than a very large GeoTiff. This is read by Mapnik. ...
1
vote
0answers
47 views

different data range of raster file in mapserver and qgis

a NetCDF file with only one band with a data range from [0,1.1]. IDV & QGIS identifies it correctly, however gdalinfo finds the data range [1,26505]. Or could it be possible that gdalinfo is right ...
1
vote
3answers
102 views

Lossless conversion from shapefiles to XML/JSON/Text

I'm a software developer and I need to create a custom tool that imports ESRI shapefiles (shp + shx + dbf), manipulates the content and then exports new shapefiles. To avoid a complex binary reading ...
1
vote
1answer
86 views

Remove values below a given value with GDAL for DEM geo tif files

I have a DEM in GeoTiff format. The processing have made areas that have really low values that I would like to remove. Are there any GDAL tools that could set all values below -30 to NODATA or -30.
2
votes
0answers
118 views

Compression in GeoTiff files with GDAL

I have an output from 3th party program in tif format where gdalinfo gives the following information: Driver: GTiff/GeoTIFF Files: J:\ascend\Lufthavne - ...
2
votes
1answer
90 views

Can ogr2ogr write data to a file geodatabase feature dataset?

I am dumping postgis tables into a file geodatabase and it works great--the result is a file geodatabase full of featureclasses. But is it possible to write these featureclasses into feature ...