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

Join CSV using gdal/ogr

I'm trying to use gdal/ogr to join a shapefile called myshape.shp with a csv called mycsv.csv. The join field in the shape in joinshp and the join field in the csv is joincsv. I've found this command ...
0
votes
1answer
47 views

GDAL: Update spreadsheet from polygons in a shapefile

I am a beginner with GIS data, and have been looking for an answer but have not yet found it. I have a spreadsheet with lat/long points. I also have a shapefile with several numbered polygons over a ...
1
vote
0answers
550 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 ...
0
votes
1answer
102 views

Adding a Vertex To Exisiting Feature(linear) Using OGR Python

I am working on a project that involves finding the intersection of linears in a shapefile, if that vertex(the found intersection) doesn't already exist for each(or one) of the features, then that ...
0
votes
2answers
462 views

gdal2Tiles output differs from MapTiler

I'm in need to process a large batch of Landsat images with MapTiler and since the version I'm using (1.0 Beta 2) does not support multiple input files, I created a simple python script that calls ...
2
votes
1answer
98 views

Invert gdalwrap cutline

I have a big raster image with medium details, but for some regions I do have far more details as a seperate raster image. Now I'd like to display both of them in my map that I render with Tilemill. ...
3
votes
1answer
279 views

Holes between imagery tiles projected using gdalwarp

I am having the following issue when reprojecting imagery tiles using gdalwarp. These are 4096x4096 imagery tiles, which I'm gdalwarp'ing from WGS84 to a State Plane projection. When reprojected, ...
0
votes
0answers
163 views

Polygonize only pixel values of 1 from a raster

I have a raster mask that I have created in QGIS (2.2.0) that consist of 0, 1 and nodata values. I would like to end up with a polygon dataset that consist of only the value of one pixels from this ...
4
votes
1answer
285 views

Projection mismatch when plotting precipitation with ggplot

I'm trying to plot accumulated precipitation in Brazil for a given period, see code below: Load packages library(raincpc) library(SDMTools) library(raster) library(ggplot2) library(rgdal) ...
2
votes
1answer
277 views

Reproject MODIS level 3 chlorophyll data using gdal

I am trying to reproject MODIS Level 3 chlorophyll data (for instance, this file) using gdal in R. I used gdal_translate and gdalwarp functions. There was no crs information associated with ...
1
vote
1answer
1k views

Using OGR Intersection Method Python

I'm currently working to find the intersections between linears in a shapefile. I'm a newbie to the GDAL/OGR library. I've done some googling and came up with what seems to be a useful example(see ...
7
votes
2answers
1k views

Conflict between ArcPy and GDAL/OGR?

TL;DR I can't get an installation free version of GDAL/OGR to work from within ArcMap 10.0. I get an ImportError: DLL load failed on _mod = imp.load_module('_gdal', fp, pathname, description). ...
1
vote
1answer
90 views

How to connect to CouchBase with ogr2ogr?

I am trying to load a Shapefile to CouchBase (version: 2.5.0 DEV edition, installed on Ubuntu), but with no success. According to this documentation it is possible like this: ogr2ogr -f couchdb ...
5
votes
2answers
425 views

Optimize Tilemill export of project with contour-lines (and fix out of memory problems)

I have a project in Tilemill (based on OSM-Bright) that covers Europe and also displays Contour-Lines. My problem is that the performance of Tilemill is pretty bad when displaying the contour lines ...
0
votes
1answer
346 views

Combine overlapping rasters in gdal_calc.py, taking the maximum value from each pixel?

I want to combine 3 overlapping rasters to create 1 raster. Each have the same dimensions, resolution and CRS. I want to populate the new raster with the maximum value from each of the layers. ...
0
votes
0answers
258 views

GDAL Installation help in windows

I'm trying to build a GIS application in Django Framework. I completed the GeoDjango installtion (postresql, spatial libraries). I installed OSgeo4W in my C drive. I also setup environment ...
0
votes
1answer
155 views

gdal2tiles.py generating dark tiles from .img file

Dark tiles are generated when I run gdal2tiles.py on a .IMG file(a scanned and georeferenced topographical map). This question advised to apply a color relief map using gdaldem before generating ...
0
votes
1answer
230 views

problem with feeding input files into gdal_merge.py

While merging the files using gdal_merge.py the following code works: subprocess.call([sys.executable,gmerge,'-o','C:\\r.tif','-of','GTiff','D:\\a.tif','D:\\b.tif'],shell=True) However, when the ...
0
votes
1answer
112 views

gdal_rasterize : getting area not in or touching polygon (options -at and -i)

I would like to capture the area that is not in a polygon (option -i) and that doesn't touch that polygon. gdal_rasterize -at will return the polygon + edges, but gdal_rasterize -at -i will not ...
5
votes
1answer
1k views

Why the result of a merge of multiple raster is so big?

I try to merge 14 geotiff like this : Each geotiff is about 50Mb. I need a geotiff at the output My workflow : gdalbuildvrt -input_file_list list.txt test.vrt (where my list contains the name ...
0
votes
0answers
318 views

I cannot use Gdal Python but Gdal binary, how to fix?

**UPDATE: I can use GDAL with Python now by copy file from Gdal 1.11 compiled folder to Python package folder. But the problem with GDAL binary and GDAL library that they cannot read the HDF5 file. ...
2
votes
0answers
249 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 ...
3
votes
1answer
3k views

How to add a color table to a one band tiff using GDAL?

I have 1-band GeoTIFF file. The color interpretation for band is gray. I want a color table in order to transform my Col.tiff to Col.jpeg or Col.png. I tried putting a color table in this tiff via a ...
2
votes
0answers
431 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
1answer
376 views

gdal_grid with CSV file and VRT “No point geometry found on layer”

Given the following .vrt file which I manually created: <OGRVRTDataSource> <OGRVRTLayer name="DGM5_BE"> <SrcDataSource>DGM5_BE.txt</SrcDataSource> ...
0
votes
1answer
72 views

QGIS processing does not support jpg

I would like to work with the gdal functions provided in the processing toolbox. It is easy to call: import processing processing.alglist() ...
3
votes
2answers
375 views

ENVI raster file with nodata change to 0

I have a geotiff that I converted to a ENVI raster file with the following command: gdal_translate -of ENVI srtm_05_08.tif data.bil The geotiff has water values supposedly set to nodata. I am ...
3
votes
2answers
692 views

How do I convert to EPSG:3857 with GDAL?

I have some TIFs in EPSG4326 the GeoTiff info looks like: GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTHORITY["EPSG","7030"]], ...
1
vote
1answer
231 views

gdal_proximity distance to meters

I followed the answer given in this topic for searching the distance to the nearest coastline given a lat/lon coordinate. The distance output obtained is in degree (I think) but I need the distance ...
2
votes
0answers
337 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 ...
0
votes
0answers
95 views

gdal_translate from within pyqgis [duplicate]

I want to use gdal_translate from within a python script in my plugin qgis2leaf. I've tried with ...
1
vote
1answer
715 views

How do I convert a .bil (and associated .clr, .blw, .hdr, and .stx) files into georeferenced rasters in R? OR to geotiffs using gdal_translate?

The title says it all. I've got a bunch of .bil files, which I can open in R using the raster command pointed to the path to the .bil file. But I don't know how to efficiently incorporate the ...
-3
votes
1answer
713 views

tif tfw to GeoTiff [duplicate]

How do I combine a .tif file with a .tfw file to create a GeoTiff? There are many answers telling me to use gdal, but I have no clue. So can someone provide a step by step example of how to do this ...
1
vote
1answer
120 views

Extracting z values from shapefile using OGR Api

Is there an example of how we can extract the z value from a shp file given the x,y coordinate value (namely utm) using gdal/ogr c++ lib? alternatively, would it be easier if we convert the shp file ...
0
votes
0answers
134 views

What is GDAL equivalent to ArcGIS Lookup (Spatial Analyst) tool?

I use the Lookup (Spatial Analyst) tool to create several raster layers from a combined raster layer. I know ArcPy can do the same thing. However, I would like to know whether GDAL (or another ...
1
vote
1answer
396 views

Map project a raster having separate latitude and longitude raster bands

I have a satellite image where one of the channels is a raster having the longitude values, another channel a raster containing the latitude values. What would be the best way for map projecting such ...
0
votes
1answer
132 views

Not able to find some of the gdal drivers when using the java bindings as opposed to gdal by itself

I am new to the gdal scene so there might be a simple explanation for my problem. When executing ogrinfo from the command line there is a large list of supported formats. ...
2
votes
1answer
76 views

Zoom option error on gdal_tiler.py

I'm using the tiler-tools library from: https://code.google.com/p/tilers-tools/wiki/QuickStart And specifically this command: ./gdal_tiler.py --zoom=9-15 my_map.vrt To generate tiles. When I use ...
2
votes
2answers
996 views

Is there a way to use gdal functions from OSGeo4W out of their shell?

I would like to use some gdal functions from OSGeo4W (gdal2xyz, and gdal_merge) for my work. I can use them in the OsGeo4W shell provided after installing the package. But I want to use them in ...
4
votes
3answers
311 views

How to compute the outline of multiple tracks?

I have multiple gps tracks with depth data gathered from the lake. The tracks are overlapping/intersecting. I want to get an outline polygon. I have a multi-line shp. I've tried to diffuse it in qgis ...
0
votes
1answer
686 views

Creating records in output shapefile using Python

I have a list (output from a previous loop) containing attributes from an input shapefile. What I am stuck with is inputting those attributes to a new shapefile. I have plotted the point data using ...
1
vote
0answers
108 views

GDAL Color distortion around nodata value stitching

I am trying to extract JPEGs from SID images using GDAL. I ran into a an issue where black lines were present and figured out that it was a nodata issue, so I added -src/dst/vrt nodata arguments and ...
0
votes
0answers
36 views

Copy queried attributes of a shapefile to a new shapefile [duplicate]

I have a polyline-shapefile, which I have successfully read into python using pyshp module. My problem is to save the 2nd node of every line feature into a separate, single point file. I have ...
5
votes
1answer
429 views

ogr SimplifyPreserveTopology does not keep the Topology

I use ogr2ogr to simplify a layer of polyline, from the document: -simplify tolerance: (starting with GDAL 1.9.0) distance tolerance for simplification. This method will preserve topology, in ...
2
votes
1answer
2k views

How to change GDAL version in QGIS to get ECW import to work

I’ve previously been successful in importing ecw into QGIS 2.0 but the same methods (for example How to add *.ecw files to QGIS 2.0 on Windows?) have failed with a fresh install of QGIS 2.2 with ...
0
votes
1answer
203 views

Intersecting great circle arcs and polygons

I am trying to find the point of intersection between a great cirle arc and a polygon. I am working in long/lat, where I am given the endpoints of the great circle arc, A and B, and the vertices, C, ...
4
votes
2answers
140 views

“Levelling” elevation data coinciding with polygon?

I have a shapefile containing building geometries (polygons) and a elevation dataset. The elevation data doesn't perfectly match the buildings. I want to "flatten" the elevation data so that ...
2
votes
0answers
111 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 ...
0
votes
0answers
247 views

gdal read whole geotif band into vector

i have another c++ GDAL question. I try to read alle elevation data out of an GeoTiff file into a c++ vector so that i can get elevation at x,y like this: vector[x][y]. I also like to use this data to ...
2
votes
2answers
576 views

Import xyz data from csv into gdal

I have many CSV file in the form of lon, lat, Z that in would like to import into gdal for further processing. I need to do this on the fly in a python script so any step involving gdal_grid or ...