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
65 views

Convert a KML ImageOverlay to World Files with GDAL

I have many raster outputs from an Open Source radio transmission topographic path tool called SPLAT, georeferenced in KML image overlay format. I have figured out an easy way to use tools to ...
1
vote
1answer
113 views

PostGIS Raster in MapServer

I'm having trouble getting a PostGIS raster to load and display in MapServer. I'm not sure, but suspect it's a version problem. My mapfile layer and the error message are shown below. LAYER ...
2
votes
2answers
659 views

How to convert a .tif with .tfw to a GeoTiff?

I have a tif image file, and a .tfw world file (created from OSM data with Mapertive). How can I (using linux command line) convert that to a GeoTIFF file? (i.e. I want to have one file, not 2). I'm ...
0
votes
1answer
127 views

Create Tiles for a raster(georeferenced) to use in leaflet

I have raster, with each pixel containing some data. What i want is to split a large image with each pixel's color is determined by its value, to show on leaflet layer.How to go about this ? raster ...
1
vote
0answers
59 views

Extra row and column in a raster after GDAL grid import

I'm exporting an ESRI grid to TIFF format and after that I'm reading this exported raster file into SAGA GIS. However when Saga reads this file in there appears to be an additional row and column in ...
0
votes
0answers
77 views

GEOS assertion failed

The following script fails, if I import the 'rasterstats' library. import rasterstats import ogr geom1 = ogr.CreateGeometryFromWkt('POLYGON ((415924.806155724043492 ...
2
votes
1answer
67 views

How to use GDAL to import lat/lon points from ESRI File Geodatabas into MySQL

I have a file geodatabase from here: ftp://ftp.csc.noaa.gov/temp/MarineCadastre/AIS.SampleData.zip Using ogrinfo, I can see that it contains a layer called "Broadcast_" that has lat/lon points: ...
3
votes
3answers
132 views

how to digitize areas of a raster and snap to the raster cells?

I need to digitize the extent of an ecw raster image where the RGB value is not black.(see attatchment)ie. not the envelope, just the area where the image is shown without the border areas. The black ...
1
vote
1answer
78 views

WKT for local mine grid

I have been trying to create a WKT entry to use in gdalwarp etc. I need an affine transform from coordinates in MGA94 Zone 54 to our local mine grid. All efforts have so far failed, including using ...
0
votes
0answers
28 views

How do i install liblas using OSGeoW4

I have downloaded and installed the OSGeo4W program on my system. I recently downloaded gdal file and liblas. When i try opening the osgeow4 shell, it keeps telling me the gdal file is missing. when i ...
0
votes
1answer
40 views

Missing Georeferencer

I installed Qgis Master 64 bits on win7 64 16gb of ram and I am unable to locate the georeferencer. All other Gdal application are there, Gdal is active, but no georeferencer. I tried yesterday, ...
1
vote
1answer
419 views

Clip raster layer with *raster* mask layer in QGIS

Given two input layers 1. The clip layer - an arbitrary layer (though in my case it's a float 1 band geotiff) 2. The mask layer - a paletted raster in which every pixel is either NO DATA or a single ...
1
vote
2answers
32 views

Comparing projections of two datasets when SRID not available

Similar to this question: How to test if two GDAL datasets are in the same projection? However, the above does not cover cases where an SRID is not available. It can often be the case that an SRID is ...
1
vote
2answers
458 views

os.system /subprocess.call error handling when error in command line function?

In my Python script I call a GDAL utility using: subprocess.call('gdal_rasterize -a ICE_TYPE -where \"ICE_TYPE=\'Open Water\'\" -b 1 -burn 5 -l ' + shapefileshortname +' ' + shapefile + ' ' + ...
0
votes
1answer
203 views

GDAL WMS Question

I'm new to GDAL and am currently trying to hit a WMS using gdal. I need a way to save the image to my desktop. I currently have an xml file with the WMS description syntax that gdal wants -> ...
3
votes
2answers
152 views

How to check which version of GDAL do I have and how many?

I have installed in my computer and running the following software: PosgresSQL 9.3.2.3 PostGIS 2.1.1.1 QGIS 2.0.1 Dufour (this one installed through the OSGeo4w package) QGIS 1.7.0 Wroclaw I am ...
0
votes
0answers
40 views

GDAL - least cost path analysis height data

I would like to use Pythons GDAL bindings to find the least cost path through an array and then show the route on top of the original array. I am using DEM data and would like the slope to be taken ...
1
vote
2answers
235 views

Processing HDF4 data in R

I am trying to process MODIS files in R, which are in HDF4 format (link to data). I have experience reading in HDF5 and netCDF files in R, but HDF4 seems to be a common problem in R. I am using a MAC ...
3
votes
2answers
646 views

How to convert a NetCDF file to GeoTIFF using GDAL

I have a NetCDF file which contains data for the number of rain days in the UK. I'd like to convert it to a GeoTIFF, but can't get GDAL correctly georeference it. See dumps from gdalinfo below. I ...
4
votes
2answers
234 views

how would one get the end points of a polyline?

I am attempting to write a script to get the end points of a series polylines, and then use those end points to create a polygon. I'm not sure which tools to begin to use on this project. I would ...
2
votes
3answers
110 views

Geotiff reprojection

I am new for GIS programming. I try to change geotiff coordinate system from epsg 32635 to epsg 4326. I used the following gdal command: gdalwarp -rc -s_srs epsg:32635 -t_srs epsg:4326 -dstnodata 0 ...
3
votes
2answers
126 views

Python, GDAL, VRT: help wth code efficiency and processing time

I'm going to attempt to make this question as succinct as possible while still including relevant information, as I have trouble with being concise and am trying to avoid writing unnecessary stuff. ...
1
vote
1answer
82 views

How can I use the proj string in java code?

I want to do projection from UTM/oblique mercator CRS to WGS84. And I know the solution to use a proj.4 string in commandline like+proj=omerc +lat_0=51.4 +lonc=7 +alpha=-10 +k=1 +x_0=0 +y_0=0 +gamma=0 ...
0
votes
0answers
30 views

Installing GarminCustomMap: “missing _gdal_array module” error

I've used the GarminCustomMap plugin in QGIS 1.8 without troubles. But now, with QGIS 2.0, when I try to install it, I get the following message "missing _gdal_array module". Any idea how to fix this? ...
1
vote
1answer
144 views

how to create CF-compliant NetCDF from ESRI Grid?

been trying to make a CF-compliant netCDF file from an ESRI Grid. I want to serve it w/the THREDDS server. Been using gdal_translate. I know the netcdf output option is definitely not CF compliant, so ...
1
vote
1answer
94 views

Trying to build GDAL with ECW SDK 5.0 on Windows

I'm using Visual Studio 2010 to open and build the GDAL 1.10 source (using the "makegdal10.sln"). This builds without an issue. Now I want to include the ECW support by including the ECW v5.0 SDK ...
0
votes
0answers
60 views

Update a NetCDF file using GDAL

I have to edit a NetCDF file using GDAL in C++ code. When I try to update a subset I have this error message: ERROR 6: The NETCDF driver does not support update access to existing datasets. How can ...
3
votes
2answers
167 views

Merge two rasters with no data value

I have two rasters: Raster A and B. Raster A has has only values of 0. The white area in Raster B are No Data values. I want to merge them to one raster (Raster C). So basically I want all pixels in ...
1
vote
1answer
161 views

Why does OGR driver.Open() for File Geodatabase give AttributeError: NoneType object has no attribute Open?

There are a number of existing Questions that hint at a possible Answer to this Question, but none seem to be providing me with the solution to my problem: How to access feature classes in file ...
2
votes
1answer
186 views

“ImportError: No module named gdalsrsinfo” What's the issue?

in the python console I can call import gdal, and gdal is recognized as a module; however, most gdal programs, such as gdalinfo and gdalsrsinfo are not recognized. I checked the GDAL source files to ...
2
votes
0answers
88 views

Zonal Statistics using Gdal and Python

I recently came across a loop_zonal_stats function by @ustroetz found here: Issue Trying to create Zonal Statistics using Gdal and Python.. I am using a (single) polygon shapefile and a raster file to ...
1
vote
1answer
65 views

gdal_grid has some problem with reading vrt file

I am using gdal_grid module(from gdal 1.10) and it has some problem with reading the vrt file. I checked this link but it does not work. gdal_grid problem reading vrt file The command is ...
0
votes
0answers
130 views

Add labels for printed map (GeoTools)

I have some problems with labels for printing map with GeoTools. Next code (with some cuts, code on Scala language) can help me for printing map: val map = new MapContent() val parser = new ...
2
votes
2answers
90 views

How to show only raster attribute table using GDAL?

I have a binary ArcGRID raster (AIG/Arc/Info Binary Grid) with an attribute table. I would like to get only this attribute table but gdalinfo outputs also projection information and bounds and then ...
1
vote
0answers
62 views

Write raster with its attribute table and symbology

I am classifying lots of rasters and I'd like to automate the legend creation and symbology displaying for easy manipulation. Here is the kind of raster I'm dealing with: library(raster) r <- ...
0
votes
0answers
21 views

Is there a way of converting polygon vertices to points using GDAL python? [duplicate]

I have an ArcGIS standard license and I want to convert polygon vertexes to points in a script. Since I don't have a good enough ArcGIS license I've been using ET Geowizards to get the points. The ...
0
votes
2answers
69 views

Is it possible to alter raster values without Arcpy? [closed]

I'd like to be able to alter raster values, through fairly simple equations, without using the module Arcpy. Is that possible? It seems like gdal, numpy, and rasterlang might make this possible but ...
1
vote
2answers
186 views

how to georeference my png using gdal

Hello I have created my own fantasy map, I have the map in a full size png and the map broken down into tiles ( with no geo data ) I installed MS4W and tried to use it to get the georeference data ...
1
vote
1answer
65 views

gdal_merge Sum/Agregation files into single output band

I need to create a single tiff files from several others (each have the same extent and resolution) and the output file need to have a single band containing the aggregated/sum of the values in the ...
0
votes
0answers
45 views

How to assign no value to pixels with (255,255,255)

I have Imagery where I want to replace white edges (255,255,255) with no value. I have used gdal translate within QGIS and put in no value 255 but this give no value to other pixels eg 255,240,200. ...
0
votes
1answer
24 views

Break Mosaic into Large Tiles

I have 607 aerial images of Massachusetts in .jp2 format. They are in a mosaic however I need to break down the 607 images into groupings of 16(4x4) spatially so that I can more easily process ...
4
votes
1answer
191 views

Why is tiff different in epsg 3857 and epsg 900913 in gdal

When I transform tiff from 4326 to 3857 and from 4326 to 900913 in gdal, I get different result(when display in Qgis) eventhough these two projections(3857 and 900913) should be the same. What is ...
3
votes
3answers
321 views

What format and settings to use for aerials photos in QGIS?

Following question which was more about handling Aerials with ArcGIS: Most Effective Format to Manage Aerial Photography for Viewing Only It seems there are 2 major options for storing / resampling ...
0
votes
0answers
76 views

Choose what zoomlevel to read from for mbtiles file in GDAL

GDAL is supposed to read MBTiles files and to do so, it selects the highest zoomlevel available. However, is there a way to choose another zoomlevel? i.e. override the default and go for a lower one? ...
1
vote
1answer
53 views

Is there any way to add an area column to new shapefile created from raster with gdal's polygonize in the same step?

I've been going through a fairly lengthy process of polygonizing very intricate raster images, then creating a new layer with the extra area attribute, so that I can then remove the millions of tiny ...
3
votes
0answers
376 views

Reprojecting WGS 1984 Web Mercator (EPSG:3857) in Python with GDAL

I am reprojecting rasters in python using GDAL. I need to project several tiffs from geographic WGS 84 coordinates to WGS 1984 Web Mercator (Auxiliary Sphere), in order to use them later in Openlayers ...
2
votes
2answers
58 views

Using gdal_merge to create a grayscale mosaic

My gdal fu is weak and I need help! ;) I'm creating a mosaic in geotiff from hundreds of ecw. I created a list of all my rasters and used this command: gdal_merge -o c:\temp\output_image.tif -q -v ...
0
votes
1answer
75 views

Need help installing QGIS 2.0 Need GDAL then SQLite3 3.7.16

OSX 10.7.5 I need help installing QGIS 2.0 It says I then need GDAL. When I try to install GDAL it says SQLite3 3.7.16 or newer framework is required. I can't find SQLite3 3.7.16.
1
vote
1answer
97 views

How do I find out the resolution in (centimeters per pixel) of a TIF?

I tried using GDAL and got: # gdalinfo nicholas_map_2.tif Driver: GTiff/GeoTIFF Files: nicholas_map_2.tif nicholas_map_2.tfw Size is 28799, 27875 Coordinate System is: GEOGCS["WGS 84", ...
1
vote
1answer
254 views

Convert polygons in shape file to a numpy array?

I have a shape file containing a number of polygons and all attributes being text. I'm trying to import the polygons into numpy as an array where each polygon is represented as unique values. I ...