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

1
vote
1answer
100 views

Load the raster image which is in TIF format using QGIS API

I am quite new inthe field of GIS. I am interested in is to load the raster image which is in TIF format using QGIS API I have followed the steps given the in the follwing link to load shape file ...
1
vote
2answers
100 views

How do I convert a GeoTiff's neatline metadata from vector format to pixels using GDAL?

I have a GeoTIFF that produces the following when run through gdalinfo: Driver: GTiff/GeoTIFF Files: PAA_1.tif Size is 14400, 6600 Coordinate System is: ...
1
vote
0answers
18 views

gdal_polygonize floating values [duplicate]

gdal_polygonize converts whatever rasters values into vector integers as it is based on the GDALPolygonize() function (http://gdal.org/gdal__alg_8h.html#a3f522a9035d3512b5d414fb4752671b1). It is ...
1
vote
0answers
42 views

Using OGR2OGR with LIBKML (Examples please)

I am trying to figure out how to use GDAL (ogr2ogr) to export spatial data from SQLServer to KML. The standard KML driver works, but it is not very robust. LIBKML has so much promise (layers, styling, ...
0
votes
0answers
81 views

How to Read GDB from c#

I am trying to read .gdb file using osgeo api i go for following code but cant read. suggest me. string File= @"D:\GDB File\Demo.gdb"; public IConnection OpenFDOSourceConnection(string File) ...
0
votes
2answers
78 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
83 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
42 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
19 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
18 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
64 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
546 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
92 views

Offset problem with GeoTIFF and Google Maps

I have several GeoTIFFs generated with GDAL, and just some of them present a weird offset shift of a range of ~5 km in average. The range of these problematic rasters is as much as -85/-70/-20/0. ...
0
votes
0answers
27 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
162 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
194 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
20 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 ...
3
votes
1answer
73 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
48 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
30 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
632 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
99 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
46 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
65 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
148 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
45 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, ...
2
votes
0answers
75 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 ...
3
votes
0answers
42 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
213 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
64 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
26 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
46 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
52 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
44 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
40 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
190 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
206 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
71 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
147 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
62 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
68 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
45 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
144 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
105 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
168 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
116 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
51 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
1answer
108 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
29 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 ...