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
196 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
138 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
256 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
314 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
168 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
442 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
140 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
348 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
181 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
375 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
47 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
258 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
207 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
183 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 ...
0
votes
0answers
49 views

Getting through the jungle of gdal2tiles.py versions and patches

After trying out different methods of greating tiles from satellite images, I have decided to use gdal2tiles.py which seems fastest of the ones I looked at. Then I realise that there are several ...
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
303 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
956 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
332 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
776 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. ...
0
votes
0answers
90 views

GDAL Create( ) for a netCDF file segfaults

I have a problem with the creation of NetCDF dataset using GDAL. When I use "gdalinfo" and "gdal_translate" utilities on a NetCDF format, they work fine, but when I create a new dataset using the ...
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
233 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
289 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
120 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
251 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
78 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
602 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
147 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
140 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
147 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 ...
11
votes
1answer
569 views

How to create a KML that supports animated ground overlays?

I have 12 images that represent weather conditions at day 1 through day 12. I am using gdal2tiles to produce map tiles from each of these images. gdal2tiles also produces a KML file that can be used ...
4
votes
1answer
67 views

Raster Calculator Bug at start

I use ArcGIS 10.1 on a Windows XP. ArcGIS completely shuts down when I attempt to launch the "Raster Calculator" tool. Has this happened to someone else? Any ideas on how to solve this annoying issue? ...
0
votes
2answers
599 views

JP2 fails to open in QGIS and GDAL

I am trying to open a JP2 file in QGIS and it fails, and I'd appreciate any advice to try and track down the cause of the issue, be it with my system or with the file itself. The file is a HiRISE ...
1
vote
0answers
160 views

GDAL in Java fails on ReadRaster()

To summarize, trying to read a raster keeps giving me an UnsatisfiedLinkError that makes sense, but looking into the compiled GDAL jar, I don't know how to fix. gdal.AllRegister(); TheFile = ...
2
votes
1answer
121 views

gdal2tiles.py generating black tiles from .img file

I have a source HFA Erdas Imagine (.img) file which I'm trying to convert into tiles using gdal2tiles.py. I'm running: gdal2tiles.py -w google -g GOOGLEKEY input_file.img The HTML is generated ...
2
votes
0answers
123 views

Why is GDALPolygonize so much slower than ArcGIS Raster to Polygon?

I am attempting to polygonize a raster using GDALPolygonize() in a Python script. The script began polygonizing yesterday at 5pm and is still polygonizing now at 9:30 am. I have no clue how far ...