* GDAL替换影像0值为NAN

由于ArcMap中的影像背景值需要为NAN才能设置为透明，而有的经过预处理的影像数据的背景值为0，可通过以下方式修改:

from osgeo import gdal
import numpy as np

dataset = gdal.Open("G:\mypython3\gdal\data\mosaic.tif")
cols = dataset.RasterXSize
rows = dataset.RasterYSize
proj=dataset.GetProjection()
im_bands = dataset.RasterCount

datas = dataset.ReadAsArray(0,0,cols,rows)
dimen = np.array(datas).shape
t_datas = np.array(datas)

print(dimen)

#修改为NAN值
t_datas=t_datas.astype(np.float32)
t_datas[np.where(t_datas == 0)]=np.nan