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

eliminate small polygons of a PNG with gdal_sieve

I want to eliminate polygons wich consist of less than 3550 px in this PNG . I tried using gdal_sieve but there seems to be an error. gdal_sieve -st 3550 E8.png -of PNG E8sieve3550_4.png ERROR 6: ...
1
vote
3answers
1k views

convert csv to kml or geojson, etc using ogr2ogr or equivalent command line tools (use “latitude” and “longitude” attributes for coordinates)

I've just unsuccessfully tried to use ogr2ogr to convert CSV into KML like so: ogr2ogr -f "KML" output.kml input.csv The csv has "latitude" and "longitude" columns. The attributes are all loaded ...
1
vote
0answers
518 views

How to fix “Unable to load GdalTools plugin. The required ”osgeo [python-gdal]“ module is missing.”

When I start QGIS, it display error message: Unable to load GdalTools plugin. The required "osgeo [python-gdal]" module is missing. Install it and try again. Couldn't load plugin 'processing' from ...
0
votes
0answers
316 views

Set coordinate origin to the top left corner

I have a BMP-picture and I know that it's top left corner has the EPSG:3068 coordinates 35500,24390. I would like to use gdal_polygonize and than calculate the EPSG:3068 coordinates of the produced ...
4
votes
2answers
474 views

How to create a mosaic in QGIS with cutline and feathering for Landsat-8 imagery

1 - I have two scenes Landsat-8 RGB 6/5/4 and I need to mosaic them: 2- This is my cutline in polygon format that will smooth the overlapping area: 3 - So, I need to create a mosaic with ...
2
votes
0answers
84 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 ...
0
votes
1answer
104 views

I can't configure authentication with Google Maps Engine driver from GDAL

I'm using last GDAL from trunk, and trying to read from Google Maps Engine (GME) vía GME driver I already exported GME_API_KEY: My API Key, for the app I registered GOA2_CLIENT_ID: Generated with ...
0
votes
1answer
204 views

gdal_polygonize ignores some polygons

I want to polygonize this PNG. I try gdal_polygonize pic.png -f "GeoJSON" polygonized.geojson but I get only 11 polygons instead of more than 20. What is the cause or how can I fix it?
0
votes
1answer
835 views

gdal_rasterize error: “attempt to create 0x0 dataset is illegal..”

I'm trying to use gdal_rasterize to convert a polygon shapefile into a raster, where the value comes from the 'COVTYPE' attribute. I try calling gdal rasterize like so (I'm calling from within a ...
2
votes
1answer
386 views

Install QGIS on Ubuntu with GDAL built from source?

I am trying to install QGIS on Ubuntu 14.04. While I can easily install it via the package manager or apt-get, I find that QGIS, when installed, does not have Python support as expected. The symptoms ...
5
votes
1answer
376 views

merging GeoTiffs with GDAL causes grid-like artifacts

I'm processing ASTER DEM tiles, merging them into larger blocks (to minimize visible tile edges) and importing the contours into a postgis database. Workflow: gdal_merge.py -of GTiff -o merged.tif ...
5
votes
1answer
923 views

Determine if shapefile and raster overlap in python

I am building a script in python using osgeo/gdal. I have a set of shapefiles and a set of GeoTiff raster files. I would like to have my script ignore shapefiles if they do not intersect with the ...
0
votes
0answers
173 views

Linear interpolation on raster data, save result (geoserver)

I am trying to perform a linear interpolation on a coverage layer (Source data is in EHdr Coverage format, file extension: .flt). Ultimately, I want to use the interpolated layer as input to a WPS ...
5
votes
2answers
344 views

GDAL (or other open-source) equivalent of Arc's Point To Raster Tool

I'm trying to convert a python script that relies on arcpy into one that only relies on open source modules, but I'm having some trouble. The bulk of my processing lies in one Arc tool "Point to ...
1
vote
1answer
63 views

How to use gdal's “forceToLineString()” with Python

I'm working with Python and GDAL and I'm quite new to this ^^ For my script I'm looking for a method to convert a MultiLineString into a simple LineString. Example: MULTILINESTRING ((0 0,0 1),(0 1,0 ...
3
votes
2answers
2k 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
236 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
128 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 ...
4
votes
1answer
472 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
729 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
152 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
443 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
87 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
4
votes
1answer
200 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
447 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
234 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
1answer
103 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
187 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 ...
1
vote
0answers
276 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
110 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
34 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
14 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
91 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
325 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
105 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
142 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
235 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
103 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
44 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
246 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 ...
1
vote
1answer
273 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
223 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
132 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 ...
1
vote
1answer
396 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
211 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. ...
3
votes
2answers
354 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
551 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
155 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 ...
5
votes
2answers
496 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
128 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), ...