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
2answers
145 views

System call gdal function in R on Mac

Hopefully someone can help me figure this issue out. I am trying to system call gdal functions inside of R using the system / system2 functions. Here is an example of the code: > system("gdalinfo ...
2
votes
2answers
112 views

GDAL reads NEX Climate projection datasets as netCDF on Mac, HDF5 on Ubuntu

I'm trying to work with the climate projection data that NASA released, and is hosted publicly on S3: http://aws.amazon.com/nasa/nex/. The data is available at s3://nasanex/NEX-DCP30. The code I'm ...
2
votes
1answer
64 views

Why does GDAL Polygonize flip features vertically?

I am trying to convert a 1 band jpg raster into polygons. It is converting the data and I get a polygon shape file. However, the output data is vertically flipped. Input data have no reference ...
1
vote
0answers
35 views

How to set -srcwin with python gdalbinding and driver.createcopy

I am using gdal_translate currently to cut a 60000x60000 geotiff into smaller 1000x1000 geotiff files. This works but it takes about 24 hours(and I have 44 of those big images to compute). Im not ...
0
votes
1answer
22 views

gdaladdo adding stripes to data

I created some image overviews for a large (1.2 Gb) topo relief image. gdaladdo -r cubic -ro --config COMPRESS_OVERVIEW JPEG --config PHOTOMETRIC_OVERVIEW YCBCR --config BIGTIFF_OVERVIEW YES ...
0
votes
1answer
95 views

Help interpret GeoTif pixel value python gdal: connection between pixel value, RGB values and physical concentration?

I was running python script that queries a GeoTiff file at the coordinates of a centroid of blocks and as a result I got a list of pixel values, which I have trouble interpreting. When I run gdalinfo ...
2
votes
1answer
1k views

Install GDAL on CentOS 7

I'm having issues with installing GDAL on a fresh install of Centos7 using the ELGIS repo. Here's my commands and output: sudo rpm -Uvh http://elgis.argeo.org/repos/6/elgis-release-6-6_0.noarch.rpm ...
0
votes
0answers
35 views

GOES 13 EAST How to Georeferencing

I have a many image for Georeferencing, Goes 13 EAST this is metadata: GOES 13 Imager frame 78 at UTC 14:45:22.485 day 009 of 2015 Vis pixels : 4009, 24840 lines : 2621, 13440 Vis Lon : ...
1
vote
1answer
316 views

Changing color of raster images based on their data values - gdal

Let's say I have 1 band raster images with Gray Color Interpretation and I want to change their color based on their data values. Data values -1, will be Blue, -0.5 will be Light Blue, 0 will be ...
3
votes
2answers
316 views

Downloading raster data into python from postgis using psycopg2

I have raster data in a postgres table that I want to get into python as a numpy array. I am using psycopg2 to connect to the db. I am able to download the data but it comes back as a string (probably ...
1
vote
0answers
24 views

GDAL Warp cuttline not working

Has anyone ever got the cutline to work in C++? I have it working off the command line no problem, but it doesn't seem to do anything when I populate that argument in C++ as below: *//Create a ...
4
votes
1answer
105 views

How do I get GDAL_translate to output better resolution

I have a set of GeoPDFs that I converted to GeoTiffs. However, the end results are horrible. They are very pixelated. The command I passed was gdal_translate -of GTiff inFile outFile -co ...
0
votes
0answers
65 views

Mosaic several geotiff files

how to mosaic several geotiff file to submap that does not contain ALL the source maps I need GDAL C++ API?
0
votes
1answer
33 views

Code of _gdal.pyd

I want a deeper understanding of how gdal_polygonize.py works. It uses gdal.Polygonize and gdal.Polygonize refers to _gdal.pyd. Where can I see the actual algorithm where the actual polygonization ...
4
votes
2answers
1k views

Obtain coordinates and corresponding pixel values from GeoTiff using python gdal and save them as numpy array

How to obtain projected coordinates as well as the actual pixel values at those coordinates from GeoTiff file and then save them into a numpy array? I have arsenci020l.tif file, and its coordinates ...
0
votes
1answer
131 views

How do I create a projection file using GDAL

I am trying to obtain a buffer from a polyline file using GDAL with the following function def createBuffer(inputfn, outputBufferfn, bufferDist): inputds = ogr.Open(inputfn) inputlyr = ...
1
vote
0answers
47 views

How to change values in a raster file at particular positions?

I have a collection of coordinates which comes from a QGSGeometry object. The shape is a polygon. example: qgsGeometry [[(75.2172,17.6566), (75.8499,17.9613), (75.5218,17.2113), (75.2172,17.4574), ...
1
vote
0answers
82 views

QGIS Mac OSX version cannot batch write output files after rasterization

Happy new year everyone, I have been trying to batch rasterize a series of .shp files using QGIS for Mac OS X, and keep coming up against a strange problem. I am using the QGIS 2.6.1 user interface, ...
0
votes
2answers
211 views

How to subtract rasters using the gdal_calculate command

I just want to know how to subtract two raster using gdal_calculate command? Or what is the proper way of subtracting using this? Any idea would be greatly appreciated. Note: gdal_calc is different ...
1
vote
0answers
65 views

CentOS 6.5 - After compiling gdal.so: Unable to load dynamic library

CentOS release 6.5 (Final) Linux version 2.6.18-028stab101.1 (root@rhel5-build-x64) (gcc version 4.1.2 20080704 (Red Hat 4.1.2-46) I'm trying to Install Cartaro CMS (http://cartaro.org/). GeoServer, ...
5
votes
0answers
128 views

Working with rasters in file geodatabase (.gdb) with GDAL

I'm working on a tool that converts raster layers to arrays for processing with NumPy, and ideally I would like to be able to work with rasters that come packaged in a .gdb without exporting them all ...
0
votes
1answer
17 views

Rastercalc.py and Rastercalcengine.py

Is there any examples on how to use the rastercalc.py or rastercalcengine.py? And what are the purpose of these python module? Are these modules the same with gdal utilities which can be run in ...
0
votes
0answers
25 views

Error in warping the merged hdf4 files

I've tried to merge a set of hdf4 files(unknown projection) and it run well. But when I tried to warp the output of the merge hdf4 files an error like this appears: ERROR 1: Unable to compute a ...
3
votes
0answers
56 views

GDALRasterizeGeometries doesn't burn correct Values

What I have is a vector shape file with a string field. What I've done is written a CSV file with the field ID, field_value and now I would like to rasterize the vector layer as a byte raster with the ...
4
votes
2answers
341 views

How to convert float raster to vector with python GDAL

I have float raster and now I want to convert it to vector. How is it possible with python gdal library? I have tried with gdal_polygonize.py of gdal utilities in command line and it worked excellent. ...
2
votes
1answer
95 views

Creating Spatialite DB from Swedish OSM

Good afternoon, I have just started to use GDAL/OGR and I am trying to create a spatialite DB from OSM data. I am using a cut of the Swedish OSM data downloaded from Geofabrik in .osm format. I have ...
0
votes
0answers
30 views

increase gdal_polygonize's capacity

I have a very big PNG (82500x40900) with more than 5000 polygons. I want to use gdal_polygonize but the geojson output shows that gdal_polygonizer stopped after 4770 polygons and the geojson isn't ...
0
votes
1answer
57 views

Georeferencing jpg via Gdal_translate output a GIS tif 15.5 times bigger [duplicate]

Given as input a light 874kb jpg, I use gdal_translate -a_ullr (source) to reinject GIS georeferencing : -a_ullr ulx uly lrx lry : Assign/override the georeferenced bounds of the output file. ...
1
vote
0answers
70 views

GDAL “Permission Denied” [closed]

I'm trying to convert a file using gdal-translate and keep getting a "Permission Denied" message. Looking online (including here) the solution is usually to make sure you have write access to the ...
0
votes
1answer
49 views

Is gdal_polygonizer creating a wrong geojson?

I have this PNG and I use gdal_polygonize 0VxOhek.png -f "geojson" 0VxOhek.geojson This is the output. There is a "features"-array which contains many objects with "type", "properties" and "geometry" ...
0
votes
1answer
46 views

eliminate small polygons of a PNG with gdal_sieve

I want to eliminate polygons wich consist of less than 3550 px in this PNG . I tried using gdal_sieve but there seems to be an error. gdal_sieve -st 3550 E8.png -of PNG E8sieve3550_4.png ERROR 6: ...
1
vote
3answers
394 views

convert csv to kml or geojson, etc using ogr2ogr or equivalent command line tools (use “latitude” and “longitude” attributes for coordinates)

I've just unsuccessfully tried to use ogr2ogr to convert CSV into KML like so: ogr2ogr -f "KML" output.kml input.csv The csv has "latitude" and "longitude" columns. The attributes are all loaded ...
0
votes
0answers
259 views

How to fix “Unable to load GdalTools plugin. The required ”osgeo [python-gdal]“ module is missing.”

When I start QGIS, it display error message: Unable to load GdalTools plugin. The required "osgeo [python-gdal]" module is missing. Install it and try again. Couldn't load plugin 'processing' from ...
0
votes
0answers
113 views

Set coordinate origin to the top left corner

I have a BMP-picture and I know that it's top left corner has the EPSG:3068 coordinates 35500,24390. I would like to use gdal_polygonize and than calculate the EPSG:3068 coordinates of the produced ...
1
vote
1answer
219 views

How to create a mosaic in QGIS with cutline and feathering for Landsat-8 imagery

1 - I have two scenes Landsat-8 RGB 6/5/4 and I need to mosaic them: 2- This is my cutline in polygon format that will smooth the overlapping area: 3 - So, I need to create a mosaic with ...
2
votes
0answers
67 views

Help with Python/GDAL documentation

The instructions from the url http://gdal.org/gdal_tutorial.html give the following example: dataset = gdal.Open( filename, GA_ReadOnly) with the following description: It's interpretation is ...
0
votes
1answer
72 views

I can't configure authentication with Google Maps Engine driver from GDAL

I'm using last GDAL from trunk, and trying to read from Google Maps Engine (GME) vía GME driver I already exported GME_API_KEY: My API Key, for the app I registered GOA2_CLIENT_ID: Generated with ...
0
votes
1answer
78 views

gdal_polygonize ignores some polygons

I want to polygonize this PNG. I try gdal_polygonize pic.png -f "GeoJSON" polygonized.geojson but I get only 11 polygons instead of more than 20. What is the cause or how can I fix it?
0
votes
1answer
287 views

gdal_rasterize error: “attempt to create 0x0 dataset is illegal..”

I'm trying to use gdal_rasterize to convert a polygon shapefile into a raster, where the value comes from the 'COVTYPE' attribute. I try calling gdal rasterize like so (I'm calling from within a ...
2
votes
1answer
175 views

Install QGIS on Ubuntu with GDAL built from source?

I am trying to install QGIS on Ubuntu 14.04. While I can easily install it via the package manager or apt-get, I find that QGIS, when installed, does not have Python support as expected. The symptoms ...
3
votes
1answer
220 views

merging GeoTiffs with GDAL causes grid-like artifacts

I'm processing ASTER DEM tiles, merging them into larger blocks (to minimize visible tile edges) and importing the contours into a postgis database. Workflow: gdal_merge.py -of GTiff -o merged.tif ...
4
votes
1answer
207 views

Determine if shapefile and raster overlap in python

I am building a script in python using osgeo/gdal. I have a set of shapefiles and a set of GeoTiff raster files. I would like to have my script ignore shapefiles if they do not intersect with the ...
0
votes
0answers
73 views

Linear interpolation on raster data, save result (geoserver)

I am trying to perform a linear interpolation on a coverage layer (Source data is in EHdr Coverage format, file extension: .flt). Ultimately, I want to use the interpolated layer as input to a WPS ...
5
votes
2answers
176 views

GDAL (or other open-source) equivalent of Arc's Point To Raster Tool

I'm trying to convert a python script that relies on arcpy into one that only relies on open source modules, but I'm having some trouble. The bulk of my processing lies in one Arc tool "Point to ...
1
vote
1answer
34 views

How to use gdal's “forceToLineString()” with Python

I'm working with Python and GDAL and I'm quite new to this ^^ For my script I'm looking for a method to convert a MultiLineString into a simple LineString. Example: MULTILINESTRING ((0 0,0 1),(0 1,0 ...
3
votes
2answers
363 views

Gdal: clipping a raster with another raster

I'm writing a simple utility to crop batches of multi-band geotiff raster files to the same (smaller) area. Using gdalwarp, I can easily crop a file using a single-polygon clipping shapefile: ...
1
vote
0answers
97 views

Unable to call GDAL functionalities from Python

I am trying to call GDAL functionalities from the following script I made in the past. import os, fnmatch from subprocess import call call("dir", shell=True) # define the variables inFolder= ...
2
votes
0answers
86 views

QGIS polygonize GDAL DLL load failed

I'm using QGIS 2.6.1 and i've been having problems using the polygonize function. I keep getting a "DLL load failed" import error. I've searched older threads which unfortunately weren't able to ...
4
votes
1answer
228 views

Merge tif files with python

I'm also having problems with running gdal_merge.py within another python script and I cannot run it. I already checked a similar problem here. My example is as follows: First I set my workspace and ...
2
votes
1answer
336 views

Orfeo tool box (OTB) in Ubuntu 14.04 and SAGA incompatibility (gdal versions)

I need to test segmentation of a 12000 x 12000 pixels scene. I successfully compiled OTB in Ubuntu 14.04 and ran the process. But when I tried to start SAGA, I got a startup error. I found that gdal ...