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

2
votes
1answer
60 views

GDAL set negative values to nodata or NaN

How can I set all negative values in a raster to NaN or no data? I've made a virtual raster to link data and coordinates from a HDF5 then using something like this: gdalwarp -geoloc -tr 3000 3000 -...
0
votes
1answer
49 views

Converting Raster Orthography from MrSID to GeoTIFF Creates Colourful Artifacts

I'm using geoserver 2.9.0 to display our orthography. I've managed to get the MrSID file displaying properly however it's very slow. From what i've read to increase the performance I should convert to ...
1
vote
0answers
56 views

Merging features by admin or sovereingt

I followed Mike Bostock's 'Let's Make a Map' tutorial to make a GeoJSON map of the EU starting from a 1:10M cultural admin0 shapefile from Natural Earth. The problem that I'm facing if that some of ...
2
votes
1answer
74 views

Cannot update GDAL on Ubuntu

I am trying to use GDAL_Rasterize and found a link via Google that suggested that there were bugs in older versions. I am running Ubuntu 14.04 LTS and 'gdalinfo --version' gives me 'GDAL 1.11.2, ...
1
vote
0answers
27 views

Create Diference TIFF image from two TIFFs

I am trying to use gdal_calc to calculate the differences between two almost identical images. The problem is that the images are tiffs with three bands (RGB) and the result is difficult to interpret. ...
0
votes
0answers
23 views

gdalwarp -t_srs EPSG:XXXX failed on data overlaping “Antimeridian”

I have program which process some image from satellite (crop and projection). For data overlaping "antimeridian" my process failed. This is the FootPrint of my image: Corner Coordinates: Upper Left ...
0
votes
0answers
18 views

What is the difference between ways of merging rasters? [duplicate]

Does anyone know the difference between the ways of merging rasters using QGIS and gdal_merge.py of (.tif) raster?
1
vote
2answers
27 views

Remove Color Attribute Table - GDAL - geoTiff

How can I remove the GDALRasterAttributeTable with GDAL? It interferes with styling in Geoserver. GDALInfo output below: Driver: GTiff/GeoTIFF Files: gcev3-product1-cropland-extent-australia.tif ...
4
votes
0answers
53 views

Cannot Open VRT of Sentinel2 image created by AutoCreateWarpedVRT

I am trying to create a warped VirtualRaster of a Sentinel2 image. Here the Python code I tried: from osgeo import gdal, osr srs = osr.SpatialReference() srs.SetFromUserInput('EPSG:32632') gdal....
2
votes
1answer
72 views

Converting a CSV to DXF using ogr2ogr from the command line

I am trying to convert a .csv file to .dxf using ogr2org from the command line using the following command: ogr2ogr -f DXF control.dxf opua.vrt This is the .vrt file I am using: <...
1
vote
1answer
20 views

How to reuse memory pointer of GDAL Memory driver

I am doing the same as in [GDAL API]: can't save image in some formats That is, saving a JPEG image with GDAL. I know I'd have to create a MEM buffer and then CreateCopy to JPEG. However, I'd ...
2
votes
1answer
34 views

What's the difference between gdal_cachemax and setting the cache size with -wm in gdalwarp?

I'm trying to get better performance out of gdalwarp, and am wondering if there is any difference between setting gdal_cachemax and using the -wm option in gdalwarp?
1
vote
1answer
27 views

Making a shapefile or kml tiles from a file that says it has no GCP

Learning my way through the GIS world, i was able to create a nice script to handle some geotiffed nitf files. Could create png, thumbnails and even kmls from these files with GDAL. Then I have ...
0
votes
0answers
8 views

Rasterize all layer attributes with gdal.RasterizeLayer

I want to rasterize a layer and to have a 2 bands raster in output where the first band is the first feature of the layer and the second one represent the second feature. I tried that but i have an ...
1
vote
2answers
30 views

Gdal vsicurl not working

I'm currently using gdal version 1.11.2 (also tested this on 2.0.0) to open AWS Landsat-8 data directly from its url. If I manually test the link in the browser, it works. If I use /vsicurl/ the ...
4
votes
1answer
173 views

PostGIS Out of db VRT containing JPEG2000 GDAL 2.1 OpenJPEG

Build GDAL 2.1 with OpenJPEG Build PostGIS with above GDAL build Created VRT file with SENTINEL2 Driver / JP2 files. gdal_translate SENTINEL2_L2A:...
0
votes
0answers
30 views

Creating superoverlay from geotiff covered whole globe, shouldn't of?

I am new to this GIS world and found this blog post: https://developers.google.com/kml/articles/raster#step-1-download-and-install-gdal Following its instructions I used my own geotiff (ntif) image ...
0
votes
0answers
11 views

Include partial grid cells with gdalwarp -cutline

I am using gdalwarp (v 1.10) to crop a large raster to the size required by my application using the command line utility with the following command: gdalwarp -cutline border.shp -crop_to_cutline ...
0
votes
0answers
29 views

How to rasterize points in python?

I have a problem rasterizing a vector layer base from Count values. I am getting a black image. Below is my sample code base fromhttps://pcjericks.github.io/py-gdalogr-cookbook/raster_layers.html. Any ...
1
vote
0answers
51 views

ogr2ogr not running on Ubuntu 15.10

I'm trying to convert an osm.pbf file to a spatialite db but I can't get ogr2ogr to run. I have qgis 2.14.2 and ubuntu 15.10. I've reinstalled gdal-bin and I'm getting the following message from the ...
0
votes
0answers
29 views

Out of memory when writing file on /vsimem with GDAL

I would like to write out some intermediate results to the memory using the GDAL library, to reduce the number of I/O operations required by my program. Simplified code: GDALAllRegister(); const ...
2
votes
3answers
112 views

GDAL error in R path.expand(path): invalid 'path' argument

Using windows 10. I would like to stick to using R since that's where a lot of my pre-processing and analysis will take place for a bigger project. I have downloaded a number of MODIS files for the ...
0
votes
0answers
128 views

“ QFSFileEngine::open: No file name specified ” installing QGIS Server on Windows 2012R2

I am attempting to get QGIS Server installed on a Windows 2012 R2 Server. After following the information at the https://anitagraser.com/2012/04/06/qgis-server-on-windows7-step-by-step/, and looking ...
0
votes
1answer
56 views

kml to shape file (ogr2ogr)

I am trying to convert kml to shape file. I'm using ogr2ogr (gdal version = 2.0.0). Here is the command : /home/vineet/Project/gdal-2.0.0/bin/ogr2ogr --config SHAPE_ENCODING UTF-8 -update -append ...
1
vote
2answers
64 views

Extraction of raster with shapefile python

I have a raster tiff file containing vegetation index & also a shapefile of fields. I want to extract the field wise values from that raster file using GDAL. How can I do this? I can read the ...
1
vote
1answer
42 views

Merge Natural Earth SQLite country polygons using ogr2ogr and ST_Union

I'm trying to create a map with the world continents (with a few modifications). As a source I use a sqlite file from Natural Earth. I already introduced my modification in the definition of the ...
0
votes
1answer
42 views

GDAL Warp giving garbage data on edges when transforming to WGS 84

I want to change one tiff file to USGS DEM file. For that I am doing following steps: Transform the tiff file to WGS84 coordiantes using Gdal Wrap utility gdalwarp -t_srs EPSG:4326 -dstnodata -...
1
vote
1answer
44 views

gdal_polygonize returns with a strange error

i'm trying to execute gdal_polygoinize.py command. I have a main python method that takes in input a raster file and a output shapefile, and then calls a method named polygonize.py that builds and try ...
0
votes
0answers
17 views

Land-Water mask for MODIS 0.05deg grid

I need to create a worldwide land-water mask for a series of MODIS LST/VI images (0.05 deg). I know that each MODIS hdf file contains a land-water mask. However, this mask changes between different ...
0
votes
1answer
25 views

NED (NAVD88) to WGS84 heights

I need to convert the 1/3 arc second DEM files from the NED into WGS84 ellipsoid heights. The documentation says that the NED is referenced to NAVD88. I've tried using gdalwarp as follows: gdalwarp ...
1
vote
0answers
63 views

How to Convert from GMT Line files to Shapfile for Qgis?

I have several line text-files of (used in GMT) saved in the format (this is a small sample): -111.29 42.28 -111.29 42.31 -111.29 42.33 > -111.37 41.99 -111.38 42.02 I am trying to convert this ...
4
votes
2answers
111 views

Convert MODIS HDF file in sinusoidal projection into GeoTiff with

I'm looking for a way to convert MODIS LST data from HDF to Tiff format, and reproject from sinusoidal to WGS84 - all using python GDAL. An example file is here (I only care about band 0) : http://...
3
votes
1answer
65 views

Changing the origin of a raster?

I have a raster file that is projected not detailed enough. Every pixel in the raster should be replaced, just a bit. With ArcGIS I calculated the offset per pixel. I couldn't find a way to do this ...
-1
votes
1answer
38 views

Gdalbuildvrt dosen't output the vrt file

I am trying to use gdalbuildvrt to create a time series of landsat images. I used the following code but it is not outputting the vrt file. In the code inlist is the txt file with path to input ...
0
votes
0answers
11 views

Resampling using gdalwarp gives me no data

I have an image with the following information. Driver: GTiff/GeoTIFF Files:test_DEM.tiff Size is 91, 73 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", ...
0
votes
0answers
40 views

Gdal ogr simplify/intersects

Can anyone give an example: How to use the simplification? For example // read the coordinates of the polygon char * wgsRef = NULL; OGRPolygon *poPol = (OGRPolygon *)...
1
vote
1answer
164 views

Installing GDAL2.1

Here will be my sequence of command lines while trying to install gdal2.1 in a UBUNTU virtual machine. My virtual machine is a UBUNTU 16.04 LTS(64bit) I would need gdal2.1 and especially the Python ...
1
vote
1answer
62 views

Unable to add MrSID extension to GeoServer

I'm trying to serve some MrSID files with GeoServer following the geoserver.org guide, however, I can't get the MrSID Raster Data Source to show up when attempting to add a new GeoServer Store. This ...
1
vote
0answers
54 views

Converting TIFF (DIMAP) to .xyz file format with GDAL

We are attempting to convert a .tif file into .xyz file. We have been using gdal_translate on a linux box and this is the line we've used: gdal_translate -strict -of XYZ ABC.tif ABC.xyz -mo ABC.dim ...
0
votes
1answer
31 views

Set PATH variables for ogr2ogr without admin privileges

In the office I use ogr2ogr to convert several file formats such as KML or ESRI shapefiles. Since the commands are the same almost every time I'd like to run ogr2ogr with a batch. I'm using the GDAL ...
1
vote
1answer
55 views

tiff to shapefile using GDAL - slow speed

I am converting a 500MB .tif file to .shp file using GDAL. I am using the following: gdal_polygonize.py .tiff file -f “ESRI Shapefile” out.shp However, I abandoned the first attempt as in 12hrs it ...
3
votes
0answers
35 views

GDAL gdal_polygonize.py TypeError

I am using GDAL to convert a .tif to .shp file. gdal_polygonize.py input.tif -f “ESRI Shapefile” out.shp However this does not work and results in: Traceback (most recent call last): File "/usr/...
1
vote
0answers
50 views

Is there any way to connect features (entities) of dxf to database in python?

I'm using gdal, ezdxf, sqlite python packages and sqlite ODBC driver. They work fine, and in AutoCAD I can easily manually connect to my dbase.sqlite (read only) using _DBCONFIGURE and _DBCONNECT ...
3
votes
0answers
37 views

How to blend edges of aerial imagery with gdalwarp

I need to fade/blend the edges of a series of aerial images with gdal. The closest thing that I've been able to accomplish is using the -cutline and -cblend switches of gdalwarp to get the edges to ...
4
votes
2answers
181 views

How to add support for FileGDB (Esri file gdb API) driver in fiona?

There has been a tons of questions on that but none of them solved my problem. I have a machine with: Windows 7 x64 Python 3.4.3 x64 gdal installed: gdalinfo --version --> GDAL 1.11.4, released ...
3
votes
0answers
96 views

GDAL RasterIO performance in C++ vs ReadAsArray in python/numpy

I'm currently writing optimisations for a raster viewer program which uses gdal as it's base. It's currently written purely in python, and has some major speed issues which cause problems when we are ...
2
votes
0answers
27 views

Qgis gdalogr:contour is not generating output

I am trying to use gdal:contour script in Qgis console and then script. Here is code which I tried write to console. raster = QgsRasterLayer('C:/Users/Bulva.user/Desktop/test_raster.asc', "DEM") ...
0
votes
1answer
45 views

Can gdalwarp reproject from ESPG:4326 (WGS84) to UTM?

I've a GeoTIFF image in ESPG:4326 (WGS84) and when I try with gdalwarp to reproject to UTM, it has problem gdalwarp -s_srs "EPSG:4326" -t_srs '+proj=utm +zone=33 +datum=WGS84 +units=m +no_defs' -of ...
2
votes
0answers
36 views

Non-zero error with ImportFromEPSG() in GDAL/OSR/Python

I'm currently getting an error when testing ImportFromEPSG() on my GDAL. import osr testSR = osr.SpatialReference() testSR.ImportFromEPSG(4326) #this returns value of 6 print testSR.ExportToPrettyWkt(...
6
votes
2answers
232 views

Transforming geostationary satellite image to lon/lat

The ICARE Data and Services Center distributes the SEV_AERUS-AEROSOL-D3 product in HDF5 format (http://www.icare.univ-lille1.fr/archive/?dir=GEO/MSG+0000/SEV_AERUS-AEROSOL-D3.v1.03/). The main ...