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

Transformation of OGR style data

I have been using GDAL/OGR in a project and it is going quite well. I have been re-projecting data from various source formats into a single CRS (using the C# bindings), but one issue I have come ...
2
votes
2answers
431 views

Mask image based on no-data pixels in another image with different projection - GDAL or Python

I've got a land cover image which covers the whole globe, and a satellite image which covers a very small part of the globe. These two images have different projections, co-ordinate systems and pixel ...
1
vote
0answers
185 views

Broken PostGIS support GDAL

I'm having some problems with Postgres/PostGIS after rebuilding my GDAL library:( I wanted to add MrSID support in GDAL so rebuilt it with the MrSID SDK. Everything compiled ok and I am now able to ...
2
votes
1answer
85 views

how do i get a GDAL(using gdalwarp) clipped layer projected in pseudo mercator ESPG:3857 to not be shifted?

how do i get a GDAL rasterized layer projected in pseudo mercator ESPG:3857 to not be shifted ? i have a vector layer projected in ESPG:3857,after rasterizing using this command: gdal_rasterize ...
1
vote
0answers
189 views

Workflow with raster in Qgis 2.2

I'm trying to build a workflow with the Qgis 2.2 processing toolbox. I want to clip a raster layer with a vector layer using the GDAL tool "clip with a layer mask". But it doesn't work. The error ...
3
votes
1answer
1k views

Import shp to Postgis using Python and ogr

I just export a Postgis Table to shp using this tips but I'm not able to import a shp to Postgis using the same library(ogr). Any idea? Thanks a lot f.
1
vote
2answers
162 views

Remove unneeded columns from geographic data file?

I have some GML data that has some unneeded columns in it. I'd like to remove these columns in order to reduce the file size. (To clarify: I mean that the file has many points in it, and each of these ...
1
vote
0answers
31 views

ogr2ogr giving “ERROR 1: Failed connect to 192.168.4.18:8181; Operation timed out” error?

This is my command, which was working OK yesterday: ogr2ogr -f "GeoJSON" ../geojson/myfile.geojson ../data/gml/myfile.gml -s_srs EPSG:27700 -t_srs EPSG:4326 The output is this: ERROR 1: Failed ...
0
votes
0answers
94 views

Fantasy map, how do i slice large png and georefrence it the right way?

Im and IOS developer working on a map based ios game, and new to GIS stuff. I've read a lot of info and i started to get a hint about what i need to do in order to meet my goals. I have a png map i ...
1
vote
1answer
81 views

why is map tiling computationally intensive?

I recently starting using gdal2tiles.py (and its parallel implementation) to generate tiled images from raster datasets for kml interpretation and I'm shocked by how computationally intensive the ...
0
votes
1answer
1k views

AttributeError: 'module' object has no attribute 'Open'

Hy, I am running the following code which I found on: GDAL - Perform Simple Least Cost Path Analysis which I adapted to my case, but I get the following error when running the code in the python ...
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
256 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
163 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
111 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
175 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
43 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
934 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
147 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
637 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
578 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
325 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
382 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
41 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
79 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
182 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
95 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
437 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
439 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
558 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 ...
1
vote
2answers
523 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
112 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
444 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
217 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
230 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
401 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 ...
3
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
289 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
103 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
115 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
114 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
253 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
210 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
57 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
64 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 + ' ' + ...