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

Reclassify rasters using GDAL and Python

Using @SS_Rebelious very useful answer, I am trying to reclassify rasters using GDAL and Python. I have a grayscale raster (with values ranging from 0 to 255) and I want to reclassify it to a binary ...
3
votes
2answers
92 views

gdalwarp leaves horizontal artifacts regridding from EASE-Grid (laea) to Polarstero (stere)

I'm trying to understand how to use gdalwarp to warp an image. I believe I'm doing the basics correctly, but I might be missing some gdalwarp options? Any help would be greatly appreciated. The ...
6
votes
2answers
539 views

How to install GDAL and QGIS on a Raspberry Pi?

I would like to run GDAL QGIS on Raspberry Pi to act as an affordable GIS processing engine. The plan is to buy about 50 Pis. Does anybody have any links or information on installing and running GDAL ...
3
votes
2answers
73 views

Export Raster Attribute Table using Open Source software

I've got 120 or so geotiff files that I'd like to build and export the raster attribute table to something a little more friendly (e.g. csv). I'm curious if there is an open source way to do this. ...
1
vote
0answers
69 views

Using gdal ReadAsArray or Arcpy.RasterToNumPyArray on python 2.6 32bit without it crashing

I am trying to run ReadAsArray on an entire raster and on a single band, python will crash without an error message when trying to do either. Is there a way I can create a 3d array without having to ...
1
vote
1answer
52 views

how to open shp file which created by using python-ogr in qgis?

I have done following steps, so how I open that created shp file in Qgis driver=ogr.GetDriverByName('ESRI Shapefile') out=driver.CreateDataSource('pts.shp') ...
0
votes
1answer
35 views

how to beginning with GDAL, how can i use it for image format conversions

I like to work with GDAL. is it a software with or without graphical user interface? or the commands are run in command prompt.
9
votes
0answers
236 views

Geotransformation for polar stereographic?

I'm currently working to import CANGRID climate data (provided as Surfer Grid ascii, ".grd" files) into ArcGIS. The grid is 95 rows by 125 columns in size.Metadata provides lat/lon of origin (lower ...
2
votes
0answers
44 views

How to Use GetGeomtype Gdal Python

I am just wondering what I am doing wrong when trying to use the GetGeomType. I need to focus on separating layers from a pix file. Two categories are the bitmaps(which are vectors) and then the ...
0
votes
0answers
28 views

Ignore empty tiles in gdal2tiles.py when generating tiles

I am in the process of generating a large set of tiles from a .vrt file composed of many images. The image is following a road section, being very long and small strip of images and due the shape of ...
0
votes
0answers
38 views

gdalbuildvrt not working with star in path on windows

In my directory C:\dev\Jpeg-geo\del1 I have a set of .jpg and .jgw files that I am trying to put together with gdalbuildvrt but its giving me a weird error. ...
4
votes
1answer
117 views

Removing some bands or layers from a tiff using gdal

I have been working on a multiband raster image. It is a pix image. I need to just work on the 26 bands that are raster files and not the 4 bitmaps that are are also included in the pix image. I could ...
3
votes
1answer
222 views

Using GDAL in python to stack Landsat bands

I am processing a large amount of landsat data. Landsat scenes ship as a collection of tiffs- one file for each spectral band, plus a bunch of masks. I need to stack some of the bands into a single ...
1
vote
1answer
196 views

How to get the mapserver WFS OGR driver to encode unicode in geojson

Running mapserver 6.4.1, gdal 1.14.3 I try to output some WFS features using outputformat=geojson. The relevant bits of the mapfile are appended below. Everything works fine, but unicode characters ...
1
vote
0answers
48 views

How to run cmd prompt in the background and show progress bar

I am warping files using gdal. Since I have numerous files to warp how will I run the cmd in background and just have a progress bar to indicate the percentage in warping the files? Because when I ...
1
vote
2answers
216 views

gdalwarp batch process

I have downloaded numerous hdf files from an ftp site. Half of them were level 1 and the rest were level 2. Lvl1 format is Afile.A1_AC.sample.hdf and the lvl 2 format is: Afile.A1_AC.sample.hdf. How ...
3
votes
3answers
212 views

How to find all distinct values from a single column of an attribute table with GDAL (Python)?

I am working with GDAL and modify ESRI shapefiles with it. I am searching a possibility how to find all distince values from a single column of an attribute table with GDAL (Python). In SQL it would ...
1
vote
1answer
64 views

Raster created by python GDAL interpreted upside down

I am using CRS which has axis in left bottom corner. Code below produce Geotiff which is shown allright in QGIS but when I open it in other software it is upside down. Also gdal_calc takes it as ...
2
votes
1answer
165 views

Read DGN file in java or c#

I'm looking solution to read DGN file in java or c# without any success. Gdal makes many problem and doesn't support all element eg. complex element. Does anyone know of a library that supports the ...
0
votes
0answers
28 views

Adding a column to several feature classes in a file gdb using ogr2ogr

I'm trying to add the same column to 25 feature classes in an ESRI FileGDB. I managed to do this for one feature class by using ogr2ogr: ogrinfo C:\temp\test.gdb -sql "ALTER TABLE table01 ADD COLUMN ...
2
votes
2answers
67 views

Can GDAL fill elevation below a certain value?

I have several GeoTIFFs that I'm working with to create contours using gdal_contour. One of them crashes the utility. gdalinfo tells me that the computed min values for band 1 and 2 are ...
2
votes
2answers
104 views

Compression Artifacts and GDAL

I am working with large tiled (no overlap) datasets of ortho imagery that has been mosaicked in ArcGIS and exported as a single, merged TIF for the entire AOI. These images need to be compressed so ...
2
votes
1answer
73 views

GDAL Python Not Working Correctly - Raster Calculations

I am sure this question has been asked and answered before. I am working with a x64 bit installation of OSGeo4W. I am trying to perform raster calculations using the gdal through the OSGeo4W shell, ...
2
votes
2answers
400 views

Exporting 16 bit single-band to 8 bit in SAGA/QGIS

Everytime with every function (e.g convert data storage type or gdal-translat -ot Byte) the output is a completely white map with a black border. only the "save as" export works, but the TIFF also ...
0
votes
0answers
40 views

georeference a grid by using a script

I've got an image and a set of coordinates, which represent a grid, that looks similar to this: What I'm looking for is a script, that would transform my image, by using such a grid. Edit: My ...
0
votes
0answers
51 views

Precision loss when transforming coordinates systems in GDAL (vector layers)

Anybody knows the order of magnitude of the error introduced when transforming coordinates using GDAL? More specifically, I would like to know if significant error is introduced when the ...
0
votes
1answer
80 views

Mapserver: “Premature end of script headers: mapserv,”

Some of requests to the Mapserver return "Internal server error". When I checked the Apache log I saw that there are several error messages: GDAL: In GDALDestroy - unloading GDAL shared library., ...
2
votes
1answer
96 views

ogr2ogr update failing to update shapefile int PostGIS

I have a bunch of shapefiles that I am looping in with a bash script to load and update files but the ogr2ogr fails. I am running gdal 1.92 on ubuntu 14.04. for file in `ls *.shp`;do ogr2ogr ...
1
vote
1answer
57 views

Why does the commands listgeo and gdalinfo don't show the same Proj4 Strings?

As the title says, why doesn't the unix tools listgeo and gdalinfo show the same proj4 strings? Check this file: ...
1
vote
1answer
207 views

Using GDAL to read dted files

How do I get the elevation at a certain point using gdal to read in a dted file? I don't see much useful information out on the web. I'm sorry I should have been more specific. How would I ...
0
votes
1answer
108 views

Create Georeferenced Raster For JPEG file with Mapnik Python

I am able to use Mapnik generate_image.py to generate image with a BBOX and projection. Also, I have added a python script to generate a JGW file. But to completely georeference the generated image. ...
0
votes
3answers
224 views

Georeferencing window won't open QGIS 2.4

I had the georeferencer window open in QGIS 2.4 while project was running, I accidentally closed the Qgs document with the georeferencing window open (worked absolutely fine till then). When I opened ...
1
vote
1answer
53 views

Aerial Imagery Showing in Different Colour while pan/zoom geoserver

I have img/tif files of a small aerial imagery. I published it under ImageMosaic plugin (GDAL plugin installed), but when I pan or zoom the map in OpenLayer preview I am getting different colour for ...
1
vote
1answer
102 views

Problem installing gdal - works in build directory, not in /usr/bin

I've been trying to install gdal with the FileGDB driver. I installed it fine on a linux mint machine. I am having problems installing on a CentOS machine though. It builds fine, I can run the apps in ...
1
vote
0answers
83 views

NASA SRTM: precipitous contour lines in TileMill

I was trying to make a nice looking print map of Bali (Indonesia) in TileMill. After creating a terrain layer (using this walkthrough) I noticed that the island's coastline was precipitous. Is ...
2
votes
1answer
65 views

Python and gdalwarp: change resolution of raster

I am new to GDAL and have been having some hard time figuring out how to change a raster's resolution and how to align it to another existing raster. I have come to the point where I have a raster ...
1
vote
0answers
81 views

Batch GDAL_Calc multiple band calculation

I have a question concerning gdal_calc. I´m running windows 7 64b and am working from the OSGeo4W shell. I have multispectral images that need to be prepared for analysis but they are not saved in a ...
3
votes
0answers
464 views

How to correctly import GeoJSON in PostGIS? [closed]

I'm trying to migrate GeoJSON data to PostGIS database and I thought that ogr2ogr is the tool: ogr2ogr -f PostgreSQL PG:"host={} user={} port={} dbname={} password={}" test.geojson -nln ogr_test ...
1
vote
1answer
115 views

Using gdalwarp with multiple elevation models

I'm trying to use gdal_warp to slice and convert some SRTM data. The data I'm interested in are in S30E152.hgt and S30E153.hgt. What I'd like to do is: Convert these DEM files from EPSG:4326 to ...
1
vote
1answer
56 views

How can I fix a broken gdal-tools plugin in QGIS?

I have QGIS installed on my 2011 macbook pro and I recently upgraded my OS to 10.9. After the upgrade, I noticed that QGIS gives me this error (see attachment 1) saying that the plugin is not ...
0
votes
2answers
356 views

GDAL for Python: extracting subdomains from NetCDF file?

I need to make a script that processes a NetCDF file that contains 3 days of hourly forecast data from the norwegian meteorological office. The NetCDF file contains various data I need ...
1
vote
1answer
434 views

qgis polygonize error 4 “file does not exist in the file system”

When I try to convert a raster DEM (in asc or geotif format) to polygons I get the error below. ERROR 4: `Y:/GIS/oversvmmelser/tmp/out.tif' does not exist in the file system, and is not recognised ...
2
votes
1answer
120 views

Reprojecting MODIS Ice Sea Temperature products with gdalwarp

I have problems projecting MODIS Ice Sea Temperature (IST) products to EPSG:3031 (Antartic Polar Stereographic) using gdalwarp. The product I'm working with is MOD29E1D, documented here, p33: The ...
2
votes
1answer
147 views

Add band to GDAL VRT dataset and convert to TIFF

I'm trying to add a band to a TIFF dataset by converting it to a VRT format and using the AddBand() function, then convert it back to a TIFF to access the band and write the data. Here is an example ...
0
votes
0answers
42 views

Dxf to Shp - Coordinates wrong

I've converted a dxf file to a shape file using the folowing ogr2ogr command: ogr2ogr -s_srs EPSG:31464 -where "LAYER='0'" -f "ESRI Shapefile" dgm.shp DGM1.dxf Or ogr2ogr -where "LAYER='0'" -f ...
2
votes
1answer
136 views

QGIS show gdal_polygonized file in EPSG:54004 instead of EPSG:3857

I am polygonizing a TIFF with the gdal_polygonize.py and the resulting Shape is shown in QGIS in the wrong projection - EPSG:54004 instead of EPSG:3857. I have a TIFF in EPSG:3857 and polygonize it ...
0
votes
1answer
95 views

Gdal: How to reproject from equirectangular to mercator?

I work with raster ETOPO1, made a crop of my area of interest. curl -L -C - ...
1
vote
2answers
111 views

GDAL Application crashes when loading GDAL drivers

I'm developing on QT an application that uses GDAL library (which I compiled using MinGW 4.8). When I use the application on my windows 8 laptop, it works fine, but when I tried using it on several ...
1
vote
1answer
147 views

Problem with Converting XY Coordinates to LatLon

Here is the function that I use to convert from pixel x,y to lat,lon coordinate: public static Coordinate imageToWorld(Dataset gds, int x, int y) { double[] adfGeoTransform = new ...
2
votes
1answer
72 views

Is this a known issue in gaps between masked raster and spatial polygon in R with mask command?

This is my codes raster file: RAD polygon(shape): ws_sfc RAD.sub <- crop(RAD, extent(ws_sfc)) plot(RAD.sub) RAD.sub <- mask(RAD.sub, ws_sfc) plot(RAD.sub) plot(ws_sfc, add = TRUE) ...