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
0answers
8 views

GDAL: defining nodata value for gdal_calc.py

This old post exactly describes my question: http://lists.osgeo.org/pipermail/gdal-dev/2011-September/030281.html I want to copy the numpy ma.masked_outside functionality to mask values in my ...
1
vote
1answer
134 views

geotiff conversion to coordinate system with gdal

I would like to get the coordinates X Y Z from a geotiff file. I am using the gdal translate for the conversion. Which output format is the more adequate for containing a 3 numpay array X Y Z ?
1
vote
0answers
16 views

RuntimeError opening rasterio array in GDAL

I've opened a GeoTiff with rasterio: import rasterio with rasterio.open('Africa.tif') as src: transform = src.meta['transform'] array = src.read(1) array([[203, 203, 203, ..., 0, 0, ...
2
votes
2answers
479 views

How to create contour from NetCDF?

I have downloaded Bathymetry data from GEBCO. Its in NetCDF format. I want to create contour from it. I tried so many tools to convert the file into raster but couldn't do so. How can I create raster ...
1
vote
1answer
28 views

Installing GDAL with HDF4 driver on Debian system

I am working on a Debian 7 Wheezy system using Anaconda and working from Spyder. I've had a hard time installing a version of GDAL that reads HDF4 drivers. Could someone walk me through the ...
1
vote
3answers
112 views

GDAL HDF4 Driver on Linux System

I'm running code that works with MODIS HDF4 datasets on a Debian Linux distribution. I am also running this code on a Windows machine. On both the Windows and the Linux machines, I have the most ...
0
votes
0answers
27 views

Vector to Raster using GDAL (gdal.RasterizeLayer) : Error in Output Raster(generating output raster with NAN values)

For rasterizing a vector layer, I have tried gdal.RasterizeLayer method for converting vector(shapefile) to raster(tiff). But its giving raster output with NAN values(complete black image). I need ...
0
votes
2answers
214 views

Conflicts while installing postgresql9.3-postgis2.1,QGIS and gdal-bin

I want to use postgresql9.3-postgis2.1 and gdal-bin parallel . I have refer instruction here for installing postgresql-9.3 on ubuntu 12.04. I have install qgis from here While installing gdal-bin it ...
5
votes
2answers
270 views

Raster: mosaic first, or project first?

I have a thousand raster tiles in un-projected decimal degrees from which to build a seamless mosaic in a metres-based coordinate system. The final image will be about 70,000px by 50,000px. Which is ...
0
votes
1answer
47 views

Facing a problem when re-projecting a file?

I am trying to re-project from Equi-rectangular. The file is provided as HDF5. I wonder why half of the map is not projected correctly. gdal_translate -of "ENVI" -gcp 0 0 0 90 -gcp 1440 0 360 90 -gcp ...
2
votes
1answer
482 views

How to re-project the EASE (Equal Area Scalable Earth) grid with a ~25 km cylindrical projection to WGS84 0.25 degree?

I have nc files. from the metadata ,the projection is cylindrical and the resolution is 25 km: indrical NC_GLOBAL#ease_resolution=25 ...
1
vote
1answer
40 views

gdal RasterizeLayer gives different extent to arcpy PolygonToRaster

I have an existing arcpy script which uses PolygonToRaster_conversion to convert a polygon shapefile to raster (tiff). I am now exploring GDAL/OGR and have used the cookbook approach to carry out the ...
2
votes
2answers
119 views

How to make a movie from a ascii grid files?

Lets say I have several ascii grid files for each time step for a particular location. I would like to spit out the image for each of these ascii grid files and legend on the side bar. I also would ...
1
vote
1answer
44 views

ERROR 1: None of ECW_ENCODE_KEY and ECW_ENCODE_COMPANY were provided

I downloaded ERDAS-ECW_JPEG_2000_SDK-5.1.1 and tried to use it with GDAL to create ECW files from georeferenced PNG files, but they are not created, it basically complains about: ERROR 1: None of ...
1
vote
1answer
356 views

Using GDAL to read dted files

How do I get the elevation at a certain point using gdal to read in a dted file? I don't see much useful information out on the web. I'm sorry I should have been more specific. How would I ...
1
vote
1answer
36 views

GDAL doesnot support HDF4 e [duplicate]

I am trying to import the pymodis library within python and I keep getting error that "GDAL installation has no support for HDF4, please update GDAL" I am running the latest version of GDAL ...
1
vote
0answers
209 views

GeoPDF Support in GDAL

I installed OSGeo4W using a Windows installer... it contains GDAL 1.7.3. I can't get GDAL to recognize PDF as a file type using commands in MSYS shell. Does this version of GDAL not support PDFs? ...
2
votes
2answers
555 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
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, ...
4
votes
2answers
676 views

How to add Georeferenced PDF as layer to QGIS 2.0?

I have some georeferenced pdf's that I would like to open in QGIS. I have read that GDAL version 1.8 and higher will support pdf, and I have GDAL 1.10 installed with QGIS 2.0 on Windows 7, but when I ...
0
votes
1answer
12 views

creating a pyramide tiles of a geotiff using gdal

I have a big geotiff file (about 300MB), and I saw an example of a big data organized in folders : folder 0 contains one folder 0 which contains one .png file named 0.png folder 1 contains two ...
0
votes
3answers
91 views

gdal_translate : projection error?

I am using GDAL library routines to build a tiled hierarchy of images for use with Google Earth. The images are just PNG files created using Matplotlib. To tile the images, we are using ...
4
votes
1answer
169 views

Merge tif files with python

I'm also having problems with running gdal_merge.py within another python script and I cannot run it. I already checked a similar problem here. My example is as follows: First I set my workspace and ...
17
votes
2answers
397 views

How to create transparent hillshade?

I get elegant transparent hillshades via a combination gdal and convert. Compare to grey-based hillshade commonly used, such transparent hillshades are very cool because they can be placed between the ...
0
votes
2answers
127 views

How to crop Geotiff image in C++ (GDAL or other possibilities)?

I am developing a code to process a geotiff image. But firsly I need to crop a geotiff image preserving its metadata (projection...) in C++. My idea it's to provide the corners and the geotiff image ...
0
votes
1answer
36 views

Can't get the correct version of GDAL installed to use ReadAsArray without getting error in Python

Using the instructions here . . . https://pythongisandstuff.wordpress.com/2011/07/07/installing-gdal-and-ogr-for-python-on-windows/ I downloaded the following files: GDAL-1.11.1.win-amd64-py2.7.msi ...
1
vote
2answers
34 views

How to define PROJ.4 string for GDAL projection

I'm doing some zonal computations in GDAL using the code posted in the solution here: Issue Trying to create Zonal Statistics using Gdal and Python. It runs fine for me using the example data posted ...
0
votes
1answer
120 views

Combine overlapping rasters in gdal_calc.py, taking the maximum value from each pixel?

I want to combine 3 overlapping rasters to create 1 raster. Each have the same dimensions, resolution and CRS. I want to populate the new raster with the maximum value from each of the layers. ...
0
votes
1answer
61 views

OSGeo4W and GDAL - gdal2tiles.py error

Wondering why gdal2tiles is throwing errors, I've added the gdaldata and gdalplugins environmental variables and my python is at Python27/ArcGIS10.3/, and am wondering if that is messing with the ...
0
votes
1answer
272 views

gdalnumeric in python virtualenv [closed]

I'm having a problem importing gdalnumeric when working on virtualenvs. I'm getting the following error: from osgeo import gdalnumeric Traceback (most recent call last): File "<stdin>", line 1, ...
0
votes
0answers
14 views

Which command is used to check the coordinates and resolution of each pixel for raster binary data(.dat) format?

I downloaded the modis composite data(8day) of MOD09A1 product mainly band 4 and 6 in HDF-EOS format, did the conversion from HDF-EOS to header format(hdf) and then merged all the hdf into one single ...
4
votes
1answer
535 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 ...
1
vote
0answers
25 views

GDAL equivalent to ArcGIS Union?

I have a single shapefile full of overlapping polygons. The ArcGIS Union tool will generate individual, non-overlapping polygons with that geometry. Is there a way to do this with GDAL/OGR or another ...
0
votes
1answer
28 views

where to download gdal_polygonize.py?

I am trying to help out a friend of mine. She learnt GIS in school in the past using ArcGIS but has recently installed QGIS package on her Mac OSX. When she tries to render raster into shape file, ...
0
votes
0answers
13 views

Create a list of affected tiles from polygon / Gridding a polygon

Given a projection, a zoom level, a grid definition and a coverage polygon - what is the easiest/fastest way to obtain a list of map tile indexes touched by the polygon? Are there any existing tools, ...
0
votes
1answer
35 views

gdal_array.SaveArray() leaves dataset open in Python

When I use gdal_array.SaveArray() to create a raster, the newly created dataset appears to stay open in Python, preventing other processes from working with it. For instance, consider the following ...
1
vote
1answer
28 views

Getting Feature names from WFS in C++

to start off with: Coding in C++ Using WFS Version 1.0.0 Using QT 5.4.1 Windows 8 Latest versions I'm trying to read in a WFS in code via Gdal compiled with Xerces - I'm able to read the GML ...
1
vote
1answer
254 views

Changing color of raster images based on their data values - gdal

Let's say I have 1 band raster images with Gray Color Interpretation and I want to change their color based on their data values. Data values -1, will be Blue, -0.5 will be Light Blue, 0 will be ...
2
votes
1answer
277 views

How to remove strange gdaldem hillshade artifacts?

I've downloaded 9 chunks of 1 arc-second NED in IMG format from the National Map. I'm trying to create a hillshade from them. I'm using: GDAL 1.10.0, released 2013/04/24 First, I combine them and ...
1
vote
2answers
36 views

Raster Caculator Qgis 2.6 v/s Gdal_calc.py by Terminal

I need calculated monthly mean of Snowcover data by MODIS product, and then modify this script gdal_calc.py in this section # create alphabetic list for storing input layers ...
3
votes
1answer
56 views

GDAL can i try to create a workflow?

When i try to compose a workflow im confused to the funcionality of that. I want to create a workflow to execute commands gdal. what syntax I can use if I wanted to create a .bat file? For example ...
1
vote
2answers
3k views

Get extent of Georeferenced Rasters in Python and output to polygon shapefile

I would like to create a shape file containing the extents of each of the rasters in a directory. Is it possible to capture the extent of a raster using Python? I have tried ...
0
votes
1answer
15 views

How to update an existing GeoJSON with Fiona

With the code below I can create a new GeoJSON and add data to it: data =[(LINESTRING (-98.98234 19.38234, -98.990341 19.342345, -98.992449 19.378621, -98.991299 19.378233, -98.992349999 19.377381), ...
9
votes
3answers
179 views

How to create composite hillshade?

The concept of composite hillshade is to combine multiple hillshades with different sunlight orientation to avoid blind spot. It requires to merge 3 hillshades images with 315⁰ (NW sunlight, default ...
0
votes
1answer
21 views

How to change band type?

I have an hillshade tiff for which the current band type=Bytemay be problematic. How to change the band type to type=float ?
2
votes
2answers
264 views

Driver Error : DTM to Raster using GDAL_GRID

I am looking to interpolate a DTM (eventually over 4,000 DTMs) that is in a personal geodatabase format. I would like to do this in as few step as possible, eventually writing a batch script to do ...
0
votes
1answer
43 views

Import BigTiff files in Grass

I am trying to import a raster (geotiff) in Grass with size of 64 GB . It makes me an error like this: ERROR 4: This is a BigTIFF file. BigTIFF is not supported by this version of GDAL and libtiff. ...
1
vote
1answer
21 views

ogr2ogr does not set osm_id for all multipolygons

I used ogr2ogr to convert from osm data (pbf) to spatialite : ogr2ogr -progress -f SQLite -dsco spatialite=yes france.simple.spatialite france-latest.osm.pbf then, a large part of the imported data ...
1
vote
0answers
20 views

Enable LWGEOM for OGR tools (QGIS 2.8)

I am doing some scripted operations on a Spatialite DB using OGRINFO, OGR2OGR etc. in the QGIS 2.8 OSGEO4W shell. Unfortunately, LWGEOM functions such as ST_Split() are unavailable and ...
0
votes
0answers
40 views

Incorrect input type geoserver, gdal

I have been working with gdal, mrsid. I have successfully built gdal 1.11.2 with MrSid plugin and everything else and it works without problem when i deploy my webapp on tomcat. But when i deploy ...