Branched from: Help: Batch - composite bands Arc-python Hi, I would like use this script, but it fails. Some hint? import arcpy , os arcpy . env . workspace = r 'C:\Teste_Auto_CLEIA\Ima_compactadas' # list all folders in a directory folders = arcpy . ListWorkspaces ( "L8*" , "Folder" ) # Inside main directory, already have other folders, but the bands are always inside folders that begin with "L8", e.g. "L8227069". #The output files place outws = r 'C:\Teste_Auto_CLEIA\Imagem' # Line 9 for folder in folders : arcpy . env . workspace = folder rasters = arcpy . ListRasters ( "*.tif" ) # I already erased the bands which I won't use chk4 = [ 'band{0}.tif' . format ( i ) for i in range ( 0 , 7 ) ] # I changed the range rasters = [ a for a in rasters if a [ - 9 : ] in chk4 ] name = os . path . join ( outws , rasters [ 1 ] . split ( "_" ) [ 0 ] + ".tif" ) # I wrote "outws" arcpy . CompositeBands_management ( rasters , name ) # I put it print ( rasters ) I'm using ArcMap 10.1 and the bands names aren't the same, for instance, "LO82270682017183CUB00_B1", ... , "LO82270682017183CUB00_B7". Thanks!

nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Batch composite bands using arcpy

BrunoDeus1

Branched from: Help: Batch - composite bands Arc-python

Hi,

I would like use this script, but it fails. Some hint?

<SPAN class="keyword token">import</SPAN> arcpy<SPAN class="punctuation token">,</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">'C:\Teste_Auto_CLEIA\Ima_compactadas'</SPAN>

<SPAN class="comment token"># list all folders in a directory</SPAN>
folders <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListWorkspaces<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"L8*"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"Folder"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># Inside main directory, already have other folders, but the bands are always inside folders that begin with "L8", e.g. "L8227069".</SPAN>

<SPAN class="comment token">#The output files place</SPAN>
outws <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\Teste_Auto_CLEIA\Imagem'</SPAN> <SPAN class="comment token"># Line 9</SPAN>

<SPAN class="keyword token">for</SPAN> folder <SPAN class="keyword token">in</SPAN> folders<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> folder
&nbsp;&nbsp;&nbsp; rasters <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*.tif"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># I already erased the bands which I won't use</SPAN>
&nbsp;&nbsp;&nbsp; chk4 <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">'band{0}.tif'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>i<SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">for</SPAN> i <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">7</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="comment token"># I changed the range</SPAN>
&nbsp;&nbsp;&nbsp; rasters <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>a <SPAN class="keyword token">for</SPAN> a <SPAN class="keyword token">in</SPAN> rasters <SPAN class="keyword token">if</SPAN> a<SPAN class="punctuation token">[</SPAN><SPAN class="operator token">-</SPAN><SPAN class="number token">9</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="keyword token">in</SPAN> chk4<SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; name <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outws<SPAN class="punctuation token">,</SPAN> rasters<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>split<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"_"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># I wrote "outws"</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CompositeBands_management<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">,</SPAN> name<SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># I put it</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>rasters<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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

I'm using ArcMap 10.1 and the bands names aren't the same, for instance, "LO82270682017183CUB00_B1", ... , "LO82270682017183CUB00_B7".

Thanks!

Trouver plus de publications taguées avec

Accepted answers

XanderBakker

Can you try this?

<SPAN class="keyword token">import</SPAN> arcpy<SPAN class="punctuation token">,</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">'C:\Teste_Auto_CLEIA\Ima_compactadas'</SPAN>

<SPAN class="comment token"># list all folders in a directory</SPAN>
folders <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListWorkspaces<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"L8*"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"Folder"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># Inside main directory, already have other folders, but the bands are always inside folders that begin with "L8", e.g. "L8227069".</SPAN>
<SPAN class="keyword token">print</SPAN> <SPAN class="string token">"folders: "</SPAN><SPAN class="operator token">+</SPAN>str<SPAN class="punctuation token">(</SPAN>folders<SPAN class="punctuation token">)</SPAN> <SPAN class="comment token">#To see if the search above finds a result</SPAN>

<SPAN class="comment token">#The output files place</SPAN>
outws <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\Teste_Auto_CLEIA\Imagem'</SPAN> <SPAN class="comment token"># Line 10</SPAN>

<SPAN class="keyword token">for</SPAN> folder <SPAN class="keyword token">in</SPAN> folders<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> folder
&nbsp;&nbsp;&nbsp; images <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"L*2*0*_B*.tif"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token">#I changed the user word for an exclusive. This line keeps the bands that I will use</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"images: "</SPAN><SPAN class="operator token">+</SPAN>str<SPAN class="punctuation token">(</SPAN>images<SPAN class="punctuation token">)</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; chk4 = ['band{0}.tif'.format(i) for i in range(0, 7)] # I changed the range</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; print "chk4: "+str(chk4)</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; rasters = [a for a in images if a[-9:] in chk4]</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; print "rasters: "+str(rasters)</SPAN>
&nbsp;&nbsp;&nbsp; name <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outws<SPAN class="punctuation token">,</SPAN> images<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>split<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"_"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># I wrote "outws"</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CompositeBands_management<SPAN class="punctuation token">(</SPAN>images<SPAN class="punctuation token">,</SPAN> name<SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># I put it</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Finished"</SPAN> <SPAN class="comment token">#23‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍</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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Since you already filtered the folder and removed the images that you don't want to use, there is no need to filter them using lines 16 to 19. The error you were receiving was caused due to trying to extract an item with index 1 from an empty list. The list was empty since you were trying to filter the list using a name "band*" when the bands have a different name (there will be no match).

All comments

NihalAlbuquerque

This is the code I used for my Landsat work, put all the Landsat folders into a working folder. You can alter the code to suit your purposes.

import arcpy
import glob
import os

#working drirectory with landsat data in individual folder
os.chdir('E:\GIS\landsat2019')
#list folders
folders = glob.glob('*/')
print(folders) 

#below code, makes a composite file for landsat bands in each folder    
for folder in folders:
    #list file name for each band
    filelist = glob.glob(folder +'/L*_B*.tif')
    print(filelist)
    #makes the folder the working directory
    arcpy.env.workspace = folder
    #creates file name based on folder name
    name = str(folder)
    names = str(name[0:-1]+'.tif')
    print (names)
    ##Compose multi types of single band raster datasets to a TIFF format raster and dataset
    arcpy.CompositeBands_management(filelist,names)
    print 'complete'
    #puts the file with the folder name, in each landsat folder

BrunoDeus

It works fine!

And the first (and original) script didn't worked to "n" scene(s) just because one line

<SPAN class="keyword token">import</SPAN> arcpy<SPAN class="punctuation token">,</SPAN> os
<SPAN class="comment token">#Define the workspace that contains all of the folders with the Landsat imagery</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\Teste_Auto_CLEIA\Ima_compactadas'</SPAN>

<SPAN class="comment token">#Specify where you would like the output to go</SPAN>
outws <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\Teste_Auto_CLEIA\Imagem'</SPAN>

<SPAN class="comment token">#List all of the workspaces in the previously defined workspace</SPAN>
folders <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListWorkspaces<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"L8*"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"Folder"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token">#Filter to retrieve just the folder(s) with Landsat's bands</SPAN>

<SPAN class="comment token">#Iterate through this list of workspaces and create a new list within each iteration of unstacked raster bands e.g. ['LC80260272014159LGN00_B1', 'LC80260272014159LGN00_B2',...]</SPAN>
<SPAN class="comment token">#Landsat files follows this form: LC80260272014159LGN00_B1.tif, so we need to strip off anything after the "_" and use the first basename as the output name. You can do this with various slicing methods and string manipulation in Python.</SPAN>
<SPAN class="keyword token">for</SPAN> folder <SPAN class="keyword token">in</SPAN> folders<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> folder
&nbsp;&nbsp;&nbsp; rasters <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"L*2*0*_B*.tif"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="comment token"># Filter and does that just the Landsat scenes with path/row from my ROI will be used</SPAN>
&nbsp;&nbsp;&nbsp; name <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outws<SPAN class="punctuation token">,</SPAN> rasters<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>split<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"_"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CompositeBands_management<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">,</SPAN> name<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">)</SPAN><SPAN class="comment token">#Line which allows composite bands for "n" scene(s)</SPAN>
<SPAN class="comment token"># Ao completar, por a frase:</SPAN>
<SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Processing complete"</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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Thank you.

XanderBakker

Can you try this?

<SPAN class="keyword token">import</SPAN> arcpy<SPAN class="punctuation token">,</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">'C:\Teste_Auto_CLEIA\Ima_compactadas'</SPAN>

<SPAN class="comment token"># list all folders in a directory</SPAN>
folders <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListWorkspaces<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"L8*"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"Folder"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># Inside main directory, already have other folders, but the bands are always inside folders that begin with "L8", e.g. "L8227069".</SPAN>
<SPAN class="keyword token">print</SPAN> <SPAN class="string token">"folders: "</SPAN><SPAN class="operator token">+</SPAN>str<SPAN class="punctuation token">(</SPAN>folders<SPAN class="punctuation token">)</SPAN> <SPAN class="comment token">#To see if the search above finds a result</SPAN>

<SPAN class="comment token">#The output files place</SPAN>
outws <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\Teste_Auto_CLEIA\Imagem'</SPAN> <SPAN class="comment token"># Line 10</SPAN>

<SPAN class="keyword token">for</SPAN> folder <SPAN class="keyword token">in</SPAN> folders<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> folder
&nbsp;&nbsp;&nbsp; images <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"L*2*0*_B*.tif"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token">#I changed the user word for an exclusive. This line keeps the bands that I will use</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"images: "</SPAN><SPAN class="operator token">+</SPAN>str<SPAN class="punctuation token">(</SPAN>images<SPAN class="punctuation token">)</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; chk4 = ['band{0}.tif'.format(i) for i in range(0, 7)] # I changed the range</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; print "chk4: "+str(chk4)</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; rasters = [a for a in images if a[-9:] in chk4]</SPAN>
<SPAN class="comment token">##&nbsp;&nbsp;&nbsp; print "rasters: "+str(rasters)</SPAN>
&nbsp;&nbsp;&nbsp; name <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outws<SPAN class="punctuation token">,</SPAN> images<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>split<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"_"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># I wrote "outws"</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CompositeBands_management<SPAN class="punctuation token">(</SPAN>images<SPAN class="punctuation token">,</SPAN> name<SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># I put it</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Finished"</SPAN> <SPAN class="comment token">#23‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍</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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

BrunoDeus

Can you give more instructions?

I inserted some "prints lines" to show the steps executed correctly and I saw the script finds the first folder and him bands. Besides, I don't understand how the original script is fine and this one isn't.

<SPAN class="keyword token">import</SPAN> arcpy<SPAN class="punctuation token">,</SPAN> os
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> 
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\Teste_Auto_CLEIA\Ima_compactadas'</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> 
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> <SPAN class="comment token"># list all folders in a directory</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> folders <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListWorkspaces<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"L8*"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"Folder"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># Inside main directory, already have other folders, but the bands are always inside folders that begin with "L8", e.g. "L8227069".</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"folders: "</SPAN><SPAN class="operator token">+</SPAN>str<SPAN class="punctuation token">(</SPAN>folders<SPAN class="punctuation token">)</SPAN> <SPAN class="comment token">#To see if the search above finds a result </SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> 
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> <SPAN class="comment token">#The output files place</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> outws <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\Teste_Auto_CLEIA\Imagem'</SPAN> <SPAN class="comment token"># Line 10</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> 
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> <SPAN class="keyword token">for</SPAN> folder <SPAN class="keyword token">in</SPAN> folders<SPAN class="punctuation token">:</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> folder
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; images <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListRasters<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"L*2*0*_B*.tif"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token">#I changed the user word for an exclusive. This line keeps the bands that I will use</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"images: "</SPAN><SPAN class="operator token">+</SPAN>str<SPAN class="punctuation token">(</SPAN>images<SPAN class="punctuation token">)</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; chk4 <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">'band{0}.tif'</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>i<SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">for</SPAN> i <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">7</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="comment token"># I changed the range</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"chk4: "</SPAN><SPAN class="operator token">+</SPAN>str<SPAN class="punctuation token">(</SPAN>chk4<SPAN class="punctuation token">)</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; rasters <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>a <SPAN class="keyword token">for</SPAN> a <SPAN class="keyword token">in</SPAN> images <SPAN class="keyword token">if</SPAN> a<SPAN class="punctuation token">[</SPAN><SPAN class="operator token">-</SPAN><SPAN class="number token">9</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="keyword token">in</SPAN> chk4<SPAN class="punctuation token">]</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"rasters: "</SPAN><SPAN class="operator token">+</SPAN>str<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">)</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; name <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>outws<SPAN class="punctuation token">,</SPAN> rasters<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">.</SPAN>split<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"_"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">".tif"</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># I wrote "outws"</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CompositeBands_management<SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">,</SPAN> name<SPAN class="punctuation token">)</SPAN> <SPAN class="comment token"># I put it</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">)</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Finished"</SPAN> <SPAN class="comment token">#23</SPAN>
<SPAN class="punctuation token">.</SPAN><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>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

folders: [u'C:\\Teste_Auto_CLEIA\\Ima_compactadas\\L8227068', u'C:\\Teste_Auto_CLEIA\\Ima_compactadas\\L8227069']
images: [u'LO82270682017183CUB00_B1.TIF', u'LO82270682017183CUB00_B2.TIF', u'LO82270682017183CUB00_B3.TIF', u'LO82270682017183CUB00_B4.TIF', u'LO82270682017183CUB00_B5.TIF', u'LO82270682017183CUB00_B6.TIF', u'LO82270682017183CUB00_B7.TIF']
chk4: ['band0.tif', 'band1.tif', 'band2.tif', 'band3.tif', 'band4.tif', 'band5.tif', 'band6.tif']
rasters: []‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Runtime error
Traceback (most recent call last):
File "<string>", line 20, in <module>
IndexError: list index out of range

Kind regards.

DanPatterson_Retired

List index out of range means that python was expecting a list of something and the 'index' items in the list are out of range... short translation, list can't list stuff in a folder if you are already in the folder, hence there is nothing to list (ie index is out of range since there is nothing). If you fed in a folder path into a script then it would look in the folder, retriebve the items and cycle through them one by one (ie index from 0 to X)

BrunoDeus

I tried it. I saw that my line #18 was inside "for folder in folders", while the original wasn't, but none change happens when I fix it.

<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>rasters<SPAN class="punctuation token">)</SPAN>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ <SPAN class="comment token">#18‍‍</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN></SPAN>

The message error still the same.

Runtime error
Traceback (most recent call last):
File "<string>", line 16, in <module>
IndexError: list index out of range

Thanks!

BrunoDeus

Xander Bakker,

I appreciate your help. Right now, I've 2 scenes (some times are more), the folder from a scene (227/068) is showed below.

*I didn't receive the e-mail about your comments, I will set my account to fix it.

Thankful

XanderBakker

Hi BrunodeDeus , I branched your post to a new question and I have a question if you can provided some additional information before I look at the code, You mention that your band are not named the same, but in the example you included the bands seem to have the same name (apart from 1 and 7). Can you list the rasters that you have and post that list here?

Thanks

DanPatterson_Retired

failure is relative... did you get an error message? and if so, what was it?

indentation on what you have posted is completely wrong, so perhaps you could format the code in case there is an issue in the processing because of it

Rejoindre la conversation

Les membres connectés peuvent publier, suivre les mises à jour, et plus encore. Nouveau ici ? Inscrivez-vous gratuitement.

Liens rapides

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.