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
17 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 ...
3
votes
4answers
219 views

List of GDAL raster file extensions

I am trying to compile a list of file extensions that are accepted by GDAL, something like ['.tif', '.png', '.jpg', '.jpeg', ...] I found a list of GDAL Raster file formats, but the values in the ...
1
vote
1answer
46 views

GDAL “Simple Tiling” (Or yet another - how do I tile question)

Before I start, if anyone considers this a duplicate, then that's fair enough, I'll not object, I do happen to think however that it's a question that I don't see a clear answer for either here or ...
2
votes
1answer
66 views

Using Python to convert line shapefile to raster, value=total length of lines within cell

I have already read the similar question Convert line shapefile to raster, value=total length of lines within cell using R. And I want to realize the function using Python. ===== My attempt ==== ...
2
votes
1answer
42 views

How to create multiple external overviews for a geotiff with gdal?

I want to have mutiple overview files (one for each overview level) for a GeoTiff, so I'm able to edit this files in a raster editor. How can this be done? Unfortunately, gdaladdo -ro creates only ...
1
vote
1answer
34 views

Bad corner coordinates transforming SRTM *.hgt files to GRD format using gdal_translate

I'm transforming some SRTM *hgt files (1 second resolution version) to GRD format using dal_translate. For example, one of my orders is: gdal_translate -of AAIGrid N43W001.hgt N43W001.grd The ...
0
votes
0answers
63 views

How to access ArcSDE metadata using OGR?

I'm trying to get a list of all layers in an ArcSDE database hosted in oracle here is the command I'm using: ogrinfo -al SDE:localhost,5151,sde,sde,sde Note: I've replace the database, username ...
3
votes
1answer
61 views

Geo-Referencing an image from google maps with gdal_translate

I am trying to geo-reference an image using gdal_translate in a command prompt. Currently I have the following gdal_translate.exe -of GTiff -a_nodata 0 --config GDAL_DATA "Config location" -a_ullr ...
1
vote
0answers
9 views

Georeferencing image from google maps with gdal_translate? [duplicate]

I am trying to georeference an image using gdal_translate in a command prompt. Currently I have the following gdal_translate.exe -of GTiff -a_nodata 0 --config GDAL_DATA "Config location" -a_ullr ...
1
vote
1answer
135 views

Transform an image with projection EPSG:4326 to EPSG:3857

I have an image with known boundary and projection EPSG:4326 and want to overlay it to openstreetmap using Leaflet.js. I have tried using function imageoverlay, but the image is not aligned with the ...
1
vote
2answers
66 views

How to access Shapefile metadata and metatypes information using OGR?

I'm trying to use ogrinfo to get some details on a shapefile. Now, the only way that I know is using orginfo: ogrinfo -ro -so -al BRA_001.shp But its display all information about the shape: ...
1
vote
0answers
34 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 ...
3
votes
0answers
48 views

GDAL. Interpolation with breaks

I have a csv file dem.csv with data: LON, LAT, SO2(concentration) with NODATA=NAN. I need to interpolate my data from *.csv with breaks. I used gdal_grid with the command: gdal_grid -a ...
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
36 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
1answer
77 views

Gdal ogrinfo CSV unable to open datasource

When I try to import a csv file into a PostgreSQL table with the Gdal ogr2ogr tool I get the "Unable to open datasource" Failure. So I tried to open it with the tool ogrinfo and I get the same ...
1
vote
2answers
34 views

get first levels of Raster Data

I have a huge raster tif file 40GB, when I use this file in ArcMap it is very slow, I want to use this file only for the first levels, so I don't care about resolution. I don't want a compression I ...
2
votes
0answers
96 views

display of raster in QGIS - cells are not rotated

when viewing rotated rasters in QGIS - the overall raster is rotated (as designed), though each grid cell is not. (the nodata value doesn't appear to pick up as the default in QGIS either - ...
0
votes
2answers
207 views

extract scientific layers from MODIS HDF Dataeset using Python GDAL

I have been trying to extract scientific data layers from MODIS HDF files but without success. Before I have been doing this type of task using Arcpy, but I want the script that I am currently working ...
0
votes
0answers
30 views

How to generate raster “GeoTIFF” for given lat long bounds at particular zoom using GDAL

How to convert Lat Long bounds to geotransform values for a particular zoom in EPSG:4326 I can generated a raster with width, height and the geotransform values. But for an Lat Long bounds at a ...
1
vote
0answers
46 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
2answers
326 views

Install GDAL 2.0 on Windows

Ho do I go about getting GDAL 2.0 on Windows? I can get the source from the OSGEO site http://www.osgeo.org/node/1591 but I have no idea how to compile it. Is there not a built version of GDAL 2.0 ...
1
vote
1answer
52 views

Select points situated on zero degree slope

I need to choose points from SHP file, which are positioned on raster cells with 0 degree slope. Creating a slope raster from DME is not a problem, but I don’t know how to write a script that chooses ...
0
votes
0answers
28 views

GDAL/OGR in C++: Opening a vector file

So I'm making a vector file this way: OGRSFDriverH hDriver=OGRGetDriverByName("KML"); OGRDataSourceH hDS=NULL; hDS = OGR_Dr_CreateDataSource (hDriver, fname_dst, NULL); //I add layers ...
1
vote
1answer
77 views

GDAL/OGR in C++: Merge multiple vector datasets to single [closed]

I need to merge multiple kml files into a single one. I found that I could use "ogr2ogr -append" but I have to do it through c++ and my Google-fu hasn't been successful so far. So, any ideas? Also, ...
2
votes
1answer
116 views

Conditional statement in gdal_calc with result as float, not boolean

Many examples on gis.se that refer to gdal_calc give a boolean example. example 1 example 2 If I need to calculate a conditional statement on a float: If Raster_Value >= 1024, then 1024; else ...
0
votes
1answer
117 views

Converting multi-geometry GML with ogr2ogr not working

I am having difficulties with GML loading in PostGIS with ogr2ogr. It is a GML of buildingparts containing both polygons and (label)-points. In this example (below) I have a selection of 9 objects. ...
3
votes
1answer
99 views

How to ‘invert’ or ‘reverse’ a vector cutline when gdalwarp'ing a raster?

I have a raster of North America and I'd like to remove the larger lakes from it, replacing them with the same no-data value as the original raster. So I downloaded Natural Earth's 10 meter lakes ...
1
vote
1answer
13 views

Create single map image from two images with different color tables

Is it possible to merge the color tables of two images when using gdal to create a single image? Or is there another way to merge said images not using a color table?
1
vote
1answer
59 views

Getting elevation readings at any lat,lon from shapefile with scattered irregular points

I have a set of x,y coordinates, each of which I'm trying to get an elevation reading for. The elevation data source, however, is not a regular grid, but consists of a shapefile containing scattered, ...
1
vote
0answers
112 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
41 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 ...
0
votes
0answers
46 views

I want to make a small TIFF file from a big TIFF file

I just study on GIS system. When I create the Tiff file, there is no any data in Tiff file I want to make 4 channel TIFF file My code is below: void crop(const char *inputPath, const char ...
1
vote
0answers
59 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 ...
0
votes
2answers
47 views

Skip a specified input layer in ogr2ogr

I am using ogr2ogr to load data from multiple GML files into PostGIS. As shown in the ogrinfo output below, each GML file has multiple layers. Had to open data source read-only. INFO: Open of ...
3
votes
2answers
165 views

Can GDAL extract a raster window from disk without loading the whole file to memory?

I have a 500MB GeoTIFF with 3 bands on disk. I can easily read a very small, 25x25-pixel subwindow using GDAL like so: gdal_translate -srcwin 600 550 25 25 hugeimage.tif smallwindow.tif Does this ...
1
vote
2answers
143 views

Determining Min and Max elevation in Python using GDAL? [closed]

I have a DEM file and shape file of specific area ,and i want to Create a module in Python that writes name of each area which is shows in the picture, with the greatest differences in altitude on ...
4
votes
1answer
110 views

Optimizing Python GDAL ReadAsArray

I am using the GDAL ReadAsArray method to work with raster data using numpy (specifically reclassification). As my rasters are large, I process the arrays in blocks, iterating though each block and ...
1
vote
1answer
229 views

Python Gdal ERROR 6: GEOS support not enabled

I installed GDAL from here (https://pypi.python.org/pypi/GDAL/1.9.1) Before attempting the python GDAL I installed the prerequisite gdal libraries using Kyngchaos's description some time ago ...
1
vote
2answers
66 views

gdal_translate crashes on huge jp2

I have 44GB jp2 file, I want to translate it to tif file, but once it reaches 10% it crashes. No messages, just stopped. This is my command: gdal_translate.exe -of GTiff --config GDAL_CACHMAX 8000 ...
2
votes
1answer
36 views

Libraries for creating surfaces output in geotiff

I am looking for some libraries that can assist me in creating geotiffs from 3d polygons that becomes surfaces. Given a set of polygons where each node is [x,y,z] coordinates, i would like to create ...
6
votes
1answer
154 views

Using a geo-referenced image as an overlay to display markers on floorplan

I've created a tool based on the Google overlay tiler. This allows our users to upload a floor plan, adjust the position etc. using Google maps. Nothing fancy. We have a server-side job which creates ...
1
vote
0answers
77 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
1answer
110 views

GDAL Python OGRFeature .GetGeometryRef vs .Clone().GetGeometryRef

Just recently I was trying to get the geometry of a polygon stored in an ESRI Shapefile. In the script I was trying to read the polygon as follows: import ogr path = "shape.shp" ds = ogr.Open(path, ...
1
vote
0answers
42 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 ...
2
votes
0answers
141 views

NetCDF to Shapefile with Python and GDAL?

I need to save my outputs, which are netcdf files or panda python tables, to shapefiles. They cannot be rasters because the software I am loading them into accepts shapefiles only. I have googled ...
0
votes
1answer
121 views

How to open SMAP's Surface Moisture data files

I downloaded surface moisture data from Nasa's Soil Moisture Active Passive (SMAP) sensor website (https://nsidc.org/data/smap/smap-data.html). For each processing level (1-4), three types of files ...
0
votes
1answer
37 views

gdaldem hillshade using command to calculate hillshade of whole folder [closed]

I need to calculate hillshade of folder using gdaldem hillshade which has 60 files in it. How can i do it using command prompt window.
0
votes
0answers
29 views

USGS Raster Conversion with GDAL Returns Blank .TIF File

I am new to GIS systems. I installed GDAL to convert .ADF files to .TIF files to obtain image heightmaps. It outputs a GeoTIFF .TFW file that I can't open, and a .TIF file that is blank (all white). ...
1
vote
1answer
48 views

raster2pgsql size problem

I am about to import a raster (.tif with 50mb, size 3301x21676, pixel size 10x10) into postgres with raster2pgsql. With small extracts of the image it works but ... raster2pgsql [myfilename].tif -F ...