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

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
73 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
145 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
100 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
91 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
184 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
148 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
186 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
218 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
79 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
533 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
144 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
375 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
174 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
200 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
749 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
285 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
99 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
96 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
1answer
92 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
467 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
66 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
91 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
110 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
217 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 ...
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
158 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 ...
1
vote
2answers
168 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
77 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
127 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
194 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
235 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
82 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
171 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
51 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
106 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
344 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 ...
19
votes
3answers
1k 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 ...
3
votes
1answer
43 views

Reprojecting to same srs shifting the layer

I am using ogr2ogr to reproject a shapefile. Here is what I have done. I have a shapefile A.shp. I used the following command to create B.shp: ogr2ogr -t_srs EPSG:3857 B.shp A.shp Now I reprojected ...
0
votes
1answer
68 views

Merging SHP lines with attribute data combined?

I have one shapefile which contains hiking routes and another shapefile which has ski routes. Sometimes they follow the exact same path. In QGIS if the lines mentioned above are in the same SHP I ...
1
vote
2answers
520 views

Generate PNG/JPG maps from multiple ASCI/GeoTIFF in Python/GDAL

I need to generate maps (saved in PNG, JPG or PDF) from raster files (for each file - one map) which meets following criteria; * each map has exactly the same color pallette, so that it will be easy ...
0
votes
0answers
184 views

Creating geoTiff from MBTiles

I have mbtiles file with various zoom levels (8-14) and I want to create geotiff file from it. I think that this could be acomplished with gdal but I dont know how. Any suggestions on how to complete ...
2
votes
3answers
569 views

Using GDAL/C++ to calculate distance in meters

I'm using GDAL in C++, with WGS84. I want to calculate distance between two points in meters. The code snippet below tries to find the distance between 25N/75E and 25N/76E. The code, as is, returns ...
1
vote
0answers
130 views

Suppress verbose output when something goes wrong with ogr2ogr

Whenever ogr2ogr is not able to read a file, it seems to produce pages of output: Unable to open datasource `file.geojson' with the following drivers. -> ESRI Shapefile -> MapInfo File ...
0
votes
1answer
576 views

How to setup GDAL setting for python using OSGEO4W?

I have installed OSGEO4W installer. In that I have installed Express Desktop install option. QGIS GDAL uDig and GRASS GIS were installed in C:/OSGEO4W64. I have added 'C:/OSGEO4W64/bin' in the path. ...
1
vote
1answer
631 views

How Does Ogr2Ogr Simplify Work?

Can someone explain how the ogr2ogr simplify works? I am trying to use raster data which has a spectral resolution of 30m. The vector file that I am trying to import into Google Fusion is large in ...
0
votes
1answer
56 views

Error in using gdal java binding

I am using the java binding of gdal to project a set of coordinate points from UTM to lat/log using the CoordinateTransformation Class of gdal. I have written the code as follows: ...
1
vote
2answers
440 views

gdalwarp changing value in nodata pixels

I'm using gdalwarp to convert MODIS NDVI mosaics from Sinusoidal to epsg:4326. To save space I'm using VRT files until I reach the final product. Here is my workflow. Mosaic NDVI Modis tiles using ...
1
vote
1answer
35 views

Most efficient way to retrieve pixel distributions inside many polygones (custom zonal stats)

I want to get the distribution of pixels from a 1 band raster that fall in each polygon of a shapefile to proceed to some statistics. I guess a zonal statistics does it in some kind of way but I need ...