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

ogr2ogr failing to open/convert GML

I'm trying to convert a GML file into GeoJSON using ogr2ogr, I've compiled GDAL, Expat and Xerces from source, but when I run my command I get an error. My convert command ./ogr2ogr -f "GeoJSON" ...
1
vote
0answers
18 views

QGIS Georeferencer adds padding to image after save

I am attempting to georeference an aerial image using satellite imagery. I'm successful in getting an output, however I notice there is an issue where QGIS's georeferencer is adding padding to all the ...
6
votes
2answers
120 views

How to convert the geographical projection of an extremely large file (GDAL)

My Task I need to convert 9 x 14Gb GeoTiff files to WGS84 geo projection. The rough size in pixels is about 200k x 100k. The data-type is byte and it is not a subject to interpolation, it is some ...
0
votes
2answers
40 views

It's possible to use gdal_calc.py with an entire path of rasters?

I try to use gdal_calc.py to scale and transform a time-series of MODIS images. I've done that when i worked with a single image in the -A variable and a mask in -B, but i need to work with a ...
1
vote
1answer
39 views

Read raster at request size through Python

I'm trying to read the raster (tif) file using Python. What I want is to read the raster at a specific request size rather than its original resolution. Example: If a raster has a cellsize of 5 x ...
1
vote
1answer
23 views

Grass r.out.gdal, wrong values when importing to QGis

I have a problem with export/import functionalities between Grass and QGis. I create a raster file in Grass (with r.wind.sun, but that most likely does not matter) It has the following properties: ...
0
votes
0answers
13 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
25 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 ...
2
votes
0answers
49 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
50 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
25 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
12 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
13 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
28 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
48 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
53 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
26 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
55 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
156 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 ...
2
votes
2answers
103 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
22 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
51 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
75 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
18 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
41 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
53 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
21 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
54 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
34 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
25 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
64 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
58 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
25 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
47 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
43 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
37 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
50 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 ...
1
vote
1answer
47 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
37 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
45 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
52 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
40 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 ...
3
votes
0answers
109 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
47 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
67 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
26 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
65 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
53 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 ...
2
votes
1answer
128 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
25 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: ...