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

5
votes
1answer
424 views

How to read Geotiff tags and pixel per pixel values using GDAL C#

I'm new here and I'm new at programing with GDAL C#. My question is how I can get the value of SRTM GEOTIFF pixel to compare with it's neighbors, and how can I get the values of Geotiff tags.
5
votes
1answer
399 views

GDAL.RasterizeLayer not ported to the C# bindings?

I am trying to use the OSGeo.GDAL.Gdal.RasterizeLayer(…) from GDAL 1.9.2 and it requires some SWIGTYPE parameters: OSGeo.GDAL.Gdal.RasterizeLayer( OSGeo.GDAL.Dataset dataset, int bands, ...
4
votes
1answer
52 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 ...
4
votes
1answer
571 views

How to create tiles with points from MongoDB?

I have a (potentially) large data set stored in MongoDB. Each record is just lat, lon and some float-point value. Each pair of lat and lon is unique, but the whole set is not covering the Earth, there ...
3
votes
1answer
220 views

install gdal extention in GeoServer

I am trying to install the GDAL extension (I want to use ecw files in GeoServer). So here is what I did so far: 1) downloaded the GDAL extension and extracted to: C:\Program Files ...
0
votes
0answers
33 views

Using GDAL and Python to replace the Spatial Analyst TabulateArea method

Has anyone tried to phase out the use of ESRI's Spatial Analyst extension by replacing it with GDAL in Python? I have a script that only uses the TabulateArea and Statistics analysis but not much else ...
0
votes
0answers
28 views

gdal_grid data from 2 different utm zone

I'm developing a component which grid data over Vietnam region using gdal_grid from gdal. Follow gdal example I've successfully griding data within 1 specific UTM zone. However, for data over the ...
0
votes
0answers
98 views

MODIS products reprojection problem

I am working with Modis Data (MOD02, MOD11, MOD03, MOD35). In particular I am looking for a solution to reproject the subdatasets below, into albers equal. Because I am using NLCD dataset which uses ...
0
votes
0answers
256 views

GDAL HDF not recognized as supported file

I am having the same issues as in this thread; however, I'm new to programming and need more advice on how to get gdal to open hdf files. I downloaded hdf from MODIS, which I believe are hdf 4. I have ...
0
votes
0answers
68 views

need help creating layer with geom type wkbGeometryCollection with gdal python

I am having trouble creating a layer with geom type as geometry collection. Thus I am unable to output to a shapefile. I have attached code below. dstLayer is null. This does not happen if geom type ...
0
votes
0answers
47 views

Import GRD and count days with temperature below zero, using python

Minimum temperature of a year is given in GRD file. I have to read that file and then calculate the frost day (less than zero degrees celcius) from it in python.
0
votes
0answers
105 views

aggregate counts on a grid and determine grid adjacency

So here is the problem I am working on. I have a bunch of lat, long coordinates for event occurrences. I want to create a grid and then count how many points fall into each grid pixel. Once I have ...
0
votes
0answers
104 views

Convert Shp to KML error org2org Error4

I am attempting to convert Shp file to Kml using command line org2org application. I have used this in the past to convert more simple shp files and was sucessful. However I get an error message when ...
0
votes
0answers
83 views

Problems rendering tilestache served mbtiles in QGIS

I came across Geoffery Hing's walk through on using Tilestache to serve up mbtiles to QGIS. I followed the configuration and I was able to view tiles in the browser without any problems. However, when ...
0
votes
0answers
129 views

Make a raster with pixel scale 100 by 100 meter with gdal_rasterize

I'm asking here, because, I've got great results last time. Here it is the next of my adventure with gdal_rasterize, I want to make a raster picture where each pixel correspond to a square of 100 by ...
0
votes
0answers
83 views

How to rasterize a point layer? Cell value should be 1 if feature attribute is 1

I’m a fairly experienced GIS user, but I’m new to QGIS. I’m using version 2.2, and I’m trying to convert a point data layer into a raster data layer using the Rasterize tool. The point data layer has ...
0
votes
0answers
53 views

Python Error in recognising ENVI file

I am trying to read a .bip AVIRIS file using GDAL in Python. here is the code: from osgeo import gdal from osgeo import ogr from osgeo import osr from osgeo import gdal_array from osgeo.gdalconst ...
0
votes
0answers
107 views

No data appears when using gdal_calc.py but source had data values

I'm in the process of calculating topographic exposure using a mix of GDAL, grass, R and python in Windows. I need to process a very large tiff file using the following calculation: ...
0
votes
0answers
75 views

Fantasy map, how do i slice large png and georefrence it the right way?

Im and IOS developer working on a map based ios game, and new to GIS stuff. I've read a lot of info and i started to get a hint about what i need to do in order to meet my goals. I have a png map i ...
0
votes
0answers
96 views

Adding vector layers on top of gdal rasterized image

I used gdal2tile.py on a big image (actually its screenshots of image map) and now I have a nice openlayers-powered map I can pan and zoom around. But now I want to add a vector layer on top so I can ...
0
votes
0answers
76 views

How to use ogr2ogr to remove features from a shapefile using sql statements (Unable to open datasource `action/North' with the following drivers.)

OK, I've spent the obligatory 15-25 minutes on this and haven't gotten anywhere so, here's the command I'm trying to get to work right now: ogr2ogr -f "ESRI Shapefile" -where "area > 4000" ...
0
votes
0answers
18 views

save a WMS image on disk and apply a classitem (gdal commands?)

I need to save an image on disk trough WMS coloring the pixels applying a CLASSITEM (e.g. Mapserver CLASSITEM). Is it possible using gdal, if so, how? Thanks in advance, D
0
votes
0answers
221 views

ECW files to geoserver 2.4

I'm trying to add .ecw files (20 orthophotos) to geoserver 2.4. Geoserver is running on an ubuntu 12.4. Vserver. I'v already installed the GDAL plugin. Is this the right format for the ...
0
votes
0answers
13 views

UnIntersected part of gdal-geometry?

i have a rectangle and polygon geometries . polygon.Instersection(rectangle) gives me the polygon as it is inside the rectangle. But , what i want is the unintersected(non-intersecting) part of both ...
0
votes
0answers
105 views

command line alternates for softwares like tilemill, maptiler

tilemill works just fine for my needs but , i want to do the whole process on command line . since new geotiff files are going to be added in my project, and i need to automate this whole process . ...
0
votes
0answers
73 views

gdal2tiles colorpallete error

merged few files using gdal_merge gdal_merge.py -pct -o 713.tif sa713.tif na713.tif wa713.tif ea713.tif When i try to make tiles it says gdal2tiles.py: error: Please convert this file to RGB/RGBA ...
0
votes
0answers
96 views

GEOS assertion failed

The following script fails, if I import the 'rasterstats' library. import rasterstats import ogr geom1 = ogr.CreateGeometryFromWkt('POLYGON ((415924.806155724043492 ...
0
votes
0answers
45 views

How do i install liblas using OSGeoW4

I have downloaded and installed the OSGeo4W program on my system. I recently downloaded gdal file and liblas. When i try opening the osgeow4 shell, it keeps telling me the gdal file is missing. when i ...
0
votes
0answers
57 views

GDAL - least cost path analysis height data

I would like to use Pythons GDAL bindings to find the least cost path through an array and then show the route on top of the original array. I am using DEM data and would like the slope to be taken ...
0
votes
0answers
31 views

Installing GarminCustomMap: “missing _gdal_array module” error

I've used the GarminCustomMap plugin in QGIS 1.8 without troubles. But now, with QGIS 2.0, when I try to install it, I get the following message "missing _gdal_array module". Any idea how to fix this? ...
0
votes
0answers
68 views

Update a NetCDF file using GDAL

I have to edit a NetCDF file using GDAL in C++ code. When I try to update a subset I have this error message: ERROR 6: The NETCDF driver does not support update access to existing datasets. How can ...
0
votes
0answers
236 views

Add labels for printed map (GeoTools)

I have some problems with labels for printing map with GeoTools. Next code (with some cuts, code on Scala language) can help me for printing map: val map = new MapContent() val parser = new ...
0
votes
0answers
46 views

How to assign no value to pixels with (255,255,255)

I have Imagery where I want to replace white edges (255,255,255) with no value. I have used gdal translate within QGIS and put in no value 255 but this give no value to other pixels eg 255,240,200. ...
0
votes
0answers
125 views

Choose what zoomlevel to read from for mbtiles file in GDAL

GDAL is supposed to read MBTiles files and to do so, it selects the highest zoomlevel available. However, is there a way to choose another zoomlevel? i.e. override the default and go for a lower one? ...
0
votes
0answers
88 views

gdalbuildvrt ------— Problems with binary rasters used by grass

r.in.gdal is used for importing tif images in grass environment and they were changed to binary by default. I am trying to make a stack from them but unfortunately grass does not see my rasters ...
0
votes
0answers
66 views

How do I assign parameters for gdal proximity in QGIS?

I used to use the Proximity (Raster distance) option from the QGIS Raster menu a lot, which evidently corresponds to the commandline command, gdal_proximity.py. I'd like to be able to now use this ...
0
votes
0answers
255 views

gdal python set projection of a raster - does not work

Gdal 1.9.2, ECW driver: ECW/ERDAS Compressed Wavelets (SDK 3.x) I am having problems with the script that supposed to set projection of the ECW to EPSG:27700, i.e. the British_National_Grid. Please ...
0
votes
0answers
94 views

Strange black pixel on edge when generating with gdaldem -hillshade

I'm getting a strange result when I'm generating an hillshade with the gdaldem hillshade command. Here is my command (see the compute_edges parameter): ...
0
votes
0answers
328 views

What is meaning of GDAL on Ubuntu install [ERROR] Cannot make qgis.db private copy?

I am so lost and confused. I have spent days on this and not had any luck, so I'm really hoping someone can help! I am a new opensource user with an Ubuntu 13.04 OS. I used the software GUI (Ubuntu ...
0
votes
0answers
163 views

why readOGR() is not reading projection information from .spj - r

0 down vote favorite I successfully import a polygon shapefile with readOGR, however the information in the .prj file is not being imported. test class : SpatialPolygonsDataFrame nfeatures : 19407 ...
0
votes
0answers
231 views

How to read vector datasets using GDAL library

Up to now, I have used GDAL in order to read raster information and geotiffs and show them with the help of OpenGL. Now I have a vector dataset, and by that, I mean I have the files with types: DBF ...
0
votes
0answers
306 views

GDAL/OGR to extract pixel values

I want to use GDAL/OGR to extract pixel values from a raster file. The extraction should be based on polygons, stored in a shapefile. My thought was to iterate over the features in the shapefile, ...
0
votes
0answers
63 views

Why does ERDAS not like my GDAL output?

I am creating a raster using GDAL and the IMAGINE file doesn't work in ERDAS because it doesn't recognize the projection. Here is the code to generate a sample: from osgeo import gdal projectionWkt ...
0
votes
0answers
148 views

How to check are there any OGRGeometry (polygons) within an OGRGeometry (polygon) ?

Please see my image: I used the below statements to get Geometry: OGRGeometry *poGeometry; poGeometry = poFeature->GetGeometryRef(); Do you know how to check whether the Polygon 1 ...
0
votes
0answers
51 views

How to retrieve PROJECTION value for a map?

I am trying to read the Projection from Geotiff/shape file, using gdal 1.9.0 by the following command in the command line... gdalinfo.exe -proj4 map-path but its retuning empty string for ...
0
votes
0answers
387 views

Reprojection problem with gdalwarp

I have a GeoTiff file with WGS84 datum projection and I'm trying to reproject it using gdalwarp. Using this command gdalwarp -t_srs '+proj=tpers +azi=19 +tilt=30 +lon_0=-109.5 +lat_0=23 +h=1274199.4 ...
0
votes
0answers
611 views

How do I add elevation data to a GeoJSON line from a GPX file?

I understand that elevation can be added as an optional third value in Coordinates (lat, long, elev) in a geojson file, but I can't find a way to actually do that from a gpx file. I've used ogr2ogr to ...
0
votes
0answers
148 views

QGIS GDAL polygonize - how to set scale

I've got some old contour plots, paper copy only, that I'd like to be able to put into QGIS. I get as far as a successful polygonize into shp and it looks fine, but the scale is wrong. How might I ...
0
votes
0answers
48 views

How to select a subwindow from a source image in EPSG 3857?

I need to select a subwindow from an image projected to epsg 3857. If I use, for example: gdal_translate -projwin ${UL_x_3857} ${UL_y_3857} ${LR_x_3857} ${LR_y_3857} source_3857.tif ...
0
votes
0answers
157 views

SAGA grid VRT, file how to apply a LUT?

I've created a heap of saga grid files for some hydro survey work. I'm looking at a way to read these into a mapguide raster/grid layer, I'm able to load up the grids using a vrt file, I can't figure ...