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

5
votes
1answer
286 views

Generate hillshading with GDAL, and respect nodata

I'm still trying to nut out this problem of merging overlapping geotiffs. From the original dataset (here, a 20m resolution digital elevation model of the state of Victoria), I have reprojected it ...
-1
votes
1answer
219 views

How to Geocode use GeoPy API combined with GDAL Python bindings?

My question is related to the 'How to Geocode 300,000 addresses on the fly?' topic that has been asked by Mehul Nahar. I'm interested in running the processes explained by capooti in response to mehul ...
0
votes
1answer
305 views

ogr2ogr DXF to PostgreSQL some LINESTRING coordinates mirrored to the negative longitude

I have been using ogr2ogr GDAL tool to migrate geospatial data from DXF file to PotgresQL database without much problems. But now there is a special problem I can not fix. I have a DXF file with ...
0
votes
1answer
174 views

Why do Shapefiles and GeoJSON behave differently in GDAL Python?

I have a test.shp. I loop over the file, buffer the features, and write the buffers to a new file. I tested writing the buffers to a Shapefile and to a GeoJSON. If I try to use the new buffer layer ...
1
vote
2answers
204 views

Calculate logarithm of a raster?

I'm trying to calculate the logarithm of each cell in a raster. Is there a way to do this with F/OSS GIS? gdal_calc.py only has basic operators and the QGIS raster calculator feature doesn't seem to ...
1
vote
0answers
140 views

Trying to build GDAL's ECW plugin in 64 bit Linux

I'm trying to build the GDAL's ECW plugin in Linux. I'm using the latest version of the SDK from Intergraph's site (5.0). It fails because it tries to compile against the 32 bit version of the ...
1
vote
0answers
223 views

Merge overlapping DEMs with irregular boundaries

I have a 10m DEM covering part of a state, a 20m DEM covering the whole state, and a 90m DEM covering the whole country. How can I merge them together, getting the best quality available in each area? ...
1
vote
1answer
105 views

Load Raster Values into Shapefile [duplicate]

I have a shapefile containing many points. Some of these points fall within an area for which I have raster data. I would like to load the values of the raster which underlie each point into the ...
1
vote
0answers
201 views

Can GDAL Retile run multi-threaded?

gdal_retile is currently using only one CPU of the available eight and it is taking a while (48 hours already) to retile 5TB of data.
1
vote
1answer
132 views

How to get WGS84 lat lon bounds from a VRT file

I have vrt files generated using gdalbuildvrt that reference GeoTiffs referenced using the NZTM coordinate system. Is there an easy way to determine the bounds as WGS84 latitude longitude ...
1
vote
0answers
139 views

GeoDjango with MongoDB

Has anyone had experience successfully deploying GeoDjango with MongoDB? Is it possible? If so, how is it as simple as as installing the non-relational fork of Django and setting MongoDB as the ...
2
votes
1answer
257 views

gdal_translate produces tif with different colors as original

I am using SDE as my raster data source with the ultimate goal to tile this dataset: I am running the following gdal_translate command to generate the initial VRT: gdal_translate -projwin 390605 ...
1
vote
3answers
317 views

KML to Shapefile using GDAL (LIBKML) - MultiGeometry data not transforming

I am trying to convert a KML file to Shapefile using GDAL (ogr2ogr). Regular KML files convert fine but KML files containing MultiGeometry features (eg: Geometry containing both Point and Polygon) do ...
2
votes
0answers
82 views

C# GDAL DEM Speed

My 64-bit, Win 7 C# application using GDAL creates a number of constant-bearing paths with 1000+ points along each path. At each point, I want to look up the 4 closest elevations for interpolation. ...
3
votes
1answer
288 views

gdal_translate- why jpg compressed tif is 2 times greater than jpg file?

Using the following two commands yields completely different size files: gdal_translate test.tif -of GTiff -co "COMPRESS=JPEG" -co "JPEG_QUALITY=100" output.tif and gdal_translate test.tif -of JPEG ...
0
votes
2answers
437 views

GDAL and Python: how to create new files for (e.g.) clipping with gdalwarp

I'd like to learn using GDAL, and I want to use it in Python. Following the tutorial (http://www.gdal.org/gdal_tutorial.html), I succeded in opening vector and raster files using ogr.Open and ...
2
votes
2answers
170 views

display GeoReferencing image through MapServer

I have a tiff image and got the world file which contains the needed information for GeoReferenceing, how to start add this image on Mapserver ? if I use C# and use GMAP API to call the google server ...
1
vote
2answers
389 views

How to remove z & m coordinates from shapes with QGIS? (Or arcgis, grass, gdal, Python…)

I've recently received some parcel data in shapefile format that has z & m coordinates (all = 0.0). The WKT for a typical parcel looks like this: "MULTIPOLYGON (((2672464.560700002 ...
5
votes
2answers
448 views

In Python, reading a GDAL raster from memory instead of a file

I'd like to analyze some raster data in Python using NumPy. This data doesn't necessarily exist as a file on the filesystem, it could be the result of a query against PostGIS, or something ...
0
votes
1answer
141 views

GDAL command lines on windows [duplicate]

I am trying to get started with GDAL on a windows machine using Python2.7. I have tried commands such as: gdalinfo RasterFileName.TIF and gdal_translate.py -projwin ulx uly lrx lry -of GTiff ...
0
votes
1answer
99 views

How to consider more than one no data value in utility gdalwarp of GDAL

I have to do a bilinear interpolation on an image with two no data values, using the utility gdalwarp of GDAL (http://www.gdal.org/gdalwarp.html), so these two values don't must be considered in ...
0
votes
1answer
89 views

GeoReference use wordfile in arcgis [duplicate]

I have a tiff image and have the required information to Geo-reference it on the MAP. which called world file like this. using GMAP.NEt API there are many limited functions not like the java script ...
4
votes
1answer
349 views

How to apply “Band” settings using gdal Python bindings?

I am trying to replicate the gdal_translate command, but using pure Python: gdal_translate infile.tif outfile.tif -co COMPRESS=JPEG -co PHOTOMETRIC=YCBCR -co TFW=YES -b 1 -b 2 -b 3 I have ...
1
vote
1answer
182 views

gdal2tiles doesn't use resampling method correct?

I'm using gdal2tiles to create tiles from MODIS satellite imagery (250 m resolution) and I have noticed that I get the same bad result if I use -r near or -r bilinear (or lanczos for that matter). See ...
1
vote
1answer
182 views

Can I edit pixels values of an image with GDAL?

There is a method with GDAL (used in C/C++) to edit pixel values of a georeferenced image, maintaining the original file? Or I have to create a new dataset? My goal is to maintain ALL originals ...
3
votes
0answers
377 views

Fortifying re-projected shapefile for ggplot2 in R

I have a shapefile of US county boundaries which I got from the NHGIS. I can load the file using the rgdal package and see information about it: library(rgdal) shp_1860 <- readOGR("shp", ...
0
votes
1answer
48 views

OGR FeatureCount returns '-1' with OSM data

I am trying to read in an OSM file with OGR. import ogr driverOSM = ogr.GetDriverByName('OSM') roadsDs = driverOSM.Open('data/map.osm', 0) roadsLayer = roadsDs.GetLayer() The feature count returns ...
2
votes
1answer
265 views

How can I dynamically trim the informational border off of FAA sectional raster charts?

I'm trying to write a script that can automatically use GDAL to crop the informational border from the FAA's published sectional charts (GeoTIFF) found here: I know that I can use gdal_translate ...
3
votes
2answers
209 views

Homebrew OGR with FGDB support on OSX 10.9

This answer worked beautifully in 10.8, but isn't working immediately in 10.9. The link to download the FileGDB API still works, but there are a bunch of options available I don't remember from ...
3
votes
2answers
184 views

How to work with OGR and GDAL in GP Service in ArcGIS for Server 10.1

I have got a lot troubles with GDAL and OGR in the context of ArcGIS Server on 64 bit server machine. I can run script that use GDAL in ArcGIS as Tool and in the debug mode in PyScripter, then ...
1
vote
0answers
29 views

How to reproduce the gdal utility programs in Java by GDAL Java bindings [duplicate]

Can anyone tell me how can I reproduce the gdal utility programs gdal_translate, gdalwarp and gdal2tiles in java by GDAL Java bindings? Thanks!
1
vote
1answer
2k views

Any experience with OSX and Mac OSX Maverics using OSGeo? [closed]

There is a new Apple Release of OSX called Mavericks and my questions are: Will GDAL run? Will QGIS run? Will Tilemill run? Will Mapnik run? Also are there any known issues? A similar question ...
0
votes
0answers
305 views

What is meaning of GDAL on Ubuntu install [ERROR] Cannot make qgis.db private copy?

I am so lost and confused. I have spent days on this and not had any luck, so I'm really hoping someone can help! I am a new opensource user with an Ubuntu 13.04 OS. I used the software GUI (Ubuntu ...
2
votes
2answers
961 views

how to for loop a folder to batch clip rasters by polygon using python and QGIS

I am using python and QGIS 2.0. I am trying to clip rasters in a folder by one polygon feature. It's the first time for me using (let's say) "PyQGIS", I was used to arcpy before. Anyway, I don't get ...
1
vote
2answers
339 views

How to open OSM file in Python with GDAL

I am trying to open a *.osm file in Python with GDAL. But I doesn't work. What am I doing wrong? import ogr driver = ogr.GetDriverByName('OSM') datasource = driver.Open('map.osm') layer = ...
1
vote
2answers
789 views

Installing python GDAL with pip in local (user) space on RHEL6 server: failure due to missing headers

There are similar themed questions for other platforms, but the solutions out there do not work for me. (e.g. ...
1
vote
1answer
118 views

How to work with 16 GB orthophoto raster

I would like to use one GeoTIFF raster in a custom ArcGIS-Engine Application, but my provider gave me 4 orthophotos GeoTIFF 4 GB each, because GeoTIFF is limited to 4 GB. Is there any other way to ...
0
votes
3answers
240 views

How to find the attributes of a shapefile from the command line?

I'm following Mike Bostock's tutorial "Let's Make A Map." He explains how to filter a shapefile using ogr2ogr based on the attributes of the features. I can open a shapefile in QGIS and look at the ...
2
votes
1answer
108 views

Problems in gdalwarp

I was trying to use gdalwarp inside a qgis 2.0 model via processing but I think that the implementation is not working properly. So I would like to know if anyone has encountered this problem?
3
votes
1answer
292 views

How to change the field value of a Shapefile using GDAL/OGR?

I am trying to change the field value of a shapefile. However it doesn't change the value. What's wrong with my code? import ogr driver = ogr.GetDriverByName('ESRI Shapefile') fn = 'dist.shp' ...
1
vote
0answers
122 views

GDAL and Python: How to “pad the extent” or otherwise handle band.ReadRaster() for area falling slightly outside the image bounds?

In a nutshell, I'm trying to create a bare-metal tiling utility using GDAL and Python. As you can see in the graphic below, which shows the perimeter of a raster and the OSM-structured tile index, ...
2
votes
1answer
254 views

How to run OGR in ArcGIS Scripts?

When I run my script outside ArcGIS it works because during execution script uses a GDAL and OGR modules from osgeo. When I use script as a tool in ArcGIS it does not work. In PySripter I can import ...
1
vote
1answer
79 views

OGRFeature::DestroyFeature fails

I am using OGR to read features from a shape file; when the feature is destroyed I am getting an error 0xC0000005: Access violation reading location 0xffffffffffffffff. This secton of the code looks ...
6
votes
2answers
1k views

Best way to clip a large raster ECW file ?

I am trying to clip a large ECW (details below) but the raster file is too big to be processed entirely. Some details of the ECW below Driver: ECW/ERDAS Compressed Wavelets (SDK 5.0) file size: ...
5
votes
5answers
609 views

Impossible to sum rasters because of nodata pixels

I have a set of rasters (.vrt) with daily soil moisture data. I want to sum the pixels of all rasters in order to have a measure by month. However, the whole world is not covered each day, which ...
2
votes
1answer
1k views

QGIS 2.0 without MySQL connector

On a mac (Mountain Lion), after a major upgrade to GDAL (1.10) and related frameworks from KyngChaos page, QGIS 2.0 doesn't show any MySQL connector: Add Vector > Database > ... The driver ...
0
votes
0answers
148 views

why readOGR() is not reading projection information from .spj - r

0 down vote favorite I successfully import a polygon shapefile with readOGR, however the information in the .prj file is not being imported. test class : SpatialPolygonsDataFrame nfeatures : 19407 ...
1
vote
0answers
98 views

Efficient computation of path profile with Numpy and GDAL

I have numpy array with SRTM elevation data and I need calculate path profiles - vectors with elevation data from one point (transmitter) to every points of my map. Here's my current solution: def ...
1
vote
1answer
141 views

Writing array to vrt

I would like to write an array from a single band image to a 3 band vrt where each band has the ColorInterpretation() tag of red, green, blue respectively. My code below; from osgeo import gdal ...
1
vote
1answer
148 views

SetColorInterpretation() for a GeoTiff using GDAL not working

I have a single band image which I would like to copy to a 3 bands of a new image and set the ColorInterpretation to red, green and blue respectively. I have tried two approaches, one using the GTiff ...