Histograma Combinando usando Python para melhorar o processo de mosaico de várias rasters sobrepostas?

Eu estou tentando fazer a correspondência de histograma usando Python para melhorar o processo de mosaico de várias rasters sobrepostas. Estou baseando meu código no encontrado em: http://www.idlcoyote.com/ip_tips/histomatch.html Até o momento, consegui recortar a área sobreposta de duas rasters...