机载多光谱扫描仪用于工业温排水监测的技术方法

1.一种机载多光谱扫描仪用于工业温排水监测的技术方法，其特征在于它包括下列步骤：
(1)获取扫描型机载传感器影像和DGPS/IMU组合导航数据；
(2)结合实验室辐射定标结果利用高低温黑体辐射源测量数据进行辐射校正；
(3)提取与影像数据行对应的位置和姿态数据，根据由POS系统的导航角，计算像空间坐标系到地面坐标系的旋转矩阵MIMatrx：
MIMatrx＝MEMatrix×EGMatrix×GIMatrix×ICMatrix×CIMatrix.
其中，CIMatrix为像空间坐标系旋转到传感器坐标系的旋转矩阵，ICMatrix为传感器坐标系旋转到IMU坐标系的旋转矩阵，GIMatrix为IMU坐标系旋转到局部地理坐标系的旋转矩阵，EGMatrix为局部地理坐标系旋转到地心坐标系的旋转矩阵，MEMatrix为地心坐标系旋转到地面坐标系的旋转矩阵；
(4)以沿轨方向为y轴，以穿轨方向为x轴，以扫描行中心为原点建立像空间坐标系，根据扫描型传感器影像中心的共线方程进行几何粗纠正，得到地面点坐标，共线方程式为：
其中，(Xs，Ys，Zs)是投影中心在地面坐标系中的坐标，λ为摄影比例尺，θ为主光轴方向与像点摄影方向的夹角，fsinθ，0，-fcosθ是像点的像空间坐标，M为像空间坐标系到地面坐标系的旋转矩阵，对于扫描型传感器y＝0，(X，Y，Z)是像点对应的地面点坐标；
(5)对原始影像采用直接插值法，利用计算的地面点坐标，构建纠正影 像的范围，遍历原图像每一个像素点(i，j)，其地面坐标为(x，y)，对应纠正影像的位置(x-xmin，y-ymin)，令m＝x-xmin，n＝y-ymin，其中，m0≤m≤m0+1，n0≤n≤n0+1，将各点灰度值按权分配给周围四个像素点，权值取距离的倒数，则
∑DN(m0，n0)＝∑DN(m0，n0)+(m0-m+1)×(n0-n+1)×DN(i，j)
∑P(m0，n0)＝∑P(m0，n0)+(m0-m+1)×(n0-n+1)，
∑DN(m0，n0+1)＝∑DN(m0，n0+1)+(n-n0)×(m0-m+1)×DN(i，j)
∑P(m0，n0+1)＝∑P(m0，n0+1)+(n-n0)×(m0-m+1)，
∑DN(m0+1，n0)＝∑DN(m0+1，n0)+(m-m0)×(n0-n+1)×DN(i，j)
∑P(m0+1，n0)＝∑P(m0+1，n0)+(m-m0)×(n0-n+1)，
∑DN(m0+1，n0+1)＝∑DN(m0+1，n0+1)+(m-m0)×(n-n0)×DN(i，j)
∑P(m0+1，n0+1)＝∑P(m0+1，n0+1)+(m-m0)×(n-n0)，
判断新图像每一个像素，若P(i，j)≠0，则该点像素值DN′(i，j)＝∑DN(i，j)/∑P(i，j)；
(6)基于地面点坐标的影像拼接；
(7)用海面实测温度对校正后数据进行标定，以外围均匀海温为本底，以能反映潮汐变化规律为标准，以温度升高0.5度为区间制作温度场分布图。

机载多光谱扫描仪用于工业温排水监测的技术方法\n技术领域\n[0001] 本发明涉及热电厂、核电站等沿海工业基地运行过程中排放冷却水的空间分布信息的自动获取与分析技术领域，特别是基于实现航空图像、实验室辐射定标结果、机上黑体定标数据、POS系统数据等多元数据集成的高分辨率遥感技术进行温排水量化监测方法，并应用于工程方案的选址及环境保护措施的制定。\n背景技术\n[0002] 我国能源消费量居世界第二位，能源问题更加严峻。核动力也将再次成为化石燃料强有力的替代能源，未来3年内中国将开建9个核电站，16台核电机组，装机容量在1000万千瓦以上，这意味着未来3年的核电建设总量将超过过去23年的总和。电力企业运行需要大量的冷却水，同时也要排放大量的温排水，势必对工业基地附近海域温度场造成影响，不同季节、不同潮汐状况下影响情况不同。这种变化反过来也可能影响电力企业的运行。这些都需要客观真实的海面温度场数据来说明。\n[0003] 当前，卫星遥感技术在我国水环境监测中发挥极其重要的作用，但是受空间分辨率和重访周期的制约对于电力企业温排水的定量化监测分析难以有效进行。而航空遥感技术要满足这种应用需求，还需要解决成图质量和测量精度、稳定性方面的问题，而且要尽可能提高作业效率。因此，需要研究针对温排水监测的扫描型遥感成像数据处理方法。\n发明内容\n[0004] 本发明的目的是，针对现有的温排水监测技术和方法中存在的缺陷或不足，提供一种新的高分辨率的温排水遥感定量化监测技术方法，并用于取排水口工程方案的选址及环保措施的制定。该发明能实现温排水的识别、信息提取与量化分析，对温排水给附近海域造成的影响程度、范围以及是否降低电力企业冷却水循环效率等进行评估。\n[0005] 为了实现上述任务，本发明采用以下技术措施：\n[0006] 一种工业温排水的高分辨率遥感量化监测技术方法。其特征在于：通过实验室试验测试标定遥感器黑体源，基于机上高低温黑体进行遥感数据的辐射校正，利用POS系统数据，构建扫描型遥感成像的几何粗校正模型，进行图像几何畸变纠正、定位和航带拼接，以现场实测水温标定校正后的影像，根据能够反映潮汐变化规律温度升高0.5度为区间，从而实现一种工业温排水高分辨率遥感监测方法。\n[0007] 本技术方法包含下列步骤：\n[0008] 1)制定飞行计划，获取监测区的遥感数据。\n[0009] 2)热红外通道传感器的机上定标，采用内部温度参考源法，热遥感器附有内部温度参考源，采用在旋转扫描镜角视场的两侧放置两个黑体辐射源的形式，这两个黑体辐射源的温度被精确控制，并设置为地面监测目标的“最冷”与“最热”，对于每一条扫描线，扫描器先记录冷参考源的辐射温度，然后扫描地面，最后记录热参考源的辐射温度，所有的信号均记录下来，两个温度源也随图像数据记录，用来推算整幅热图像的辐射温度，也可由一个绝对辐射值作参考与其他热扫描仪输出值比对。转换公式基于这样一个假设，即普朗克函数在一定温度范围内近似为线性的。通过实验室辐射定标确定机上定标黑体的辐射亮度与亮温的转换关系。\n[0010] 3)相对卫星平台，航空平台受到气流影响较大，姿态(侧滚角，俯仰角和偏转角)变化剧烈；此外，飞机自身状态(如飞行速度变化、航高变化、航偏)以及摆扫式的工作模式等影响，导致机载多光谱影像几何畸变严重。DGPS/IMU组合导航与定位系统能以较高的精度和频率提供传感器的位置和姿态，可直接测定影像的外方位元素，进行扫描型机载多光谱扫描仪影像的几何校正。\n[0011] 4)测量获取GPS天线相位中心、IMU及机载传感器之间空间关系，直接对POS系统获取的GPS天线相位中心的空间坐标(X，Y，Z)及IMU系统获取的侧滚角、俯仰角、航迹角进行数据处理，获取航空影像像元成像瞬间的扫描行中心点三维空间坐标(XS，YS，ZS)及其三个姿态角，实现无地面控制点辅助条件下直接恢扫描型复航空遥感的成像过程。计算插值影像的四个角点坐标，然后基于地面点坐标作影像的拼接。\n[0012] 5)根据完成辐射校正、几何校正、航带拼接的航空遥感数据，以外围均匀海水表面水温为本底，以温升0.5度为区间，制作温度场分布图。\n[0013] 本发明与现有技术相比，具有以下优点：\n[0014] 第一，采用结合POS系统数据的方法，可以实现一种基于摆扫成像的高分辨率遥感方法手段，能够查清工业温排水的扩散规律。\n[0015] 第二，采用实验室定标和机上定标结合的方法，提供了一种定量化数据分析的手段，解决多时相观测比较的难题。\n[0016] 第三，采用本发明的技术方法，实现航空红外传感器的数据获取和处理流程，有利于进行依据潮汐规律的遥感探测分析。\n[0017] 第四，采用本发明的技术方法，实现扫描型航空红外传感器的数据获取和处理流程，即实现高分辨率要求又兼顾了大视场扫描，从而提高效率。\n附图说明\n[0018] 图1扫描型航空传感器的成像方式\n[0019] 图2对地面和高低温黑体的扫描成像示意\n[0020] 图3为本发明的技术流程\n[0021] 图4为扫描型传感器的原始影像(左)和经过几何校正后影像(右)\n[0022] 图5为利用本发明输出的海面温度场分布图\n具体实施方式\n[0023] 下面结合附图，对本发明作进一步详细描述：\n[0024] 第一步，遥感影像和参数测量数据获取。\n[0025] 在执行一组飞行任务前，对遥感器进行实验室辐射定标，获得定标数据。制定飞行计划，获取监测区遥感数据和同步POS系统数据，其中遥感数据中包含高低温黑体测量数据。\n[0026] 第二步，热红外遥感影像辐射校正。\n[0027] 首先通过实验室辐射定标过程，确定高低温黑体的辐射亮度与亮温的关系。机载热红外测量数据的辐射定标，通过传感器内的高低温黑体来完成，采用下面的转换公式：\n[0028] L(i，j) ＝ [(LH(i)-LL(i))×DN(i，j)-LH(i)×DNL+LL(i)× DNH(i)]/[DNH(i)-DNL(i)]\n[0029] 其中：L为辐亮度值，(i，j)代表第i行第j列的像元，DN为图像测量值，下标H和L分别代表对应高温黑体和低温黑体的测量值。\n[0030] 第三步，POS系统同步测量数据处理。\n[0031] 利用POSPac软件，对获取的原始DGPS/IMU数据进行预处理。处理POS数据文件要解决两个问题，一是时间匹配问题和延迟处理；二是坐标系转换。\n[0032] 根据触发文件(EVENT)，采用一次线性内插算法，提取POS数据中的传感器姿态信息。首先对EVENT进行分割，判断若相邻触发事件时间间隔大于10倍的单位时间间隔，则认为这是一条航带的结尾，若相邻EVENT事件间隔大于1.5倍的单位时间间隔，则认为这里存在漏行，需要填补。由于EVENT触发的滞后性，需要判断外方位元素(EO)的偏移行数。\n[0033] 若传感器与POS之间没有触发信息，则使用扫描影像记录的世界协调时间(UTC)与POS记录的GPS时间进行对比，根据影像记录的起始和中止时间，内插出每条扫描线的记录时间，然后从POS数据中提取相应的传感器姿态信息，其中需要考虑GPS时间与UTC时间之间的误差。\n[0034] 第四步，计算像空间到地面坐标系的旋转矩阵。\n[0035] 像空间坐标系旋转到地面坐标系可分解为以下步骤：像空间坐标系旋转到传感器坐标系(旋转矩阵CIMatrix)；传感器坐标系旋转到IMU坐标系(旋转矩阵ICMatrix)；IMU坐标系旋转到局部地理坐标系(旋转矩阵GIMatrix)；局部地理坐标系旋转到地心坐标系(旋转矩阵EGMatrix)；地心坐标系旋转到地面坐标系(旋转矩阵MEMatrix)。\n[0036] 则像空间坐标系旋转到地面坐标系的旋转矩阵：\n[0037] MIMatrx＝MEMatrix×EGMatrix×GIMatrix×ICMatrix×CIMatriEG\n[0038] 第五步，计算扫描行的外方位元素。\n[0039] 线元素：设瞬时摄影中心在WGS84下的坐标(X，Y，Z)，转为大地经纬度(B，L，H)，再转为地面坐标系坐标(x，y)，以飞行区域的中心(B0，L0)为地面坐标系原点。\n[0040] 角元素：建立以X轴为主轴的 系统，ω为旁向倾角，为航向倾角，κ为像片旋角，旋转矩阵的表示形式为：\n[0041] \n[0042] 则角元素 为：\n[0043] \n[0044] 第六步，计算像点的地面点坐标。\n[0045] 对于线阵摆扫式成像方式，根据成像特点可采用共线方程：\n[0046] \n[0047] 令：\n[0048] \n[0049] 则对应地面点坐标为：\n[0050] \n[0051] 其中(Xs，Ys，Zs)是投影中心在地面坐标系中的坐标，λ为摄影比例尺，θ为主光轴方向与像点摄影方向的夹角，(fsinθ，0，-fcosθ)是像点的像空间坐标，M为像空间坐标系到地面坐标系的旋转矩阵，对于扫描型传感器y＝0，(X，Y，Z)是像点对应的地面点坐标。\n[0052] 第七步，影像插值。\n[0053] 对原始影像采用直接插值法，利用计算的地面点坐标，构建纠正影像的范围，遍历原图像每一个像素点(i，j)，其地面坐标为(x，y)，对应纠正影像的位置(x-xmin，y-ymin)，令m＝x-xmin，n＝y-ymin(m0≤m≤m0+1，n0≤n≤n0+1)，将各点灰度值按权分配给周围四个像素点，权值取距离的倒数。则\n[0054] ∑DN(m0，n0)＝∑DN(m0，n0)+(m0-m+1)×(n0-n+1)×DN(i，j)\n[0055] ∑P(m0，n0)＝∑P(m0，n0)+(m0-m+1)×(n0-n+1)，\n[0056] ∑DN(m0，n0+1)＝∑DN(m0，n0+1)+(n-n0)×(m0-m+1)×DN(i，j)\n[0057] ∑P(m0，n0+1)＝∑P(m0，n0+1)+(n-n0)×(m0-m+1)，\n[0058] ∑DN(m0+1，n0)＝∑DN(m0+1，n0)+(m-m0)×(n0-n+1)×DN(i，j)\n[0059] ∑P(m0+1，n0)＝∑P(m0+1，n0)+(m-m0)×(n0-n+1)，\n[0060] ∑DN(m0+1，n0+1)＝∑DN(m0+1，n0+1)+(m-m0)×(n-n0)×DN(i，j)[0061] ∑P(m0+1，n0+1)＝∑P(m0+1，n0+1)+(m-m0)×(n-n0)\n[0062] 判断新图像每一个像素，若P(i，j)≠0，则该点像素值DN′(i，j)＝∑DN(i，j)/∑P(i，j)。\n[0063] 第八步，基于地面点坐标的影像拼接。\n[0064] 根据影像插值的范围，计算插值影像的四个角点坐标，然后基于地面点坐标作影像的拼接；还需在影像的重叠部分引入少量同名点，对拼接的影像做局部的调整。\n[0065] 第九步，海面温度场制图\n[0066] 根据完成辐射校正、几何校正、航带拼接的航空遥感数据，用海面实测温度进行标定。以最外围均匀海温为本底，以能反映潮汐变化规律为标准，以温度升高0.5度为区间制作温度场分布图。