修改gdal python 程序 ,增加无值区域的判断

在海洋区域的温度、盐度、风场等信息处理中,由于陆地和云层遮挡部分无值,需要剔除这些无值区域(-999)进行计算。本文介绍了如何通过掩膜方法在GDAL Python脚本中增加无值区域的判断,将无值区域在结果影像中标记为-999,以避免影响栅格运算。
摘要由CSDN通过智能技术生成

是这样的,在栅格运算过程中,要对每一块像元处理,因为有的像元不存在值,项目背景是,获取的为海洋区域的温度、盐度、风场等信息,因此陆地部分无值,海洋中被云层遮挡部分无值,所以,要将所有的无值区域不参与计算。


ENVI中打开,看到海洋区域中明显的黑点,显示为-999.

栅格运算中,-999需要剔除掉,才能和其他像元点计算,栅格计算重要用到了一个python脚本(需要用到的可以去我的资源页下载)。

下面介绍掩膜方法增加无值区域的判定,并在最终结果影像中显示无值区域为-999,这一部分是原生gdal python没有的。

            # create empty buffer to mark where nodata occurs
            myNDVs=numpy.zeros(myBufSize)
            myNDVs.shape=(nYValid,nXValid)

	    my999=numpy.zeros(myBufSize)
            my999.shape=(nYValid,nXValid)

myNDVs是原始无值区域掩膜,其用处是将所有无值的区域标注出来,用来以后无值区域的判定。

在这里增加了my999数组,初始值为0,大小为栅格图像的大小。

            for i,Alpha i
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值