I am doing land use classification using multiple scenes from the Landsat satellite for 12 months. In each scene, I have removed the cloudy pixels and replaced them with pixels with no data value. Now I composited the 7 bands from 12 scenes together, but most of the bands have missing values, and the classified image also has missing values from those bands used.

Is there a way I could interpolate the missing values in each band from an average of the corresponding pixel from all other bands with values? I would like to do change detection afterward.

I am doing the analysis in ArcGIS Pro.