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
1answer
83 views

Read shapefile as a numpy array

Is it possible to read a shapefile with ogr as a numpy array, in order to use it as a mask to filter out a geotiff image? What I would like to obtain is a numpy array from the layer below: from ...
1
vote
0answers
204 views

Reading, modifying and writing a geotiff with GDAL in python

I'm trying the learn the ropes of Remote Sensing image processing using Python GDAL bindings and numpy. As a first attempt, I'm reading a Landsat8 geotiff file, do a simple manipulation and write the ...
1
vote
0answers
38 views

how to reproject an image programatically in c# using gdal

I have to reproject an image programmatically in c# using gdal library. For the moment i'm a little bit lost in gdal and i don't really now which methods i have to use. I think i will have to ...
1
vote
0answers
22 views

gdal: Pixelvalues to nodata

I have a GRD-File where all pixelvalues which are less than 0.05 should be converted to nodata. I know how to do this with ArcMap, but the process should be done automatically, for example with ...
1
vote
1answer
51 views

Which coordinate system does MapTiler use for its “Photoshop Overlay for Google Maps” page?

I am attempting to use this MapTiler webpage to create a background image that I can position my image over and then use GDAL to generate tiles. My process is: Position the map on the MapTiler page ...
1
vote
0answers
90 views

What if 'geometry_columns' already exists in MSSQL Server Database?

When trying to load a .SHP into MSSQLSpatial using OGR2OGR, I get the following error: ERROR 1: Error initializing the metadata tables : [Microsoft][ODBC SQL Server Dr iver][SQL Server]There is ...
5
votes
1answer
190 views

Python GDAL documentation

Is there an official documentation for the python GDAL module? There are various cookbooks, but scant information on each function or arguments for the most basic functions. halp!?
0
votes
1answer
73 views

Use lamber 93 imageOverlay with leaflet on top of OSM tiles layer

I'm trying to use a lambert 93 image overlay on top of an Open Street Map title layer. var map = L.map('map'); L.tileLayer('http://{s}.tile.osm.org/{z}/{x}/{y}.png', { }).addTo(map); var coords = { ...
4
votes
1answer
128 views

How to correctly import a raster file with these metadatas into Tilemill?

I have a GeoTIFF file for which gdalsrsinfo reports this: PROJ.4 : '+proj=utm +zone=32 +datum=WGS84 +units=m +no_defs ' OGC WKT : PROJCS["WGS 84 / UTM zone 32N", GEOGCS["WGS 84", ...
0
votes
1answer
50 views

how do i create OGRSpatialReference from raster file's georeference C++

I want to create .shp file using ogr with georeference based on my raster file, so i create an ogr layer using CreateLayer( "point_out", NULL, wkbPoint, NULL ); i want to set the OGRSpatialReference ...
0
votes
2answers
203 views

How can I read MrSID files in Python (using gdal)?

I have to read a large number of MrSID images and process them in Python. I usually access GIS file formats using the gdal application and the gdal Python library. However, the gdal package for Python ...
4
votes
1answer
67 views

gdal_grid out of memory allocating 64 bytes

Python 2.7.4 32-bit GDAL 1.10.1 I have a Python script which takes a point file with x,y,z coordinates, makes a vrt of it and then tries to create a Geotiff from it with gdal_grid. During testing ...
0
votes
0answers
45 views

Define band specific metadata of a GeoTIFF file using GDAL (gdalwarp/gdal_translate)

Does anybody know how to change or set the "Description" option/tag (band specific metadata) of a GeoTIFF file using GDAL? To specify what I mean, this is an example of gdalinfo return from a GeoTIFF ...
2
votes
3answers
140 views

How to create a Projection from a CRS string using pyproj?

I'm trying to convert between two coordinate systems using pyproj: import pyproj p1 = pyproj.Proj("+init=EPSG:4326") p2 = pyproj.Proj("""PROJCS["NAD27 / California zone ...
0
votes
0answers
56 views

GDAL and WMTS supprt

I have researching possibility to fetch raster file from WMTS server via GDAL. I have found some old post about that GDAL doesn't support WMTS (WMS only). However, I also find this information: ...
2
votes
3answers
144 views

What is the altitude difference between WGS84 and NAD27 Zone II?

I have a point cloud in a piece of software and have exported the point in both WGS84 and NAD27 State Plane California Zone II (in US survey feet). The WGS data looks like this, with elevation in ...
1
vote
1answer
40 views

Add Color To GTiff output raster band

I am using GDAL in Python to create a tiff raster of data that is originally from a shapefile. I have an auxillary data channel, which I would like to use to color the pixels of the tif. I have made ...
1
vote
1answer
38 views

Connecting points located in opposite side of international date meridian

Given a list of coordinates, I need to create a line shapefile joining each of the point with a Python GDAL/OGR script. The procedure described here Convert XY points to a line? works fine in almost ...
1
vote
1answer
32 views

How to re-project only one band from a hdf?

I have this file which contains: gdalinfo product.hdf Coordinate System is `' Metadata: calendar=standard center=gsfc Subdatasets: Corner Coordinates: Upper Left ( 0.0, 0.0) Lower Left ( ...
2
votes
1answer
127 views

How to get pixel values of indvidual pixels falling under a poygon to shapefile attribute table or CSV?

I want to be able to write pixel values of each pixel falling under a polygon onto a polygon attribute table or a CSV. I have about 1200 polygons and a georefernced GeoTiff file with NDVI values. ...
0
votes
1answer
79 views

How to install OSM .pbf support in GDAL on Debian?

I'm trying to use Anita Graser's post on loading OSM data from .PBF files into a SQLite .db file using ogr2ogr but it stopped working. I'm trying to run the following command: ogr2ogr -f "SQLite" ...
1
vote
1answer
64 views

Adding fileGDB driver to gdal 1.11.1

I need to use the fileGDB driver in gdal, however I can't figure out how to install the driver. This is the only kind of documentation I found but it doesn't make any sense to me: ...
2
votes
1answer
133 views

Reclassify a raster value to -9999 and set it to the nodata value using python and or gdal

I modified a script which reclassifies a raster. It is based on the script by Roger in his blog (the second script) and it works pretty well.The problem I am facing is as follows: I would like to ...
0
votes
0answers
98 views

stack raster using GDAL and numpy dstack

I am testing Savitzky-Golay filter with my MODIS dataset. Currently I am testing for 10 rasters (but later on will use my whole dataset) only just to test how should I create my script. So my initial ...
1
vote
1answer
82 views

How to create .jgw file out of jpg picture with exif header?

If i have picture with this metadata in exif header: Size is 4000, 3000 Coordinate System is `' Metadata: EXIF_ApertureValue=(3.625) EXIF_ColorSpace=1 EXIF_ComponentsConfiguration=0x1 0x2 0x3 ...
1
vote
2answers
94 views

Convert jpg with exif header into geotiff

I have bunch of jpgs with exif headers like this: Size is 4000, 3000 Coordinate System is `' Metadata: EXIF_ApertureValue=(3.625) EXIF_ColorSpace=1 EXIF_ComponentsConfiguration=0x1 0x2 0x3 ...
0
votes
1answer
50 views

Reverse gdal_rasterize (-i option) does not work

When I use: gdal_rasterize -b 1 -b 2 -b 3 -burn 255 -burn 255 -burn 255 polygon.shp tiff.tiff I get a new tiff with a white hole in it with the shape of the polygon shp. However, when I use the ...
1
vote
1answer
57 views

Clipping a coloured raster with polygon shp in gdalwarp changes original values, why?

I need to clip my raster based on a polygonal shapefile, basically leave raster values inside the polygon shape (inside red perimeter in the figure below) as shown here: Therefore, I use: gdalwarp ...
3
votes
1answer
63 views

How can the Planet Labs TMS be used with GDAL?

How can Planet Labs' mosaic Tile Map Service (TMS) be used with GDAL like in the Open Layers example? The API requires authentication and other specific settings.
1
vote
0answers
34 views

gdal2tiles issues with color scale between tiles

Background: I have a series of about 40 small GeoTIFFs that need to be merged together to make one large image. Each GeoTIFF has 2000 bands. I've tried using both gdal_merge and gdal_retile to pull ...
1
vote
2answers
297 views

Is gdal 2.1 is available for windows?

I would try Gdal 2.1, does someone knows if it is available for windows ? For example I would try those WMTS features : http://www.gdal.org/frmt_wmts.html
0
votes
0answers
65 views

ogr2ogr conversion error - Can't create fields of type Binary on shapefile layers

I'm having trouble using the ogr2ogr tool on the command line. My data is stored in SQL as binary, type is mixed polygons and multipolygons. I am using the following command for conversion from SQL to ...
4
votes
3answers
177 views

Python GDAL/OGR: How to Erase features?

I have a big round polygon feature (5km diameter) and a small round polygon feature (3km diameter) that I would like to clip with Python GDAL/OGR to obtain the "donut". Once this is done, I will apply ...
2
votes
2answers
200 views

Raster reclassify using python, gdal and numpy

I would like to reclassify a raster file from a raster with 10 classes to a raster with 8 classes using pyhton, gdal and/or numpy. The classes are represented as integers. I have tried following the ...
0
votes
0answers
78 views

Tile generation (gdal2tiles): 'merge' old and new (overlapping nodata) tiles?

I have a task that can be done by brute force - if a specific map region has changed at a given scale, upload new geoTIFF at that scale, then rebuild the tile pyramid for that scale. (I don't create ...
2
votes
1answer
50 views

UnsatisfiedLinkError caused by org.gdal.ogr.ogrJNI.GetDriverCount()

I am developing in Java on Windows 7 with Eclipse (Luna). I am using gdal for some ogr processes and running the application from eclipse is fine, but running from command line on the same machine ...
0
votes
1answer
96 views

Problem using raster2pgsql to upload a raster to PostGIS

I am trying to upload a large number of 3 band rasters to a PostGIS database using raster2pgsql on the command line in Windows 7. Here is my command: raster2pgsql -I -C -s 27700 -F -t 100x100 -M ...
2
votes
3answers
100 views

gdal2tiles and qgis

I have made a tile pyramid of a number of Dutch topographical maps using gdal2tiles. This gives apart from the tiles openlayers.html and tilemapresource.xml. The openlayers shows perfectly and ...
1
vote
0answers
45 views

gdallocationinfo used in javascript with node-gdal?

I've been scratching my head for a white trying to figure it out if it's possible to get pixel information from my rasters using gdallocationinfo in javascript. I found node-gdal, but after reading ...
0
votes
1answer
24 views

Using ArcMap georeferencing link file with in GDAL or another alternative?

I've a lot of georeferencing link files from ArcMap 10 that are like so: 5.693887 17.605231 37.000000 -17.000000 5.896141 5.922077 37.000000 -19.000000 16.877354 5.957769 ...
0
votes
1answer
177 views

How to extract style data from DXF when converting to geoJSON

I am using ogr2ogr (GDAL) to convert DXF files (from autoCAD) to geoJSONS for usage inside Leaflet.js The issue is that a lot of information, like line stroke size and colors are not kept during the ...
3
votes
0answers
73 views

gdal_retile.py leaving gaps in mosaic

I'm trying to create a seamless mosaic from the FAA's VFR aeronautical charts using gdal_retile.py. The end goal is to serve the mosaic using the ImagePyramid extension in GeoServer. However, I'm ...
3
votes
0answers
28 views

Building GEOS / GDAL halts on embedding manifest files [closed]

I'm trying to build GEOS and GDAL from source on a Windows machine. I have Visual C++ for Python 9.0 installed with Lib / Include / PATH environment variables established (as far as I am aware). But ...
0
votes
1answer
66 views

GeoTiff image not aligned after using gdal_translate and gdalwarp

I don't see what I am doing wrong. The coordinate system was assigned to the GeoTiff using the following gdal_translate command: (Usage: gdal_translate -of GTiff -a_srs EPSG#### -a_ullr upper left ...
0
votes
1answer
79 views

Building GDAL with libcurl support on Debian Linux

Does anyone know how to build GDAL with libcurl support on Debian Linux (specifically LinuxMint and/or Ubuntu)? After following build instructions here and tip from @EvilGenius below, I'm specifying ...
1
vote
0answers
24 views

How Do I Split a Raster into parts using GDAL? [duplicate]

I have a large Raster with High resolution DEM of a large area, and I need to split it into smaller squares of 500mx500m, so that I can run other raster functions on it. Basically I'm looking for the ...
4
votes
0answers
122 views

How to remove nodata values from XYZ output with GDAL?

I have a buffered line route that is overlayed on a DEM: I've clipped the DEM to the buffered line, and exported it as an ADF ESRI grid file. I then used gdal_translate to convert it from the ...
4
votes
2answers
490 views

Convert multi-layer DWG into several KML or geoJSON files (one per layer)

I have some big DWG files, with 5 to 50 layers, which I draw and geolocated using AutoCAD (with some external modules). I need to convert them to KML or geoJSON so that the JS library can load them ...
0
votes
2answers
97 views

Gdal driver.Create TypeError

I have been writing a python script to convert an ascii grid into a geotiff file and when running the Create method in gdal i am receiving the following error message: ...
1
vote
0answers
81 views

WKID issue when adding raster layer from ArcGIS “mapserver” to QGIS via REST

QGIS docs here describe how to do this. I am trying to add a layer from this MapServer: gps.digimap.gg/arcgis/rest/services/JerseyBaseMap/MapServer When I try the following command from the Python ...