全色和多光谱融合pansharpen的尝试

一个1波段全色地理tiff图像文件的文件名(.tif扩展)。
一个3或4波段的多光谱geotiff图像文件(RGB，近红外波段，按此顺序排列)。
使用GDAL工具
将多光谱Geotiff图像文件重采样到相同的更高维度
作为全色图像的地理tiff文件采用双三次插值。
这个文件被写到磁盘上。然后，对重采样的多光谱和
全色geotiff切片成块(使用gdal_retile.py)。为
每一对，更小的泛锐多光谱geotiff创建使用
4种泛锐化方法(每片产生4个新的MS geotiff):
(1)Brovey,
(2)主成分分析(PCA)，
(3) FIHS(快速强度色相饱和度)，
(4)小波
稍后，对于每种泛锐化方法，所有的瓦片都合并在一起(马赛克)
用四种方法进行泛锐化，形成4个较大的马赛克(使用gdal_merge.py)。

GitHub上找到代码，clone。
打开后发现为python2程序，进行改写为python3程序
运行，发现pywt.ravel_coeffs之类的方法显示不存在，
1、打开cmd终端。

2. pip uninstall pywt 。
3. pip install PyWavelets 。
4. 重装之后正常

运行，发现生成csv失败

  procTileMS  = subprocess.Popen([tileCmdMS],shell=True,stdout=subprocess.PIPE)
  procTileMS.wait()

  procTilePan = subprocess.Popen([tileCmdPan],shell=True,stdout=subprocess.PIPE)
  procTilePan.wait()

cmd命令中运行可以成功生成csv，考虑windows和linux的差别,去掉列表符号,即中括号后可以运行
procTilePan = subprocess.Popen(tileCmdPan,shell=True,stdout=subprocess.PIPE)

融合函数createMosaic,将小图通过gdal_merge进行融合

   def createMosaic( outname, outputTileNames ):
    '''function createMosaic( outname, outputTileNames ): 
    This static (class) method calls the gdal_merge.py Python script (which should be
    installed onto your system) to stitch-together a list of Geotiff image
    files into a single Geotiff. Hence, this function takes the list of Geotiff image
    files in outputTileNames, and passes them to gdal_merge.py such that they may
    be merged into a single Geotiff by a mosaic process. 

    Args:
      outname (str): Output string name of Geotiff mosaic.
      outputTileNames (list): List of strings for Geotiffs that should be mosaic together.

    mosaicCmd = gdalMergeLocation+' -o '+outname+' -of GTiff '
    mosaicCmdList = []
    mosaicCmdList.extend([mosaicCmd])
    mosaicCmdList.extend(outputTileNames)
    mosaicCmd = ' '.join(mosaicCmdList)
    proc = subprocess.Popen([mosaicCmd],shell=True,stdout=subprocess.PIPE)
    proc.wait()

运行到最后的时候,可以看到文件夹里的4种小图片都已经成功生成,但程序运行完成后却没有想象中的4种融合图出现,
推测是融合阶段出了问题,找到融合函数,发现cmd命令中竟然没有加上python解释器....
mosaicCmd = pythonInterpreter + ' ' + mosaicCmd
加上后成功.

生成的4种图分别为
(1) a 3 or 4 band multispectral Geotiff image created using Brovey pan-sharpening
(2) a 3 or 4 band multispectral Geotiff image created using FIHS pan-sharpening
(3) a 3 or 4 band multispectral Geotiff image created using Wavlet pan-sharpening
(4) a 3 or 4 band multispectral Geotiff image created using PCA pan-sharpening

其中小波的图Wavlet 是黑的,不知道为什么...