Use arcpy to extract bands from ECW image

AnthonyCheesman1

Hi folks

This may be more of an arcpy question than an SA question, but I'll ask anyway.

I'm looking for an automated way using arcpy to extract bands from ECW images. The only way I've been able to do this so far is to do it manually by loading the individual band to an ArcMap session, and then right-click > Export Data. I haven't been able to work out how to do this with any ArcToolbox tools, or any arcpy functions.

Have you been able to do this with arcpy or ArcToolbox in the past? How did you do it?

Thanks

Anthony

MehdiPira1

Anthony Cheesman‌

Let's say your image contains 3 bands, you can use the following model:

I also wrote a python code for this:

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">import</SPAN> os
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">" "</SPAN>
rasterfile <SPAN class="operator token">=</SPAN> Input Raster
desc <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>rasterfile<SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">for</SPAN> band <SPAN class="keyword token">in</SPAN> range <SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">,</SPAN> desc<SPAN class="punctuation token">.</SPAN>bandCount<SPAN class="operator token">+</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
       outputRaster <SPAN class="operator token">=</SPAN> <SPAN class="string token">"rdlayer_"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>band<SPAN class="punctuation token">)</SPAN>
       <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>outputRaster<SPAN class="punctuation token">)</SPAN>
       arcpy<SPAN class="punctuation token">.</SPAN>management<SPAN class="punctuation token">.</SPAN>MakeRasterLayer<SPAN class="punctuation token">(</SPAN>rasterfile<SPAN class="punctuation token">,</SPAN> outputRaster<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">''</SPAN><SPAN class="punctuation token">,</SPAN> envelope<SPAN class="operator token">=</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> band_index<SPAN class="operator token">=</SPAN>str<SPAN class="punctuation token">(</SPAN>band<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Hope that helps.