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
0answers
121 views

Is Django's LayerMapping capable of putting a static value into all features?

I'm using GeoDjango's LayerMapping utility to try and import a set of features. The catch is that aside from the attributes in the source layer, I want all the destination features to have the same ...
3
votes
1answer
122 views

rgb2pct.py increases size of file 20 fold!

I wanted to reduce the number of bands of a tif file from 3 to 1... I used this command rgb2pct.py -n 256 fnp_braunschweig.tif fnp_bs.tif and it increased the size from 31mb to 778mb I actually ...
7
votes
3answers
1k views

Poor image quality of raster pyramid files in ArcGIS (no anti-aliasing)

There are several GeoTIFF files (size ~7000x7000 px) which I try to open in ArcGIS Desktop 10.0. I have already created raster pyramids for them (.ovr files) using geoprocessing tool but their quality ...
1
vote
1answer
384 views

Gdal tiling error for EPSG: 3857

i have lots of data which was coordinated with gdal some programatic way at below: gdal_translate -of GTiff -a_ullr %s %s %s %s -a_srs EPSG:3857 %s/%s %s/%s.tif then for tiling them all i am ...
2
votes
2answers
571 views

How to install GDAL with TIFF support on Linux?

I am trying to install GDAL on a linux machine and tested some sample commands like gdalinfo --version This worked fine, but when I started working with images like .tif I got the following error . ...
1
vote
2answers
1k views

How to convert GeoTIFF to grayscale and add gaussian blur?

I have a giant GeoTIFF I'd like to convert to grayscale and add a gaussian blur to. As far as I know, GDAL can strip two of the three color bands and call it grayscale, which might work, but is there ...
2
votes
1answer
326 views

Converting x, y, zoom for Calculating tile bounds

i am working on calculating tile bounds lat lon from x,y,z tile coordinate. i have looked gdal2tiles code structure and i have found this good function: def TileLatLonBounds( tx, ty, zoom ): "Returns ...
12
votes
2answers
3k views

API documentation for Gdal/Ogr with C#

I want to assign projection to shapefile in my web C# application, for this i am using Gdal/OGR/OSR C# bindings and add osr_csharp.dll and ogr_csharp.dll references. But i am facing some difficulty ...
2
votes
0answers
55 views

Tiled image transforming to new profile

i have a set of tile which was tiled profile of mercator and proper to tms scheme mechanism. unfortunately i dont have any tif or georeference table or anything for this. all was tiled with mercator ...
2
votes
2answers
356 views

gdal_polygonize error with Arc/Info binary grid raster file

I am trying to convert an Arc/Info binary grid raster file to a shapefile using the gdal_polygonize.py utility. But the operation fails with the following error Warning 1: Unsupported Arc/Info ...
0
votes
1answer
353 views

Gdal2Tiles getting tif with dtype = uint8 but want int16

I have a gdal2tiles code for getting tiff file. In my code I am getting true tiff format but it makes my dype uint8. I want int16 for dtype. how can I overcome it? I know I have to add gdal.GDT_Int16 ...
2
votes
1answer
3k views

How to process GRIB2 Data for Visualization?

My goal is to extract wind information from data provided to me in the GRIB2 format as it has been done in the attached image. I was curious if anybody had any experience with this format. I'm not ...
5
votes
1answer
358 views

What is the best way to prepare digital orthophotos for the web?

My system: Linux Ubuntu 10.04 server (100gb: 4gb RAM: 32bit) /tomcat6/geoserver/geowebcache/FWTools My Objective: Create Aerial Photo coverage in ETRS89 N32 My Data: 1300 3 band Orthophotos in Tif ...
1
vote
1answer
500 views

Convert MGRS using gdaltransform

I'm using gdaltransform to convert coordinates. What's the best way to convert MGRS using gdaltransform given that it doesn't have an EPSG code?
2
votes
0answers
441 views

GDAL tiling DTED (dt2) dataset as gdal2tiles for 3D Maps

i am searching some webgl technologies to create 3d map with digital elevation data. for web tech. i have found ReadyMap SDK to create 3d map from TMS... in ReadyMap for elevation data it uses sth ...
4
votes
2answers
331 views

How to import GDAL in Jython?

I'm trying to use jython with GDAL on Windows. Python and GDAL are installed through MS4W and work perfectly together. Jython is installed manually, but I can not import GDAL. I tried to update ...
1
vote
1answer
225 views

GDAL Tiling, Missing Tiles

i have used gdal2tiles.py for tiling about 100 tiff tiles which i have generated virtual dataset with gdalbuildvrt for my 1:250.000 scaled maps... everything goes well from 15 to 11, but when it ...
0
votes
1answer
940 views

gdal2tiles format of argument for the '--srcnodata' flag

When you run gdal2tiles, by default it sets pixels outside your map area to black-transparent. I would like to set them to white opaque. Looking at the documentation... ...
3
votes
2answers
1k views

How to fix for the GDAL clipper tool in QGIS?

I would like to clip a DEM with an admin boundary using qgis, in order to calculate a measure of variation. However, when I use the Raster clipper tools and select mask, I get the error message: "An ...
21
votes
1answer
2k views

Can ESRI's File Geodatabase API be Redistributed?

ESRI's file geodatabase API can now be used with GDAL/OGR - http://www.gdal.org/ogr/drv_filegdb.html however each individual seems to have to login to ESRI and download the API files individually. For ...
2
votes
1answer
171 views

Define gdalwarp tile size

gdalwarp has -co TILES=YES (256x256 as default) Can I set different tile size? Thanks
2
votes
1answer
618 views

overlap two GeoTiffs to fill missing information (transparency)

I have several hundered geotiff files. Some of them (almost all of them) have black areas (irregular) on the edges. However, these black areas are covered by some other file in the file set. I read ...
3
votes
2answers
814 views

ImageMagick's convert program changes PNG file format so that color/alpha info doesn't survive GDAL

It appears that ImageMagick's convert program changes something in my PNG file that causes GDAL processing to mess up the color and/or transparency bands. My original PNG file gets processed ...
4
votes
2answers
1k views

How to use NASA Digital Elevation model data?

I've heard it is possible to use NASA elevation model to determine approximate elevation at (lat;lng) point. This might be really useful for my upcoming project since Google Elevation API has weird ...
1
vote
0answers
291 views

ValueError when reading a raster band as a numpy array

Here's the problem shown in Python shell: In [32]: from osgeo import gdal In [33]: geotiff = '/Users/calvin/work/luc/proj_public/media/GLO ...
7
votes
2answers
219 views

How to handle interpolation of multiple (quality) tiers of irregular data?

I'm going to try to provide a high level overview of the issue - realizing that more details/info will be necessary - but I don't want to write a novel that seems too overwhelming for people to ...
6
votes
3answers
763 views

How to convert factfinder2.census.gov map data to KML?

I've downloaded a map from http://factfinder2.census.gov/faces/tableservices/jsf/pages/productview.xhtml?pid=DEC_10_SF1_P1 (limited to San Francisco census tracts), and I'd like to convert it to KML ...
13
votes
2answers
2k views

How to smooth/interpolate a raster in python using GDAL?

This is my first post so please feel free to critique where necessary. Background: I am developing in Python and using GDAL from OSGEO to manipulate and interact with rasters and shapefiles. ...
4
votes
1answer
4k views

Transforming EPSG:3857 to EPSG:4326

I am new to gdal and have a simple question. I want to transform an image from OpenStreetMap projection to Nasa Blue Marble projection. What I understand is that OpenStreetMap projection is EPSG:3857 ...
4
votes
1answer
3k views

Calculate lat lon bounds for individual tile generated from gdal2tiles

I have multiple tile sources generated using gdal2tiles I'd like to present on the same map. So when presenting a tile I need to determine which source to serve it up from, checking against the ...
4
votes
2answers
446 views

How to fix a raster with wrong Eastings?

I get an error when trying to warp a tif file in QGIS: ERROR 1: Too many points (35721 out of 35721) failed to transform, unable to compute output bounds. I am using QGIS 1.7.0 Wroclaw with gdal ...
6
votes
2answers
2k views

How can I turn a shapefile into a mask and calculate the mean?

I have been using python to process huge arrays stored as netcdf files. I would like to calculate the mean of an area defined by a shapefile. I have just installed GDAL but if there are other tools I ...
8
votes
2answers
407 views

How to add different sized rasters in GDAL so the result is only in the intersected region

I'm writing a Python method that adds two rasters and generates a single raster output. For reasons beyond my control, the extents of the input rasters are different, but they do overlap. Is it ...
2
votes
2answers
156 views

Tilecache positions my gdal_translated map 20km south of where I want it

I have created a geotiff file through gdal_translate: gdal_translate -a_ullr 547812.2384341886 6834734.336439242 547967.5011478929 6834622.35081697 -a_srs EPSG:900913 pdz3.jpg pdz.tif When I load ...
11
votes
3answers
4k views

Clipping raster with vector layer using GDAL

I have installed GDAL using the Osgeo installer. How can I clip a raster layer with a vector layer programmatically? Is there any GDAL API available which can help me with this? I am using Python.
5
votes
1answer
3k views

How to import a raster into PostGIS?

I'm trying to follow the instructions here to load a raster into a PostGIS database: python raster2pgsql.py -s 4269 -I -r *.tif -F myschema.demelevation -o elev.sql I understand that I replace * ...
9
votes
2answers
664 views

GeoTIFFs made in GDAL don't have min/max range when loaded in Arc/QGIS

I'm creating some floating point rasters in GDAL with the GeoTIFF driver. When I load the resulting images into QGIS or Arc the default symbology puts the min at -3.40282e38 and max to 3.40282e38, so ...
2
votes
1answer
473 views

Simple band arithmetic on GDAL virtual rasters (VRT)

I want to calculate the difference of two bands. I have the individual bands stored somewhere, but I was thinking that maybe there's a way of producing a VRT file that is the result of this ...
2
votes
1answer
155 views

ArcGIS Geoprocessing Script - Calling External EXE Works in Desktop - Causes GP Service to Quit Abruptly

I have an interesting problem on my hands... I have a generic python script that does a bunch of geoprocessing functions using the standard arcgisscripting module (Yes, I am still at 9.3.1). At one ...
1
vote
1answer
513 views

Seeding Tilecache with gdal2tiles output

Very simple question. I've executed gdal2tiles against a GeoTIFF coverage, composed by several TIFF files. I want to serve the resulting tiles using TileCache. In the TileCache documentation, I see ...
2
votes
1answer
1k views

Tiling rgba geotiffs using vrt rather than gdal_warp

I've carefully cut the collar off a set of USGS geotifs and I'd like to run through them with gdal2tiles but producing a vrt that undersands the alpha layer, or having gdal2tiles overlap two source ...
6
votes
2answers
681 views

Can I change the GDAL/OGR library that QGIS uses?

I have been adding a new format to OGR, all seems well and OGRInfo works. Now I would like to test my code out using QGIS. Is there an easy way to switch which GDAL lib QGIS uses on a windows box? ...
0
votes
1answer
1k views

Import Error while Importing osgeo package [duplicate]

Possible Duplicate: How to install GDAL with Python on windows? I have downloaded Fwtools for python2.5 from http://fwtools.maptools.org(windows installer) after installing it, it was ...
7
votes
2answers
560 views

TRI, TPI, or Roughness

I am working on a site selection projection where I will be unable to perform any in situ observations. A key component of the selection process is the ruggedness of the underlying terrain. I plan ...
13
votes
2answers
1k views

Good practice for generating contours from a DEM

We have just received a large set of DEMs at work and I would like to generate contours from them. The DEMs have a resolution of 1m and a size of 1kmx1km. Output from gdalinfo: Driver: ...
1
vote
0answers
136 views

Geoserver Custom CRS IAU_Moon_2000

I am trying to setup a custom basemap using a GTiff in the IAU's Moon2000 projection (30100 or 30101). The basemap image was originally created through gdal_merge and gdalinfo returns the expected, ...
5
votes
2answers
1k views

Converting Ordnance Survey raster maps to WGS84 webmap tiles

I would like to generate Google/Bing-like map tiles from Ordnance Survey's 1:25000 and 1:50000 raster maps, which are provided as 4000x4000 UK national grid tiles (GeoTIFF). A sample tile can be ...
2
votes
1answer
588 views

gdal2tiles problem

Im using the following to try to produce a tms tileset: gdalbuildvrt test.vrt *.tif gdal2tiles -p raster --zoom=0-6 test.vrt out Two tifs work fine, more than that and it seems to miscompute the ...
4
votes
1answer
437 views

Is it possible to create color tables for floating point raster maps in GDAL?

I've got a script that dumps a floating point raster that I'd like to automatically color. Currently I load this raster in Arc or QGIS and manually create a color table every time I re-run my model. ...
2
votes
1answer
445 views

How to determine area of pixel in meters with GDAL and Python?

I have a raster loaded in Python using GDAL bindings. I'd like to be able to determine the area of a pixel in meters although I'm not confident that the projection is in meters. On this page I see a ...