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
96 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
403 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
211 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
206 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
438 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
52 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
293 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
221 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
202 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
328 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 ...
3
votes
2answers
1k 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
393 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
1k 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
128 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
305 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
114 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
354 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
137 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
295 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
86 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
763 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
164 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
104 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
163 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
166 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
666 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
69 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
716 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
176 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
130 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
150 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 ...
2
votes
0answers
31 views

Workflow to assign and reproject VRT [duplicate]

I have the following workflow. generate vrt from imagelist file(with paths to tifs/grids) check if coordinate transformation is required, if so convert to target SRS Create tile using gdal_retile ...
0
votes
1answer
502 views

GDAL RasterIO grabbing memory

I have written an application in c++ to read LiDAR files and convert to various rasters using GDAL; this worked very well on small scale tests (47 LAS files) however, when used on a live task the ...
1
vote
2answers
776 views

Python - GDAL : Polygon to point shape file

I want to read all the points in a polygon shape-file and use their coordinates to make a new point shape file. then find the min and max of the X coordinates. So here is the code I wrote for the ...
3
votes
0answers
80 views

Tabulate Area using GDAL and Numpy

How could I perform a Tabulate Area using GDAL and Numpy? I am trying to summarize how many pixels of each class in a Landcover raster are within each polygon of a Census Blocks shapefile. I've ...
1
vote
1answer
730 views

Delineate polygonal catchment areas from a raster using open-source GIS

I'm 90% sure this basic question has been asked before but looking through the other posts I simply do not see a useful answer. I need to calculate polygonal(vector) catchment areas from a raster ...
1
vote
1answer
96 views

Need assistance setting center coordinate with GDALWARP

I'm having some issues converting my project into an MBTiles format. I'm trying to use gdalwarp to convert my main file to geotiff format, but I can't seem to get the center coordinate to how it is ...
3
votes
0answers
187 views

Can GDAL Perform Vertical Transformations?

Can GDAL perform vertical transformations? If the answer is yes, please indicate which amongst these listed pages can explain how exactly to accomplish vertical transformations? If there is some ...
0
votes
2answers
228 views

GDAL and Python - Don't print GDAL Error messages

According to the GDAL Python Documentation, GDAL functions do not raise exceptions or errors but do print errors to the console via stdout and stderr. I have some code in my script that calls ...
3
votes
2answers
253 views

How practical is this workflow, tiling from a complex, multi-step GDAL VRT?

This is probably a cw-candidate, but here goes.. I've got 1002 ~100 MB (5000'x5000') .tif files in a state plane (transverse merc) projection comprising a county-wide, 1-ft aerial dataset. After two ...
4
votes
0answers
245 views

How to convert a raster (array) to a NetworkX graph and do a * search?

I want to use a slope-raster for a a_star path analysis in NetworkX. I am using Python for this project. My test raster looks like this First I read in the raster with GDAL input_raster = ...
1
vote
1answer
62 views

Can I create a virtual mosaic using other virtual mosaics?

I am following this example - http://linfiniti.com/2009/09/image-mosaicking-with-gdal/ I have 6TB of aerial imagery in 6 folders 1TB each and each folder has many subdirectories including the ...
1
vote
1answer
167 views

Python/GDAL: Best way to read in raster for use in A* search algortihm

I want to create a least cost path program with Python and GDAL since I couldn't find any open source cost path additions for use in Python. As recommended here the best approach is to use A* ...
5
votes
2answers
3k views

How to add *.ecw files to QGIS 2.0 on Windows?

How do I add and view *.ecw raster images in Qgis 2.0 ?
1
vote
1answer
122 views

How to duplicate this QGIS workflow with gdal/GRASS? (reproject and crop a raster)?

I have a number of rasters (rgb images) that I want to re-project and crop. I have established a workflow for doing this in QGIS, but I'd like to do it using gdal modules, because the output is being ...
3
votes
2answers
555 views

How to set output resolution for gdalwarp to that of source file?

I have a large raster, which I want to re-project and crop for a new location that I'm creating in GRASS. I seem to have gotten the re-projection and cropping right with this call to gdalwarp: ...