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
2answers
941 views

How to overwrite metadata in a tif file with GDAL?

I have a set of tif files generated from a process with multiple steps. Some metadata tags are added at different steps. I would like to update the metadata of those file before distributing the files....
0
votes
2answers
87 views

gdalcompare.py not found in GDAL 1.10.1

I cannot find the script gdalcompare.py within the osGeo distribution of GDAL by http://www.osgeo.org/. It is the 1.10.1 though. Is there any way of downloading that script separately from any other ...
2
votes
0answers
210 views

Improve gdaldem hillshade output (lighten up flat parts)

I have a hill-shading that I made using gdaldem hillshade which is already pretty good. I used the option -alt 60 to get nice contrast, but the problem is, that with lowering the light-source, also ...
3
votes
1answer
348 views

Conflicts while installing postgresql9.3-postgis2.1,QGIS and gdal-bin

I want to use postgresql9.3-postgis2.1 and gdal-bin parallel . I have refer instruction here for installing postgresql-9.3 on ubuntu 12.04. I have install qgis from here While installing gdal-bin it ...
2
votes
1answer
2k views

How to correctly resize raster GIS images to a given px width?

I work on ETOPO1, a 466MB file, and want to move to SRTM (~20GB Earth). My workflow involve several time-consuming processing steps so I strongly want to work on a minimum raster. It usually want to ...
1
vote
1answer
99 views

why gdaldem is not producing any result?

I have a geotiff file with class values (0, 1, 2, 3, and 4). I want to show each class as different rgb color. So, I created the text file (color_relief.txt) which contans the following: 0 55 78 94 ...
0
votes
2answers
265 views

gdal_grid performance issue — alternatives?

I use gdal_grid within a Python script, but for my large files the calculation time is enormous. Apparently this is a gdal_grid issue http://trac.osgeo.org/gdal/ticket/2411 Does anyone have an advice ...
0
votes
1answer
157 views

going from hillshade to percent shaded (gdaldem)

I assumed hillshades output percent shaded (scaled to 255) at the given azimuth and altitude but this is incorrect. A flat raster at 80 degrees altitude has a maximum value of 251 instead of 255. How ...
0
votes
1answer
233 views

Gdal_sieve not working?

I'm using the Sieve command within QGIS (Raster > Analysis > Sieve), but when I run the command with any variation of parameters, it completely changes the histogram of my raster data from between -0....
1
vote
0answers
33 views

Sending commands into FGDB

I have a large .gdb of all blocks in the United States, available from Census. I would like to create a point shape file of all census tract centroids in the US. I have QGIS with FileGDB enabled, as ...
8
votes
2answers
190 views

How can I use data from PDF map?

I have a map over Sweden with plenty of data represented as circles of various sizes. The map is saved as PDF from Illustrator, fully editable. I wonder if there is any way I could get these circles ...
4
votes
4answers
795 views

Fast and robust open source portable script or tool for interpolation - triangulation or similar

I am looking for some fast and robust open source interpolation library/script/tool. The goal is to re-interpolate relatively big (1M+ points) unstructured (xyz) grids into high resolution cartesian ...
1
vote
3answers
613 views

Using very large image formats with GDAL

I would like to use gdalwarp and gdal_translate with some very large raster images (313,040 x 313,040 px) and I'm wondering what the best format or method is to do this, if it's even possible. I have ...
1
vote
0answers
167 views

GDAL WMS in QGIS 2.4

I'm trying to add in QGIS a layer shown here : propluvia When I watch in my console I see this request corresponding to the layer I need. The server is not able to respond to a getCapabilities ...
1
vote
0answers
105 views

Converting MDB File With Unicode Characters

I'm using C# wrapper of GDAL and I'm having a problem converting MDB file that contains unicode characters. Lets take for example conversion to KML file, non-English characters are replaced by '?' ...
2
votes
0answers
55 views

Error in polygonising raster data with QGIS on Mac

I'm using QGis Chugiak on a Mac and I'm trying to converting a raster layer in a vector data by: Raster/Conversion/Polygonize. So, I select the layer needed from the drop-down menu, I give a name to ...
4
votes
2answers
1k views

GDAL (Python) Aggregate Raster Into Lower Resolution

I have a (continent-scale) raster with 30m resolution in an Albers equal area projection. Each 30m cell has a value 1-34. I would like to aggregate information from this raster into a 1km grid such ...
4
votes
3answers
1k views

Install gdal in CentOS without Root

I have been struggling to get python setup to do some image processing, so I really need the http://www.gdal.org/ GDAL libraries installed on CentOS 6.5. I want to install the package in a local ...
3
votes
2answers
3k views

Problem in displaying tile images using openlayers3

I have problem with displaying my TMS layer using OL3. This is my first time working with tile images this way so maybe it doesn't have to do anything with ol3, maybe I just don't understand it right. ...
3
votes
1answer
296 views

Merging h5 files from VIIRS to gtiff

I am trying to mosaic all the data from VIIRS of a day into a single gtiff. I've already done it with MODIS data using GDAL. The problem it seems is that the coordinates of each h5 files doesn't seem ...
1
vote
2answers
587 views

Is there any efficient way to create a raster with pixel value representing a unique id calculated based on pixel lat/lon?

I have 10m raster for each state in US. Current pixel value represents something. I want to have 'same' raster but pixel size representing a unique id calculated based on pixel lat and lon. I am ...
1
vote
1answer
372 views

raster2pgsql does nothing

I try to load rasters through raster2pgsql. I try to call it through pgAdmin command line, but when I call any raster2pgsql command (even raster2pgsql -G), it does nothing, I just get another line to ...
1
vote
2answers
511 views

How to restore gdalinfo metadata after ImageMagick processing?

I first print out the file info : $gdalinfo crop.origin.tif result: + gdalinfo crop.origin.tif Driver: GTiff/GeoTIFF Files: crop.origin.tif Size is 1920, 1950 Coordinate System is `' Origin = (66....
2
votes
1answer
217 views

(How) can I write string data to a raster with Python/GDAL?

In Python, I have a NumPy array of string values (categories of some variable) that I'd like to write to a raster file (.tif or .img) using GDAL. The impression I get is that it's not possible to ...
1
vote
1answer
534 views

Convert GML to DXF using ogr2ogr

I cannot make ogr2ogr.exe work. I'm using GDAL 1.11. I keep getting the following error message: ERROR 1: Did not get at least 2 values or invalid number of set of coordinates 709993.98 6181946....
1
vote
1answer
340 views

Convert list of coordinates to OGRGeometry or WKT

I am using a Visvalingam-Whyatt simplification script and the output is a list of coordinates such: [[[[-10885205.10690245 3525662.26531131] [-10885369.01690829 3525374.0741439 ] [-10885424....
2
votes
1answer
164 views

Error reading large virtual raster with rasterio: “Cannot invert degenerate transform”

I've created a virtual raster using gdalbuildvrt, which covers a large extent (England) at a high resolution. The data is fairly sparse, focused around the coastlines with no files covering inland ...
1
vote
0answers
158 views

Rendering a geotiff block with gdal and gdk-pixbuf - Some flaws in the result

I currently try to render a geotiff file with GTKmm\GDKmm (gdk-pixbuf). The file is a tiled geotiff with a block size of 256*256 and three rasterbands. The shortened gdalinfo output is: Driver: GTiff/...
3
votes
1answer
198 views

Are there examples of high-level interfaces to NetCDF

The R raster package and GDAL both have independent wrappers around NetCDF-4* that treat an array variable as a "complete" dataset, i.e. they provide objects or interfaces that seamlessly give access ...
-1
votes
2answers
467 views

Converting *.kml file to *.lyr file so Google Earth can create spatial data to use in ArcGIS for Desktop? [duplicate]

Want to convert the .kml file into .shp file so that I can create new the point, line and polygons data to work in ArcGIS. Actually I want to use Google Earth to make data for ArcGIS.
2
votes
0answers
252 views

Memory leak with GDALOpen?

I'm very new in using GDAL, and trying the tutorial given here. But with a piece of code with only GDALRegister_GTiff(); GDALDataset *poDataset = (GDALDataset *) GDALOpen(argv[1], GA_ReadOnly ); ...
2
votes
1answer
1k views

Is there any way to fast read pixel value using gdal in C#

I create with C# a series of lat/lon value pair in the code(using i, j loop) and I want to extract the pixel value based on either lat/lon pair or pixel location (i, j). I know there is one gdal ...
2
votes
0answers
1k views

Correct way to extract pixel data using GDAL. Current method gives repeated but changing value.

I am trying to open a spectral file and move bands of data from GDAL to OpenCV's Mat data structures. As of now I am just trying to make sure that I can correctly access each band using GDAL and pull ...
1
vote
0answers
102 views

ogr2ogr “Error: pointer being freed was not allocated” on mac mavericks

I've successfully installed GDAL via homebrew on mac But while I'm using command ogr2ogr, keep getting error below $ ogr2ogr -f geoJSON tl_2010_20_tract10.json tl_2010_20_tract10.shp ...
0
votes
1answer
166 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
3k 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 "+init=...
0
votes
2answers
367 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
422 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 ...
1
vote
1answer
413 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 ├── [1....
0
votes
2answers
664 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
127 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
1k 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
85 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
95 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
363 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
52 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 ...
2
votes
0answers
675 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
110 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
507 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
111 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. ...