C# GDAL 数字图像处理Part6 大气辐射矫正

        最近作业是真的多。

        应不愿意透露姓名的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:大家可能会问，怎么只讲思想不讲太多的实现？其实主要还是没时间，而且各种控件，参数在不同窗口间的传递讲起来也很费力，而且网上都有教程，所以就没怎么讲)