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
1answer
532 views

Can't get gdal_merge Python script to run

I downloaded and installed GDAL. I am using Python 27. I had to use the following commands at the IDLE Python GUI: from osgeo import gdal,ogr sys.path.append('C:\Program Files\GDAL') import ...
1
vote
0answers
56 views

Rasterstats throwing error: Cannot open as GDAL raster

I'm using the rasterstats zonal_stats function (pip installed), and running it on a .vrt file. Everything runs perfectly in my Ipython notebook, but when I try to run it in my regular python ...
0
votes
0answers
42 views

gdal tool for extracting data

Is there any gdal tool available for obtaining data from number of files. actually i need atmospheric parameters which are downloaded from modis website. some datasets are not having data i.e ...
5
votes
0answers
161 views

can I preserve a fading alpha layer when mosaicing images with gdal?

I am mosaicing some images with gdal and would like to improve the final result by using a fading / gradual alpha layer towards the edge of each image to remove the sharp edges in the middle of the ...
0
votes
1answer
65 views

Join line segments sharing a common attribute, using GDAL

I'm having trouble rendering a shapefile consisting of trail routes broken down into many short segments. So I'd like to combine them into a small number of line features, which share a NAME and touch ...
1
vote
0answers
70 views

MS Visual 2013 with GDAL and HDF 4

I'm using MS Visual Studio 2013, C++ and newest gdal downloaded from here http://www.gisinternals.com/query.html?content=filelist&file=release-1800-x64-dev.zip In readme and libraries I can see ...
1
vote
1answer
243 views

How to use gdal.FillNodata in python?

I am trying to use gdal.FillNodata in python on a ".tif" raster. The nodata value is set to -3.402823e+38. Here is my code : import gdal from gdalconst import * ET = ...
0
votes
1answer
189 views

Grid-type shapefile from raster GeoTIFF

I need to produce a shapefile from a raster. I need to be able to specify either an interval, or manually define the intervals. The input file is a GeoTIFF. I would like to use GDAL and/or Python. ...
1
vote
0answers
57 views

Extra white space in ASCII file due to Translate GDAL from QGIS Modeler

Could anyone enlighten me why do I get an extra white space in from of my first column in an ASCII file generated through a model? With ArcGIS it doesn't happen. Output raster type: Int32 Output ...
2
votes
1answer
2k views

Interpolation grid for scattered data having (latitude, longitude) coordinates

I've got some scattered data in the form of (latitude, longitude, someParameterValue). I'm using inverse distance weighting interpolation method to interpolate them in a rectangular grid of pixels. ...
0
votes
1answer
49 views

projection for high accuracy small areas

I need to convert (lat,long,height)points into a cartesian system for an area as small as 40kmx40km or even 10kmx10km. I'm currently using the UTM projection, but there are some issues regarding ...
1
vote
1answer
1k views

Sample for MS Visual Studio C++ using GDAL library

I need to create a sample console application on c++ that will use *.hdf file as an input and will read it. As I got, I need a GDAL library. I've found a GDAL tutorial on official GDAL website, but ...
1
vote
1answer
900 views

How to place an numpy array into GeoTIFF image using python-gdal

I am attempting to export data of varying scales into GeoTIFF images...a new method for me. The problem I'm running into is when using the GDT_Byte pixel type, data values < 0 and > 256 are not ...
0
votes
0answers
95 views

GDAL Create Layer WFS

So I'm trying to write a program that has the ability to read and write data from a WFS (Web Feature Service) I've managed to pull all the data I want from the WFS. Now I'm trying to write to it and ...
1
vote
0answers
215 views

How to build gdal_calc.py from source in MinGW

I was wondering how to build the gdal python binding from source using MinGW. In particular, I need to install python utility scripts such as gdal_calc.py. I've searched around and found this ...
2
votes
1answer
1k views

Properly Converting from UTM to Lat/Long with gdalwarp

I have a file in GeoTiff format in UTM coordinates. I convert the file to LatLong, then back to UTM, and the coordinates aren't what I started with. I have the printout of the process below. ...
0
votes
0answers
74 views

How to overlay dem (100m) over MODIS (500m) in square grid?

I have SRTM Dem data and the resolution 90m approx at the equator however the resolution varies as we move from equator to poles. For my case the study area is Himalayas. My plan is to get the square ...
0
votes
1answer
95 views

Changing the origin of a Geotiff and resampling

I have georeferenced data from 1 to 10M resolution in GeoTiff format. I can use gdal to resample everything to 5M, but I would like to store the data so that rather than having an origin at a decimal ...
0
votes
0answers
87 views

Proper projection name for NAD83 2011 revision

I have two PRJ files, one from ESRI and another from some custom GDAL/OGR code I wrote (pasted below). I was trying to have both work in NAD83 UTM10N with the 2011 revision. The ESRI file has 2011 in ...
0
votes
0answers
123 views

Create GeoTIFF from JPG and JFW files with gdal

I have a JPG image and the corresponding geospatial information on a JFW wolrd file in a ETRS89 geodesical system and UTM coordinates (EPSG:25830) 5.0000000000 0.0000000000 0.0000000000 -5.0000000000 ...
0
votes
1answer
81 views

How to cut contour lines at the edge of a data source

I have calculated some contour-lines using gdal-contour for two elevation model into a Postgis database. One is high res and the other lower res. Now I am seeing problems at the edge of my high ...
0
votes
0answers
184 views

ogr2ogr failing to open/convert GML

I'm trying to convert a GML file into GeoJSON using ogr2ogr, I've compiled GDAL, Expat and Xerces from source, but when I run my command I get an error. My convert command ./ogr2ogr -f "GeoJSON" ...
1
vote
0answers
30 views

QGIS Georeferencer adds padding to image after save

I am attempting to georeference an aerial image using satellite imagery. I'm successful in getting an output, however I notice there is an issue where QGIS's georeferencer is adding padding to all the ...
6
votes
2answers
228 views

How to convert the geographical projection of an extremely large file (GDAL)

My Task I need to convert 9 x 14Gb GeoTiff files to WGS84 geo projection. The rough size in pixels is about 200k x 100k. The data-type is byte and it is not a subject to interpolation, it is some ...
0
votes
2answers
217 views

It's possible to use gdal_calc.py with an entire path of rasters?

I try to use gdal_calc.py to scale and transform a time-series of MODIS images. I've done that when i worked with a single image in the -A variable and a mask in -B, but i need to work with a ...
1
vote
1answer
68 views

Read raster at request size through Python

I'm trying to read the raster (tif) file using Python. What I want is to read the raster at a specific request size rather than its original resolution. Example: If a raster has a cellsize of 5 x ...
1
vote
1answer
54 views

Grass r.out.gdal, wrong values when importing to QGis

I have a problem with export/import functionalities between Grass and QGis. I create a raster file in Grass (with r.wind.sun, but that most likely does not matter) It has the following properties: ...
0
votes
0answers
114 views

Escaping the MODIS Swath Projection with GDAL: MOD06_L2: MRT Swath Tool not working

Stackexchange, I am looking at ftp://ladsftp.nascom.nasa.gov/allData/6/MOD06_L2/2012/222/MOD06_L2.A2012222.0800.006.2015062054108.hdf Trying to get the SDS (in WGS84 Projection and as GeoTiff): ...
0
votes
1answer
84 views

gdal_translate cannot convert MSG3 channels to TIFF

following this guide http://www.gdal.org/frmt_msg.html I want to convert my MSG3 HRIT format files in TIFF but I get this error gdal_translate H-00 ...
2
votes
0answers
233 views

Zonal computations: computing raster pixel counts within polygons using GDAL versus rasterstats

I'm computing some zonal statistics, and was comparing the results obtained using the GDAL cookbook example code here with the rasterstats zonal_stats function. I want to compute the mean value per ...
1
vote
0answers
133 views

Use Gdal from Python IDE

I am just starting to experiment with writing python scripts. I have been using QGIS and the Osgeo4w shell for sometime now. I would like to be able to call Gdal from the python IDE i have installed ...
4
votes
1answer
112 views

Does gdal_calc only support 26 input raster files at a time?

In gdal_calc.py, input files are specified with letters A-Z. Here are two files, input.tif and input2.tif, being averaged into result.tif: gdal_calc.py -A input.tif -B input2.tif ...
0
votes
1answer
33 views

Gdal/OGR LayerMetaData not parsing its value

To start off with: Using Gdal/Ogr compiled with Xerces to read a WFS Reading from a WFS version 1.0.0 Compiler is QT version 5.4.1 ( C++ ) So I'm trying to get all the capabilities of a WFS via ...
1
vote
1answer
181 views

Accessing bands in .pix (PCIDSK) file in python/arcpy

I am trying to perform zonal statistics on a single band from a multi-band raster in .pix format. I have things working great with a tiff file. I am able to access the bands using this code: desc = ...
1
vote
2answers
451 views

GDAL raster shifting with python error: 'NoneType' object has no attribute 'GetGeoTransform'

I'm trying to shift a georeferenced raster by 0.5m N and 0.5m W in QGIS using the python console: from osgeo import gdal # Open in read/write mode rast_src = gdal.Open('test_raster.tif', 1) # Get ...
0
votes
0answers
75 views

RuntimeError opening rasterio array in GDAL

I've opened a GeoTiff with rasterio: import rasterio with rasterio.open('Africa.tif') as src: transform = src.meta['transform'] array = src.read(1) array([[203, 203, 203, ..., 0, 0, ...
0
votes
1answer
310 views

Vector to Raster using GDAL (gdal.RasterizeLayer) : Error in Output Raster(generating output raster with NAN values)

For rasterizing a vector layer, I have tried gdal.RasterizeLayer method for converting vector(shapefile) to raster(tiff). But its giving raster output with NAN values(complete black image). I need ...
2
votes
2answers
314 views

How to make a movie from a ascii grid files?

Lets say I have several ascii grid files for each time step for a particular location. I would like to spit out the image for each of these ascii grid files and legend on the side bar. I also would ...
2
votes
2answers
187 views

ERROR 1: None of ECW_ENCODE_KEY and ECW_ENCODE_COMPANY were provided

I downloaded ERDAS-ECW_JPEG_2000_SDK-5.1.1 and tried to use it with GDAL to create ECW files from georeferenced PNG files, but they are not created, it basically complains about: ERROR 1: None of ...
0
votes
1answer
151 views

creating a pyramide tiles of a geotiff using gdal

I have a big geotiff file (about 300MB), and I saw an example of a big data organized in folders : folder 0 contains one folder 0 which contains one .png file named 0.png folder 1 contains two ...
1
vote
1answer
136 views

GDAL doesnot support HDF4 e [duplicate]

I am trying to import the pymodis library within python and I keep getting error that "GDAL installation has no support for HDF4, please update GDAL" I am running the latest version of GDAL ...
1
vote
2answers
374 views

How to define PROJ.4 string for GDAL projection

I'm doing some zonal computations in GDAL using the code posted in the solution here: Issue Trying to create Zonal Statistics using Gdal and Python. It runs fine for me using the example data posted ...
2
votes
0answers
119 views

GDAL equivalent to ArcGIS Union?

I have a single shapefile full of overlapping polygons. The ArcGIS Union tool will generate individual, non-overlapping polygons with that geometry. Is there a way to do this with GDAL/OGR or another ...
1
vote
1answer
312 views

where to download gdal_polygonize.py?

I am trying to help out a friend of mine. She learnt GIS in school in the past using ArcGIS but has recently installed QGIS package on her Mac OSX. When she tries to render raster into shape file, ...
0
votes
1answer
226 views

Can't get the correct version of GDAL installed to use ReadAsArray without getting error in Python

Using the instructions here . . . https://pythongisandstuff.wordpress.com/2011/07/07/installing-gdal-and-ogr-for-python-on-windows/ I downloaded the following files: GDAL-1.11.1.win-amd64-py2.7.msi ...
1
vote
1answer
78 views

Getting Feature names from WFS in C++

to start off with: Coding in C++ Using WFS Version 1.0.0 Using QT 5.4.1 Windows 8 Latest versions I'm trying to read in a WFS in code via Gdal compiled with Xerces - I'm able to read the GML ...
1
vote
1answer
56 views

How to change band type?

I have an hillshade tiff for which the current band type=Bytemay be problematic. How to change the band type to type=float ?
3
votes
2answers
161 views

Raster Caculator Qgis 2.6 v/s Gdal_calc.py by Terminal

I need calculated monthly mean of Snowcover data by MODIS product, and then modify this script gdal_calc.py in this section # create alphabetic list for storing input layers ...
3
votes
1answer
74 views

GDAL can i try to create a workflow?

When i try to compose a workflow im confused to the funcionality of that. I want to create a workflow to execute commands gdal. what syntax I can use if I wanted to create a .bat file? For example ...
0
votes
1answer
73 views

How to update an existing GeoJSON with Fiona

With the code below I can create a new GeoJSON and add data to it: data =[(LINESTRING (-98.98234 19.38234, -98.990341 19.342345, -98.992449 19.378621, -98.991299 19.378233, -98.992349999 19.377381), ...