[GIS] Using gdal merge correctly

gdallandsatqgis

Currently trying to use gdal_merge within qgis to stack three raster bands, the command it's showing is like this:

gdal_merge.bat -separate -of GTiff -o out.tif "LC80190302014126LGN00_B4.TIF" "LC80190302014126LGN00_B3.TIF" "LC80190302014126LGN00_B2.TIF"

however it keeps giving me ERROR 1: TIFFOpen:out.tif: Permission denied

I realize this should be a fairly simple procedure, is there something obvious I'm missing here?

Best Answer

"TIFFOpen:out.tif: Permission denied" message means that you do not have enough rights to write a new file into the target directory. Edit the -o parameter to point to some directory where you have write access.

Related Solutions

[GIS] Using GDAL in python to stack Landsat bands

Although it would require another library (which currently is only available on OS X and Linux) you could use RSGISLib (http://rsgislib.org), which is built on top of GDAL to do this. There is a function to stack bands, as an example:

#!/usr/bin/env python
import rsgislib
from rsgislib import imageutils
# Create list of images
imageList = ['LC80430342013291LGN00_B1.TIF',
             'LC80430342013291LGN00_B2.TIF',
             'LC80430342013291LGN00_B3.TIF',
             'LC80430342013291LGN00_B4.TIF',
             'LC80430342013291LGN00_B5.TIF',
             'LC80430342013291LGN00_B6.TIF',
             'LC80430342013291LGN00_B7.TIF']
# Create list of band names           
bandNamesList = ['Coastal','Blue','Green','Red','NIR','SWIR1','SWIR2']
# Set output image
outputImage = 'LC80430342013291LGN00_stack.tif'
# Set format and type
gdalFormat = 'GTiff'
dataType = rsgislib.TYPE_16UINT
# Stack
imageutils.stackImageBands(imageList, bandNamesList, outputImage, None, 0, gdalFormat, dataType)

It also allows passing in a list of band names.

However, if you were going to install RSGISLib to stack the Landsat bands you could just install ARCSI, which would stack the bands and convert to radiance or reflectance (if you were looking to do this). More details are here: http://spectraldifferences.wordpress.com/2014/05/27/arcsi/

Best Answer

Related Solutions

[GIS] Using GDAL in python to stack Landsat bands

Related Question