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

3
votes
1answer
1k views

Georeferencing an image with GDAL, based on QGIS Georeferencer Plugin

I would like to georeference an image with GDAL. The image already has a spatial reference set (EPSG:31276). For the georeferencing, I have a list of 4 ground control points (also in EPSG:31276). ...
6
votes
1answer
262 views

GDAL equivalent of ArcGIS Mosaic

The word 'mosaic' is often used for quite different things within the worlds of ArcGIS and GDAL. With GDAL and it's downstream applications like Qgis it is common to see mosaic used in the context ...
2
votes
1answer
175 views

File size inflation normal with gdalwarp?

After using gdalwarp to project and align-to-grid (via -tap) a number of rasters I noticed that the output rasters were significantly larger than the original rasters. A fairly thorough web search ...
1
vote
0answers
115 views

How can I draw features on GeoTIFF in Python?

I have a set of data (coords + value) and want to turn it into a heatmap in GeoTIFF. Should I calculate pixels and set them as array, or there are shortcuts? This is what I do now: driver = ...
1
vote
1answer
191 views

Set a RANGE of values as nodata in a .vrt?

Is it possible to set a range of values to nodata in a vrt file? How? I have some elevation models with incorrect values, down in the -15,000s. I could just clip the rasters along that edge, or apply ...
1
vote
0answers
32 views

Problem defining the tiling starting point

We are tiling a image datasets. However, we want to specifically define custom windows. To do this, we are using -te option during warping. The warped output is then tiled. The problem though is that ...
0
votes
1answer
48 views

Missing Layer in OpenLayers

When I try to include the "fwp layer" from geoserver 2.4 (Ubuntu 12.4 VServer) into my code it doesn't show up. But the geoserver layer preview (OL) shows the right image. The GeoExplorer preview ...
3
votes
3answers
977 views

Per-pixel (statistical) calculations on a raster stack using GDAL

In R, it's relatively trivial to perform per pixel calculations based on a raster stack (e.g., get std.dev for each pixel on a 12 layer GeoTIFF). Unfortunately, the speed is less than desirable when ...
0
votes
1answer
152 views

gdal's polygonize function (via qgis) is not turning monochrome rasters into readable shapefiles regardless of inversion of color scale

I've tried polygonizing this monochrome bitmap as this: and this: and either way it ends up looking like this when (zoomed in) after gdal's polygonize operation: Is there any way I can tell ...
2
votes
2answers
652 views

GDAL Conversion of US Topo GeoPDFs

I've used GDAL_translate with "--config GDAL_PDF_LAYERS" option to convert US Topo GeoPDFs to GeoTIFF and it works great for creating collarless topos based on the new format. However, I can't get it ...
5
votes
1answer
602 views

ArcGIS Raster Mosaic Dataset and GDAL VRT, when to use?

ArcGIS Desktop has a the Raster Mosaic Dataset, QGIS and friends have Virtual Raster Tables. A skim of the documentation reveals much similarity. They both create or enable an alternate ...
0
votes
0answers
338 views

ECW files to geoserver 2.4

I'm trying to add .ecw files (20 orthophotos) to geoserver 2.4. Geoserver is running on an ubuntu 12.4. Vserver. I'v already installed the GDAL plugin. Is this the right format for the ...
0
votes
1answer
411 views

How can I correct a .pgw file for a PNG image?

I have a PNG file and I need its .pgw file to be set with NAD 83 as the projection. I have created a .pgw file: 799 0 0 -799 -122.562 37.1946 However when I open the PNG in QGIS or ArcMap 10.2 it ...
2
votes
0answers
42 views

GeospatialPDF javascript access

Using javascript in GeospatialPDF, is it possible to zoom in to an extent, say when the document is opened the first time ? To what extent can I use javascript to access the geospatial data? pan, ...
0
votes
1answer
83 views

Is it possible to get all projections form OGR?

I'm working on a project where the user can give my program so GIS data which may or may not have which projection it is using. So is it possible to get a list of all the projections from OGR or ...
2
votes
1answer
195 views

gdal2tiles generates small picture inside 256x256 tile

I've used gdal2tiles for generating tiles in number of different polar projections and always experienced the same problem — strange small picture inside normal 256x256 png with the rest of space ...
0
votes
1answer
99 views

Cannot display Oracle raster images using QGis

I am dealing with QGIS 2.2 and the plug-in Oracle spatial georaster. Apparently the connection to the database can be made, as the raster images are listed once you press 'Connect', but when try to ...
1
vote
1answer
494 views

ogrinfo sql example

I'm having trouble learning how to query my shapefiles for certain features. I have a shapefile that looks something like this when I summarize it using the following command: ogrinfo -al -geom=OFF ...
1
vote
1answer
453 views

datum conversion with GDAL

I run gdaltransform like so: gdaltransform -s_srs nad27 -t_srs wgs84 -104.5 39.5 -104.5 39.5 0 gdaltransform -s_srs wgs84 -t_srs nad83 -104.5 39.5 -104.5 39.5 0 gdaltransform -s_srs nad27 ...
5
votes
2answers
1k views

How to access Shapefile metadata using OGR?

I'm trying to use ogrinfo to get some details on a shapefile I downloaded. Currently, the only way I know how to do this is to load it into QGIS and manually click around to find any information on ...
0
votes
0answers
15 views

UnIntersected part of gdal-geometry?

i have a rectangle and polygon geometries . polygon.Instersection(rectangle) gives me the polygon as it is inside the rectangle. But , what i want is the unintersected(non-intersecting) part of both ...
3
votes
2answers
616 views

How to batch “gdal_merge” on Win7 in one command line command?

I want to create a mosaic of all tif files in a directory on windows 7. How could I create a one-line command in the windows command line to loop through all tifs and merge them all together using ...
2
votes
2answers
553 views

I need to extract a specific layer from a Geospatial PDF using GDAL/OGR

I have GDAL/OGR 1.10 installed on a 64bit Windows 7 machine. I have some Geospatial PDFs that I want to crop to the outline of the map and convert to Rasters. However I cannot clip to the neatline. ...
2
votes
1answer
117 views

GDAL not accepting *.tif to list through files in directory

I am using gdal_merge.py to mosaic tif. If I explicitly write the file names it works, but if I use *.tif it fails with the following error. Traceback (most recent call last): File "C:\Program ...
4
votes
1answer
528 views

install gdal extention in GeoServer

I am trying to install the GDAL extension (I want to use ecw files in GeoServer). So here is what I did so far: 1) downloaded the GDAL extension and extracted to: C:\Program Files ...
0
votes
1answer
228 views

How to get the pixels from a geotiff file in gdal python for a given bbox?

I have a GeoTIFF file in WSG 84 Projection reference (it's a GeoTIFF of SRTM). I can use python & gdal to read the whole lot of data into an array (with ReadAsArray()), but I want to only have a ...
2
votes
1answer
241 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 ...
2
votes
1answer
420 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 ...
4
votes
2answers
3k 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
298 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
112 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
119 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
119 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
262 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
218 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
1answer
58 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
2k 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 ...
2
votes
2answers
65 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
2k 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
338 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 -> ...
8
votes
2answers
1k 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 ...
1
vote
2answers
1k 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 ...
4
votes
2answers
2k 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 ...
5
votes
2answers
668 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
206 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
319 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
173 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 ...
2
votes
1answer
301 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 ...
2
votes
1answer
214 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
85 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 ...