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
0answers
29 views

Confuse about how to make a tiles to support retina device

I have a map picture(15772*13796) so I break them to tiles using gdal2tiles(512*512px not 256*265px) for support retina device but is this correct? Because I a had read this issue and this. They ...
0
votes
0answers
104 views

Merge/mosaic categorical rasters in QGIS/GDAL

I am attempting to make yearly bins of some monthly MODIS data for all of the US. The cells of the raster are defined with Julian dates indicating the date at which a burn event took place. After ...
0
votes
0answers
56 views

Read esri shape point data and spline interpolation using python shell scripts?

I am quite new in python scripting. I am a dedicated user of arcgis but this time I have to work without arcgis. I have a point shapefile with 10 different attribute fields of different peat layer ...
0
votes
0answers
82 views

Read large raster into lat/lon value pair to csv file using python and/or gdal?

I have rasters for all of the states in US and they are in 10m resolution. So they can be very large, containing a large number of pixels if you think of TX or AK. I want to create a csv file with ...
0
votes
0answers
78 views

Polygonize only pixel values of 1 from a raster

I have a raster mask that I have created in QGIS (2.2.0) that consist of 0, 1 and nodata values. I would like to end up with a polygon dataset that consist of only the value of one pixels from this ...
0
votes
0answers
35 views

Gdal_Calc.py random pixel error

I'm using the gdal_calc.py in python in order to do the same thing as the Con tool in ArcGIS. I'm calling the .py using os.system: os.system('gdal_calc.py -A%s -B%s --outfile=%s ...
0
votes
0answers
91 views

Combine overlapping rasters in gdal_calc.py, taking the maximum value from each pixel?

I want to combine 3 overlapping rasters to create 1 raster. Each have the same dimensions, resolution and CRS. I want to populate the new raster with the maximum value from each of the layers. ...
0
votes
0answers
81 views

GDAL Installation help in windows

I'm trying to build a GIS application in Django Framework. I completed the GeoDjango installtion (postresql, spatial libraries). I installed OSgeo4W in my C drive. I also setup environment ...
0
votes
0answers
23 views

Opening Spot 5 DIMAP imagery in GDAL

Would like to know if someone is familiar of opening SPOT 5 DIMAP imagery using GDAL. Though the GDAL provides DIMAP driver, when I tried to open, it can not show the geotransform and geoprojection.
0
votes
0answers
150 views

I cannot use Gdal Python but Gdal binary, how to fix?

**UPDATE: I can use GDAL with Python now by copy file from Gdal 1.11 compiled folder to Python package folder. But the problem with GDAL binary and GDAL library that they cannot read the HDF5 file. ...
0
votes
0answers
27 views

Gamma correcting ECW

we would like to modify an existing ECW file with a gamma correction in order to make the image less colourful. Is there a way to do this using e.g. GDAL?
0
votes
0answers
76 views

What is GDAL equivalent to ArcGIS Lookup (Spatial Analyst) tool?

I use the Lookup (Spatial Analyst) tool to create several raster layers from a combined raster layer. I know ArcPy can do the same thing. However, I would like to know whether GDAL (or another ...
0
votes
0answers
77 views

gdal read whole geotif band into vector

i have another c++ GDAL question. I try to read alle elevation data out of an GeoTiff file into a c++ vector so that i can get elevation at x,y like this: vector[x][y]. I also like to use this data to ...
0
votes
0answers
41 views

How do I convert a VMAP to a shapefile using GDAL?

This is my first time trying to convert a VMAP to a shapefile and the examples online aren't as straightforward to me. Any guidance is appreciated.
0
votes
0answers
57 views

SST process using just GDAL and python

I am trying to process AVHRR 4km SST using a module that uses only GDAL and Python. OSGEO4W, GDAL, Python 3.4.1. I borrowed the download code from http://grasswiki.osgeo.org/wiki/AVHRR. The rest I ...
0
votes
0answers
47 views

Comparing datums with GDAL/OGR

I have a bit of code in Java that takes in a SpatialReference object and checks that it is one of the acceptable datums (WGS84, NAD83, NAD83 HARN, NAD27), but because of the lack of standardization in ...
0
votes
0answers
155 views

How to import SRTM data for the whole planet into TileMill

I need to get SRTM elevation data for the whole planet into TileMill. I already walked through the instructions on Steve Bennett's blog and it worked fine for small areas. The problem is now that for ...
0
votes
0answers
65 views

Snag loading Gdal 1.11 for Qgis 2.2, on a Mac

Can't load Q2.2 on a Mac without Gdal 1.11. but when I try to load 1.11 'Complete' version, it fails and gives me an error message to contact the developer. Anyone else had this kind of trouble? ...
0
votes
0answers
69 views

gdal_retile: to tile or pyramid only

I'm having trouble understanding when to use the -pyramidOnly option in gdal_retile. If I am trying to tile 10 GeoTiffs, am I right in assuming that if the images contain overlays and or inner tiles ...
0
votes
0answers
97 views

gdal commands cannot import gdal — possibly ArcGis conflict with separate Python?

I installed ArcGIS 10.1 yesterday (Windows 7 64-bit), but have a separate Python 2.7 installation. Now my gdal does not work correctly anymore. It looks like all commands without the ending "*.py" ...
0
votes
0answers
35 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
31 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
112 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
339 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
76 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
57 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
120 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
122 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
147 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
104 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
58 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
118 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
83 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
100 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
90 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
19 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
285 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
14 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
108 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
79 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
134 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
93 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
308 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
103 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
353 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
186 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
270 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
348 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
65 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
155 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 ...