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
1answer
355 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
325 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
197 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
160 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
142 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
400 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 ...
2
votes
0answers
127 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 ...
1
vote
1answer
339 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
1answer
237 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
80 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
59 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 ?
3
votes
2answers
165 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 AlphaList=["A","B","C","...
3
votes
1answer
74 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
78 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), ...
3
votes
0answers
156 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 lwgeom_version(...
1
vote
1answer
103 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
94 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
1answer
201 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
163 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 ...
2
votes
1answer
211 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 ...
1
vote
1answer
228 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
83 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
575 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] OPTIONS:...
0
votes
1answer
154 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
0answers
393 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
1answer
179 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
215 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 ...
3
votes
1answer
822 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 ...
1
vote
0answers
293 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
104 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
103 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 they'...
0
votes
1answer
98 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 ...
10
votes
3answers
486 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 ...
2
votes
1answer
71 views

How to average gdal_hillshades?

For some unknow reason, my average of 3 hillshades is very dark : gdaldem hillshade input.tif hillshades_C.tmp.tif -s 111120 -z 5 -az 315 -alt 60 -compute_edges gdaldem hillshade input.tif ...
0
votes
1answer
95 views

How to not process when query is outside the raster extent ? (Error when Completely Outside)

Given a West North East South bounding box, I want to crop this data from a hundreds of inputXX.tif. Nearly every times but once, the query first print out that it "falls completely outside raster ...
1
vote
1answer
112 views

In a GeoTiff, what unit are extents in?

I was looking at the GDAL metadata of a USGS GeoTiff quad (o48118e8) with the GdalToTilesWin program. It dumped the following info: Name o48118e8.tif Size X/Y 5315/6780 ResolutionX/Y 2.4384/2.4384 ...
2
votes
1answer
248 views

gdal_calc works but I get a Python error at the end of each process

I have a problem while using gdal_calc by a .bat file. Here is the code : call "C:\OSGeo4w64\bin\o4w_env.bat" gdal_calc -A"D:\tuilage_script_1\tuilage\ndvi\4_810000_6370000_820000_6360000_ndvi....
3
votes
3answers
1k views

How to create a compressed GeoTIFF with the QGIS Clipper tool?

I have a geo-tiff created by Pix4D which is about 375 Mb in file size. I have clipped it to a smaller area using Qgis 2.8(using Raster ‣ Extraction ‣ Clipper), and the output file is about 900 MB in ...
0
votes
0answers
171 views

GDAL/OGR - Ogr_wrap.dll entry point not found

I am using the latest stable compiled binaries of GDAL 1.11 from gisinternals.com with a C# project in Visual Studio 2013 and im having some problems with the SetSpatialFilterRect/SetSpatialFilter-...
1
vote
2answers
173 views

Problem with contour lines thickness in bigger zoom levels

I'm a beginner user of GDAL tools. I'd like to rasterize contour lines into Tiff file (for example such region: N47E006), everything works for quite small zoom levels e.g. 10... but the problem is ...
0
votes
0answers
84 views

gdal2tiles only rendering half of tif

I am using gdal2tiles.py to render map tiles of some data. The file I use is generated by combining a .tif and .tfw using gdal_translate XX.tif XXn.tif Whenever I open up the file in QGIS, I see ...
2
votes
1answer
129 views

Perform affine transformation on shapefile with OGR, GDAL, or other command line tools?

Transformations with WKT or EPSG code are made easy with ogr2ogr, but are there tools for me to define my own affine transformation to apply to the coordinates of a shapefile? Or, in general, is ...
0
votes
1answer
207 views

How to get projection data (.prj) of raster image using GDAL/OGR?

I start with a boundary shape-file (think U.S. states) in lat/long that I convert to some projection, and then rasterize the result as a .tif. What I need now are the raster coordinates of the ...
0
votes
1answer
264 views

“Too many points” when using gdal2tiles.py in EPSG:4326

I want to put a custom map onto cartodb, which, I believes, requires EPSG 4326. When I try making tiles with gdal2tiles, using a 4326 georeferenced image, I get the following error: ERROR 1: Too many ...
1
vote
0answers
83 views

GRASS GIS r.slope.aspect write errors

I am running into an error running r.slope.aspect I have never seen before. When I run it, I get the following output: r.in.gdal complete. Raster map <dem> created. WARNING: map [slope] - ...
0
votes
0answers
188 views

gdal_merge.py import gdal error

normal Gdal utilities work. gdal utilities with .py extension do not. they had worked very recently, within the last two months. checked the path variable and python version called from the command ...
1
vote
2answers
56 views

gdal_calc : fails with value 0?

I working to invert values of a one-band tif files (download sample). If i use 254 it works: $ gdal_calc.py -A ./shadedrelief.tmp.tif --outfile=whited.tmp.tif --calc="254*(A>180)" 0 .. 10 .. ...
0
votes
1answer
54 views

How to pass value of a band to and other band and change its value?

Using gdaldem hillshade, I get a one-band hillshade (download here, 2MB) : gdaldem hillshade crop_xl.tmp.tif Shadedrelief.tmp.tif -s 111120 -z 5 -az 315 -alt 60 -compute_edges I currently can ...
1
vote
1answer
112 views

gdal_contour is missing lines below a certain elevation

I have a very strange problem with gdal_contour. I have a very large GeoTif of Europe and would like to generate contour-lines for that region. Everything worked fine when I had lines every 25m but ...
7
votes
2answers
345 views

Migrating from ArcGIS to gdal [closed]

I would like to collect a table of equivalent or close enough gdal based functions to replace ArcGIS functions (esp in spatial analyst). Are there functions available using gdal to perform the ...