十几张影像图分别逐像元分波段怎么做比如mse评价指标

有一个参考影像，怎么将每张图的mse逐像元做出来并可视化，python或者idl