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
35 views

Gdal_grid gives blank .tiff file after conversion

After using gdal_grid I get blank .tiff file. I use the following .vrt file : <OGRVRTDataSource> <OGRVRTLayer name="p"> <SrcDataSource>p.csv</SrcDataSource> ...
2
votes
1answer
253 views

how to use gdal_rasterize? how to know which attribute to burn with ogrinfo?

I have a vector layer I need to rasterize. I guess this is a basic operation, but I haven't found the information I need yet. First, I have reprojected the vector shapefile: ogr2ogr -s_srs ...
1
vote
1answer
79 views

geotiff conversion to coordinate system with gdal

I would like to get the coordinates X Y Z from a geotiff file. I am using the gdal translate for the conversion. Which output format is the more adequate for containing a 3 numpay array X Y Z ?
2
votes
1answer
80 views

How to convert ASCII in x,y,z format using gdal_grid?

I have a map saved in (AFAIK) Arc/Info ASCII format and I would like to convert it with gdal grid but I am getting errors. My data (p.csv) is in the following format: 319053.25,319301.12,544.28 ...
0
votes
0answers
80 views

Merge/mosaic categorical rasters in QGIS/GDAL

I am attempting to make yearly bins of some monthly MODIS data for all of the US. The cells of the raster are defined with Julian dates indicating the date at which a burn event took place. After ...
0
votes
0answers
45 views

Read esri shape point data and spline interpolation using python shell scripts?

I am quite new in python scripting. I am a dedicated user of arcgis but this time I have to work without arcgis. I have a point shapefile with 10 different attribute fields of different peat layer ...
1
vote
0answers
65 views

How to read .adf (TIN) map files?

I have a set of maps in few formats : TIN ├── [4.0K] info │   └── [2.6K] arc.dir ├── [ 708] log ├── [374K] M-33-45-D-d-4.asc └── [4.0K] m-33-45-d-d-4tin ├── [ 96] tdenv.adf ├── ...
0
votes
1answer
89 views

Does gdal_calc.py not calculate arithmetic with in the same file?

I have one GeoTIFF file with 2 bands. The data type is Byte. I want to calculate (Band2-Band1)/(Band2+Band1). But the following command provided a dark image: gdal_calc.py -A infile --A_band 1 -B ...
0
votes
0answers
63 views

Read large raster into lat/lon value pair to csv file using python and/or gdal?

I have rasters for all of the states in US and they are in 10m resolution. So they can be very large, containing a large number of pixels if you think of TX or AK. I want to create a csv file with ...
2
votes
3answers
129 views

Use GDAL command line to copy projections

I have two images, a georeferenced Tiff and a second Tiff without a reference. The second image is a copy of the first, but it lost its reference. Now, I want to copy the reference from the first Tiff ...
1
vote
0answers
41 views

Clip Aerial Image Pair by Overlapping Area?

I have a set of aerial images which I am trying to clip by their overlapping areas, for use in a more involved program I am writing in Python with Python Imaging Library (PIL). Is there a command ...
0
votes
1answer
47 views

Trouble with r.out.gdal in GRASS GIS

I'm processing some landsat 8 files. I'm using this scene LC82260982013325LGN00 and I have to apply i.landsat.toar and export the result as a GTiff. However when I use r.out.gdal it takes a very long ...
3
votes
1answer
113 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
31 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
97 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
42 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
116 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
42 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. ...
2
votes
1answer
79 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
63 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 ...
3
votes
1answer
105 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) ...
0
votes
1answer
126 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
117 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
365 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
41 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 ...
0
votes
0answers
31 views

Gdal_Calc.py random pixel error

I'm using the gdal_calc.py in python in order to do the same thing as the Con tool in ArcGIS. I'm calling the .py using os.system: os.system('gdal_calc.py -A%s -B%s --outfile=%s ...
4
votes
2answers
166 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
0answers
82 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
59 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
61 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
91 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
0answers
21 views

Opening Spot 5 DIMAP imagery in GDAL

Would like to know if someone is familiar of opening SPOT 5 DIMAP imagery using GDAL. Though the GDAL provides DIMAP driver, when I tried to open, it can not show the geotransform and geoprojection.
0
votes
1answer
34 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 ...
4
votes
1answer
205 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
135 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. ...
1
vote
0answers
115 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 ...
0
votes
1answer
326 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
120 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
151 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
49 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() ...
2
votes
1answer
68 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 ...
2
votes
2answers
126 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
57 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 ...
1
vote
0answers
72 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
56 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
111 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 ...
0
votes
0answers
20 views

Gamma correcting ECW

we would like to modify an existing ECW file with a gamma correction in order to make the image less colourful. Is there a way to do this using e.g. GDAL?
-3
votes
1answer
68 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
38 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
69 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 ...