专利名称：一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法
技术领域：
本发明涉及一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法，具体地说，基于溢油模拟实验(如附图1所示)，获取溢油油膜厚度变化标准光谱反射率曲线，利用一种改进的油膜厚度估算模型并结合基于溢油模拟实验构建的参数查找表，来计算卫星高光谱影像数据中溢油油膜厚度的方法(流程如附图2所示)。
背景技术：
海洋石油开采、加工、运输等诸多过程中产生的溢油污染，在风、浪、流等因素的影响下，对海水、海表大气、海底、海岸带等都产生诸多不良影响，严重损害了海洋生态环境，并产生了巨大的经济损失。如2010年美国墨西哥湾的深海油井溢油事件、2010年中国大连新港溢油事件、2011年蓬莱19-3油田溢油事件，都产生了巨大的环境、经济损失。海洋溢油污染事故发生后，有效的、及时的空间监测与评估技术，能确定海洋溢油污染的位置、圈定污染的面积、查明溢油污染形成的类型，将有助于海洋溢油污染后续清理工作的开展；于此同时，也将有助于海洋溢油污染造成的损失评估，为海洋生态环境损失索赔提供依据。遥感技术因其大面积同步观测、时效性、数据的综合性、可比性及经济性等特点，成为海洋溢油污染监测与评估的有效支撑。多光谱、热红外、微波雷达、激光荧光等诸多技术手段在海洋溢油污染遥感研究中具有不同的技术特点。目前，海面溢油面积监测尚不能完全满足应用部门的需求，还需要深入开展溢油量估算研究，准确评估溢油量是溢油污染遥感应用研究需关注的重点。海面油膜是海洋溢油污染的最直接形式，也是海洋溢油污染转化的重要承载体之一。海洋溢油油膜随着厚度的不同，其光学特性也发生了显著的变化，同时海水背景影响逐渐增强，因此海面油膜探测的难度也逐渐增大，成为海洋溢油污染监测中的一个难点。高光谱遥感技术是20世纪90年代后期才趋于成熟的新型光学遥感前沿技术，该技术可提取出相应像元(像素)大小地表物质的、精度如同实验室光谱仪测量得到的精细光谱曲线，从而进行直接识别地物类型，并对其组成成分含量进行定量计算、分析与空间分布制图。在准确探测油膜污染面积的基础上，开展油膜厚度的定量遥感估算，得出对溢油量的估算将具有重要的现实意义。
发明内容
本发明的目的是提供一套利用地物光谱仪测量的标准油膜厚度变化的高光谱数据，根据星载传感器的性能指标，建立油膜厚度估算模型及其参数查找表，并以此计算卫星高光谱遥感影像中海洋溢油油膜厚度及分布的方法。本发明的原理如下利用较新的高光谱遥感技术手段，针对海洋溢油污染油膜，获取溢油油品与海水样品，模拟并测量标准油膜厚度连续变化的高光谱反射率数据，获取该次溢油海域的高光谱卫星数据，根据卫星传感器的参数指标，将标准实验数据转换成卫星传感器的模拟数据，利用模拟数据结合油膜厚度计算模型去建立计算模型的参数查找表，卫星影像数据预处理后，通过归一化处理，消除不同水色背景差异的影响，利用掩膜技术，识别溢油污染区，再结合计算模型与参数查找表，计算海洋溢油污染区油膜厚度。本发明的技术方案主要包括以下步骤(I)首先利用地物光谱仪，设计油膜厚度连续变化模拟实验，获取标准油膜厚度变·化的光谱反射率数据(图3)，其中包含了油膜厚度指标、反射率大小和范围。根据标准油膜厚度变化光谱，分析溢油油膜厚度估算的可用波段和油膜厚度的有效探测范围；海面溢油油膜厚度估计依靠对背景反射率的大小，考虑水体的光学辐射特性，因此选定油膜厚度参数表与计算模型的可用光谱区间为380nm-760nm,油膜的最大探测厚度为54 μ m。(2)根据卫星传感器的性能指标与参数，进行模拟数据的光谱重采样，地物光谱仪所测量的可用光谱数据范围为380um-760nm，光谱分辨率为Inm ;根据星载Hyperion高光谱传感器的性能指标(表I)，确定可用探测范围为426. 8nm-762. 6nm，光谱分辨率约为10nm，将实验模拟数据转换成星载传感器的模拟光谱数据(图4)。(3)模拟光谱数据的归一化处理，消除水色背景差异，得到归一化反射率值NR，NR=(Rt5il-RminZ(Rs--Rmin)，其中Rt5il为某一个波段油膜反射率值、Rmin为该波段最小油膜反射率值、Rsrawate为该波段海水背景反射率值，归一化反射率值(NR)的变化区间为
(图5)。(4)油膜厚度估算模型，根据Beer定律，光强随着穿透介质长度的增加而逐步的衰减，因此油膜厚度光谱反射率与油膜的厚度具有负指数变化的关系，根据这一理论基础，可以建立油膜厚度与油膜归一化反射率的关系为NR = e(_a<M)，贝帳据步骤3中的归一化反射率公式，可以将该模型描述为(Rt5il-RminV(Rseawate-Rmin) = e(_aW)，d为油膜厚度，Rtjil为所测量的标准油膜反射光谱，则可以根据这些参数估算各波段衰减系数a。(5)海洋溢油高光谱卫星数据获取与处理，获取海洋溢油事故发生地的卫星高光谱影像数据，该数据经过大气校正、辐射校正等预处理后，得到卫星遥感反射率数据；对该数据进行取样分析，此过程为人工过程，需要均匀选择清洁海水样点、海洋溢油样点与海面其它目标样点(图6a)，分析这些样点的光谱特征(图6b)，可以发现反射率的大小排列为海水反射率>溢油油膜反射率，对这些目标取平均，并结合影像统计，获取平均海水背景反射率值、油膜最低反射率值等信息(图6c);平均反射率曲线也有助于最优波段的选择，根据海水背景的反射率大小，选择背景海水反射率最高值所在波段为最优监测波段。(6)对该影像反射率数据进行归一化处理，处理方法为NR= (Roil-Rfflin)/(Rseawate-Rmin)，考虑到卫星数据量较大，且只需要一个波段进行计算即可，因此选取其中一个最佳波段，其中心波长为548. 9nm，进行影像数据反射率归一化处理，其中Rtjil为该波段油膜反射率值、Rmin为最小油膜反射率值、Rseawate为该波段海水背景反射率值，这些数据可采用图像统计的方法获得，计算得到该波段的归一化反射率影像(图7)。(7)溢油区掩膜处理与油膜厚度计算，经过归一化处理后的反射率影像，其归一化值的范围为
，其中D > I为非溢油区，因此先对该区做掩膜处理，保留溢油区信息，利用D = (-l/a548.9)*LN(NR548.9)公式计算溢油区油膜厚度，其中NR548.9为掩膜处理后的卫星影像数据548. 9nm波段的归一化反射率，a548.9为步骤3中计算得到的548. 9nm波段衰减系数，可以计算得到该次溢油事故中溢油油膜厚度及其分布情况(图8)。与现有技术相比，本发明优势体现在高光谱遥感技术为海面溢油油膜厚度估算提供了有效的技术手段，与传统监测手段相比，该技术方法无需进入溢油污染现场，只需要在室内模拟测量溢油油膜厚度连续变化光谱，根据获取溢油信息的遥感传感器性能参数进行转换，建立相应的参数查找表，能客观高效的估算海洋溢油油膜厚度及其范围分布，该方法能满足海洋溢油应急监测的需求，提高溢油量应急监测与估算的准确度。具体有益效果如下本发明利用地物光谱仪，模拟并获取了标准油膜厚度连续变化过程中的精细光谱反射率数据，基于溢油油膜厚度计算模型建立了相应的参数查找表，明确了利用该方法进行卫星影像中溢油油膜厚度估算的计算流程，利用该方法成功计算了海面溢油的油膜厚度。这种基于参数查找表计算卫星影像中油膜厚度的方法，能为今后海洋卫星溢油污染应急监测提供应用基础与相关信息，能降低溢油区监测的时间、人力和物力成本，保证监测精度，提高实时监测效率。实际应用也表明，本发明提供了一套有效的方法，对今后的溢油污染监测与评估，只需要通过该方法建立一定的参数数据库，针对不同的卫星遥感影像的性能指标，均可以同过该发明进行参数的转换，并进一步用于海洋溢油油膜厚度的定量估算。该方法克服了传统方法需花费大量时间和大量人力物力，且降低了对溢油污染区的难进入的实际困难。该方法也能进一步增强高光谱遥感技术在海洋溢油油膜污染监测应用的普及性和有效性，可以更好地服务于海洋环境监测、海洋溢油污染评估、海洋溢油量估算与溢油污染定损索赔等海洋应用行业。
四
图1为本发明的模拟实验2为本发明的流程示意图；图3为实验室测量的标准油膜厚度变化光谱反射率曲线；图4为模拟卫星传感器获取的油膜厚度变化光谱反射率曲线；图5为归一化反射率模拟光谱曲线；图6为溢油区典型目标及其影像光谱a.溢油区典型目标及其空间位置；b.典型目标的采样光谱曲线；c.典型目标的平均光谱曲线；图7为548. 9nm波段的归一化反射率影像；图8为海洋溢油油膜厚度估算；
五具体实施例方式以下通过实例对本发明作进一步解释以2007年中国辽东湾海域一次典型溢油污染监测为例，发生溢油事件后，获取了溢油发生区域的海水与油品样品。使用ASD地物光谱仪(具体参数如表I所示)，利用2盏50瓦卤素灯作为入射光源，入射光角度与上午10时的太阳高度角相近，在实验室内获取海水背景下不同油膜厚度的高光谱反射率数据(如图1)，该数据的光谱范围为350-2500nm，光谱分辨率为lnm，选取可见光波段380-760nm作为原始数据(如图3)。表IASD光谱仪的技术参数
权利要求
1.一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法，至少基于美国ASD 公司生产的FieldSpec 3Hr便携式地物光谱仪I台、50瓦卤素灯2盏和25°视场的探测镜头I个等设备开展实验模拟，其特征在于包括以下步骤(1)首先利用地物光谱仪，通过油膜厚度连续变化模拟实验，获取标准油膜厚度变化的光谱反射率数据，其中包含了油膜厚度指标、反射率大小和范围等参数，根据标准油膜厚度变化光谱，分析溢油油膜厚度估算的可用波段和油膜厚度的有效探测范围；海面溢油油膜厚度估计依靠对背景反射率的大小，考虑水体的光学辐射特性，选定油膜厚度参数表与计算模型的可用光谱区间为380nm-760nm,油膜的最大探测厚度为54 μ m。(2)根据卫星传感器的性能指标与参数，进行模拟数据的光谱重采样；地物光谱仪所测量的可用光谱数据范围为380nm-760nm，光谱分辨率为Inm ;根据星载高光谱传感器的性能指标，可用探测范围为426. 8nm-762. 6nm，光谱分辨率约为10nm，将实验模拟数据转换成星载传感器的模拟光谱数据。(3)模拟光谱数据的归一化处理，消除水色背景差异，得到归一化反射率值NR= (Roil-RminV(Rs--Rmin)，其中Rt5il为某一个波段油膜反射率值、Rmin为可探测最厚油膜该波段最小油膜反射率值、RSMWatOT为该波段海水背景反射率值，归一化反射率值的变化区间为O 到I。(4)油膜厚度估算模型，根据Beer定律，光强随着穿透介质长度的增加而逐步的衰减， 因此油膜厚度光谱反射率与油膜的厚度具有负指数变化的关系，根据这一理论基础，可以建立油膜厚度与油膜归一化反射率的关系为NR = e(_a<w)，则根据步骤3中的归一化反射率公式，可以将该模型描述为(Rt5il-RminVOUrate-Rmin) = e(_aW)，d为油膜厚度，Rtjil为所测量的标准油膜反射光谱，可以根据这些参数估算各波段衰减系数a。(5)海洋溢油高光谱卫星数据获取与处理，获取某次海洋溢油事故发生地的卫星高光谱影像数据，该数据经过大气校正、辐射校正等预处理后，得到卫星遥感反射率数据；对该数据进行取样分析，此过程为人工过程，需要均匀选择清洁海水样点、海洋溢油样点与海面其它目标样点，分析这些样点的光谱特征，可以发现反射率的大小排列为海水反射率>溢油油膜反射率，对这些目标取平均，并结合影像统计，获取平均海水背景反射率值RSMWate、 油膜最低反射率值等信息Rmin ;根据海水背景的反射率大小，选择背景海水反射率最高值所在波段为最优监测波段。(6)对该影像反射率数据进行归一化处理，处理方法为NR= (Roil-Rmin) / (Rseawater-Rmin)， 考虑到卫星数据量较大，只需要一个波段进行计算即可，因此选取其中一个最佳波段进行影像数据反射率归一化处理，其中Roil为该波段油膜反射率值、Rfflin为最小油膜反射率值、 Rseawater为该波段海水背景反射率值，这些数据可采用图像统计的方法获得，计算得到该波段的归一化反射率影像。(7)溢油区掩膜处理与油膜厚度计算，经过归一化处理后的反射率影像，其归一化值的范围为从O到D，其中D > 1，因此先对归一化反射率影像中归一化反射率值小于I的区域做掩膜处理，保留溢油区信息，利用D = (_l/a548.9)*LN(NR548.9)公式计算溢油区油膜厚度， 其中殿548.9为掩膜处理后的卫星影像数据548. 9nm波段的归一化反射率，a548.9为步骤3中计算得到的548. 9nm波段衰减系数，可以计算得到该次溢油事故中溢油油膜厚度及其分布情况。
2.根据权利要求1所述的一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法，其特征在步骤1，所述油膜光谱反射率数据是有地物光谱仪获取的油膜厚度连续变化过程中的光谱曲线，其中包括了油膜厚度指标，光谱分辨率指标为lnm，连续光谱指标 350nm-2500nm，对获取的光谱反射率曲线进行平滑处理，作为建立参数查找表的基础信息源。
3.根据权利要求1所述的一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法，其特征在步骤2，对实验室获取的标准反射率数据进行光谱重采样，根据需要计算的卫星光学传感器性能指标、探测波段选择标准油膜反射率光谱的波段范围，根据卫星传感器的波段宽度、中心波长、半波高宽将标准油膜反射率光谱进行转换，将标准油膜反射率光谱数据模拟成卫星传感器所获取的光谱数据。
4.根据权利要求3所述的一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法，其特征在步骤3，对于获取的基础数据，将油膜光谱反射率曲线进行归一化处理，将油膜反射率光谱大小转换到从O到I的数值区间，消除实验背景的差异，使归一化油膜反射率光谱可以应用到不同海洋环境背景、不同光学卫星影像数据中去。
5.根据权利要求1-4之一所述的一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法，其特征在步骤4，基于Beer定律构建的油膜厚度光学模型能用来描述透明或半透明介质的光衰减行为，通过计算衰减系数作为油膜厚度计算模型的参数查找表，查找表数据用于所模拟卫星光学传感器获取的影像中油膜厚度计算。
6.根据权利要求5所述的一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法，其特征在步骤5，卫星传感器获取的海洋溢油区影像，需要经过数据预处理，将其转换反射率影像，在影像上手动提取溢油油膜和背景海水的反射率光谱曲线，确认海水背景反射率 > 油膜反射率，对海水背景反射率光谱求平均，确认背景海水最高反射峰所在位置，以此作为计算的最优波段；统计最优波段海水背景反射率值和最低油膜反射率值，作为影像数据归一化处理的基础数据。
7.根据权利要求6所述的一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法，其特征在步骤6，以步骤5中获取的最优波段背景海水反射率值和最低油膜反射率值，结合步骤4中计算归一化油膜反射率公式，计算最优卫星遥感影像数据，获得归一化影像数据，归一化遥感影像反射率值< 1，可以视为影像中油膜覆盖区；归一化遥感影像反射率值> I，可以视为闻反射率目标物，为闻悬浮物海水、航迹线水体与船等。
8.根据权利要求7所述的一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法，其特征在步骤7，对归一化卫星影像数据进行掩膜处理，保留归一化反射率值小于I 的海域，结合步骤4中建立的参数查找表与油膜厚度计算模型，对最优波段进行计算，获取该海域溢油油膜厚度及其分布数据。
全文摘要
本发明涉及一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法，属于海洋环境监测的研究领域。其具体步骤为通过设计模拟实验，获取标准油膜厚度连续变化高光谱数据并进行预处理；海洋溢油卫星遥感数据的获取与预处理；基于卫星传感器性能指标的标准油膜光谱重采样；光谱反射率归一化；油膜厚度光学计算模型；衰减系数参数查找表；卫星数据的反射率归一化与掩膜处理；卫星数据的最优波段及其相关参数查询建立；利用参数查找表及光学模型计算卫星数据中的油膜厚度。本发明与传统观测手段相比，野外工作量小，能量化溢油油膜的厚度，无需直接进入海洋溢油污染区就能计算溢油区油膜的厚度，能满足海洋溢油应急监测与评估的需求。
文档编号G01B11/06GK102997856SQ20121053183
公开日2013年3月27日 申请日期2012年12月12日 优先权日2012年12月12日
发明者陆应诚, 李想 申请人:南京大学