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

Troubling using MODIS package in R (mac). Error with gdal

I'm new to using MODIS data, and have been trying to set up a a way to use the HDF files produced by MODIS on R. I had trouble with gdal (could only download the package through system dialog), and ...
1
vote
0answers
48 views

How do I calculate the percentage of progress using raster data using Python?

How do I calculate the progress amount of two or more raster data being processed? I have seen examples of vector data but it becomes difficult to calculate the progress when no vector data are ...
1
vote
0answers
49 views

ogr2ogr (or other) simplification of polygons with shared boundaries

I am trying to simplify a collection of 100,000 or so polygons containing nearly 4 million vertices. The data is stored in a 2GB shapefile so there is no explicit representation of the topology ...
1
vote
0answers
22 views

How can I read the TIFFTAG_DATETIME using GDAL in Python

I am currently using GDAL from Python to read and transform coordinates and that works well. I would now like to be able to read some of the none geographic tags, for example, TIFFTAG_DATETIME using ...
1
vote
0answers
62 views

gdal_merge.py is not working with NODATA values

I am trying to merge a couple of GeoTiff files with gdal_merge.py. Example input: In the input data the nodata value is -32678. Hence I tried: gdal_merge.py -n -32678 -o DEM.tif -of GTiff na_DEM....
1
vote
0answers
54 views

Storing non-orthorectified images in GeoTIFF

I'm reading GeoTIFF and GDAL documentations and I'm trying to understand what type of transformations are supported. GeoTIFF documentation says you can use ModelTransformationTag to define a ...
1
vote
0answers
18 views

Local Maxima (or Pour Point) Pixel Value extract in Raster

I am doing watershed & pourpoint using python GDAL. I have calculated Fill, Flow-Direction & Flow-Accumulation using Python & GDAL. The next is to get Pour Point in certain window ...
1
vote
0answers
62 views

error using gdal_translate with -projwin option

I have to crop geotiff image, I'm using gdal assembly and I'm trying to use the gdal_translate method If I crop my image using the -srcwin option all is fine gdal_translate -srcwin 0 0 400 400 ...
1
vote
0answers
31 views

Flipped (Y-axis) Raster from point cloud

I have been working on a project where I need to create rasters from irregularly spaced point clouds. I am resampling my point clouds to a regular 0.25 cm grid using a nearest neighbor method ...
1
vote
0answers
43 views

How to correctly rasterize a polygon which crosses 180E longitude?

I want to rasterize a polygon wich crosses 180E lon with gdal_rasterize. The polygon is in the WGS84 projection. I used the following command: gdal_rasterize -burn 1 -of GTiff -ot Byte -ts 1500 ...
1
vote
0answers
29 views

How delete feature from layer using Gdal binding java

I am trying delete one feature from layer using GDAL in java https://trac.osgeo.org/gdal/wiki/GdalOgrInJava But nothing happens. ogr.RegisterAll(); Driver driver = ogr.GetDriverByName("ESRI ...
1
vote
0answers
25 views

gdal_retile — force full tiles?

I have a set of large geotifs (DEM + orthoimagery) that I need to generate image pyramids from. So far gdal_retile from osgeo4w is working very well for my use case with one exception: I need all my ...
1
vote
0answers
69 views

undefined reference to `TIFFDefaultStripSize@LIBTIFF_4.0'

I am trying for the first time to use gdal in Ubuntu 15.10 with gcc. I call gdal.h in my program. At first I got the error "undefined reference to `GDALOpen'", so I check online and I found that I ...
1
vote
0answers
29 views

Creating tiling scheme to get multiple extents before gridding elevation data

Looking to create 1 degree by 1 degree tiles with pixel overlap to get extents before I start gridding. So basically I would like to have a shapefile of the tiling scheme with extents before having ...
1
vote
0answers
42 views

Generating continuous aspect surfaces from DEM with QGIS

I'm trying to extract aspect surfaces from different DEM sources at different resolutions in QGIS. I've had success working with SRTM-derived grids (3-arcseconds) and am trying to perform the same ...
1
vote
0answers
22 views

Access MEM data from GDAL?

I don't have any answers (yet) from from previous question on GDAL memory raster + OpenLayers, so i'll try to ask an easier question :) Anyone knows how to access GDAL memory with command lines when ...
1
vote
0answers
120 views

Write stacks of GeoTIFF image to NetCDF by chunk size

I converted GeoTIFF image stacks to NetCDF file with writing data by Index order (layer by layer) (this was a main solution). The problem is, it writes slowly (for 25 images - 1Gb takes 90 seconds). ...
1
vote
0answers
40 views

Retaining resolution of MODIS data while reprojecting using gdalwarp

Can anyone help me to use gdal(gdalwarp specifically) for reprojection? I am using it to reproject MODIS data : I originally have an hdf file - MCD43A4.A2005001.h08v05.005.2007351192333.hdf which has ...
1
vote
0answers
19 views

Check if raster (part of) lies within a pixel - python/gdal

I have a list of the corner-coordinates of different grid (c1-c4). First i sorted these coordinates to ensure that c1 is the upper left, c2 is the upper right, c3 is the lower left and c4 is the lower ...
1
vote
0answers
36 views

getting started designing satellite projection raster

How do I get started designing a map with a satellite projection like this:http://bl.ocks.org/joyrexus/6619251. So far using d3 I've created a map of North Africa w/ topojson and tiles from Mapbox but ...
1
vote
0answers
37 views

ReadAsArray leaves Qgs in unstable state

I have a QGIS plugin that reads a raster to get the elevation at vertex points in a polyline. The issue is that even though a single call to ReadAsArray(...) appears to work, if I immediately exit ...
1
vote
0answers
39 views

CF convention and GDAL projection

I'm trying to create some NETCDF4 CF compliant files from some code that generates GDAL compatible files. I was wondering whether there is a simple translation from GDAL (or proj4) projection ...
1
vote
0answers
54 views

How to add Geo PDF header to a PDF rendered with mapnik?

I rendered a PDF with mapnik and I am now trying to add geospatial header information to the PDF. I tried the method add_geospatial_pdf_header() from mapnik.printing.PDFPrinter. I do get some Geo PDF ...
1
vote
0answers
153 views

Error in opening Raster dataset in qgis “Cannot open GDAL Dataset”

Error in opening Raster dataset in qgis (see image below). Is there any suggestion as to why this is happening? and how do i correct it?
1
vote
0answers
45 views

Changing pixel values of GDAL VRT

Is it possible to manipulate the pixel values of a GDAL VRT and, if so, how? The original data must remain unchanged. I have a dataset which might be a bit rotated so that there are no data values on ...
1
vote
0answers
74 views

Upper limit to number of GeoTiff tiles?

I'm using the GDAL SWIG python bindings to create a new GeoTiff raster for use in a blockwise exponential-decay function. kernel_size = 177764054 driver = gdal.GetDriverByName('GTiff') kernel_dataset ...
1
vote
0answers
94 views

gdalwarp UTM Zone 1 to WGS84 - fails to transform

I have a global dataset consisting of tiles in UTM projection that I want to mosaic and therefore need to reproject to a common projection. Can I solve my reprojection error with options to gdalwarp ...
1
vote
0answers
51 views

Handling warning about DX and DY in gdal_translate from .Gtiff to .ASC?

When I try to convert Gtiff to .asc file, QGIS shows the warning: Producing Golden Surfer style file with DX and DY instead of CELL SIZE since the input pixels are non-square. Use the FORCE CELL SIZE ...
1
vote
0answers
28 views

Raster statistics change after re-projection — which one is correct

I'm looking into some global DEMs (like ETOPO1 - http://www.ngdc.noaa.gov/mgg/global/global.html), and noticed that the statistics for the elevation values can change quite a bit after reprojecting ...
1
vote
0answers
32 views

Theory notions on exporting bufr data in geotiff format

I'm trying to export some bufr data from polar satellites/sensors like DMSP, ASCAT, ATMS etc in geotiff image. When my program rasterizes satellite data (brigthness temperature for example) by ...
1
vote
0answers
73 views

WMS layer throws “IReadBlock failed at X offset 0, Y offset 1” error. Possible XML issues

I am working with a few WMS layers to try and see which one would work best with an application I'm developing out of QGIS 2.8. I need the CRS to be in EPSG 4326 but, as many of you know, most webmaps ...
1
vote
0answers
22 views

Should I set the projection if the data has already been projected?

I need to plot some data which was produced on a native grid with a lambert conformal conic projection. I need to create a geotiff (using python and gdal) and was wondering how to go about and do ...
1
vote
0answers
176 views

how to set the geotransform using GDAL in python

I am trying to create an image from a numpy array and georefenrece it so that I can import it into qgis. I have the latitude and longitude in degrees for the top left corner of the array as well as ...
1
vote
0answers
37 views

Projected Coordinate System issue using GDAL and qgis when defining x and y

I have a numpy array which I would like to save as GeoTiff using GDAL in python. The array was created using the following PCS in mind: +proj=lcc +lat_1=48 +lat_2=48 +lat_0=48 +lon_0=0 +x_0=0 +...
1
vote
0answers
115 views

How to handle GDAL dataset error when accessing Arcgis Online world imagery data?

I would like to use the world imagery data set for loading base maps to my QGIS projects. Up until now, I have had no problem by typing the following command in the Python console: qgis.utils.iface....
1
vote
0answers
67 views

GDAL MSSQL with numeric table names

I have many SHP tables organised by year, for example 2015_xx_data.shp that I have loaded into a MS SQL Server database using OGR2OGR. I have also loaded them into MySQL for a comparison. There are no ...
1
vote
0answers
573 views

Error While installing GDAL in Ubuntu 12.04

I am trying to set CartoDB in Ubuntu 12.04. I am trying to install Python dependencies as a part of my installation. I am installing with below commands export CPLUS_INCLUDE_PATH=/usr/include/gdal ...
1
vote
0answers
90 views

GDAL - replacing rgb value (0,0,0) with (nodata,nodata,nodata)

I have a set of rasters in which I want to replace any rgb value of (0,0,0) with (nodata,nodata,nodata) I've tried using gdal_translate -a_nodata "0,0,0" in.tif out.tif But this seems to convert ...
1
vote
0answers
25 views

Transform stereographic projection using GDAL

How can I read/wrap a JPG image without any embedded projection/coordinate system to an autoreferenced GEOTIFF with Mercator or similar coordinate system using GDAL? The image is 500x500 px and I ...
1
vote
0answers
45 views

How to use ST_fromGDALRaster correctly with a given GeoTIFF

I am trying to import gdal raster data (a GeoTIFF) to a Postgis database, using the ST_FromGDALRaster method (http://postgis.net/docs/manual-2.1/RT_ST_FromGDALRaster.html). The question is that I don'...
1
vote
0answers
34 views

Changes in MOD09A1 quality mask as processed with GDAL

I'm analyzing MOD09A1 spectral data using GDAL/Python. The quality mask is a uint32 dataset and want to pull out the value of the first bit to mask the spectral data (1 == poor data; 0 == good data). ...
1
vote
0answers
948 views

Reading, modifying and writing a geotiff with GDAL in python

I'm trying the learn the ropes of Remote Sensing image processing using Python GDAL bindings and numpy. As a first attempt, I'm reading a Landsat8 geotiff file, do a simple manipulation and write the ...
1
vote
0answers
105 views

how to reproject an image programatically in c# using gdal

I have to reproject an image programmatically in c# using gdal library. For the moment i'm a little bit lost in gdal and i don't really now which methods i have to use. I think i will have to ...
1
vote
0answers
25 views

gdal: Pixelvalues to nodata

I have a GRD-File where all pixelvalues which are less than 0.05 should be converted to nodata. I know how to do this with ArcMap, but the process should be done automatically, for example with ...
1
vote
0answers
176 views

What if 'geometry_columns' already exists in MSSQL Server Database?

When trying to load a .SHP into MSSQLSpatial using OGR2OGR, I get the following error: ERROR 1: Error initializing the metadata tables : [Microsoft][ODBC SQL Server Dr iver][SQL Server]There is ...
1
vote
0answers
108 views

GDAL and WMTS supprt

I have researching possibility to fetch raster file from WMTS server via GDAL. I have found some old post about that GDAL doesn't support WMTS (WMS only). However, I also find this information: http://...
1
vote
0answers
49 views

gdal2tiles issues with color scale between tiles

Background: I have a series of about 40 small GeoTIFFs that need to be merged together to make one large image. Each GeoTIFF has 2000 bands. I've tried using both gdal_merge and gdal_retile to pull ...
1
vote
0answers
116 views

ogr2ogr conversion error - Can't create fields of type Binary on shapefile layers

I'm having trouble using the ogr2ogr tool on the command line. My data is stored in SQL as binary, type is mixed polygons and multipolygons. I am using the following command for conversion from SQL to ...
1
vote
0answers
79 views

gdallocationinfo used in javascript with node-gdal?

I've been scratching my head for a white trying to figure it out if it's possible to get pixel information from my rasters using gdallocationinfo in javascript. I found node-gdal, but after reading ...
1
vote
0answers
111 views

WKID issue when adding raster layer from ArcGIS “mapserver” to QGIS via REST

QGIS docs here describe how to do this. I am trying to add a layer from this MapServer: gps.digimap.gg/arcgis/rest/services/JerseyBaseMap/MapServer When I try the following command from the Python ...