最近作业是真的多。
应不愿意透露姓名的GUI同学的呼吁,今天我们来讲一下我自己使用的、和老师讨论过了但是并不太严谨的、还算过得去的大气辐射矫正方法。
为什么要进行大气辐射校正?说实话我自己也是一知半解,按我的理解是,因为大气中对遥感影像成像有干扰,所以我们得进行大气辐射矫正。
那大气辐射校正应该使用什么方法?本人不知道什么时候异想天开,相处了如下一套方法:
由于红外波段反应的水体不受大气影响,所以理应(此公式纯属个人YY,并不严谨):
思想其实很简单:我们先将水体点与波段极值差值的比值称为水体相对值。
由于红外波段水体反射率不受影响,所以红外波段的水体相对值是一个可以参考的值,我们暂且称之为参考相对值。
然后在待矫正波段,若不受影响的情况下该波段的 水体相对值 应该与 参考相对值 相近,我们求出 该波段水体相对值与参考相对值的差值 ,我们将其视为大气影响导致的数值差异, 最后每一个像素减去这个差异,就完成了一个矫正。
再说一遍,此公式纯属个人YY,不推荐大家用,本贴就是分享个人的方法而已!
按照这个思想,我把代码给敲了出来,给大家搭配示例讲一讲:
在课件中我们有着一个4波段的多光谱图像“2006多光谱Mul01”,在 只进行简单的0-255线性拉伸且Band123分别对应RGB 的情况下获得的图像为:
整体色调时偏向红色,那么可以推断出R通道,也就Band1的数值偏大。
RGB通道的均值:115.8、92.8、65.9
接下来我们进行矫正,在Band4,红外波段,选取9个(数量自选)水体上的点:
在代码中,分别计算: 点的DN值 ÷ 最大最小值之差 ,有九个点获得九个数,对这九个数进行求均值,这样便获得了一个相对数,在这计算得出为 0.29 。
之后,我们用 0.29 x (待矫正波段最大值 - 待矫正波段最小值),就获得了 水体在该波段的期待值(本应该为的值)。
然后获取 待矫正波段下 与上面选取的九个点相同位置 的点的DN值 ,计算这些DN值与期待值的差,共有九个,之后再对这九个差求平均获得一个值,这里计算得10.2。
最后将待矫正波段的每一个像素的值减去这个10.2就得到了校正后的值。
这样就矫正成功了。
至此。
(PS:大家可能会问,怎么只讲思想不讲太多的实现?其实主要还是没时间,而且各种控件,参数在不同窗口间的传递讲起来也很费力,而且网上都有教程,所以就没怎么讲)