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

0
votes
1answer
215 views

Better way to handle HDF and GeoTIFF files

I am using GDAL and MODIS reprojection tool in my python project and want to know better way for math new indices of HDF data (MODIS data). Also I want to visualize them in web. Which format is more ...
2
votes
0answers
251 views

Buffering a Raster with GDAL and Numpy

I am using GDAL to convert a raster dataset into a numpy.array and then using numpy to buffer the data in the raster. For example, I need to know all the regions within 100ft of a school point ...
1
vote
1answer
381 views

Problem with GDAL_WMS QGIS connecting to an ArcGIS Server?

I've been working, trying to get QGIS (v2.0 32bit) to connect to an ArcGIS Server service using the GDAL_WMS xml file. This is the service I'm trying to connect to: ...
8
votes
2answers
768 views

Building FGDB support for GDAL 1.10 on Ubuntu

Because data.gc.ca is now distributing some of its data in FGDB format (example: Plant Hardiness Zones), I would like to be able to use it with gdal/QGIS. The compile/install is not going as planned, ...
2
votes
1answer
855 views

GDAL loads in Python, but not IPython (Linux)

I am trying my hardest to utilize the work of the good folks at the GeoDa center, but I cannot seem to get all of the components to load. GDAL, in particular, has been giving me some difficulty. The ...
1
vote
1answer
1k views

QGIS 2.0 Plugin error: Unable to load GdalTools plugin?

I just used the osgeo4w exe to install QGIS and all the other software included in the basic setup. When I open QGIS, I get the following error: Plugin error Unable to load GdalTools plugin. The ...
1
vote
1answer
115 views

Reading non-spatial Access tables with ogr2ogr

I am trying to read a whole Access database into Postgis using the ogr2ogr (Version 1.9.1) command. For some reason however, it only reads spatial tables and no non-spatial tables. All my attempts to ...
1
vote
0answers
89 views

How to join attributes by location with command line

I would like to joint the attributes of 2 shapefiles by location but with a command line (linux). This feature is easy in QGIS but it is VERY slow. Is there a way of doing this with ogr or pyqgis ? ...
0
votes
1answer
178 views

how to project a flat image to spherical mercator EPGS:900013 to use in a Image Layer in openlayers

I have a flat 400x400 pixels image which represents an area of 400x400 kilometers. I need to draw this image as an overlay on my OpenLayers mapping system. I use OpenStreetMap/Google Maps as base ...
2
votes
1answer
274 views

Problem with georeferencer in Qgis 2 64bit installation

I tried to use the georeferencer in the new Qgis 2 64 bits for windows. I tried with the standallone installer and with the OSgeo4w installer. The problem is when I load the raster that I want to ...
10
votes
6answers
916 views

How do I iterate through every cell in a continuous raster?

See this link for more details. The Problem: I want to loop through a continuous raster (one that has no attribute table), cell by cell, and get the value of the cell. I want to take those values ...
0
votes
2answers
736 views

How to extract subdataset from HDF raster?

I am trying to extract a subdataset from HDF raster. Later on, I intend to make operations such as mosaicking and extracting raster averages from the extracted subdataset. The following code is ...
14
votes
1answer
339 views

Why do people pronounce GDAL (Geospatial Data Abstraction Library) like goo-dal?

Ok, a bit of a silly question, but: I've heard people pronouncing the library GDAL (Geospatial Data Abstraction Library) like goooo-dal. My natural inclination was to pronounce it geeee-dal, as in ...
1
vote
1answer
141 views

GDAL reports no bands for a 12-band raster and no projection either

When opened in ArcGIS, the raster correctly unfolds 12 bands, but when I try to open it with GDAL through Python, it appears as a raster that has zero bands. Here is the code: import gdal from ...
1
vote
1answer
73 views

Which kind of grid format is this?

in a project I got some raster files with the following content (header and data are in one txt file): HOUR 17 MINUTES 30 YEAR 2012 MONTH 5 DAY 4 CITY Emilia;Bologna NROWS 99 NCOLS 156 ULXMAP ...
0
votes
1answer
113 views

gdallocationinfo is not recognized

I tried to run gdallocationinfo from my FWTools shell and got a message saying: "gdallocationinfo is not recognized as an internal or external command, operable program or batch file". So I checked ...
3
votes
1answer
406 views

Clip and copy files to same directory structure using OGR/GDAL with Python

We need to clip all our data (many TB) by a new boundary and save it while maintaining the same data structure. I want to use OGR/GDAL as the data is in shp/tab/kml/geoTif/ecw etc. I can find the ...
2
votes
2answers
252 views

MODIS calculations with non cloud contaminated pixels

I have 11 files showing 8day daytime LST temps in .tif format from MODIS. Is there any way to calculate the average 8day daytime temperature for march 2013 for example but using in the calculations ...
1
vote
2answers
180 views

will ogr get coordinates in geographic coordinate system or projected coordinate system when using getPoint(int,OGRPoint*)

Consider the code that I have written in the stackoverflow question How to extract vertexes of geometries in ESRI shapefiles using ogr library with c++, I have used the code to show shapefiles on an ...
0
votes
1answer
74 views

Where do I find the code for gd_calc.py?

I need to create a raster that is the average of a large set of other rasters (many thousands). I've learned from another question that the script gd_calc.py can do a limited version of this, limited ...
0
votes
2answers
170 views

How to average a set of n rasters using gdal or other GRASS module? [duplicate]

I need to average a set of n rasters using gdal or some other GRASS module. What command exists for this purpose?
1
vote
1answer
327 views

How to convert 4-channel PNG files to grayscale geoTIFFs?

I suspect that there exists a gdal tool (gdal_translate?) that will convert a 4-channel PNG (+ world file) to a 1-channel (gray-scale) geoTIFF. From the gdal_translate man page, it appears simple to ...
2
votes
1answer
111 views

GDAL local files VS ftp files

I can successfully get raster information when the raster file is located on my local file system: import gdal from gdalconst import * filename= "MOD11A1.A2012193.h11v11.005.2012196013534.hdf" ...
-3
votes
1answer
79 views

How can I use ogrinfo to reach information anout a .shp?

How can I use ogrinfo? where should I write it? In the visual studio's command prompt or somewhere else? Can you make some examples? when I use ogrinfo in the windows or visual studio's command ...
2
votes
1answer
67 views

What does the predictor switch do while compressing geoTIFFs in GDAL

I am merging several geoTIFFs into one and am using the -co COMPRESS=DEFLATE create option for geoTIFFs. An additional switch is the predictor switch. What does it do? The GDAL website is not clear ...
1
vote
0answers
81 views

Contours generated from merged TIFs huge - much bigger than individually

I'm generating contours from four 5x5 degree SRTM DEMs in southeastern Australia. I'm using gdalwarp to reproject the DEM, then gdal_contour to generate contours. First, I generated contours for each ...
2
votes
2answers
368 views

How do I know the vertical/horizontal units of my DEM? I'm creating a slope using gdaldem

I'm using gdaldem to create a slope map. I can apply a scale. The default is 1. The notes say: Ratio of vertical units to horizontal. If horizontal unit of source DEM is degrees (e.g WGS84 ...
2
votes
0answers
96 views

Displaying DEM raster in QGIS

Using Raster>Analysis>DEM to display hillshade produces the message "gdaldem.exe had stopped working" so I asked the community. I'm new to QGIS, and I now realise after a little experimentation that I ...
0
votes
0answers
231 views

How to read vector datasets using GDAL library

Up to now, I have used GDAL in order to read raster information and geotiffs and show them with the help of OpenGL. Now I have a vector dataset, and by that, I mean I have the files with types: DBF ...
1
vote
0answers
507 views

How to generate JPEG tiles using gdal2tiles.py

I am using gdal2tiles.py to generate tiles and packing it into MBTiles using mb-util. gdal2tiles by default generates PNG tiles. When I convert PNG tiles to JPEG large black area around tiles are left ...
2
votes
2answers
116 views

Invert gdaldem aspect

I am trying to identify lakes and other water bodies in a DEM file. The source file is located here. I am using gdaldem aspect to identify contiguous pixels with zero slope change. gdaldem aspect ...
2
votes
1answer
258 views

KML to shape QGIS conversion [duplicate]

I am trying to import a KML file to QGIS. The KML source file is about 15mb is comprised of polygons with attribute data. The KML imports fine and polygons display but the attribute data is unreadable ...
7
votes
2answers
1k views

Convert timeseries stack of GTiff raster to single NetCDF

Moving from gdal-dev mailing list: On Mon, Sep 2, 2013 at 7:09 PM, David Shean wrote: Hi list, I'm trying to package a timeseries of GTiff rasters with identical projection/extent/resolution as a ...
1
vote
0answers
910 views

How do I change from lat/long to a metric coordinate system in GRASS?

Yes, it's yet another transformation question. I am working with ASTER DEM data, which is LatLong WGS84. I want to process the DEM map into an STL file for use in a CFD simulation. The script I am ...
8
votes
1answer
274 views

Python with GDAL and PostGIS as a replacement of ArcGIS

Few times ago, I worked with daily land surface temperature satellite images of an area with the size of Alaska. I had to process around one thousand images. I downloaded them, mosaiced, projected, ...
1
vote
3answers
474 views

How to use gdal module in python

I am trying to import gdal modules from osgeo installed in my win64x machine but I get the following error: Traceback (most recent call last): File "", line 1, in from osgeo import ogr ...
1
vote
0answers
462 views

How to properly convert EPSG:3857 to EPSG:4326 tif?

I tried using the answer from "Transforming EPSG:3857 to EPSG:4326" to do this but gdalinfo tells me that this didn't really work properly. Driver: GTiff/GeoTIFF Files: geo-map.tif Size is 1024, 512 ...
1
vote
2answers
2k views

No ogr module despite having gdal installed?

I am trying to use the mapping component of the wonderful PySAL package, but for some reason I do not have the ogr module. I was under the impression that it came packaged with gdal. Am I missing ...
2
votes
0answers
388 views

GDAL: Geotransform - rotated rasters (and ArcGis)

I am working on a GDAL driver for a proprietary grid-format. In some cases, the grids are rotated in order to get smaller file sizes. In my GDAL geotransform array, I do the following (showing a ...
0
votes
1answer
53 views

USGS FLT Hdr file interpretation

I have a legacy program that attempts to read USGS 1/3 arcsec data. When I compare results to GDAL, I find that the lat/lon to pixel conversion is off by 1 line. The reason (I think) is because the ...
2
votes
2answers
168 views

OSM linestring shown as point

I have loaded the following file http://openstreetmapdata.com/data/coastlines to postgreSQL. In visual studio with gdal/ogr and geos enabled: SELECT geom FROM waterpoly WHERE geom IS NOT NULL LIMIT ...
1
vote
0answers
184 views

GeoPDF Support in GDAL

I installed OSGeo4W using a Windows installer... it contains GDAL 1.7.3. I can't get GDAL to recognize PDF as a file type using commands in MSYS shell. Does this version of GDAL not support PDFs? ...
1
vote
2answers
517 views

gdal_calc.py outputs huge files

I think this is a general enough question to apply to more than just gdal_calc, but perhaps not. When I run gdal_calc.py -A map.tif --outfile=deforestation_00-10.tif --NoDataValue=0 --calc="A >= ...
3
votes
1answer
126 views

is it possible to project a shapefile to a raster's custom projection with command line gdal/ogr?

I have a raster with a user defined sinusoidal projection and a shapefile in lat/lng coordinates (I believe this is called 'unprojected'?). I'd like to use ogr/gdal's command line utilities to ...
1
vote
0answers
76 views

How can I maintain MapInfo GRC file attributes in QGIS/GDAL?

I am trying to manage a classified MapInfo grid file in QGIS - I'm wondering if anyone has any experience with maintaining the attributes of the classified grid? QGIS seems to assign them just a ...
3
votes
1answer
284 views

Best way to visualize the PostGIS Raster in OpenLayers

I have been stuck by this problem for a long time. Several ways have already been tried, but none of them works as I was expecting. Here are the description of each of these method as well as the ...
4
votes
1answer
588 views

Converting WMS getmap request to Geotiff

I'm accessing a WMS server that doesn't provide geotiff as one of the format capabilities, but does provide PNG and KMZ. ...
2
votes
1answer
816 views

Polygonize.py outputs into a single file or simple merge?

Given an unique raster topographic file at start, I want to get vector polygons for elevations thresholds 1m, 200m, 500m, 1000m, 2000m into a single shp file. I first slice my raster into 5 (or more) ...
0
votes
1answer
377 views

Polygonize.py: AttributeError: 'NoneType' object has no attribute 'CreateDataSource'

Given a topographic raster tif of an european country where elevations are simplyfied to either 0, either 200. How to create a polygon ? Comment: With the following polygonize.py command: ...
0
votes
1answer
340 views

How to conditionnally assign a new value to pixels of a raster image?

Given a topographic GIS raster of one country crop.tif: Given a know pixel value such as elevation (z) is +73 : $ gdallocationinfo crop.tif 1 1 -valonly > 73 Given an elevation threshold n = 50 ...