I want to create composite band image with sentinel 2 data through arcpy.

Question

As we know sentinel satellite data have many sub folder. I want to create composite band of each datasets with the band 2, band 3, band 4, band 8 which are located in one folder(10m) and band 11, band 12 which are located in other folder(20m).I am making a arcpy code for that but code was not run correctly. Here is the code

import arcpy, os

#arcpy.env.workspace = "L:\Arcpy\Data\Satellite_image_raw\New folder"

basepath = "L:\Arcpy\Data\Satellite_image_raw\New folder"

endswith = ("_B02_10m.jp2", "_B03_10m.jp2", "_B04_10m.jp2", "_B08_10m.jp2", "_B11_20m.jp2", "_B12_20m.jp2")

output = "L:\Arcpy\Data\Layer_stack"

rasterlist = []

#folders = arcpy.ListWorkspaces()

#for folder in folders:
#arcpy.env.workspace = folder
for root, dirs, files in os.walk(basepath):
     for rasterlayer in files:
           if rasterlayer.endswith(endswith):
              #print(os.path.join(root, rasterlayer))
              rasterlist.append(rasterlayer)
              print "image:" + str(rasterlist)
name = os.path.join(output, rasterlist[1].split("_")[0] + ".img")
arcpy.CompositeBands_management(rasterlist, name)
print(rasterlist)

DanPatterson · Answer

you changed the code

for root, dirs, files in os.walk(basepath):

I suspect you are trying to piece something together without fully understanding the implications.

I would suggest using the builtin tools to see if you can get a workflow established before continuing down the path of trying to script it.

Good luck

SiddharthaDhara · Answer

This time I am trying this way

import arcpy, os

arcpy.env.workspace = ("D:\Composite_band\New folder")

basepath = ("D:\Composite_band\New folder")

endswith = ("_B02_10m.jp2", "_B03_10m.jp2", "_B04_10m.jp2", "_B08_10m.jp2", "_B11_20m.jp2", "_B12_20m.jp2")

output = ("D:\Composite_band\Layer_stack")

folders = arcpy.ListWorkspaces()

for folder in folders:
      arcpy.env.workspace = folder
      for root, dirs, files in os.walk(basepath):
            for file in files:
                  if file.endswith(endswith):
                  print(os.path.join(root, file))
                  image = []
                  image = file
                  print(image)
                  name = os.path.join(output, image[1].split("_")[0] + ".img")
                  arcpy.CompositeBands_management(image, name)
                  print(image)

But this time code is able to print first selected band of the first folder after that it gives error like this:  D: \Composite_band\New folder\S2A_MSIL2A_20210104T044201_N0214_R033_T45QVD_20210104T063322.SAFE\GRANULE\L2A_T45QVD_A028916_20210104T044202\IMG_DATA\R10m\T45QVD_20210104T044201_B02_10m.jp2
T45QVD_20210104T044201_B02_10m.jp2
Traceback (most recent call last):
File "D:\Composite_band\composite_band_script_1.py", line 34, in <module>
arcpy.CompositeBands_management(image, name)
File "C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\management.py", line 13595, in CompositeBands
raise e
ExecuteError: ERROR 000271: Cannot open the input datasets
Failed to execute (CompositeBands).

It means after printing first band the code is not able to open the input folder.

DanPatterson · Answer

you had better throw some print statements in there because you are assuming rasters[0] is indeed a raster and not something else and you won't be able to get "name" until you are sure it is a single file name and not a raster, which I suspect it is