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

Driver error when converting MapInfo file to Shapefile with ogr2ogr

I am trying to convert a MapInfo file to an ESRI shapefile using GDAL's ogr2ogr command as described here, but I am getting the following error: FAILURE: Unable to open datasource `Iceland2pts.mif' ...
0
votes
1answer
185 views

What's the best way to add a projection to GDAL?

I have a number of files in Maryland State-Plane (US foot) that have the following wkt: "PROJCS["NAD_1983_StatePlane_Maryland_FIPS_1900_Feet", GEOGCS["NAD83", DATUM["North_American_Datum_1983", ...
3
votes
1answer
318 views

LSA-SAF data reprojection with GDAL and R

The Land Surface Analysis Satellite Applications Facility (LSA SAF) distributes several products derived from the EUMETSAT satellites data. I am working with R to read and process the data (the file ...
1
vote
2answers
192 views

Some tiff not included on a vrt index

I have some tiff files and I'm trying to generate a vrt index with gdal tool, gdalbuildvrt. I run this command: gdalbuildvrt orto_index.vrt -srcnodata "0 0 0" -vrtnodata "0 0 0" *.tif and throws ...
2
votes
2answers
110 views

gdal_grid unable to get nodata outside of original data area

I have a vector based csv file that I want to convert into a grid using gdal_grid, but I want to force null data values strictly around the area where there is data. I can get gdal_grid to work but ...
1
vote
1answer
569 views

importing osgeo module in ubuntu 10.04

I have problem with importing "osgeo" module in ubuntu 10.04; I have installed GDAL 1.7 from source ( 1)./configure 2)make) however, still in python I cant import osgeo module. I appreciate any hints. ...
2
votes
1answer
262 views

Black zones after tiling two tiff files

I have two Geotiff files. with one side overlapping. Each Geotiff is surrounded by a black frame. I have created an Image Pyramid with Geoserver 2.3.0 but when I try to see the layer on OpenLayers, ...
1
vote
3answers
265 views

Center of a bounding box

I have bounding box coordinates of my shapefile (fetched from layer.GetExtent() Python with GDAL/OGR): (-38.6102467739459, -38.017601026186576, 33.01563382506067, 33.624945228091406) So I assume ...
1
vote
1answer
39 views

Michigan Projection with GDAL or FWTools

Has anyone come across a specific projection for Michigan, specifically: NAD_1983_CORS96_StatePlane_Michigan_South_FIPS_2113_Ft_Intl When I look at spatialreference.org I do not find any instance of ...
2
votes
1answer
113 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
160 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 ...
3
votes
0answers
120 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
97 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 ...
1
vote
4answers
268 views

What vector formats are compatible between GDAL and ArcGIS (without Data Interoperability extension)?

What vector formats can GDAL write that ArcGIS for Desktop can read without the need for its Data Interoperability extension? I do not plan to use either Personal Geodatabases or Shapefiles
1
vote
1answer
286 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: ...
4
votes
1answer
303 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 ...
1
vote
0answers
115 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: ...
3
votes
1answer
213 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
1answer
143 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: ...
2
votes
1answer
1k 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
2answers
758 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
673 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 ...
2
votes
0answers
330 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 ...
1
vote
1answer
71 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 ...
0
votes
0answers
148 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
0answers
51 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 ...
0
votes
1answer
176 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 ...
1
vote
1answer
223 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 ...
1
vote
0answers
142 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 ...
1
vote
1answer
869 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
1answer
82 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 ...
5
votes
3answers
2k 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 ...
4
votes
0answers
644 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 ...
6
votes
0answers
401 views

What is the correct syntax for using gdal_rasterize with a PostGIS datasource [closed]

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 ...
3
votes
1answer
355 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 ...
0
votes
1answer
477 views

How to export only one band from an image using GDAL?

I have a multi-band GeoTiff, and I want to extract only the first band and write a new image consisting of only that band. How do I do this using GDAL?
1
vote
1answer
377 views

Python GDAL/OGR: WGS84 & meters units

I'm dealing with shapefiles projected in: GEOGCS["GCS_WGS_1984", DATUM["WGS_1984", SPHEROID["WGS_84",6378137,298.257223563]], PRIMEM["Greenwich",0], ...
2
votes
0answers
437 views

Use GDAL/Python to add rasters (tiff) to an ESRI File Geodatabase

I'm using GDAL python to create a bunch of tiff files from a netcdf file. def modelout_to_geotif(netcdfile, outdir, VarName): """ Given the path to the netcdf output, produces geotiffs for ...
1
vote
1answer
780 views

Python GDAL: degrees to meters without reprojecting

I am trying to get senseful results in my Python script. I want to measure distance between features in a shapefile. I'm working with the SRS: GCS_WGS_1984, so my units are in degrees. But for the ...
2
votes
1answer
394 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
1answer
1k views

How to transform World data from EPSG 4326 to EPSG 3857

I have tryed different opensource tools to transform a world data shapefile from 4326 to 3857 but I always have an error. Im not sure, but is it because of Antartica? I have tryed with GDAL and QGIS ...
3
votes
0answers
61 views

How do I get the domain from an OGRSpatialReference?

I would like to know what the projected bounds are for a spatial reference using OGRSpatialReference. Is there some way to use the GetProjParm method to get this information? An example of what ...
0
votes
1answer
103 views

What's regular blocking arrangement ? (GDAL)

I am using GDAL tools to retrieve PostGiS raster. The code(command line) I am using is as follows : gdal_translate -of AAIGrid "PG:host=localhost dbname='database' user='postgres' password='admin' ...
1
vote
1answer
209 views

How to set a spatial filter with Python/GDAL?

Is there a way with GDAL (Python) to set a spatial filter (bounding-box) on a raster before I process it? I want to read it as an array. BandReadAsArray(band, xoff=0, yoff=0, win_xsize=None, ...
1
vote
0answers
34 views

GDAL Plugin with a different version of GDAL aready built into QGIS

My qgis windows version is 1.8.0 Lisboa and use the gdaltools in version1.9.2. I want to use the plugin NetCDFBrowser, but this plugin need the gdal 1.10. I install 1.10 stand-alone gdal find here ...
3
votes
2answers
940 views

Getting started with gdal/ogr

I found some batch code using ogr2ogr in order to convert and reproject vector layers. The difficulty for me is to know what to install and which gdal commandline should I use, and how to configure ...
2
votes
1answer
133 views

Using the C++ GDAL API how do I set the GDAL_DATA path?

If I want to change the directory of the GDAL_DATA path or point to a different version of a csv file for EPSG numbers how do I do that using the C++ API?
1
vote
2answers
548 views

How to calculate Image boundary / footprint of satellite images using open source tools?

I'm referring to this question&discussion: How do I create a shapefile showing footprints of Rasters? I need to create the polygon outlines of several single layer raster images, not the ...
5
votes
1answer
1k views

Python - OGR: Transform coordinates from meter to decimal degrees

I have a point in a shapefile. The shapefile is in a WGS_1984_Web_Mercator Projected Coordinate System. I want to get the coordinates of that point in decimal degrees with GDAL/OGR in Python. I am ...
-5
votes
1answer
145 views

Problem installing GDAL and fiona module in arcmap10? [closed]

I am getting message like this Traceback (most recent call last): File "", line 1, in File "C:\Python26\ArcGIS10.0\lib\site-packages\gdal.py", line 2, in from osgeo.gdal import ...