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

4
votes
0answers
20 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
55 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
19 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
18 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
0answers
15 views

HOW TO: python gdal batch process project L2 hdf and the newest files that were added everyday

I have downloaded files from an ftpsite(the files uploaded are delayed one day, means the files for today aren't available so the newest files available are the hdf files yesterday and so on). These ...
1
vote
3answers
44 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
26 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
0answers
15 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
11 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
49 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
54 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
42 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
105 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 ...
1
vote
0answers
24 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
19 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
35 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
16 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
18 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
votes
0answers
17 views

NDVI code using C++

Using C++, I need code to calculate NDVI and output a TIF format. I am working on a project, I am also using GDAL for data manipulation.
1
vote
1answer
60 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
36 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
2answers
47 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
27 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
40 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
46 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
0answers
36 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
35 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 ...
1
vote
0answers
75 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
53 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
29 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
55 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
57 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
54 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
45 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
24 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
32 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
41 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
46 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
51 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
32 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) ...
0
votes
1answer
35 views

removing spurious lines in raster before or after vector conversion?

I'm looking for a feature/function in either python opencv/ndimage/gdal/etc. or qgis/saga/grass that will help clean up a skeletonized image that I'd like to vectorize into a shapefile of ...
0
votes
0answers
12 views

Confuse about how to make a tiles to support retina device

I have a map picture(15772*13796) so I break them to tiles using gdal2tiles(512*512px not 256*265px) for support retina device but is this correct? Because I a had read this issue and this. They ...
0
votes
1answer
22 views

Shifting image in a vrt

I have a set of geotif images that I want shifted by an integer number of pixels in X and Y. However, those images are quite large, so I don't want to duplicate them (and I need to keep the original ...
1
vote
0answers
30 views

Use Gdal2tiles to generate 512px tiles

I'm a newbie in gdal2tiles & python but I would like to modify gdal2tiles to create 512px tiles for retina display (screen). Can someone suggest how I can do this?
1
vote
2answers
48 views

Error getting ESRI Shapefile driver OGR in C++

I am unable to create an OGR Layer using the C++ API. I have followed the documentation at this link exactly http://www.gdal.org/ogr_apitut.html I get an error trying to get the shapefile driver. Any ...
1
vote
1answer
46 views

Crop and reproject shapefile based on GeoTIFF sample

I have GeoTIFF with these attributes: Size is 495, 309 Coordinate System is: PROJCS["unnamed", GEOGCS["unnamed ellipse", DATUM["unknown", SPHEROID["unnamed",6370000,0]], ...
2
votes
0answers
37 views

Any solution to mosaic Modis MOD13Q1 using batch file?

I have a Modis MOD13Q1 product for Pakistan, four tiles cover my complete area and I have data consisting of 14 years, after every 16-day interval. I want to mosaic them using a batch file which ...
1
vote
0answers
38 views

FWTools Python ImportError: DLL load failed:parameter is incorrect [duplicate]

Please, can anyone help me solve this issue? Please find the attached snapshot for your reference. I am working on a Windows 7 Professional, Service Pack 1, 32-bit laptop.
1
vote
2answers
93 views

Import error:No module named _gdal

I have installed python 2.5.1 and FWTools2.4.7 in my laptop(Windows 7) when i run my pythonscript in FWTools i get following error Import error:No module named _gdal can any one suggest how to solve ...
0
votes
2answers
27 views

Aggregating rasters with GDAL

I have a set of files that have been mosaiced into a global VRT in a particular projection and resolution. However, I have a need to report statistics on a (say) one degree grid, where the reported ...