CN101551275B - 机载多光谱扫描仪用于工业温排水监测的技术方法 - Google Patents
机载多光谱扫描仪用于工业温排水监测的技术方法 Download PDFInfo
- Publication number
- CN101551275B CN101551275B CN2009100503122A CN200910050312A CN101551275B CN 101551275 B CN101551275 B CN 101551275B CN 2009100503122 A CN2009100503122 A CN 2009100503122A CN 200910050312 A CN200910050312 A CN 200910050312A CN 101551275 B CN101551275 B CN 101551275B
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- image
- data
- rotation matrix
- discharge water
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Radiation Pyrometers (AREA)
Abstract
本发明公开了一种机载多光谱扫描仪用于工业温排水监测的技术方法,它是利用扫描型机载遥感影像数据和同步测量DGPS/IMU(又称为POS系统)组合导航数据量化监测沿海工业基地温排水的方法。本发明集成应用POS系统数据、实验室辐射定标数据、机上高低温黑体数据,海面遥感影像数据,完成辐射校正、几何校正、航带拼接和制图,有效发挥了航空遥感技术的高分辨率和高精度潜力,能够监测温排水根据潮汐变化的扩散规律。本发明的技术应用对于核电等电力企业及沿海工业基地的规划和发展具有重要意义。
Description
技术领域
本发明涉及热电厂、核电站等沿海工业基地运行过程中排放冷却水的空间分布信息的自动获取与分析技术领域,特别是基于实现航空图像、实验室辐射定标结果、机上黑体定标数据、POS系统数据等多元数据集成的高分辨率遥感技术进行温排水量化监测方法,并应用于工程方案的选址及环境保护措施的制定。
背景技术
我国能源消费量居世界第二位,能源问题更加严峻。核动力也将再次成为化石燃料强有力的替代能源,未来3年内中国将开建9个核电站,16台核电机组,装机容量在1000万千瓦以上,这意味着未来3年的核电建设总量将超过过去23年的总和。电力企业运行需要大量的冷却水,同时也要排放大量的温排水,势必对工业基地附近海域温度场造成影响,不同季节、不同潮汐状况下影响情况不同。这种变化反过来也可能影响电力企业的运行。这些都需要客观真实的海面温度场数据来说明。
当前,卫星遥感技术在我国水环境监测中发挥极其重要的作用,但是受空间分辨率和重访周期的制约对于电力企业温排水的定量化监测分析难以有效进行。而航空遥感技术要满足这种应用需求,还需要解决成图质量和测量精度、稳定性方面的问题,而且要尽可能提高作业效率。因此,需要研究针对温排水监测的扫描型遥感成像数据处理方法。
发明内容
本发明的目的是,针对现有的温排水监测技术和方法中存在的缺陷或不足,提供一种新的高分辨率的温排水遥感定量化监测技术方法,并用于取排水口工程方案的选址及环保措施的制定。该发明能实现温排水的识别、信息提取与量化分析,对温排水给附近海域造成的影响程度、范围以及是否降低电力企业冷却水循环效率等进行评估。
为了实现上述任务,本发明采用以下技术措施:
一种工业温排水的高分辨率遥感量化监测技术方法。其特征在于:通过实验室试验测试标定遥感器黑体源,基于机上高低温黑体进行遥感数据的辐射校正,利用POS系统数据,构建扫描型遥感成像的几何粗校正模型,进行图像几何畸变纠正、定位和航带拼接,以现场实测水温标定校正后的影像,根据能够反映潮汐变化规律温度升高0.5度为区间,从而实现一种工业温排水高分辨率遥感监测方法。
本技术方法包含下列步骤:
1)制定飞行计划,获取监测区的遥感数据。
2)热红外通道传感器的机上定标,采用内部温度参考源法,热遥感器附有内部温度参考源,采用在旋转扫描镜角视场的两侧放置两个黑体辐射源的形式,这两个黑体辐射源的温度被精确控制,并设置为地面监测目标的“最冷”与“最热”,对于每一条扫描线,扫描器先记录冷参考源的辐射温度,然后扫描地面,最后记录热参考源的辐射温度,所有的信号均记录下来,两个温度源也随图像数据记录,用来推算整幅热图像的辐射温度,也可由一个绝对辐射值作参考与其他热扫描仪输出值比对。转换公式基于这样一个假设,即普朗克函数在一定温度范围内近似为线性的。通过实验室辐射定标确定机上定标黑体的辐射亮度与亮温的转换关系。
3)相对卫星平台,航空平台受到气流影响较大,姿态(侧滚角,俯仰角和偏转角)变化剧烈;此外,飞机自身状态(如飞行速度变化、航高变化、航偏)以及摆扫式的工作模式等影响,导致机载多光谱影像几何畸变严重。DGPS/IMU组合导航与定位系统能以较高的精度和频率提供传感器的位置和姿态,可直接测定影像的外方位元素,进行扫描型机载多光谱扫描仪影像的几何校正。
4)测量获取GPS天线相位中心、IMU及机载传感器之间空间关系,直接对POS系统获取的GPS天线相位中心的空间坐标(X,Y,Z)及IMU系统获取的侧滚角、俯仰角、航迹角进行数据处理,获取航空影像像元成像瞬间的扫描行中心点三维空间坐标(XS,YS,ZS)及其三个姿态角,实现无地面控制点辅助条件下直接恢扫描型复航空遥感的成像过程。计算插值影像的四个角点坐标,然后基于地面点坐标作影像的拼接。
5)根据完成辐射校正、几何校正、航带拼接的航空遥感数据,以外围均匀海水表面水温为本底,以温升0.5度为区间,制作温度场分布图。
本发明与现有技术相比,具有以下优点:
第一,采用结合POS系统数据的方法,可以实现一种基于摆扫成像的高分辨率遥感方法手段,能够查清工业温排水的扩散规律。
第二,采用实验室定标和机上定标结合的方法,提供了一种定量化数据分析的手段,解决多时相观测比较的难题。
第三,采用本发明的技术方法,实现航空红外传感器的数据获取和处理流程,有利于进行依据潮汐规律的遥感探测分析。
第四,采用本发明的技术方法,实现扫描型航空红外传感器的数据获取和处理流程,即实现高分辨率要求又兼顾了大视场扫描,从而提高效率。
附图说明
图1扫描型航空传感器的成像方式
图2对地面和高低温黑体的扫描成像示意
图3为本发明的技术流程
图4为扫描型传感器的原始影像(左)和经过几何校正后影像(右)
图5为利用本发明输出的海面温度场分布图
具体实施方式
下面结合附图,对本发明作进一步详细描述:
第一步,遥感影像和参数测量数据获取。
在执行一组飞行任务前,对遥感器进行实验室辐射定标,获得定标数据。制定飞行计划,获取监测区遥感数据和同步POS系统数据,其中遥感数据中包含高低温黑体测量数据。
第二步,热红外遥感影像辐射校正。
首先通过实验室辐射定标过程,确定高低温黑体的辐射亮度与亮温的关系。机载热红外测量数据的辐射定标,通过传感器内的高低温黑体来完成,采用下面的转换公式:
L(i,j)=[(LH(i)-LL(i))×DN(i,j)-LH(i)×DNL+LL(i)× DNH(i)]/[DNH(i)-DNL(i)]
其中:L为辐亮度值,(i,j)代表第i行第j列的像元,DN为图像测量值,下标H和L分别代表对应高温黑体和低温黑体的测量值。
第三步,POS系统同步测量数据处理。
利用POSPac软件,对获取的原始DGPS/IMU数据进行预处理。处理POS数据文件要解决两个问题,一是时间匹配问题和延迟处理;二是坐标系转换。
根据触发文件(EVENT),采用一次线性内插算法,提取POS数据中的传感器姿态信息。首先对EVENT进行分割,判断若相邻触发事件时间间隔大于10倍的单位时间间隔,则认为这是一条航带的结尾,若相邻EVENT事件间隔大于1.5倍的单位时间间隔,则认为这里存在漏行,需要填补。由于EVENT触发的滞后性,需要判断外方位元素(EO)的偏移行数。
若传感器与POS之间没有触发信息,则使用扫描影像记录的世界协调时间(UTC)与POS记录的GPS时间进行对比,根据影像记录的起始和中止时间,内插出每条扫描线的记录时间,然后从POS数据中提取相应的传感器姿态信息,其中需要考虑GPS时间与UTC时间之间的误差。
第四步,计算像空间到地面坐标系的旋转矩阵。
像空间坐标系旋转到地面坐标系可分解为以下步骤:像空间坐标系旋转到传感器坐标系(旋转矩阵CIMatrix);传感器坐标系旋转到IMU坐标系(旋转矩阵ICMatrix);IMU坐标系旋转到局部地理坐标系(旋转矩阵GIMatrix);局部地理坐标系旋转到地心坐标系(旋转矩阵EGMatrix);地心坐标系旋转到地面坐标系(旋转矩阵MEMatrix)。
则像空间坐标系旋转到地面坐标系的旋转矩阵:
MIMatrx=MEMatrix×EGMatrix×GIMatrix×ICMatrix×CIMatriEG
第五步,计算扫描行的外方位元素。
线元素:设瞬时摄影中心在WGS84下的坐标(X,Y,Z),转为大地经纬度(B,L,H),再转为地面坐标系坐标(x,y),以飞行区域的中心(B0,L0)为地面坐标系原点。
第六步,计算像点的地面点坐标。
对于线阵摆扫式成像方式,根据成像特点可采用共线方程:
令:
则对应地面点坐标为:
其中(Xs,Ys,Zs)是投影中心在地面坐标系中的坐标,λ为摄影比例尺,θ为主光轴方向与像点摄影方向的夹角,(fsinθ,0,-fcosθ)是像点的像空间坐标,M为像空间坐标系到地面坐标系的旋转矩阵,对于扫描型传感器y=0,(X,Y,Z)是像点对应的地面点坐标。
第七步,影像插值。
对原始影像采用直接插值法,利用计算的地面点坐标,构建纠正影像的范围,遍历原图像每一个像素点(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)。
第八步,基于地面点坐标的影像拼接。
根据影像插值的范围,计算插值影像的四个角点坐标,然后基于地面点坐标作影像的拼接;还需在影像的重叠部分引入少量同名点,对拼接的影像做局部的调整。
第九步,海面温度场制图
根据完成辐射校正、几何校正、航带拼接的航空遥感数据,用海面实测温度进行标定。以最外围均匀海温为本底,以能反映潮汐变化规律为标准,以温度升高0.5度为区间制作温度场分布图。
Claims (1)
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度为区间制作温度场分布图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100503122A CN101551275B (zh) | 2009-04-30 | 2009-04-30 | 机载多光谱扫描仪用于工业温排水监测的技术方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100503122A CN101551275B (zh) | 2009-04-30 | 2009-04-30 | 机载多光谱扫描仪用于工业温排水监测的技术方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101551275A CN101551275A (zh) | 2009-10-07 |
CN101551275B true CN101551275B (zh) | 2011-07-27 |
Family
ID=41155641
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100503122A Active CN101551275B (zh) | 2009-04-30 | 2009-04-30 | 机载多光谱扫描仪用于工业温排水监测的技术方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101551275B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113643386A (zh) * | 2021-10-13 | 2021-11-12 | 深圳市商汤科技有限公司 | 标定方法及装置、电子设备及计算机可读存储介质 |
US11830635B1 (en) * | 2022-09-06 | 2023-11-28 | First Institute of Oceanography, Ministry of Natural Resources | Method and system for determining background water temperature of thermal discharge from operating nuclear power plants based on remote sensing |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101876701B (zh) * | 2010-07-02 | 2012-10-03 | 中国测绘科学研究院 | 一种侧视雷达遥感影像定位方法 |
CN102610000A (zh) * | 2012-03-14 | 2012-07-25 | 江苏钱旺网络科技有限公司 | 基于Wi-Fi技术的员工考勤定位的方法 |
CN103512661B (zh) * | 2012-12-07 | 2016-03-02 | 核工业北京地质研究院 | 一种基于遥感影像快速识别盐湖水域隐伏断裂构造的方法 |
CN103983361B (zh) * | 2014-06-06 | 2017-04-19 | 昆明北方红外技术股份有限公司 | 在线式网络测温热像仪标定方法 |
CN104050649A (zh) * | 2014-06-13 | 2014-09-17 | 北京农业信息技术研究中心 | 农业遥感系统 |
CN104361591B (zh) * | 2014-11-13 | 2017-06-30 | 环境保护部卫星环境应用中心 | 电厂温排水监测中区域替代基准温度的获取方法和系统 |
CN105241429B (zh) * | 2015-09-22 | 2017-08-25 | 中国科学院上海技术物理研究所 | 一种基于航空遥感的近海工业温排水提取方法 |
CN107633487A (zh) * | 2017-09-01 | 2018-01-26 | 天津津航技术物理研究所 | 一种航空摆扫式多光谱扫描仪图像的系统级相对辐射校正方法 |
CN107705272A (zh) * | 2017-11-21 | 2018-02-16 | 桂林航天工业学院 | 一种空间影像的高精度几何校正方法 |
CN108364264A (zh) * | 2018-02-07 | 2018-08-03 | 大连航天北斗科技有限公司 | 一种基于无人机红外遥感技术的海水温度监测方法及系统 |
CN108759903A (zh) * | 2018-04-02 | 2018-11-06 | 深圳万智联合科技有限公司 | 一种检测快速的电气设备故障检测系统 |
CN111611544B (zh) * | 2020-05-12 | 2024-03-26 | 中国科学院上海技术物理研究所 | 一种机载大视场面阵摆扫热像仪温排水监测方法 |
CN113109831A (zh) * | 2021-03-26 | 2021-07-13 | 国家电网有限公司 | 一种使用激光雷达对输电线路进行巡检的数据处理方法 |
CN115147239B (zh) * | 2022-06-10 | 2023-04-07 | 自然资源部第一海洋研究所 | 新建滨海电厂温升区环境本底水温推算及温升计算方法 |
CN115588127B (zh) * | 2022-10-18 | 2024-05-10 | 武汉际上导航科技有限公司 | 机载激光点云与多光谱影像融合的方法 |
CN115810078A (zh) * | 2022-11-22 | 2023-03-17 | 武汉际上导航科技有限公司 | 基于pos数据和机载可见光影像赋色激光点云的方法 |
-
2009
- 2009-04-30 CN CN2009100503122A patent/CN101551275B/zh active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113643386A (zh) * | 2021-10-13 | 2021-11-12 | 深圳市商汤科技有限公司 | 标定方法及装置、电子设备及计算机可读存储介质 |
US11830635B1 (en) * | 2022-09-06 | 2023-11-28 | First Institute of Oceanography, Ministry of Natural Resources | Method and system for determining background water temperature of thermal discharge from operating nuclear power plants based on remote sensing |
Also Published As
Publication number | Publication date |
---|---|
CN101551275A (zh) | 2009-10-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101551275B (zh) | 机载多光谱扫描仪用于工业温排水监测的技术方法 | |
CN103674063B (zh) | 一种光学遥感相机在轨几何定标方法 | |
CN107504981B (zh) | 一种基于激光测高数据的卫星姿态误差修正方法及设备 | |
CN103557841B (zh) | 一种提高多相机合成影像摄影测量精度的方法 | |
CN104897175B (zh) | 多相机光学推扫卫星在轨几何定标方法及系统 | |
CN105548976A (zh) | 船载雷达海上精度鉴定方法 | |
CN103129752B (zh) | 一种基于地面导航的光学遥感卫星姿态角误差动态补偿方法 | |
CN107144293A (zh) | 一种视频卫星面阵相机的几何定标方法 | |
CN102628942B (zh) | 一种雷达影像双视向信息补偿方法 | |
CN102346033B (zh) | 基于卫星观测角误差估计的直接定位方法及系统 | |
CN106127697A (zh) | 无人机机载成像高光谱几何校正方法 | |
CN106780321A (zh) | 一种cbers‑02卫星hr传感器影像整体严密定向与纠正拼接方法 | |
CN103413272A (zh) | 低空间分辨率多源遥感图像空间一致性校正方法 | |
CN102519433B (zh) | 一种利用rpc反演星载线阵传感器几何定标参数方法 | |
Poli | A rigorous model for spaceborne linear array sensors | |
CN104807477B (zh) | 一种基于靶标控制点的卫星ccd阵列影像几何检校方法 | |
CN106873004A (zh) | 基于太阳高度角自适应的高轨面阵相机在轨几何检校方法 | |
CN105444780A (zh) | 星载摆扫式光学相机图像定位验证系统及处理方法 | |
CN105444778A (zh) | 一种基于成像几何反演的星敏感器在轨定姿误差获取方法 | |
CN110793542A (zh) | 一种基于广义探元指向角的面阵光学遥感卫星在轨几何定标方法 | |
CN111524196B (zh) | 一种摆扫大幅宽光学卫星在轨几何标定方法 | |
CN114858133B (zh) | 一种恒星观测模式下姿态低频误差修正方法 | |
CN102346252B (zh) | 光学卫星遥感数据地理定位中的大气折射补偿方法及系统 | |
Tao et al. | On-orbit geometric calibration of the panchromatic/multispectral camera of the ZY-1 02C satellite based on public geographic data | |
Xu et al. | Innovations in the data processing algorithm for Chinese FY meteorological satellites |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |