实习四 水体遥感技术.doc

实习四 水体遥感技术应用专题 实习目的 星载或机载的遥感器所接收到的来自水面的辐射亮度值由下列各项组成： 为传感器所测量得到的总辐射亮度；为离水辐射亮度；为水面对大气下行辐射及太阳直射辐射的反射所构成的辐射亮度；为路径辐射亮度；大气漫反射透过率。 离水辐射是由水分子及水中悬浮物质（如：泥沙、各种藻类）对传入水中的太阳直接辐射的后向散射所构成。则包括两个部分：海洋表面对太阳直射辐射的反射和天空下行辐射经海面反射而进入传感器的那部分能量。代表太阳直射辐射未经与地表相撞，直接由大气分子及气溶胶的散射而进入传感器的那部分能量。 从上述水体光学遥感的基本辐射传输原理可以看出，离水辐射中包含了水体中某些参数的信息，因此，可以使用该项观测值来进行水中某些参数的反演（估算）。下面以叶绿素浓度与泥沙含量估算为例进行说明（本实例中，为了计算方便，未进行洋面观测值的大气校正处理，也就是没有去除传感器观测量中的与，不考虑大气透过率的影响，直接将观测值反射率作为离水辐射率来做处理）。 4.2.2 水体的反射率光谱特征 纯净的水体反射率非常低，尤其在红外波谱区域，其反射率几乎为零，因此，通常在遥感影像上，水体显示为暗色调（如下图中的河流与湖泊）。 4.2.3 热红外水面温度反演原理 热红外辐射计探测到的总辐射由四部分组成：第一部分表示的是到达卫星辐射计的海表面发射的辐射；第二部分是卫星辐射计探测到的沿观测方向的大气上行发射的辐射；第三部分是卫星辐射计探测到的大气下行辐射的、到达海面后又经海面反射的辐射；第四部分为太阳直射辐射到达海面后又经海面反射的辐射。 上式中，为传感器接受到的辐射；为海表面热辐射；为大气下行辐射；为太阳直接辐射；为大气上行辐射；为大气透过率；为海洋表层的发射率；为海洋表层的反射率。 在对海洋表层进行温度探测时，通常假设海水为黑体，因此上式中为1，为0，可简化为：，也就是仅仅需考虑大气透过率和大气自身的热辐射。现有的地表温度反演方法中，以多通道法应用最为广泛。因为大气对不同波长不同时间的红外遥感有不同的影响效应，根据大气对不同波段的电磁辐射的影响不同，可以用不同波段测量的线性组合来消除大气的影响，从而得到海表温度（SST）。 下面以双通道法（又称分裂窗法）为例，给出该模型的理论推导过程。在海面为黑体的近似条件下，星载传感器所测得的辐射亮度值为: 此处为SST，为整层大气透过率。等式右边第一项代表海面热辐射从大气削弱后到达传感器的辐射能量，第二项为大气上行辐射的总贡献，由于假定海面为黑体，所以大气下行辐射的贡献率为零。 设为层气体的透过率，对于热红外波段大气散射可以忽略，故： ； 因此，层厚为的空气层的热辐射亮度为。该层热辐射经上层大气削弱后，最后被传感器所接受的能量应为：；所以整层大气上行热辐射的总贡献应为：。 这里再计算海洋表层热辐射与传感器接收热辐射之间的差值：（由于），因此， 注： 对于上面公式中的，可以采用亮度温度的形式来表示为。最后，由于大气的实际参数不易获取，其区域差异较大，因此引入如下两个近似条件： 由于大气窗口水汽的吸收比较弱，所以取如下线性近似： 这里，为大气的水汽密度；此处为水汽的吸收系数，假定它只是波段的函数，不随高度而变（事实上还是气压和温度的弱函数）。 把黑体辐射作近似计算： （多项式展开） 此式的含义为Z高度处的黑体辐射亮度，可以用处的黑体辐射亮度的泰勒一阶近似代替（事实上差值越大，误差越大）。 有了上述近似条件，将遥感器所测得的辐射亮度转换成亮度温度： 注：，由近似条件②，，又，将该式中的分别用近似条件②展开， ，由于，上式可简化为： ，将近似条件①引入，用替换，则可以得到： 上式表明：积分值只与大气状况有关，与波长无关。因此，对于多个波段，分别计算得到： 4.2.4 水体参数遥感反演原理 叶绿素浓度反演为例 基于叶绿素的上述反射特性，一般采用上图中的强吸收（0.45～0.50）和强反射（0.52～0.58）范围内的观测值来进行叶绿素浓度的计算。在仅仅考虑单次散射作用条件下，光学观测量与水体的固有特性之间具有如下关系（可用Gordon公式来描述）： 上式中，为反射率，分别为吸收率和散射率。一般来讲，吸收率与散射率与叶绿素浓度之间存在着非线性关系（幂函数），基于这种原理，可以采用统计建模的方法来建立叶绿素浓度反演模型，形如： 4.2.5 实习内容 （1）基于水体的反射率光谱特征，采用植被指数阈值法来进行水体像元识别，并进行矢量化处理； （2）采用双通道法（分裂窗法）实现水面温度遥感反演； （3）采用统计模型实现水中叶绿素浓度的反演。 实习数据实习步骤 NDVI计算：（在本实例中，波段1与2分别对应MODIS数据的红、近红外