技术领域本发明属于海洋监测领域，具体地说是一种适用于拉曼散射海水温度激光遥感测量系统的基于频数分布原理的滤波方法。
背景技术：
海水水下温度测量对研究海洋环境和气候监测及自然灾害的早期预报等十分重要。由于蓝绿激光在水中良好的透射性，用于海水水下温度的遥感测量可以快速获取其水平和垂直分布，具有浮标等定点站位测量和微波辐射或红外辐射等表面温度测量方式不可比拟的优点。目前国内外很多研究人员对拉曼散射海水温度激光遥感测量技术进行了广泛深入的研究，但是海水拉曼散射受波浪、外界光线、水下悬浮物质、CDOM荧光散射等因素影响信号强度波动大，导致水温测量误差较大。虽然现有根据海水拉曼光谱反演海水温度的方法中已经采用光谱强度归一化的方式降低拉曼光谱信号强度波动对水温测量结果的影响，如2015年“AppliedOptics”第54卷第10期发表的“UltimatesensingresolutionofwatertemperaturebyremoteRamanspectroscopy”和2001年“JournaloftheKoreanPhysicalSociety”第38卷第6期发表的“LIDARmeasurementofwatertemperaturebyusingfrequency-shiftedRamanscattering”。但是该方法尚不能很好地消除所述外界瞬态噪声干扰信号对水温测量精度的影响，严重限制了拉曼散射海水温度激光遥感测量技术的实际应用。迄今为止，尚未见针对拉曼光谱受瞬态噪声干扰信号强度波动的相关滤波算法报道。
技术实现要素：
本发明所要解决的技术问题是针对水的拉曼散射光谱受噪声干扰信号强度波动大的问题，提供一种方法简单、合理，实现方便，滤波效果明显，能够有效排除外界噪声信号干扰，提高拉曼散射海水温度激光遥感测量系统的实时测温精度的基于频数分布原理的滤波方法。本发明解决上述技术问题采用的技术方案是：一种基于频数分布原理的滤波方法，其特征是：其包括以下步骤：步骤(1)：连续采集多组拉曼光谱数据I1(λ)、I2(λ)……，分别求得各组光谱的峰值强度Iimax，把连续N个Iimax值看成一个队列，存入数组Ix；Ix(i)＝[I1maxI2max…INmax],i＝1,2,…N步骤(2)：将数组Ix划分为M个区间，根据数组Ix中数据的最大值Ixmax和最小值Ixmin来确定分区间距ΔIx；ΔIx=Ixmax-IxminM]]>第j个分区区间的中间值表示为：Ixm(j)＝Ixmin-0.5ΔIx+j·ΔIx,j＝1,2,…M步骤(3)：按照分区对数组Ix中的元素进行分组，若Ix(i)属于指定第j个区间，则函数γj[Ix(i)]值为1，否则函数γj[Ix(i)]值为0，即：统计各分区区间的频数Fj；Fj=Σi=1Nγj[Ix(i)],j=1,2...M]]>步骤(4)：根据上述频数分布统计结果求出频数最大值Fmax所处的区间编号K，并求得该区间对应的中间值Ixm(K)，作为第N+1组拉曼光谱筛选的参照值；步骤(5)：设定相邻两组拉曼光谱强度波动允许的最大偏差值为A，若第N+1组拉曼光谱数据的峰值强度满足Ixm(K)-A＜I(N+1)max＜Ixm(K)+A，则认为这一组拉曼光谱有效，进行光谱预处理后存入寄存器用于拉曼光谱叠加，否则认为这一组拉曼光谱无效，光谱数据舍弃，数组Ix的长度固定为N，每采集一组新的拉曼光谱数据并求得其对应的峰值强度Iimax后都放在队列尾部，并扔掉原来队列首部的一个峰值强度数据，保持这N个数据始终是最近的数据，当寄存器内有效拉曼光谱数据达到L组以后，将L组预处理后的拉曼光谱进行叠加。根据本发明，将L组拉曼光谱叠加后的结果用于光谱特征值提取，进而将光谱特征值用于物质结构及成分鉴别、材料特性分析等，例如可根据所述滤波叠加后的拉曼光谱提取光谱特征值，进而根据光谱特征值计算海水盐度或温度。本发明所述一种基于频数分布原理的滤波方法具有实现方便，滤波效果明显，实时性好的优点，可有效排除瞬态噪声信号的干扰，提高拉曼散射海水温度激光遥感测量系统的抗干扰性和水温测量精度。附图说明下面结合附图对本发明作进一步的描述。图1是拉曼光谱峰值强度的频数分布原理示意图。图2是设定相邻两组拉曼光谱强度波动允许的最大偏差值为500时，灯光长时间干扰基于频数分布原理的滤波方法测试结果。图3是设定相邻两组拉曼光谱强度波动允许的最大偏差值为500时，灯光短时间干扰基于频数分布原理的滤波方法测试结果。图4是设定相邻两组拉曼光谱强度波动允许的最大偏差值为200时，灯光长时间干扰基于频数分布原理的滤波方法测试结果。图5是设定相邻两组拉曼光谱强度波动允许的最大偏差值为200时，灯光短时间干扰基于频数分布原理的滤波方法测试结果。具体实施方式本发明一种基于频数分布原理的滤波方法，其包括以下具体步骤：步骤(1)：连续采集多组拉曼光谱数据I1(λ)、I2(λ)……，分别求得各组光谱的峰值强度Iimax，把连续N个Iimax值看成一个队列，存入数组Ix；Ix(i)＝[I1maxI2max…INmax],i＝1,2,…N(1)式(1)中Ix(i)是数组Ix的第i个元素。步骤(2)：将数组Ix划分为M个区间，根据数组Ix中数据的最大值Ixmax和最小值Ixmin来确定分区间距ΔIx；ΔIx=Ixmax-IxminM---(2)]]>式(2)中Ixmax是输入数组Ix中的最大强度值，Ixmin是数组Ix中的最小强度值，每个区间的中间值可以表示为：Ixm(j)＝Ixmin-0.5ΔIx+j·ΔIx,j＝1,2,…M(3)式(3)中Ixm(j)是第j个区间的中间值。步骤(3)：按照分区对数组Ix中的每个元素进行分组，若Ix(i)属于指定第j个区间，则函数γj[Ix(i)]值为1，否则函数γj[Ix(i)]值为0，即式(4)中γj[Ix(i)]反映数组Ix中第i个元素是否属于第j个分区的情况。统计各分区区间的频数Fj；Fj=Σi=1Nγj[Ix(i)],j=1,2...M---(5)]]>式(5)中Fj是数组Ix中每个元素在第j个分区区间的频数。频数分布统计结果如图1所示，直接反映拉曼光谱峰值强度的变化情况。步骤(4)根据所述频数分布统计结果求出频数最大值Fmax所处的区间编号K，并求得该区间对应的中间值Ixm(K)，作为第N+1组拉曼光谱筛选的参照值。步骤(5)设定相邻两组拉曼光谱强度波动允许的最大偏差值为A，若第N+1组拉曼光谱数据的峰值强度满足Ixm(K)-A＜I(N+1)max＜Ixm(K)+A，则认为这一组拉曼光谱有效，进行光谱预处理后存入寄存器用于拉曼光谱叠加，否则认为这一组拉曼光谱无效，光谱数据舍弃，数组Ix的长度固定为N，每采集一组新的拉曼光谱数据并求得其对应的峰值强度Iimax后都放在队列尾部，并扔掉原来队列首部的一个峰值强度数据，保持这N个数据始终是最近的数据，当寄存器内有效拉曼光谱数据达到L组以后，将L组预处理后的拉曼光谱进行叠加。根据本发明，将L组拉曼光谱叠加后的结果用于光谱特征值提取，进而将光谱特征值用于物质结构及成分鉴别、材料特性分析等，例如可根据所述滤波叠加后的拉曼光谱提取光谱特征值，进而根据光谱特征值计算海水盐度或温度。本发明所述一种基于频数分布原理的滤波方法具有实现方便，滤波效果明显，实时性好的优点，可有效排除瞬态噪声信号的干扰，提高拉曼散射海水温度激光遥感测量系统的抗干扰性和水温测量精度。为了验证所述一种基于频数分布原理的滤波方法的抗干扰能力，在实验室采用灯光扰动的方式进行了测试。测试中首先实时采集20组拉曼光谱，求得其对应的光强极大值Iimax，存入数组Ix。然后将Ix中20个元素分为4个区间，统计拉曼光强在各个区间的分布，以频数最高的一个区间K的中间值Ixm(K)作为下一组拉曼光谱筛选的参照值。当设定相邻两组拉曼光谱强度波动允许的最大偏差值A＝500时，若第21组拉曼光谱的峰值强度满足Ixm(K)-500＜I21max＜Ixm(K)+500，则认为这一组拉曼光谱有效，进行预处理后存入寄存器用于拉曼光谱叠加；否则认为这一组拉曼光谱无效，舍弃该组拉曼光谱。数组Ix的长度固定为20，每次采集到一组新的拉曼光谱并求得其对应的峰值强度Iimax后都放在队列尾部，并扔掉原来队列首部的一个峰值强度数据，保持这20个数据始终是最新的数据。最后，当寄存器内有效拉曼光谱数据达到5组以后，将5组预处理后的拉曼光谱进行叠加，计算光谱特征值，进而计算水温并实时显示。图2是设定相邻两组拉曼光谱强度波动允许的最大偏差值为500时，灯光长时间干扰基于频数分布原理的滤波方法测试结果，其中虚线是高精度温度传感器的测温结果，实线是拉曼散射海水温度激光遥感测量系统测温结果，t1处产生的尖峰是由于系统运行过程中打开室内照明灯8s，然后关闭，这一扰动导致拉曼光谱强度增大，从而改变了测温结果，t2处的产生的尖峰对应开灯时间为7s。图3是设定相邻两组拉曼光谱强度波动允许的最大偏差值为500时，灯光短时间干扰基于频数分布原理的滤波方法测试结果，其中虚线是高精度温度传感器的测温结果，实线是拉曼散射海水温度激光遥感测量系统测温结果，测量过程中多次开关照明灯，每次照明时间缩短到5s以内，由于频数分布滤波算法的滤波作用，干扰信号对温度测量结果没有影响。图4是当设定相邻两组拉曼光谱强度波动允许的最大偏差值A＝200时灯光长时间干扰的测试结果，其中虚线是温度传感器同步测量的结果，实线是拉曼散射海水温度激光遥感测量系统测温结果，t3处产生的尖峰对应照明时间为7s，而t4处产生的尖峰对应开灯时间为10s。由于灯光扰动时间长导致拉曼光谱强度增大，从而改变了测温结果。图5是当设定相邻两组拉曼光谱强度波动允许的最大偏差值A＝200时灯光短时间干扰的测试结果，其中虚线是温度传感器同步测量的结果，实线是拉曼散射海水温度激光遥感测量系统测温结果，测量过程中多次开关照明灯，每次照明时间缩短到5s以内，由于频数分布滤波算法的滤波作用，干扰信号对温度测量结果没有影响。在上述说明中参照着附图，通过具体实施方式适当且充分地对本发明进行了说明，但是应该认识到，对本领域技术人员来说，能够容易的对上述的实时方式进行变更和/或改良。因此对于本领域技术人员所实施的变更方式或改良方式，只要是不脱离权利要求书中记载的权利要求的权力范围的层次，则都应被解释为被包括在该权利要求的权力范围内。