用python gdal库来读取nc文件进行波段计算，添加坐标参考系，导出成tif

有一个项目是要求对数据进行修正，我这边有原始的nc栅格数据文件以及一个修正文件，要做的事情就是两个文件一一相加，并导出结果。下边是我用python做的实力

'''
Created on 2018年1月11日

@author: cyemeng
'''
import gdal, osr, ogr, os
from os.path import split

driver_nc = gdal.GetDriverByName('netCDF')

def array2raster(newRasterfn,rasterOrigin,xstep,ystep,array):
    """
    newRasterfn：输出tif路径
    rasterOrigin:原始栅格数据路径
    xstep：x方向栅格数据的宽度
    ystep:y方向栅格数据的宽度
    array:计算后的栅格数据
    """
    cols = array.shape[1]
    rows = array.shape[0]
    originX = rasterOrigin[0]
    originY = rasterOrigin[1]
    driver = gdal.GetDriverByName('GTiff')
    outRaster = driver.Create(newRasterfn,cols, rows, 1,gdal.GDT_Float32)
    outRaster.SetGeoTransform((originX, xstep, 0, originY, 0, ystep))
    outband = outRaster.GetRasterBand(1)
    outband.WriteArray(array)
    outRasterSRS = osr.Spat