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

1
vote
2answers
76 views

Convert a grid file from one set of coordinate systems to another

I have a grid file in a certain coordinate system. I'd like to change it so I have the same grid in geographic lat/long based on WGS84. This seems like something that gdal can tackle quite easily ...
1
vote
0answers
20 views

Python Band.SetNoDataValue(), sets all bands at once

Can someone tell me what I am doing wrong? I want to set bands 1 and three to 255, and band 2 to 0. The code below only should only set band 2, but gdalinfo shows that all bands get set to 0. Ubuntu ...
2
votes
1answer
50 views

Determine spatial average of several tiff files

I am relatively new to GIS. I have a set of tiff files containing a raster each. I would like to average them, though their extents are slightly different. They each have their geospatial data ...
3
votes
2answers
124 views

Piping data to ogr2ogr

Is it possible to directly pipe data to ogr2ogr? I'm trying to curl a GeoJSON file and pipe it to ogr2ogr to covert it to another format. Without using a pipe, I was able to achieve this by first ...
1
vote
1answer
31 views

Mapping NetCDF CF Mercator parameters to GDAL

I have coordinates on a mercator grid projection from a NetCDF (CF) file. The metadata attributes in the file include : standard_parallel, longitude_of_projection_origin, earth_shape (spherical), ...
3
votes
2answers
105 views

How to create a raster tiff (dem / heightmap) from a mesh of xyz points?

I have to create a terrain heightmap(raster tif preferably) from mesh of xyz ungridded points(in geographical coordinates), I tried taking min and max of bbox of mesh and making width & height by ...
0
votes
1answer
19 views

Merging separate country features together into one region

I am trying to creating a map of the UK where Scotland and England are merged into one shape and not separate shapes for each country. As a minimum working example, I'm following Mike Bostock's Let's ...
0
votes
0answers
47 views

Fiona write crashes

I am trying to write a new shapefile with fiona but after creating the features, i keep getting this output after the program crashes: The program '[6668] python.exe' has exited with code -1073741819 ...
1
vote
1answer
60 views

Unsupervised raster classification using GDAL

Is there an equivalent in GDAL to the Arcpy ISO data unsupervised classification tool, or a series of methods using GDAL/python that can accomplish this?
0
votes
0answers
13 views

Change “World Size” with gdal2tiles.py

When I use the flag "-p 'geodetic'", the output is a world that is 512x256px at zoom level 0. In order for it to overlay with other datasets, I need it to be 640x320px at zoom level 0. I am looking ...
2
votes
1answer
39 views

Which algorithm does gdal_polygonize use?

I would like to know which algorithm (i.e. Sobel operator ) is used to extract the polygons in gdal_polygonizer. I searched here and here but could not find it. Any help would be appreciated.
2
votes
1answer
123 views

Convert GeoPDF with GDAL

I'm attempting to convert a USGS US Topo GeoPDF in Windows. I'm using GDAL 1.11.1. My initial attempts were with this Python script, which automates the process by feeding a command to the Windows ...
1
vote
0answers
28 views

“fill in” .vrt with missing locations in a mosaic

I have ~2000 NASA SRTM .hgt files. I wanted to join them together in a single mosaic, and to save disk space I used GDAL to create a .VRT rather than a very large GeoTiff. This is read by Mapnik. ...
1
vote
0answers
32 views

different data range of raster file in mapserver and qgis

a NetCDF file with only one band with a data range from [0,1.1]. IDV & QGIS identifies it correctly, however gdalinfo finds the data range [1,26505]. Or could it be possible that gdalinfo is right ...
1
vote
3answers
45 views

Lossless conversion from shapefiles to XML/JSON/Text

I'm a software developer and I need to create a custom tool that imports ESRI shapefiles (shp + shx + dbf), manipulates the content and then exports new shapefiles. To avoid a complex binary reading ...
1
vote
1answer
50 views

Remove values below a given value with GDAL for DEM geo tif files

I have a DEM in GeoTiff format. The processing have made areas that have really low values that I would like to remove. Are there any GDAL tools that could set all values below -30 to NODATA or -30.
2
votes
0answers
52 views

Compression in GeoTiff files with GDAL

I have an output from 3th party program in tif format where gdalinfo gives the following information: Driver: GTiff/GeoTIFF Files: J:\ascend\Lufthavne - ...
2
votes
1answer
50 views

Can ogr2ogr write data to a file geodatabase feature dataset?

I am dumping postgis tables into a file geodatabase and it works great--the result is a file geodatabase full of featureclasses. But is it possible to write these featureclasses into feature ...
2
votes
0answers
32 views

Getting a raster extent excluding nodata area

I have a GeoTiff raster in which half of its area is nodata. The problem is that the raster extent includes the nodata area. Is there a way with GDal to recalculate the extent of the raster excluding ...
0
votes
0answers
47 views

GeoTiff using GDAL

How can I project a tiff using the position of the four corners where the picture have to be placed ? The tiff file has no georeference or projection information I have got a grid with values of ...
2
votes
1answer
42 views

Open Source counterpart to Euclidean Allocation

I'm looking for a function that can perform something similiar to the Euclidean Allocation process in arcGIS ...
1
vote
1answer
23 views

Invoke OSGeo4W commands with Java

I want to use the command gdal_polygonize image.png bla in a Java for-loop, where image.png and bla changes in every iteration. How can I do this? I searched for it on google and stackexchange and I ...
0
votes
0answers
52 views

Google Tracks convert kml to shape and preserve attributes

I've been attempting to convert a Google Tracks kmz file to an ESRI shapefile using the GDAL and libkml(http://www.gdal.org/drv_libkml.html) which I built on my Ubuntu 14.04 (64-bit) system. The kml ...
2
votes
1answer
29 views

gdaldem output file size explodes

I want to create a color-relief and a hillshade with gdaldem and use them with mapnik. System: Ubuntu 14.04 LTS Server, gdal-bin 1.10. Here are my steps: Download all SRTM tiles from ...
2
votes
0answers
30 views

Diminishing performance of gdal_polygonize in a python subprocess for GRASS to PostgreSQL

I have written a python script to polygonize GRASS rasters into a PostgreSQL database using the gdal utility gdal_polygonize: sub.call('gdal_polygonize.py -8 gdb/quercus/{0}/cellhd/{1} -f ...
2
votes
1answer
81 views

Gdal/QGIS resampling asc file

I have a DEM (.asc) file with a resolution of 1000 meters per pixel. Therefore the file has a lot of columns and rows. I would like to chance the resolution per pixel from 1000 to 3000 meters per ...
0
votes
1answer
46 views

Tools/Software for QA/QC of data

I currently have written a few custom scripts that run basic checks (null values, geocoding errors, fullness of metadata) using Arcpy(10.2 with almost all licenses) and gdal. I'm interested to know ...
2
votes
1answer
57 views

Export PostGIS line to AutoCAD with graphical direction information

Some time ago I managed to export from PostgreSQL to AutoCAD DXF format a line (LINESTRING) using ogr2ogr program with a PostgreSQL Query, and it works (I can visualise it by the QuantumGIS). Now I ...
0
votes
1answer
39 views

gdal_polygonize unable to open PNG

I use OSGeo4W on Windows 7 (64 bit) and included libpng16.dll When I try C:\OSGeo4W64\bin>gdal_polygonize star.png bla I get the following error: ERROR 4: `star.png' does not exist in the file ...
4
votes
0answers
77 views

Working with MODIS in R

I am new to using MODIS data and was hoping to get some help starting out. I have been able to download the hdf file but viewing it in R has been a problem for me. What I have found so far, and ...
0
votes
2answers
79 views

building GDAL with libKML support

I'm attempting to import a kml/kmz file generated using Google Tracks into my QGIS 2.4 system (built with GDAL 1.10). KML files currently lose all attributes, and as such I'm attempting to build GDAL ...
2
votes
1answer
79 views

Re-project from Polar projection to EPSG:4326

I simply want to project a polar stereographic raster (EPSG:3411) to Mercator (EPSG:3857) or WGS84 (EPSG:4326). The data I is from NOAA. $gdalwarp -s_srs EPSG:3411 -t_srs EPSG:3857 northpsg.20141027 ...
0
votes
0answers
22 views

Script Efficiency

I have been working on a script to read pixel values and print the values out to a text file. This is currently working but takes around two hours with a 6.6gb file. This is to replace an older script ...
1
vote
0answers
43 views

Extract polygons of PNG with GdalOgrInJava

I have PNG pictures with polygons of different colours on them. What I plan to do is to extract these Polygons using the polygonize function (see here) Then I will try to encode the polygons in ...
1
vote
1answer
40 views

GDAL external overviews not recognized by ArcGIS

I am adding external overviews to a number of TIF files with GDAL as follows: gdaladdo -ro file.tif 2 4 8 16 --config COMPRESS_OVERVIEW LZW Then, when I add the tif file to ArcGis is still asking ...
1
vote
1answer
82 views

Could not open/read USGS raster file

I could not open/read the following .img file. The document below says it is a raster file. http://edc2.usgs.gov/1KM/land_sea_mask.php I tried to open it using many ways such as Q-Gis software, and ...
2
votes
0answers
41 views

What is the issue with geopdf format?

I'd like to add the geopdf format as output of a library that I am working on. After have looking around, I have not been able to find a (c, c++, python) library that meet my license requirements ...
0
votes
1answer
46 views

Writing shapefile with projection defined crashes Fiona part two

I have the same problem as this post. Writing shapefile with projection defined crashes fiona It says that was an issue with GDAL. I cleaned my site-packages of anything GDAL, fiona, six related. ...
3
votes
1answer
43 views

Gdalwarp with “out of globe” values in sinusoidal projection

I want to project a global image (source ERA INTERIM) from WGS 84 to sinusoidal, using gdal 1.11 and bilinear resampling. It work fine BUT I have values everywhere on the rectangular extent of the ...
0
votes
1answer
50 views

Geoprocessing tool Difference from the Terminal?

I have used QGIS Geoprocessing tool Difference to cut polylines (SHP) into to smaller pieces using a polygon (SHP). This works very good in QGIS. Now, to speed up the process, I'd like to be able to ...
0
votes
1answer
48 views

gdalnumeric in python virtualenv

I'm having a problem importing gdalnumeric when working on virtualenvs. I'm getting the following error: from osgeo import gdalnumeric Traceback (most recent call last): File "<stdin>", line 1, ...
1
vote
1answer
44 views

Enabling Erdas ECW files in gdal

I want to enable read support from Erdas ECW files in gdal on windows. I have GDAL binaries from http://www.gisinternals.com/sdk/ and i also have a self compiled version. when compiling myself I ...
1
vote
1answer
41 views

Use subprocess module to run GDAL processes

I am attempting to write a python script that will reproject & merge a series of .tif files then clip them to the boundaries of a .shp file. I have successfully used the subprocess module to do ...
3
votes
0answers
48 views

How to set Vertical Coordinate system when creating a 3D shapefile?

I'm working with GDAL v1.11 and am seeing the vertical CS stripped when exporting a 3D shapefile (from oracle spatial). This is the command I'm running: ogr2ogr -f "ESRI Shapefile" B_POINT.shp ...
2
votes
2answers
85 views

I've downloaded the GDAL-1.11.1 package but why can't I call on it in python?

I have a Mac and am very new to this coding stuff.. I'm using a code that includes import osgeo as gdal and so I've downloaded GDAL from https://pypi.python.org/pypi/GDAL/ which has a folder ...
1
vote
0answers
40 views

How to add CADRG RPF A.TOC Raster Data Source in Geoserver 2.6

I'm currently trying to add CADRG as a datasource in Geoserver 2.6. I've installed the necessary GDAL plugins (1.9.2) for Geoserver 2.6. I placed my RPF folder in C:\Program Files\Apache Software ...
1
vote
0answers
25 views

Python 2.6 Module for GDAL 1.9.2, ImportError

I work on an internal network(IE I cant get the latest 1.11 stuff) and have GDAL 1.9.2 setup on a Red Hat Linux Enterprise Server 6. My linux box has python 2.6.6 however the python binding I received ...
1
vote
1answer
91 views

What options are there for geolocation-aware SQL databses?

I need to implement an SQL database that can make queries using geo-location based logic. e.g. select rows within 500m of coordinate x, y What options are available to do this with sqlite, mySQL, ...
1
vote
0answers
40 views

Reprojecting and Clipping in one script using Python Subprocess module and GDAL

I am trying to reproject a .tif file then clip it to the boundaries of a shapefile using GDAL Warp. I want to do this in a single script, and have successfully used the Python subprocess module to ...
2
votes
1answer
194 views

calling gdal_merge.py into another python script running GDAL processes

Good Afternoon All, I am wanting to use gdal_merge.py to merge a series of .tif files prior to clipping the merged file to a shapefile boundary, but want to do so as part of another larger script ...