一种扫描式大气探测激光雷达实时数据处理方法与流程

1.本发明属于大气遥感探测领域，具体涉及激光雷达大气颗粒物污染监测系 统的实时、稳定数据处理。


背景技术：

2.开放空间大气颗粒物污染监测是定位环境污染源头的重要手段，传统手段 通过布设多个原位探测传感器的方式监测大气颗粒物污染，成本高且难以精准 定位污染源。大气激光雷达发射一束脉冲激光，通过测量不同时间大气后向散 射信号，实现不同距离大气颗粒物污染监测，结合水平扫描机构，可以实现大 空间范围网格化的大气颗粒物污染遥感，单台扫描式激光雷达即可实现大区域 颗粒物污染监测，在同样探测效能下，布设密度远低于原位探测传感器。但是， 由于激光雷达的信号强度同时受大气散射与衰减效应的影响，导致实际数据处 理过程依赖一些假设条件，对数据处理算法的鲁棒性造成一定的影响。
3.大气探测激光雷达的回波信号写为：
[0004][0005]
c为系统常数，s为数据点和激光雷达安装位置的距离，β为大气后向散射 系数，α为大气消光系数，对于均匀大气，可以采用klett斜率法，直接对距离 平方修正的信号取对数后微分得到大气消光系数，计算方法如下式所示：
[0006][0007]
这种方法实现简单、算法稳定，但是对非均匀大气，会存在较大的偏差， 为此，将消光系数与后向散射系数的比值设定为常数，并且假设参考距离的消 光系数，可以由fernald方法根据以下公式计算消光系数：
[0008][0009]
α(s
ref
)为参考距离消光系数，采用这种方法计算消光系数，参考距离处消 光系数的选择对精度影响较大，往往需要人工介入选择参考位置，实时处理的 鲁棒性较低。


技术实现要素：

[0010]
为克服现有技术在实时性、鲁棒性方面的不足，本发明提供了一种扫描式 大气污染监测激光雷达实时数据处理方法，此方法不依赖于参考位置大气参数 的选择，可以实时输出扫描区域消光系数二维分布结果。
[0011]
本发明提出的实时数据处理流程具体步骤为：(1)获取扫描激光雷达不同 角度和距离原始信号；(2)将激光雷达测得的信号乘以距离的平方；(3)将数 据从极坐标系转换到直角坐标系；(4)数据均匀网格化；(5)构造大气透过率 修正信号函数；(6)利用最小二乘法拟合本底大气消光系数；(7)根据拟合得 到的本底大气消光系数计算不同距离大气透过率，并对信号进行大气透过率修 正；(8)利用本底大气消光系数对大气透过率修正后的信号进行归一化；(9) 输出不同坐标的大气消光系数。
[0012]
对于步骤(3)，极坐标和直角坐标系的转换规则为：
[0013]
x＝s
·
cosθ，y＝s
·
sinθ
[0014]
其中s为数据点和激光雷达安装位置的距离；θ为数据点的方位角，正北方 向为0度，正东方向为90度，正南方向为180度，正西方向为270度；x为数 据点和原点的东西方向距离，向东为正，向西为负；y为数据点和原点的南北 方向距离，向北为正，向南为负。
[0015]
对于步骤(4)，按照最大有效扫描区域设计均匀网格，网格x方向和y方 向间距分别为δx和δy，x和y方向扫描范围分别为[x
min
,x
max
]和[y
min
,y
max
]，网 格数分别为n
x
和n
y
，网格的中心坐标为：
[0016]
x
i
＝x
min
(i
‑
1)
·
δx,i＝1,2,3,
…
,n
x
[0017]
y
j
＝y
min
(j
‑
1)
·
δy,j＝1,2,3,
…
,n
y
[0018]
i和j为网格序号，步骤(4)计算落在网格矩形区域内的所有数据的平均值， 赋值为该网格点数据值：
[0019][0020]
n(x
i
,y
i
)为落在[x
i
,y
i
]网格内所有数据点的数量。
[0021]
对于步骤(5)，构造大气透过率修正信号函数的形式为：
[0022][0023]
其中未知数α为大气消光系数。
[0024]
对于步骤(6)，利用最小二乘法拟合本底大气衰减系数α0，残差函数构造 为：
[0025][0026]
取δ(α)最小时的值α0为大气本底消光系数，m为网格内的所有有效数据。
[0027]
对于步骤(7)，根据步骤(6)计算得到的本底大气消光系数，对信号进 行大气透过率修正，得到相对大气散射系数：
[0028][0029]
对于步骤(8)，根据步骤(7)得到的相对大气散射系数，再利用本底大 气消光系数进行归一化，即可得到大气消光系数分布：
[0030][0031]
<r
′
(x
i
,y
j
)>为r
′
(x
i
,y
j
)的平均值。
[0032]
本发明公开了一种扫描式大气污染监测激光雷达实时数据处理方法，将激 光雷达扫描数据转换到直角坐标系下的均匀网格中进行处理，可以实现稳定可 靠的大气消光系数二维分布探测。相比传统的需要假设参考距离大气参数的反 演方法，该方法不需要频繁的人工介入，可以实现稳定、鲁棒的数据输出，利 于扫描式激光雷达数据实时数据处理。
附图说明
图1为扫描式大气探测激光雷达实时数据处理方法流程图。
具体实施方式
[0033]
下面将利用典型的扫描式激光雷达对本发明所提出的数据处理方法进行 阐述。
[0034]
典型的扫描式大气污染监测激光雷达将大气探测激光雷达安装在转台上， 整个激光雷达绕垂直方向扫描轴旋转，实现区域大气监测，典型的指标包括： 激光器工作波长1064nm，数据采集系统采样率20msps，原始数据距离分辨率 7.5m，激光雷达有效探测范围0.1
‑
2km，扫描范围0～360
°
，单次扫描时间10 分钟。
[0035]
针对此系统，采用本发明提出算法的系统工作流程为：
[0036]
(1)激光器开始工作，发射脉冲激光，同步开始数据采集，扫描机构开 始按照既定速度和范围开始扫描；
[0037]
(2)激光雷达结束扫描，数据处理系统读取不同角度θ、不同距离s的原 始信号p(θ,s)，对原始信号进行预处理后，乘以距离的平方，得到距离平方修 正信号r(θ,s)＝p(θ,s)
·
s2；
[0038]
(3)将所有数据的坐标由极坐标(θ,s)变换为直角坐标(x,y)，其中θ为方 位角，s为距离，x和y分别为东西和南北向距离坐标原点距离；坐标变换采用 以下方法：
[0039]
x＝s
·
cosθ，y＝s
·
sinθ
[0040]
(4)设计均匀的网格，网格的边界按照扫描区域的有效范围来设计，网 格间距可以根据数据信噪比来定，网格间距越大，数据信噪比越高，本例中激 光雷达最大有效探测距离为2km，区域x方向边界x
min
和x
max
、区域y方向边 界y
min
和y
max
、网格间隔δx和δy、网格数n
x
和n
y
、网格中心坐标x
i
和y
j
分别 为(单位为米)：
[0041]
x
min
＝
‑
2000，x
max
＝2000
[0042]
y
min
＝
‑
2000，y
max
＝2000
[0043]
δx＝20，δy＝20
[0044]
n
x
＝200，n
y
＝200
[0045]
x
i
＝x
min
(i
‑
1/2)
×
δx，i＝1,2,
…
,200
[0046]
y
i
＝y
min
(j
‑
1/2)
×
δy，j＝1,2,
…
,200
[0047]
(5)对落在每个网格内的所有数据点计算平均值，得到网格化的距离平 方修正信号，计算方法如下：
[0048][0049]
其中n(x
i
,y
j
)为落在对角坐标为和的网 格内的所有数据点的数量，当n(x
i
,y
j
)等于0时，即网格内没有有效原始数据 点时，将对应的r(x
i
,y
j
)记为无效数据。
[0050]
(6)构造大气透过率修正函数：
[0051][0052]
其中α为大气消光系数。
[0053]
(7)按照以下方法构造残差函数δ(α)：
[0054][0055]
利用最小二乘法拟合本底大气消光系数，将δ(α)最小时的消光系数定为本 底大气消光系数α0；
[0056]
(8)按照以下方法计算大气透过率修正信号：
[0057][0058]
其中α0为最小二乘法拟合得到的本底大气消光系数；
[0059]
(9)按照以下方法对信号进行归一化：
[0060][0061]
<r
′
(x
i
,y
j
)>为所有大气透过率修正信号的平均；
[0062]
(10)输出消光系数二维分布α(x
i
,y
j
)。
[0063]
按照以上方法，可以每10分钟输出一副4000m
×
4000m范围内的消光系 数分布图。