topic How i can split multi-band raster to separate rasters ? in ArcGIS Spatial Analyst Questions

How i can split multi-band raster to separate rasters ?

DimitrisPoursanidis — Mon, 01 Sep 2014 09:06:17 GMT

Hello,

Stragling to solve this "easy" task, i have use more than 8 different softwares to solve it (it was solved by using Matlab). Any idea, suggestion on how can this be done via ArcGIS ?

Re: How i can split multi-band raster to separate rasters ?

NtriankosIoannis — Mon, 01 Sep 2014 12:36:22 GMT

hi,

take a look

Kind Regards,

John

Re: How i can split multi-band raster to separate rasters ?

DimitrisPoursanidis — Mon, 01 Sep 2014 12:52:23 GMT

Thanks but this needs lots of manual work when i have to work with MODIS multiband datasets.

I need just the opposite of the "Compose bands" aka "Decompose Bands".

I find it just only in the Orfeo Toolbox via Monteverdi.

Re: How i can split multi-band raster to separate rasters ?

RiyasDeen — Mon, 01 Sep 2014 12:56:54 GMT

Hi Dimitris,

you need Extract Subdataset.

ArcGIS Help (10.2, 10.2.1, and 10.2.2)

Re: How i can split multi-band raster to separate rasters ?

DimitrisPoursanidis — Mon, 01 Sep 2014 13:01:00 GMT

Subdataset file formats can be either Hierarchical Data Format (HDF) or National Imagery Transmission Format (NITF) files

So no lack with GEOTIFF.

No other format will work (and i test it).

Re: How i can split multi-band raster to separate rasters ?

FilipKrál — Mon, 01 Sep 2014 21:24:00 GMT

Hi Dimitris,

I remember I struggled with this too, but I cannot recall if there is any tool for that. The code below is roughly what I would try. I don't have ArcGIS at hand at the moment to test it, but hopefully any glitches are only minor.

The key is that you can refer to individual bands within a raster dataset by their correct path, just like you can drill down to them in the Catalog directory tree. Usually, bands are called Band_1, Band_2, and so on. The problem is that sometimes they have different names (Layer_1, etc.) The code below deals with this by looping through whatever children the raster has and getting the correct name.

Results (individual raster bands) are stored as individual rasters in the output directory.

import arcpy, os
in_raster = r'c:\path\to\raster.tif' # input multiband raster
out_folder = r'c:\path\to\output\folder' # preferably an empty folder

desc = arcpy.Describe(in_raster)

for band in desc.children:
    bandName = band.name
    band_path = os.path.join(in_raster, bandName)
    dest_path = os.path.join(out_folder, bandName + '.tif')
    arcpy.CopyRaster_management(band_path, dest_path, "", "", "", "NONE", "NONE", "") # change parameters here

I have a feeling there might be a more elegant solution, but if you cannot find it, give this a go. Make sure you change the parameters at row 2,3, and 11, then just paste it into Python window and hit enter. And then enter again.

The code is inspired by 41145 - Add individual raster bands from a multiband raster... and Extract one band from 3 band raster

Let us know if this works for you, I'd really want to know.

Filip.

Re: How i can split multi-band raster to separate rasters ?

DimitrisPoursanidis — Tue, 02 Sep 2014 09:04:46 GMT

And yes, it work Filip.

Now i will try to put it in a toolbox.

So simple functions and yet are not included in the default tools (not only in arc but also in envi, erdas). Only QGIS via SEXTANTE and link to Orfeo Toolbox have such capability.