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

Escaping the MODIS Swath Projection with GDAL: MOD06_L2: MRT Swath Tool not working

Stackexchange, I am looking at ftp://ladsftp.nascom.nasa.gov/allData/6/MOD06_L2/2012/222/MOD06_L2.A2012222.0800.006.2015062054108.hdf Trying to get the SDS (in WGS84 Projection and as GeoTiff): ...
0
votes
1answer
14 views

gdal_translate cannot convert MSG3 channels to TIFF

following this guide http://www.gdal.org/frmt_msg.html I want to convert my MSG3 HRIT format files in TIFF but I get this error gdal_translate H-00 ...
0
votes
0answers
6 views

Zonal computations: computing raster pixel counts within polygons using GDAL versus rasterstats

I'm computing some zonal statistics, and was comparing the results obtained using the GDAL cookbook example code here with the rasterstats zonal_stats function. I want to compute the mean value per ...
1
vote
0answers
24 views

Use Gdal from Python IDE

I am just starting to experiment with writing python scripts. I have been using QGIS and the Osgeo4w shell for sometime now. I would like to be able to call Gdal from the python IDE i have installed ...
2
votes
0answers
15 views

Does gdal_calc only support 26 input raster files at a time?

In gdal_calc.py, input files are specified with letters A-Z. Here are two files, input.tif and input2.tif, being averaged into result.tif: gdal_calc.py -A input.tif -B input2.tif ...
0
votes
0answers
6 views

Rasterizing a SSURGO shape file with GDAL — Choosing right number of Pixels

I am trying to rasterize a SSURGO shape file into UTM coordinates with 5m resolution. My first step gives me a GeoTiff with 3000x3000 pixels in lat/lon. gdal_rasterize -a MUKEY -ts 3000 3000 -l ...
0
votes
1answer
12 views

Gdal/OGR LayerMetaData not parsing its value

To start off with: Using Gdal/Ogr compiled with Xerces to read a WFS Reading from a WFS version 1.0.0 Compiler is QT version 5.4.1 ( C++ ) So I'm trying to get all the capabilities of a WFS via ...
0
votes
0answers
11 views

Calclulate the shade from a DSM and from sun elevation and azimuth angle

I have a mosaic of about 5000 tiles of aerial orthophotgraphy taken at different dates and thus with different sun angles. Furthermore, I have a Digital Surface Model (1 m resolution). My goal is to ...
0
votes
1answer
28 views

GDAL raster shifting with python error: 'NoneType' object has no attribute 'GetGeoTransform'

I'm trying to shift a georeferenced raster by 0.5m N and 0.5m W in QGIS using the python console: from osgeo import gdal # Open in read/write mode rast_src = gdal.Open('test_raster.tif', 1) # Get ...
0
votes
0answers
25 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
0answers
17 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, ...
0
votes
1answer
37 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 ...
2
votes
2answers
127 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
56 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 ...
0
votes
1answer
16 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 ...
1
vote
1answer
38 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
2answers
37 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
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 ...
1
vote
0answers
29 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
29 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
38 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
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 ...
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
46 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 ...
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), ...
1
vote
0answers
21 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 ...
1
vote
1answer
33 views

Tilemill and GeoTiff not working together

Loading in a geotiff with projection of EPSG:4326 isn't fully rendering: The original: The original in Tilemill: When using gdalwarp to change the projection from EPSG:4326 to EPSG:3857, ...
1
vote
1answer
22 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 ...
0
votes
0answers
34 views

Qgis 2.8, gdal python scripts gives error at end of process

I am using Qgis 2.8.1 64-bit (Standalone installer) which also packages a number or other applications, like OSGeo4W Shell, Saga GIS et.c. In the OSGeo4W Shell I am encountering a problem with some ...
0
votes
0answers
44 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 ...
1
vote
1answer
29 views

GDAL does not recognize EPSG code of a shapefilefile regrojected by ogr2ogr

I am using ogr2ogr to reproject shapefiles. For example: ogr2ogr -s_srs EPSG:1234 -t_srs EPSG:3857 dest.shp srs.shp ogr2ogr then creates dest.shp, dest.shx, dest.dbf and dest.prj as expected. The ...
0
votes
1answer
28 views

Convert from WGS84 to 2D projection

I have a map in the form of .bil & .hdr using WGS84. I'd like to convert this into another map using a 2D projection while keeping the lat/long information. What I'd really like with this new map ...
0
votes
1answer
32 views

Error 4 in rasterizing a polygon shapefile

I'm using QGIS Wien 2.8.1 trying to convert a polygon shapefile to a raster, but I get the following error in a qgis-bin pop-up: ERROR 4: `' does not exist in the file system, and is not recognised ...
0
votes
1answer
36 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
34 views

NAD83 2011 AZ central projection

I am having trouble getting GDAL to read a file geodatabase with data projected in NAD83(2011) State Plane Arizona Central FIPS 0202. I have developed the following code: # standard imports import ...
2
votes
0answers
64 views

How does gdal_calc numpy operators work?

Gdal_calc manual Gdal_calc is introduced as : gdal_calc.py - Command line raster calculator with numpy syntax gdal_calc.py [-A <filename>] [--A_band] [-B...-Z filename] [other_options] ...
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. ...
0
votes
0answers
52 views

How to write features in file geodatabases (FileGDB) with GDAL and Python?

I am trying to copy features from a spatiallite database to a filegeodatabase. The target FileGDB already exists with the required feature classes. I am using the FileGDB-API 1.4 on Ubuntu and Python ...
0
votes
0answers
14 views

How can I merge multiple grib2 files into a single file with a time dimension that I can upload to Geoserver?

I have a large number of grib2 files representing the same grid for different times. I want to be able to select an arbitrary amount of them and merge them into a single grid with time dimension. That ...
0
votes
1answer
52 views

QGIS downsample then upsample raster

I'm trying to downsample a raster then upsample it again. I am using the gdalwarp function in QGIS. When downsampling the raster looks OK in the front end, I set the cell size to 10x the original ...
1
vote
0answers
38 views

Convert grey,green,blue to red,green,blue tiff

==== Edit ===== Issue tracked to strange behaviour of gdal_warp when -dstalpha is specified and source is RGBA. I had assumed that the grey,green,blue tiffs (see example below) were genuine and ...
1
vote
1answer
65 views

Extract data from netcdf using a shapefile polygon in python

I need to subset data from a NetCDF using a specific shapefile. The data are sea surface temperature and ocean color at 1/4 degree resolution. I have 4 polygons describing the US. Northeast ...
0
votes
0answers
14 views

Problem Installing GDAL extension for Geoserver - java.lang.UnsatisfiedLinkError

As the title states im getting an error trying to install the GDAL extension for Geoserver 2.7 (My goal here is to enable MrSID support). The complete error message is: ...
1
vote
0answers
57 views

Compression Algorithms output files larger than the originals. Why?

I'm currently attempting to write a script that will compress a few hundred .tif files. I've tried two approaches, first I wrote a script in Python using CopyRaster_management() function and the ...
0
votes
1answer
19 views

gdal2tiles tile origin differences using Windows/OSGeo4W and Linux/GDAL

I am getting different X and Y values for tiles from gdal2tiles.py when in a PC (OSGeo4W - GDAL 1.9.2) environment and Linux (GDAL 1.10.1). It seems to me that the newer 1.10.2 x's and y's are wrong ...
0
votes
1answer
15 views

Grey gdal2tiles map output possibly caused by incorrect color interpolation

I'm working on a project that involves generating tiles from several geotiffs, and I've found that several of them are producing these odd, greyish maps that look almost, but not quite like what ...
0
votes
0answers
20 views

GDAL Python and postgis raster?

I'd like to create a postgis raster using GDAL and python. I can get the driver easily: drv = gdal.GetDriverByName("PostGISRaster") But how can I create a new raster table? I can't find any ...
9
votes
3answers
184 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 ...