Tagged Questions

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

4
votes
2answers
953 views

How to use gdal utilities to subset from a raster?

How to you use gdal to subset a smaller raster from a bigger one given the bounding coordinates? EDIT: Can you subset a region from a vrt ?
1
vote
1answer
165 views

gdal_rasterize -at function still won't capture covered pixels properly

I'm currently using the gdal_rasterize tool in bash. My code looks like this: gdal_rasterize -l $(basename $file .shp) -burn 1 -a_nodata 0 -te -180 -90 180 90 -tr 0.1 0.01 -at -ot Byte $file ...
1
vote
0answers
155 views

`gdalbuildvrt` and selecting individual bands

I want to produce a set of VRTs (GDAL's internal virtual rasters) that comprise a re-arrangement of bands of a number of individual GeoTIFF files. Say I have my initial set of GeTIFF files, all having ...
2
votes
2answers
388 views

How to specify gdal_contour output format as postgresql table?

I want to generate contours with gdal_contour command. It generates an output as shapefile by default. I want this utility to output a postgresql (postgis) table. The documentation says there is a -f ...
0
votes
1answer
1k views

How current is GDAL in Java

I have used the ogr2ogr to import various formats to Postgres and have used GDAL utilities. Is there a complete Java port for these tools? Or are they simply wrappers calling JNI to C++ code? I have ...
2
votes
0answers
269 views

How do I use scipy.ndimage.gaussian_filter1d on an arbitrary axis?

I have a GDAL raster that looks like this: And I would really like to blur this raster along an arbitrary transect. I see that scipy.ndimage has a gaussian filter that allows me to blur the raster ...
1
vote
0answers
95 views

Calculating tiles extents using some level data

I have some data looking like this: <?xml version="1.0" encoding="utf-8"?> <levels> <level nr="0" tileWidth="512" tileHeight="512" envSizeX="5000" envSizeY="5000" minX="465000" ...
3
votes
1answer
531 views

ECW/MrSid plugin for GDAL from Kyngchaos not working on Mountain Lion?

Subject says it all. I have installed all of Kyngchaos packages, but gdalinfo --formats does not list ECW or MrSid formats. Could someone else with Mountain Lion try if they can get any JP2 driver ...
3
votes
0answers
764 views

How to create mask bands in GDAL and Python?

I have small program in Python that uses gdal and numpy to perform a logarithmic scale of some data into the RGB 0-255 range. Now this input data has a nodata value set to mask the valid values, and ...
4
votes
1answer
288 views

How to convert a png with known bounds from Plate Carree to Mercator with GDAL

I have a png image that shows the region defined by the geo bounds UpperLat, LeftLon, BottomLat, RightLon (in degrees & -ve for West Lon and South Lat) in Plate Carree EPSG:4326. I want to convert ...
1
vote
1answer
684 views

how to get the projection from a vectorlayer in qgis?

Hi I am trying to set the projection on a raster to match that of a vector point layer. Thus I need to findout what is the projection of a give layer, to use it in the GDAL.Dataset.SetProjection() so ...
3
votes
2answers
2k views

Adding georeferencing data into tiff image

I want to create a georeferenced .tiff file from a .png file. So I convert the .png into a normal .tiff file using the following command line : convert image.png image.tiff Here, all works fine, I ...
4
votes
2answers
1k views

FWTools ogr2ogr in Python

I am trying to convert MapINFO tabs to ESRI Shapefiles using FW Tools. I downloaded the tools and can run the FW Tools Shell and convert no problem: ogr2ogr -f "ESRI Shapefile" ...
2
votes
1answer
159 views

Why does ogr's CreateLayer() return None?

Why does ogr's driver method CreateLayer return None? I'm using python 2.7 with gdal 1.9.1. path = r'D:\_tempAuto' out_shp = '%s\test_OGR.shp' % (path,) driver = ogr.GetDriverByName('ESRI ...
6
votes
2answers
309 views

How to use GDAL to convert Int16 data to Byte

I have some data with Int16 bands in HDF format. I want to convert this to Byte bands in GeoTIFF format. The data range for the Int16 is -2000 to 12000, so I want to map this to the Byte range 0 to ...
2
votes
1answer
619 views

merge raster images: error message procedure not found

I am using QGIS 1.8.0. Whenever I try to merge raster images I get the following error message: Traceback (most recent call last): File "C:\PROGRA~1\Quantum GIS Lisboa\bin\gdal_merge.py", line 36, ...
5
votes
1answer
275 views

Python/GDAL get unique values in discrete-valued raster

I would like to get some advice on the most efficient way to return a list of unique values of a discrete-valued raster using Python and GDAL. I had thought that the most obvious way would be to ...
6
votes
3answers
1k views

How do I extract a raster's extent in python?

I am looking for an easy way to extract a raster's extent in python. Currently, I am using gdalinfo and then parsing out the corner coordinates. The problem with this approach is that gdalinfo is ...
4
votes
1answer
1k views

Qgis - How to change GDAL version

im working on Ubuntu 12.04 server; i noticed this when i installed additional support for ECW format into the GDAL 1.9.1 which worked ok. Then i tried to add ECW files to Qgis but it said they still ...
1
vote
1answer
480 views

Why do I get different results with gdal_calc within the osgeo4w shell and QGIS' raster calculator?

I'm just wanting to make a simple map algebra operation in the osgeo4w shell using gdal_calc but I get wrong results when I divide one raster with the other. When I use the same rasters in QGIS and ...
4
votes
1answer
2k views

GDAL VRT Format - do i understand the concept?

Okay, so im a little fed up with different raster-in-database solutions that I tried over the last 2 weeks which don't work currently (Postgis raster) or cant handle my volume of data (Rasterlite), so ...
8
votes
2answers
3k views

How to fully load a raster into a numpy array?

I have been trying to check my filters on DEM raster for pattern recognition and it is always resulting in missing last rows and columns(like..20). I have tried with PIL library, image load. Then with ...
11
votes
1answer
2k views

How to access feature classes in file geodatabases with Python and GDAL?

I am trying to access a vector dataset in a ESRI File Geodatabase using Python + GDAL. I have successfully compiled GDAL with the file geodatabase API. The FileGDB driver is working correctly since ...
0
votes
1answer
659 views

gdalwarp results in incorrect extent

I have a GeoTIFF that's 172800x86400 on the MODIS "250 meter" Sinusoidal grid (more like 231.7, but who's counting?). (You can find its metadata at the bottom of this post.) This bounty of resolution ...
3
votes
2answers
318 views

Generate a pyramid of tiles by 1/2 scaling a higher zoom level

Given a set of tiles covering a large area at a particular zoom level (These are online on someone's tile server and will need scraping, but let's assume these are on disk if it's easier. With TMS ...
0
votes
2answers
660 views

How can I add a loop in gdal_translate? [duplicate]

Possible Duplicate: gdalwarp in a bash shell for loop I have several modis images (.hdf) from which I want to use only the one band (NDVI) and create a ".tif" image. I have already ...
2
votes
1answer
317 views

GDAL 1.9.1 ogr2ogr KML to PostGIS not incuding ExtendedData

I have a KML file utilizing a Schema in the document and ExtendedData in each placemark. The ExtendedData does not seem to pass in conversion to PostGIS or shapefile. Any assistance would be much ...
3
votes
0answers
360 views

How to compile GDAL with HDF4? (for dummies)

I need GDAL with both Shapefile and HDF4-Support, if possible in 64 bit. Since the normal GDAL-versions from http://www.gisinternals.com/sdk/ have Shapefile, but no HDF4 and the version from FWTools ...
2
votes
1answer
447 views

How to manage null values in QGIS

I made a merge of several AsterGdem blocks near Egypt and Israel. I have portions of the Dead and Mediterranean Seas. Since near Galilea elevations are below sea level this makes the following problem ...
0
votes
0answers
149 views

SAGA grid VRT, file how to apply a LUT?

I've created a heap of saga grid files for some hydro survey work. I'm looking at a way to read these into a mapguide raster/grid layer, I'm able to load up the grids using a vrt file, I can't figure ...
1
vote
2answers
2k 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). ...
8
votes
1answer
2k views

GDAL RasterizeLayer Doesn't Burn All Polygons to Raster

I am trying to burn a shapefile to a raster using GDAL's RasterizeLayer. I pre-create an area of interest raster from a different shapefile, given a specific pixel size. This AOI then serves as a base ...
9
votes
1answer
1k views

Importing GDAL from Python window of ArcGIS for Desktop?

I am trying to run some Python code using the Python window in ArcGIS 10.1 that uses both the arcpy and gdal modules. However, when I try and import the gdal module I get an error: ImportError: No ...
0
votes
2answers
238 views

Why does gdalbuildvrt not allow these things?

I don't know why does gdalbuildvrt not allow: Unreferenced images Rotated images South up images Could someone throw some light? Many thanks in advance
2
votes
1answer
79 views

What should I expect in this case?

If while I'm executing gdalinfo against a raster layer I get a positive value for pixel size Y, what of these scenarios is possible? The raster's map coordinate system sets the origin of the data in ...
4
votes
1answer
437 views

raster algebra in python with rasters of different extents

I am trying to find out how to use GDAL & numpy modules to average a set of rasters....which are Sigma0 values from satellite passes at different times within a year. Each raster has been ...
3
votes
1answer
258 views

Projecting GCP points into translated file

I have file A.png and file B.png. I translate file A.png so i could translate it into B.png (creating multilayer file). Is there any easy way to get new GCP points coordinates so i could combine those ...
1
vote
0answers
234 views

Homebrew GDAL installation ECW or MrSID

I have experimented with installing gdal 1.9.1 via homebrew. It works wonderfully. Unfortunately, their is no real JP2000 support. I realize that JP2000 is supported by OpenJp2000, but that ...
3
votes
0answers
402 views

How to save table from PostgreSQL to MapInfo File?

I have two tables: 1 with geometry in WGS 2 with geometry in Local System I tried this: ogr2ogr -f "MapInfo File" output2.tab PG:"host=localhost user=postgres password=123456 dbname=postgis" -sql ...
2
votes
2answers
424 views

Using osgeo modules in ArcGIS10.1 and QGIS 1.8 python shells

This is related to an old question of mine - Issues with osgeo - qgis.core and arcpy modules Now in new versions of all the software QGIS Lisboa (1.8), ArcGIS 10.1 and pyScripter 2.5.3. They all use ...
0
votes
0answers
510 views

How to add ECW support to GDAL bundled with PostGIS Windows

I'm currently trying the new raster functionalities of PostGIS 2.0 and I successfully imported some GeoTIFF file into PostGIS (although it only find one band so it's all black and white for some ...
0
votes
1answer
413 views

Vector boundary generation -similar to Image Boundary

I want to automate the creation of boundary polygons from vector files. This is not footprint or extent creation (I already have scripts to do this). I want to create a tool similar to the QGIS, Image ...
1
vote
1answer
427 views

GDAL: C++ API to Convert Text file of Coordinate points to Raster Pixel line?

I am using GDAL to "Grow" areas form starting seed points. I have a text file listng the coordinates of my seed locations. The X,Y coordinates are in the coordinate system of a raster image. I want ...
1
vote
0answers
174 views

gdallocationinfo from C++?

I need a C++ API to convert x,y Image Coordinates in their native projection to output image pixel,line coordnates. I see the gdallocationinfo documentation, but no API for it. Where in online GDAL ...
3
votes
1answer
925 views

How do I split polygons that cross the dateline?

I have a dataset of country borders which I obtained from: http://thematicmapping.org/downloads/world_borders.php/ I have that dataset in both shp and SQL, but I'll put out there that I don't feel ...
1
vote
2answers
930 views

GDAL does ECW but QGIS (Linux !!) does not load ECW file

Followed instructions from various posts (each one with slight variations on the theme) and most of the time have no trouble getting the installed GDAL to accept ECW files. However, QGIS stubbornly ...
3
votes
1answer
388 views

How to stop writeOGR from abbreviating Field Names when using “ESRI Shapefile” driver

I'm currently using the following script to add some attribute data from a table to lots of individual shapefiles: library(rgdal) specieslist <- read.csv("SpeciesList1.txt", header=F) attdata ...
1
vote
3answers
551 views

Getting Field value for a given coordinate out of ESRI Shapefile using GDAL/OGR

I have a bunch of ESRI shapefiles that contain Land Cover data across Canada (provided by Geobase). I need to query the files and obtain the value of a certain field (called COVTYPE) in a given ...
6
votes
3answers
2k views

How to smooth / “reinterpolate” a raster with GDAL?

Short: I want to change the resolution of a raster and smooth the gray ramp like shown in the images bellow. The preference is to use GDAL, PIL or Numpy. Description: I'm kriging points into rasters ...
7
votes
2answers
496 views

How to catch PostgreSQL/PostGIS errors in GDAL/OGR Python bindings?

Is it possible to catch PostgreSQL output error messages in OGR after a query has failed? For example conn.ExecuteSQL('SELECT * FROM non_existing_table;') In Postgres this returns an error ...