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 ...
2
votes
1answer
22 views
Lack of Options in Gdal Tool “Build Overviews” QGis
I am using QGis 1.8 and I want to build overviews (pyramids) for satellite imagery using the Gdal Tool. I would like to use the resampling method "cubic convolution", but the option doesn't exist ...
2
votes
1answer
31 views
How can I confirm the existence of pyramids in QGis?
I have created pyramid layers for my ASTER satellite images (.img) using Gdal Tools "Build Overviews" (Pyramids). The dialog said "Processing completed". I don't have the impression that zooming in or ...
2
votes
1answer
119 views
USGS GeoTiff extent via GDAL
I have 2006 land cover GeoTiff from the USGS site (for a county in VA, USA). The projection is Albers Conical Equal Area Units. The download came with an HTML metadata file that tells me the geo ...
2
votes
1answer
44 views
How to decrease GeoTIFF file size by eliminating or negating the blackfill?
I have a 5 band, 16-bit orthorectified GeoTIFF from a satellite image that contains a significant amount of blackfill and has a file size of over 21GB. I would like to decrease the file size by ...
3
votes
0answers
32 views
Why are my raster mean and standard deviation -1?
I have a set of about 25 rasters whose descriptive statistics I have been examining. The data in the rasters were calculated in Python and the rasters were created using Python-GDAL.
In all but one ...
1
vote
1answer
40 views
gdal_retile.py and gdal_translate Requirements
I basically need to use two GDAL utilities: gdal_retile.py and gdal_translate, on a Solaris 9 machine. I have tried compiling GDAL from sources, but this seems not to be working. So, since I have ...
2
votes
1answer
35 views
gdal_retile.py tile naming scheme
Is there a way to define custom tile names during a gdal_retile.py process, other than the default _01_01. For instance, suppose I want to start naming my tiles at 100_100, is this possible?
0
votes
1answer
44 views
gdalwarp cant find libproj.dlib
Thought I had finally got Python, gdal and QGIS to play nicely together but just tried to do a batch transform using a bash script to call gdalwarp and it claimed it couldn't read the libproj.dlib ...
0
votes
0answers
22 views
Offset shapefile to start from prime meridian
I have Natural Earth shapefile (for example this file: link) covering whole globe in -180:180 range.
I want to offset this file, so that it starts from Greenwich prime meridian instead -180.
How can ...
2
votes
1answer
81 views
Export KML ExtendedData with OGR2OGR
I would like to export with ogr2ogr the data I have in a kml as ExtendedData
<Placemark>
<ExtendedData>
<SchemaData schemaUrl="#15540_punti">
...
0
votes
1answer
25 views
Error Updating Gdal with Homebrew
I recently updated my version of Gdal using the homebrew package manager. However, now whenever I use gdal from within my terminal shell, I get the following error:
dyld: Library not loaded: ...
5
votes
2answers
345 views
Tiling high resolution image using GDAL tools
I have a .tif image that I georeferenced using QGIS utility. It was originally in UTM with NAD83 as datum.
I need to make tiles out of it and get a kml file to overlay onto google earth.
Here's how I ...
2
votes
1answer
32 views
Where can I find a python GDAL 1.8.0 package?
I have python 2.7.4 installed on my RHEL 5 machine along with gdal-devel-1.8.0 and gdal-1.8.0. I tried to install the GDAL 1.9.1 package available on python's page ...
1
vote
1answer
213 views
Slow tiles generation for leaflet js
I am currently trying to use leaflet to display images (jpeg or png) instead of maps, I have followed the instructions on this blog post to create the tiles using gdal2tiles.py here is the code we ...
2
votes
1answer
51 views
How to sort features in a layer using OGR/Python?
How can I sort the features of an ogr layer by the values of a field?
I tried naive variations of
def sortlayer(ds, layer, field):
lname = layer.GetName()
lsort = ds.ExecuteSQL(b'select * ...
1
vote
0answers
23 views
gdalwarp not clipping neatline properly
I have a PDF file (10mb) which contains a NEATLINE polygon in metadata.
$ gdalinfo I5601_GeoPDF.pdf
Driver: PDF/Geospatial PDF Files: I5601_GeoPDF.pdf Size is 5906, 3543
Coordinate System is: ...
1
vote
1answer
29 views
re-project from EPSG:4326 to polar projection
I have global georeferenced DEM tiff file in EPSG:4326 and want to re-project in EPSG:3411. I used this command:
gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3411 global.tif polar.tif
which resulted in:
...
0
votes
1answer
250 views
Installing GDAL with Java Bindings (gdal.jar) on Ubuntu 12.10
I am trying to build a postgis importer for mutliple shapefiles. After geotools took orders of magnitude longer than shp2pgsql to import my data, I decided to attempt to use the java bindings to GDAL ...
1
vote
1answer
26 views
python gdal_polygonize error
I am trying to use gdal_polygonize to create a polygon from a binary raster.
My code is:
os.system('gdal_polygonize.py ' + filename + ' -f ESRI Shapefile ' + shapefile)
where filename is a string ...
2
votes
1answer
91 views
Driver Error : DTM to Raster using GDAL_GRID
I am looking to interpolate a DTM (eventually over 4,000 DTMs) that is in a personal geodatabase format. I would like to do this in as few step as possible, eventually writing a batch script to do ...
2
votes
1answer
122 views
gdal_retile question
I have been trying to get gdal_retile working with the --optfile command. At first it did not work because in my files.txt I have '-' characters in the filename. I tried adding quotes around the ...
0
votes
1answer
83 views
GeoTIFF: How to remove Internal Overview without recompression
"gdaladdo -clean" does not remove the internal overview from the file see http://trac.osgeo.org/gdal/ticket/4940
How do I remove GeoTIFF overview(s) without having to re-compress the main image.
9
votes
7answers
880 views
How to add field with filename when merging shapefiles with ogr2ogr?
I'm merging some shapefiles and I had some problems doing so inside QGIS, so i'm using ogr2ogr directly.
I'm doing this (in a batch):
ogr2ogr -overwrite %destination% %n1%
ogr2ogr -update -append ...
2
votes
1answer
18 views
GDAL2Tiles: MapTiles from BSB/KAP are Switched
As an example I have a BSB file from NOAA which is a navigational map. The map has a SRS of WGS84 and a Mercator projection.
I initially translate the file from BSB format to GTiff format using the ...
1
vote
3answers
360 views
Geoserver - ImagePyramid plugin problem
I have a huge Geotiff file , and i've tiled it ..
I've created the tiles using gdal_retile.py tool , but i'm getting the following error from Geoserver while linking the layer with the tiles folder : ...
1
vote
2answers
46 views
How to convert self intersecting polygon to multipolygon?
Working with the union and polygonize methods in GDAL/OGR, I often end up with self intersecting polygons. The case is one polygon which essentially consists of two polygons:
In the above example, ...
0
votes
2answers
90 views
How to install Shapely, but got stuck into Python language
I have a MAC OSX 10.7.5 and am trying to learn how to use QGIS.
Following a couple of tutorials I came to know I have to install shapely.
I followed the instructions by the multi-cited ...
4
votes
1answer
256 views
How to create tiles with points from MongoDB?
I have a (potentially) large data set stored in MongoDB. Each record is just lat, lon and some float-point value. Each pair of lat and lon is unique, but the whole set is not covering the Earth, there ...
0
votes
0answers
49 views
ogr.GetDriverByName('FileGDB') returns None
I installed:
gdal-110-1600-x64-core.msi
GDAL-1.10.0.win-amd64-py2.7.msi
gdal-110-1600-x64-filegdb.msi
from www.gisinternals.com
I can open my gdb on the command line:
C:\data>ogrinfo ...
3
votes
1answer
92 views
OGR Info - Check for spatial index?
If we want to get information about a vector layer (shapefile for example) we could use ogrinfo:
ogrinfo -so myshapefile.shp mylayer
This returns a lot of pertinent information about the shapefile. ...
3
votes
1answer
101 views
Python - OGR: Distance different than in ArcGIS
I want to determine the shortest distance between a point and a line. I use the following Python script to do so. If I compare the result of the script with the result of ArcMap's Near Analysis, I get ...
0
votes
0answers
20 views
Build overviews before building tiles?
Does it make sense to create overviews (with gdaladdo) for my GeoTiffs (about 12000x15000 px) before creating tiles in GeoWebCache? Does resampling affect the quality of my tiles?
1
vote
1answer
35 views
Error output with gdalwarp
Following is my raster file and mask shape file
files
I use qgis 1.8.0 to clip the raster with shape file by gdal tools
but after I run warp in gdal tools it's just give me no data in the output ...
1
vote
2answers
67 views
How can I convert a ascii file to geotiff using Python?
I have a ascii raster format file. For example:
ncols 480
nrows 450
xllcorner 378923
yllcorner 4072345
cellsize 30
nodata_value -32768
43 2 45 7 3 56 2 5 23 65 34 6 32 54 57 34 2 2 54 6
35 45 65 34 ...
1
vote
2answers
726 views
Clean way to install gdaltools plugin for QGIS on OSX
I have qgis 1.8.0 on osx mountain lion and am installing the gdaltools 1.2.29. I have run into a problem that a number of people are seeing. The plugin cannot find my gdal python (see below).
...
1
vote
1answer
69 views
How to use GDAL-commands in QGIS?
I recently used a GDAL-command working with the OSGeo4W command line, but it does not feel that comfortable. Is there a way to somehow execute a GDAL-command directly in QGIS?
Right now the ...
0
votes
0answers
21 views
How to check are there any OGRGeometry (polygons) within an OGRGeometry (polygon) ?
Please see my image:
I used the below statements to get Geometry:
OGRGeometry *poGeometry;
poGeometry = poFeature->GetGeometryRef();
Do you know how to check whether the Polygon 1 ...
0
votes
1answer
33 views
Can't convert geometry from Dgn file to Shape file?
I am using below codes to try convert geometry from a Dgn file to Shape file. But, I can't exit the Feature reading loop. Are there anyone know why and how to fix it ? Please download the dgn file ...
0
votes
1answer
24 views
GDAL clips raster at wrong position
I try to clip raster data with GDAL (in fact I am using QGis, which produces me a command line for GDAL)
The raster has been transfored to EPSG:3857, which is also the projection of the shape file ...
1
vote
0answers
30 views
Why does the Clipper tool lack the mask layer option?
Fairly new at QGIS, but everything's gone well so far.
I'm running QGIS version 1.8 on a Mac (OSX 10.8) and GDAL 1.10, and other GDAL tools (like Merge and Georeferencer) work just fine.
When I ...
0
votes
0answers
21 views
How to retrieve PROJECTION value for a map?
I am trying to read the Projection from Geotiff/shape file, using gdal 1.9.0 by the following command in the command line...
gdalinfo.exe -proj4 map-path
but its retuning empty string for ...
1
vote
0answers
154 views
ArcGIS 10.1 crashes when a button in an addin invokes GDAL DLL through SWIG CSharp wrapper
I am developing an addin with a button that tries to do coordinate transformation using GDAL OSR. When I create a new CoordinateTransformation object ArcMap crashes and close. If a debug against a ...
0
votes
2answers
52 views
Why gdalwarp cannot make transparent black nodata?
I have a MODIS ariel image geotiff and I want to remove black nodata pixel (this is Aerosol Of Thickness) 1 band image.
I've done like this C:>convert -transparent black ...
1
vote
1answer
51 views
Technique for splitting raster image into tiles
I have written a standalone application which renders a globe similar to GoogleEarth/WorldWind which is currently able to display raster imagery from WMS servers. The native geodesy of this ...
1
vote
0answers
22 views
How can I combine these two commands? (Get a clipped raster ASC-II file from PostGIS)?
My objective is to clip a polygon from PostGIS raster table, and then save this raster (draw polygon shape) as an ASC-II in a folder. What I did now are 2 steps.
1st Step: Insert the clipped raster ...
0
votes
0answers
33 views
Why is GDAL ECW compression much poorer than ECW compression from ERDAS
We use ECW compressor shipped with ERDAS Imagine to export Orthoimagery for our different departments. Because we are facing problems with the Imagine license server, we had to choose GDAL (built with ...
0
votes
1answer
128 views
How to clip raster by multiple polygons in multi rasters?
Using QGIS, I want to clip a raster by multiple polygons contained in a single shapefile. I read that gdalwarp can clip my raster by multiple polygons BUT the output is a single merged clipped raster. ...
3
votes
0answers
28 views
What is the correct syntax for using gdal_rasterize with a PostGIS datasource
I have a multipolygon layer in PostGIS and I would like to use gdal_rasterize to generate a raster. I cannot find the correct syntax to specify PostGIS as the data source.
$ gdal_rasterize -a ...
5
votes
1answer
217 views
GDAL.RasterizeLayer not ported to the C# bindings?
I am trying to use the OSGeo.GDAL.Gdal.RasterizeLayer(…) from GDAL 1.9.2 and it requires some SWIGTYPE parameters:
OSGeo.GDAL.Gdal.RasterizeLayer(
OSGeo.GDAL.Dataset dataset,
int bands,
...
1
vote
1answer
45 views
gdal_merge.py with single band raster
I'm trying to mosaic the True Marble world imagery 250M geotiffs using gdal_merge.py though one of them doesn't have RGB bands like the rest and gdal_merge fails (same number of bands is listed as a ...

