CN101976353A - 低密度人群的统计方法及装置 - Google Patents

低密度人群的统计方法及装置 Download PDF

Info

Publication number
CN101976353A
CN101976353A CN 201010529670 CN201010529670A CN101976353A CN 101976353 A CN101976353 A CN 101976353A CN 201010529670 CN201010529670 CN 201010529670 CN 201010529670 A CN201010529670 A CN 201010529670A CN 101976353 A CN101976353 A CN 101976353A
Authority
CN
China
Prior art keywords
point
class
classification
image
foreground
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.)
Granted
Application number
CN 201010529670
Other languages
English (en)
Other versions
CN101976353B (zh
Inventor
游磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Netposa Technologies Ltd
Original Assignee
Beijing Zanb Science & Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Zanb Science & Technology Co Ltd filed Critical Beijing Zanb Science & Technology Co Ltd
Priority to CN2010105296704A priority Critical patent/CN101976353B/zh
Publication of CN101976353A publication Critical patent/CN101976353A/zh
Application granted granted Critical
Publication of CN101976353B publication Critical patent/CN101976353B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Analysis (AREA)

Abstract

本发明提供了一种低密度人群的统计方法及装置。所述方法包括：更新背景图像并提取当前图像的前景点；对当前图像中的每个前景点进行检测以获取特征点；获取每个特征点的运动轨迹；根据运动轨迹的运动位移和运动方向变化率，提取出正常的运动轨迹；将连续多帧图像内同一轨迹上的特征点作为一个分类点，对所有分类点进行聚类，以获得初步的分类；和对初步的分类进行类分裂和类聚合处理，以获得聚类类别数，将该聚类类别数作为所统计人群的人数。本发明能够有效地克服由于遮掩等使得跟踪丢失或者错误的问题，简单易行并且提高了准确度；此外，可扩展性强，不仅仅适用于人群检测，稍加调整，可应用于车流量等类似其他目标数量统计。

Description

低密度人群的统计方法及装置
技术领域
本发明涉及图像处理、视频监控，特别涉及低密度人群的统计方法及装置。
背景技术
低密度人数统计是智能监控的一个典型应用，在不同场合有不同作用。在百货大楼、超市等商业地区，通过统计具体人数，可侧面了解该商业地区效益情况，为该地区的发展规划做出侧面的参考；在飞机场、汽车站等公共交通地区，可为交通设施的优化配置提供重要依据，排除可能发生的安全事故；在其他如重要景点、军事要地等需要限制人数的地点提供预警帮助。因此，低密度人数统计有着广泛的应用范围和重要的应用价值。
公开号为CN101325690A的中国专利申请介绍了一种监控视频流中人流分析与人群聚集过程的检测方法和系统，该系统根据对多运动目标的跟踪结果统计人群。但上述方法对于存在人群遮挡的场景效果不甚理想。
公开号为CN101188743A的中国专利申请公开了一种基于视频的智能数人系统及其处理方法，该系统和方法包括基于前景块检测的人员数目估计。但是该系统和方法还需要依赖于基于头部检测的人员数目估计，较为复杂。
综上所述，目前迫切需要提出一种能克服人群遮挡问题并且简单易行的低密度人群的统计方法及装置。
发明内容
有鉴于此，本发明的主要目的在于克服人群遮挡问题，简单易行并准确地统计出低密度人群的人数。
为达到上述目的，按照本发明的第一个方面，提供了一种低密度人群的统计方法，该方法包括：更新背景图像并提取当前图像的前景点；对当前图像中的每个前景点进行检测以获取特征点；获取每个特征点的运动轨迹；根据运动轨迹的运动位移和运动方向变化率，提取出正常的运动轨迹；将连续多帧图像内同一轨迹上的特征点作为一个分类点，对所有分类点进行聚类，以获得初步的分类；和对初步的分类进行类分裂和类聚合处理，以获得聚类类别数，将该聚类类别数作为所统计人群的人数。
优选地，所述更新背景图像包括：假设Ik表示第k帧图像，k为整数，Bk表示第k帧背景图像，其中背景图像的初始值为B0＝I0，则按下式更新背景图像：
B k ( x , y ) = B k - 1 ( x , y ) - 1 , B k - 1 ( x , y ) > I k ( x , y ) B k - 1 ( x , y ) + 1 , B k - 1 ( x , y ) < I k ( x , y ) B k - 1 ( x , y ) B k - 1 ( x , y ) = I k ( x , y )
其中，x、y分别表示像素点的横坐标和纵坐标。
优选地，所述提取当前图像的前景点包括：第k帧图像的前景图像Fk为：Fk(x，y)＝|Ik(x，y)-Bk(x，y)|；对前景图像进行阈值分割以提取前景点。
优选地，所述对前景图像进行阈值分割以提取前景点包括：如果Fk(x，y)≥第一阈值T1，则认为该点(x，y)是当前图像的前景点，否则认为该点(x，y)是当前图像的背景点。
优选地，所述对当前图像中的每个前景点进行检测以获取特征点是利用角点检测算法来进行的。
优选地，所述获取每个特征点的运动轨迹是采用KLT特征点跟踪算法来进行的。
优选地，所述根据运动轨迹的运动位移和运动方向变化率，提取出正常的运动轨迹包括：分别计算每条运动轨迹在连续多帧图像内的运动位移、运动方向变化率；将同时满足运动位移＞第四阈值T4和运动方向变化率＜第五阈值T5的运动轨迹作为正常的轨迹并输出该正常的轨迹。
优选地，所述聚类包括：(1)将所有分类点看作为一个初始类C1，此时当前的类别数Max＝1；(2)在所有类中选出直径最大且满足直径≥第六阈值T6的类Cm，其中m表示第m类，且1≤m≤Max，Max表示当前的类别数；(3)在直径最大的类Cm中找出一个与其他分类点相异度最大的分类点，执行Max＝Max+1并将该分类点看作一个新的类CMax；(4)分别计算直径最大的类Cm中分类点与新的类CMax中分类点的最大距离outerDist，其中
Figure BSA00000330686800032
表示直径最大的类Cm中第p个分类点，并计算直径最大的类Cm中分类点
Figure BSA00000330686800033
与该类中其他分类点的最小距离interDist，如果outerDist＜interDist则将该分类点
Figure BSA00000330686800034
放入新的类CMax中；(5)如果直径最大的类Cm中没有分类点放入新的类CMax中，则执行下一步骤，否则返回到步骤(2)；(6)如果所有类的直径＜第六阈值T6，则结束，否则返回到步骤(2)。
优选地，所述类分裂包括：在每一个类中，计算每条轨迹在当前帧图像之后的多帧图像内的位移的平均方向，根据计算的平均方向所属的八个方向类别，将每个类分裂成更小的类，其中八个方向分别为：-22.5°～22.5°、22.5°～67.5°、67.5°～112.5°、112.5°～157.5°、157.5°～202.5°、202.5°～247.5°、247.5°～292.5°、292.5°～337.5°。
优选地，所述类聚合包括：计算每个类的中心位置，如果两个类满足以下两个条件：①两个类的中心位置＜第八阈值T8；②两个类同属于八个方向中的同一方向类中，则将两个类聚合成一个类，直到没有两个类满足条件为止。
按照本发明的另一个方面，提供了一种低密度人群的统计装置，该装置包括：背景建立及前景点提取单元，用于更新背景图像并提取当前图像的前景点；特征点检测单元，用于从当前图像的前景点中检测出特征点；特征点运动轨迹获取单元，用于获取特征点运动轨迹；正常运动轨迹提取单元，用于根据运动轨迹的运动位移和运动方向变化率，提取出正常的运动轨迹；运动轨迹初步分类单元，用于将连续多帧图像内同一轨迹上的特征点作为一个分类点，对所有分类点进行聚类，以获得初步的分类；和人数获取单元，用于对初步的分类进行类分裂和类聚合处理，以获得聚类类别数，将该聚类类别数作为所统计人群的人数。
优选地，所述背景建立及前景点提取单元执行以下操作：假设Ik表示第k帧图像，k为整数，Bk表示第k帧背景图像，其中背景图像的初始值为B0＝I0，则按下式更新背景图像：
B k ( x , y ) = B k - 1 ( x , y ) - 1 , B k - 1 ( x , y ) > I k ( x , y ) B k - 1 ( x , y ) + 1 , B k - 1 ( x , y ) < I k ( x , y ) B k - 1 ( x , y ) B k - 1 ( x , y ) = I k ( x , y )
其中，x、y分别表示像素点的横坐标和纵坐标。
优选地，背景建立及前景点提取单元此外还执行以下操作：第k帧图像的前景图像Fk为：Fk(x，y)＝|Ik(x，y)-Bk(x，y)|；对前景图像进行阈值分割以提取前景点。
优选地，所述对前景图像进行阈值分割以提取前景点包括：如果Fk(x，y)≥第一阈值T1，则认为该点(x，y)是当前图像的前景点，否则认为该点(x，y)是当前图像的背景点。
优选地，所述特征点检测单元利用角点检测算法来进行检测。
优选地，所述特征点运动轨迹获取单元采用KLT特征点跟踪算法来获取运动轨迹。
优选地，所述正常运动轨迹提取单元执行以下操作：分别计算每条运动轨迹在连续多帧图像内的运动位移、运动方向变化率；将同时满足运动位移＞第四阈值T4和运动方向变化率＜第五阈值T5的运动轨迹作为正常的轨迹并输出。
优选地，所述运动轨迹初步分类单元执行以下操作：将所有分类点看作为一个初始类C1，此时当前的类别数Max＝1；在所有类中选出直径最大且满足直径≥第六阈值T6的类Cm，其中m表示第m类，且1≤m≤Max，Max表示当前的类别数，在直径最大的类Cm中找出一个与其他分类点相异度最大的分类点，执行Max＝Max+1并将该分类点看作一个新的类CMax，分别计算直径最大的类Cm中分类点
Figure BSA00000330686800051
与新的类CMax中分类点的最大距离outerDist，其中
Figure BSA00000330686800052
表示直径最大的类Cm中第p个分类点，并计算直径最大的类Cm中分类点
Figure BSA00000330686800053
与该类中其他分类点的最小距离interDist，如果outerDist＜interDist则将该分类点
Figure BSA00000330686800054
放入新的类CMax中，直到直径最大的类Cm中没有分类点放入新的类CMax中并且所有类的直径＜第六阈值T6。
优选地，所述人数获取单元执行的类分裂包括：在每一个类中，计算每条轨迹在当前帧图像之后的多帧图像内的位移的平均方向，根据计算的平均方向所属的八个方向类别，将每个类分裂成更小的类，其中八个方向分别为：-22.5°～22.5°、22.5°～67.5°、67.5°～112.5°、112.5°～157.5°、157.5°～202.5°、202.5°～247.5°、247.5°～292.5°、292.5°～337.5°。
优选地，所述人数获取单元执行的类聚合处理包括：计算每个类的中心位置，如果两个类满足以下两个条件：①两个类的中心位置＜第八阈值T8；②两个类同属于八个方向中的同一方向类中，则将两个类聚合成一个类，直到没有两个类满足条件为止。
与现有的技术相比，本发明的低密度人群的统计方法及装置，有效地克服由于遮掩等使得跟踪丢失或者错误的问题，简单易行并提高了准确度；此外，可扩展性强，不仅仅适用于人群检测，稍加调整，可应用于车流量等类似其他目标数量统计。
附图说明
图1示出了按照本发明的低密度人群的统计方法的流程图；
图2示出了本发明的角点检测算法的16个周边点的图；和
图3示出了按照本发明的低密度人群的统计装置的框架图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明进一步详细说明。
图1表示按照本发明的低密度人群的统计方法的流程图。如图1所示，按照本发明的低密度人群的统计方法包括：
第一步骤101，更新背景图像并提取当前图像的前景点；
第二步骤102，从当前图像的前景点中检测出特征点；
第三步骤103，获取特征点运动轨迹；
第四步骤104，根据运动轨迹的运动位移和运动方向变化率，提取出正常的运动轨迹；
第五步骤105，将连续多帧图像内同一轨迹上的特征点作为一个分类点，对所有分类点进行聚类，以获得初步的分类；
第六步骤106，对初步的分类进行类分裂和类聚合处理，以获得聚类类别数，将该聚类类别数作为所统计人群的人数。
第一步骤：
假设Ik表示第k(k为整数)帧图像，Bk表示第k帧背景图像(其中背景图像的初始值为B0＝I0)，按下式更新背景图像：
B k ( x , y ) = B k - 1 ( x , y ) - 1 , B k - 1 ( x , y ) > I k ( x , y ) B k - 1 ( x , y ) + 1 , B k - 1 ( x , y ) < I k ( x , y ) B k - 1 ( x , y ) B k - 1 ( x , y ) = I k ( x , y )
其中，x、y分别表示像素点的横坐标和纵坐标。相应地，第k帧图像的前景图像Fk为：Fk(x，y)＝|Ik(x，y)-Bk(x，y)|。对前景图像进行阈值分割以提取前景点，即如果Fk(x，y)≥第一阈值T1，则认为该点是当前图像的前景点，否则认为该点是当前图像的背景点。其中，T1∈[10，30]且T1为整数。
第二步骤：
利用角点检测算法，对当前图像中的每个前景点进行检测以获取特征点。图2示出了按照本发明的角点检测算法的16个周边点的图。如图2所示，对每个前景点p都有与之相应的16个周边点(即图2中的点1～16)。设与前景点p对应的16个周边点为pn(n＝1，…，16)，分别计算当前图像内每个前景点p的像素值Ik(p)(其中k表示当前图像为第k帧图像)与其16个周边点pn的像素值Ik(pn)的差值Ik(p)-Ik(pn)，统计该差值Ik(p)-Ik(pn)＞第二阈值T2的连续周边点的个数num1、该差值Ik(p)-Ik(pn)＜-第二阈值T2的连续周边点的个数num2，如果该个数num1＞第三阈值T3或者num2＞第三阈值T3，则认为该前景点为特征点。其中，T2∈[10，30]，T3∈[9，12]，且T2、T3均为整数。角点检测算法可以参考文献“E.Rosten and T.Drummond.Machine learning for high-speed corner detection.European conference on computer vision”。
第三步骤：
采用KLT特征点跟踪算法，获取每个特征点的运动轨迹。设当前图像(例如第k帧图像)内特征点(x，y)的像素值为Ik(x，y)，该特征点(x，y)在x、y方向的梯度分别为
Figure BSA00000330686800081
Figure BSA00000330686800082
初始化D＝[dx，dy]T＝[0，0]T，那么下一帧图像(即第k+1帧图像)中该特征点(x，y)在x、y方向的位移公式如下：
ZD＝e
其中：
Z = &Sigma; i = - 1 1 &Sigma; j = - 1 1 I k x ( x + j , y + i ) * I k x ( x + j , y + i ) * g I k x ( x + j , y + i ) * I k y ( x + j , y + i ) * g I k x ( x + j , y + i ) * I k y ( x + j , y + i ) * g I k y ( x + j , y + i ) * I k y ( x + j , y + i ) * g
e = &Sigma; i = - 1 1 &Sigma; j = - 1 1 I k x ( x + j , y + i ) * I k + 1 ( x + j + dx , y + i + dy ) * g I k y ( x + j , y + i ) * I k + 1 ( x + j + dx , y + i + dy ) * g
g是3×3高斯滤波系数。
解方程得到D，再带入上述公式循环计算，直到结束；
那么该特征点(x，y)在下一帧图像中的位移向量为：
D＝[dx，dy]T
第四步骤：
所述根据运动轨迹的运动位移和运动方向变化率，提取出正常的运动轨迹的方法如下：分别计算每条运动轨迹在连续多帧(此处的连续多帧为至少连续5帧，帧数优选地处于[5，20]的范围之内)图像内的运动位移、运动方向变化率；将同时满足运动位移＞第四阈值T4(T4∈[5，10]且为整数)并且运动方向变化率第五阈值T5(T5∈[30，45]且为整数)的运动轨迹作为正常的轨迹并输出。
第五步骤：
将连续多帧(如前面所提到的，此处的多帧为至少连续5帧，帧数优选地处于[5，20]的范围之内)的图像内同一轨迹上的每个特征点作为一个分类点，对这些分类点进行聚类，以获得初步的分类。按照优选实施方式，可以采用分裂聚类的方法，分裂聚类的方法是：
(1)将所有分类点看作为一个初始类C1，此时当前的类别数Max＝1；；
(2)在所有类中选出直径最大且满足直径≥第六阈值T6的类Cm(m表示第m类，且1≤m≤Max，Max表示当前的类别数)；
(3)在直径最大的类Cm中找出一个与其他分类点相异度最大的分类点，执行Max＝Max+1并将该分类点看作一个新的类CMax；
(4)分别计算直径最大的类Cm中分类点
Figure BSA00000330686800091
(
Figure BSA00000330686800092
表示直径最大的类Cm中第p个分类点)与新的类CMax中分类点的最大距离，计算直径最大的类Cm中分类点
Figure BSA00000330686800093
与该类中其他分类点的最小距离，如果outerDist＜interDist，则将该分类点
Figure BSA00000330686800094
放入新的类CMax中；
(5)如果直径最大的类Cm中没有分类点放入新的类CMax中，则执行下一步骤，否则返回到步骤(2)；
(6)如果所有类的直径＜第六阈值T6，则结束，否则返回到步骤(2)。
其中，T6为整数，根据实际场景确定，一般可以设为人的高度的一半。
对分裂聚类后获取的每个类别进行判别，若类别中的分类点的个数＜第七阈值T7，则认为该类别是干扰，删除该类别并执行Max＝Max-1。T7∈[T6/6，T6/2]且为整数，例如可以选T7＝T6/3。
第六步骤：
对第五步骤输出的初步的分类进行类分裂和类聚合处理，以获得聚类类别数，将该聚类类别数作为所统计人群的人数。按照优选实施方式，可以采用以下的类分裂方法：在每一个类中，计算每条轨迹在当前帧图像之后的多帧(如前所述，此处的多帧为至少连续5帧，帧数优选地处于[5，20]的范围之内)图像内的位移的平均方向，根据计算的平均方向所属的八个方向类别，将每个类分裂成更小的类。其中八个方向分别为：-22.5°～22.5°、22.5°～67.5°、67.5°～112.5°、112.5°～157.5°、157.5°～202.5°、202.5°～247.5°、247.5°～292.5°、292.5°～337.5°。按照优选实施方式，可以采用以下的类聚合方法：计算每个类的中心位置，如果两个类满足以下两个条件：①两个类的中心位置＜第八阈值T8(T8∈[15，20]且为整数)；②两个类同属于八个方向中的同一方向类中，则将两个类聚合成一个类，直到没有两个类满足条件为止。
本发明还提供了一种低密度人群的统计装置。图3示出了按照本发明的低密度人群的统计装置的框架图。如图3所示，低密度人群的统计装置包括：
背景建立及前景点提取单元1，用于更新背景图像并提取当前图像的前景点；
特征点检测单元2，用于从当前图像的前景点中检测出特征点；
特征点运动轨迹获取单元3，用于获取特征点运动轨迹；
正常运动轨迹提取单元4，用于根据运动轨迹的运动位移和运动方向变化率，提取出正常的运动轨迹；
运动轨迹初步分类单元5，用于将连续多帧图像内同一轨迹上的特征点作为一个分类点，对所有分类点进行聚类，以获得初步的分类；
人数获取单元6，用于对初步的分类进行类分裂和类聚合处理，以获得聚类类别数，将该聚类类别数作为所统计人群的人数。
与现有的人群统计技术相比，由于本发明的低密度人群的统计方法采用了特征点检测、跟踪及聚类的方法，因此克服了实际场景中存在的人群遮挡问题，准确地统计出低密度人群的人数。
与现有的技术相比，本发明的低密度人群的统计方法及装置，有效地克服由于遮掩等使得跟踪丢失或者错误的问题，提高了准确度；此外，可扩展性强，不仅仅适用于人群检测，稍加调整，可应用于车流量等类似其他目标数量统计。
以上所述，仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围，应当理解，本发明并不限于这里所描述的实现方案，这些实现方案描述的目的在于帮助本领域中的技术人员实践本发明。任何本领域中的技术人员很容易在不脱离本发明精神和范围的情况下进行进一步的改进和完善，因此本发明只受到本发明权利要求的内容和范围的限制，其意图涵盖所有包括在由所附权利要求所限定的本发明精神和范围内的备选方案和等同方案。

Claims (20)

1.一种低密度人群的统计方法，包括：
更新背景图像并提取当前图像的前景点；
对当前图像中的每个前景点进行检测以获取特征点；
获取每个特征点的运动轨迹；
根据运动轨迹的运动位移和运动方向变化率，提取出正常的运动轨迹；
将连续多帧图像内同一轨迹上的特征点作为一个分类点，对所有分类点进行聚类，以获得初步的分类；和
对初步的分类进行类分裂和类聚合处理，以获得聚类类别数，将该聚类类别数作为所统计人群的人数。
2.如权利要求1所述的方法，其中所述更新背景图像包括：
假设Ik表示第k帧图像，k为整数，Bk表示第k帧背景图像，其中背景图像的初始值为B0＝I0，则按下式更新背景图像：
B k ( x , y ) = B k - 1 ( x , y ) - 1 , B k - 1 ( x , y ) > I k ( x , y ) B k - 1 ( x , y ) + 1 , B k - 1 ( x , y ) < I k ( x , y ) B k - 1 ( x , y ) B k - 1 ( x , y ) = I k ( x , y )
其中，x、y分别表示像素点的横坐标和纵坐标。
3.按照权利要求2所述的方法，其中所述提取当前图像的前景点包括：第k帧图像的前景图像Fk为：Fk(x，y)＝|Ik(x，y)-Bk(x，y)|；对前景图像进行阈值分割以提取前景点。
4.按照权利要求3所述的方法，其中所述对前景图像进行阈值分割以提取前景点包括：
如果Fk(x，y)≥第一阈值T1，则认为该点(x，y)是当前图像的前景点，否则认为该改点(x，y)是当前图像的背景点。
5.按照权利要求1所述的方法，其中所述对当前图像中的每个前景点进行检测以获取特征点是利用角点检测算法来进行的。
6.按照权利要求1所述的方法，其中所述获取每个特征点的运动轨迹是采用KLT特征点跟踪算法来进行的。
7.按照权利要求1所述的方法，其中所述根据运动轨迹的运动位移和运动方向变化率，提取出正常的运动轨迹包括：分别计算每条运动轨迹在连续多帧图像内的运动位移、运动方向变化率；将同时满足运动位移＞第四阈值T4和运动方向变化率＜第五阈值T5的运动轨迹作为正常的轨迹并输出该正常的轨迹。
8.按照权利要求1所述的方法，其中所述聚类包括：
(1)将所有分类点看作为一个初始类C1，此时当前的类别数Max＝1；
(2)在所有类中选出直径最大且满足直径≥第六阈值T6的类Cm，其中m表示第m类，且1≤m≤Max，Max表示当前的类别数；
(3)在直径最大的类Cm中找出一个与其他分类点相异度最大的分类点，执行Max＝Max+1并将该分类点看作一个新的类CMax；
(4)分别计算直径最大的类Cm中分类点
Figure FSA00000330686700021
与新的类CMax中分类点的最大距离outerDist，其中
Figure FSA00000330686700022
表示直径最大的类Cm中第p个分类点，并计算直径最大的类Cm中分类点
Figure FSA00000330686700023
与该类中其他分类点的最小距离interDist，如果outerDist＜interDist则将该分类点
Figure FSA00000330686700024
放入新的类CMax中；
(5)如果直径最大的类Cm中没有分类点放入新的类CMax中，则执行下一步骤，否则返回到步骤(2)；
(6)如果所有类的直径＜第六阈值T6，则结束，否则返回到步骤(2)。
9.按照权利要求1所述的方法，其中所述类分裂包括：在每一个类中，计算每条轨迹在当前帧图像之后的多帧图像内的位移的平均方向，根据计算的平均方向所属的八个方向类别，将每个类分裂成更小的类，其中八个方向分别为：-22.5°～22.5°、22.5°～67.5°、67.5°～112.5°、112.5°～157.5°、157.5°～202.5°、202.5°～247.5°、247.5°～292.5°、292.5°～337.5°。
10.按照权利要求9所述的方法，其中所述类聚合包括：计算每个类的中心位置，如果两个类满足以下两个条件：①两个类的中心位置＜第八阈值T8；②两个类同属于八个方向中的同一方向类中，则将两个类聚合成一个类，直到没有两个类满足条件为止。
11.一种低密度人群的统计装置，其特征在于，该装置包括：
背景建立及前景点提取单元，用于更新背景图像并提取当前图像的前景点；
特征点检测单元，用于从当前图像的前景点中检测出特征点；
特征点运动轨迹获取单元，用于获取特征点运动轨迹；
正常运动轨迹提取单元，用于根据运动轨迹的运动位移和运动方向变化率，提取出正常的运动轨迹；
运动轨迹初步分类单元，用于将连续多帧图像内同一轨迹上的特征点作为一个分类点，对所有分类点进行聚类，以获得初步的分类；和
人数获取单元，用于对初步的分类进行类分裂和类聚合处理，以获得聚类类别数，将该聚类类别数作为所统计人群的人数。
12.按照权利要求11所述的装置，其中所述背景建立及前景点提取单元执行以下操作：
假设Ik表示第k帧图像，k为整数，Bk表示第k帧背景图像，其中背景图像的初始值为B0＝I0，则按下式更新背景图像：
B k ( x , y ) = B k - 1 ( x , y ) - 1 , B k - 1 ( x , y ) > I k ( x , y ) B k - 1 ( x , y ) + 1 , B k - 1 ( x , y ) < I k ( x , y ) B k - 1 ( x , y ) B k - 1 ( x , y ) = I k ( x , y )
其中，x、y分别表示像素点的横坐标和纵坐标。
13.按照权利要求11所述的装置，其中背景建立及前景点提取单元此外还执行以下操作：
第k帧图像的前景图像Fk为：Fk(x，y)＝|Ik(x，y)-Bk(x，y)|；对前景图像进行阈值分割以提取前景点。
14.按照权利要求13所述的装置，其中所述对前景图像进行阈值分割以提取前景点包括：
如果Fk(x，y)≥第一阈值T1，则认为该点(x，y)是当前图像的前景点，否则认为该点(x，y)是当前图像的背景点。
15.按照权利要求11所述的装置，其中所述特征点检测单元利用角点检测算法来进行检测。
16.按照权利要求11所述的装置，其中所述特征点运动轨迹获取单元采用KLT特征点跟踪算法来获取运动轨迹。
17.按照权利要求11所述的装置，其中所述正常运动轨迹提取单元执行以下操作：分别计算每条运动轨迹在连续多帧图像内的运动位移、运动方向变化率；将同时满足运动位移＞第四阈值T4和运动方向变化率＜第五阈值T5的运动轨迹作为正常的轨迹并输出该正常的轨迹。
18.按照权利要求11所述的装置，其中所述运动轨迹初步分类单元执行以下操作：
将所有分类点看作为一个初始类C1，此时当前的类别数Max＝1；
在所有类中选出直径最大且满足直径≥第六阈值T6的类Cm，其中m表示第m类，且1≤m≤Max，Max表示当前的类别数，在直径最大的类Cm中找出一个与其他分类点相异度最大的分类点，执行Max＝Max+1并将该分类点看作一个新的类CMax，分别计算直径最大的类Cm中分类点
Figure FSA00000330686700051
与新的类CMax中分类点的最大距离outerDist，其中
Figure FSA00000330686700052
表示直径最大的类Cm中第p个分类点，并计算直径最大的类Cm中分类点与该类中其他分类点的最小距离interDist，如果outerDist＜interDist则将该分类点
Figure FSA00000330686700054
放入新的类CMax中，直到直径最大的类Cm中没有分类点放入新的类CMax中并且所有类的直径＜第六阈值T6。
19.按照权利要求11所述的装置，其中所述人数获取单元执行的类分裂包括：在每一个类中，计算每条轨迹在当前帧图像之后的多帧图像内的位移的平均方向，根据计算的平均方向所属的八个方向类别，将每个类分裂成更小的类，其中八个方向分别为：-22.5°～22.5°、22.5°～67.5°、67.5°～112.5°、112.5°～157.5°、157.5°～202.5°、202.5°～247.5°、247.5°～292.5°、292.5°～337.5°。
20.按照权利要求19所述的装置，其中所述人数获取单元执行的类聚合处理包括：计算每个类的中心位置，如果两个类满足以下两个条件：①两个类的中心位置＜第八阈值T8；②两个类同属于八个方向中的同一方向类中，则将两个类聚合成一个类，直到没有两个类满足条件为止。
CN2010105296704A 2010-10-28 2010-10-28 低密度人群的统计方法及装置 Active CN101976353B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010105296704A CN101976353B (zh) 2010-10-28 2010-10-28 低密度人群的统计方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010105296704A CN101976353B (zh) 2010-10-28 2010-10-28 低密度人群的统计方法及装置

Publications (2)

Publication Number Publication Date
CN101976353A true CN101976353A (zh) 2011-02-16
CN101976353B CN101976353B (zh) 2012-08-22

Family

ID=43576237

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010105296704A Active CN101976353B (zh) 2010-10-28 2010-10-28 低密度人群的统计方法及装置

Country Status (1)

Country Link
CN (1) CN101976353B (zh)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521646A (zh) * 2011-11-11 2012-06-27 浙江捷尚视觉科技有限公司 基于深度信息聚类的复杂场景人数统计算法
CN102663776A (zh) * 2012-03-31 2012-09-12 北京智安邦科技有限公司 基于特征点分析的剧烈运动检测的方法及装置
CN103839038A (zh) * 2012-11-23 2014-06-04 浙江大华技术股份有限公司 一种人数统计的方法及装置
CN103854292A (zh) * 2012-11-28 2014-06-11 中兴通讯股份有限公司 一种人数及人群运动方向的计算方法及装置
CN105704434A (zh) * 2014-11-28 2016-06-22 上海新联纬讯科技发展有限公司 基于智能视频识别的体育场馆人群监控方法及系统
CN106815563A (zh) * 2016-12-27 2017-06-09 浙江大学 一种基于人体表观结构的人群数量预测方法
CN108133185A (zh) * 2017-12-21 2018-06-08 山东师范大学 基于轨迹数据判断行人关系的方法和系统
CN109460723A (zh) * 2018-10-26 2019-03-12 思百达物联网科技(北京)有限公司 对鼠情进行统计的方法、装置以及存储介质
CN111192604A (zh) * 2019-12-12 2020-05-22 秒针信息技术有限公司 录音设备控制方法和装置
CN111191524A (zh) * 2019-12-12 2020-05-22 华迪计算机集团有限公司 运动人群计数方法
CN113673392A (zh) * 2021-08-09 2021-11-19 南京理工大学 基于时空聚类的动目标检测方法
CN117911949A (zh) * 2024-01-19 2024-04-19 北京易融信达科技有限公司 社会稳定风险智能评估方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070003141A1 (en) * 2005-06-30 2007-01-04 Jens Rittscher System and method for automatic person counting and detection of specific events
US20070031005A1 (en) * 2000-09-06 2007-02-08 Nikos Paragios Real-time crowd density estimation from video
CN101321269A (zh) * 2007-06-05 2008-12-10 同济大学 基于计算机视觉的客流量检测方法及系统
CN101561929A (zh) * 2009-04-29 2009-10-21 同济大学 遥感影像的主分量和模糊聚类村镇专题信息提取方法
CN101872431A (zh) * 2010-02-10 2010-10-27 杭州海康威视软件有限公司 可适用多角度应用场景的人流量统计的方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070031005A1 (en) * 2000-09-06 2007-02-08 Nikos Paragios Real-time crowd density estimation from video
US20070003141A1 (en) * 2005-06-30 2007-01-04 Jens Rittscher System and method for automatic person counting and detection of specific events
CN101321269A (zh) * 2007-06-05 2008-12-10 同济大学 基于计算机视觉的客流量检测方法及系统
CN101561929A (zh) * 2009-04-29 2009-10-21 同济大学 遥感影像的主分量和模糊聚类村镇专题信息提取方法
CN101872431A (zh) * 2010-02-10 2010-10-27 杭州海康威视软件有限公司 可适用多角度应用场景的人流量统计的方法及系统

Cited By (18)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102521646B (zh) * 2011-11-11 2015-01-21 浙江捷尚视觉科技股份有限公司 基于深度信息聚类的复杂场景人数统计方法
CN102521646A (zh) * 2011-11-11 2012-06-27 浙江捷尚视觉科技有限公司 基于深度信息聚类的复杂场景人数统计算法
CN102663776A (zh) * 2012-03-31 2012-09-12 北京智安邦科技有限公司 基于特征点分析的剧烈运动检测的方法及装置
CN103839038A (zh) * 2012-11-23 2014-06-04 浙江大华技术股份有限公司 一种人数统计的方法及装置
CN103854292B (zh) * 2012-11-28 2018-02-23 中兴通讯股份有限公司 一种人数及人群运动方向的计算方法及装置
US9576199B2 (en) 2012-11-28 2017-02-21 Zte Corporation Method and device for calculating number and moving direction of pedestrians
CN103854292A (zh) * 2012-11-28 2014-06-11 中兴通讯股份有限公司 一种人数及人群运动方向的计算方法及装置
CN105704434A (zh) * 2014-11-28 2016-06-22 上海新联纬讯科技发展有限公司 基于智能视频识别的体育场馆人群监控方法及系统
CN106815563B (zh) * 2016-12-27 2020-06-02 浙江大学 一种基于人体表观结构的人群数量预测方法
CN106815563A (zh) * 2016-12-27 2017-06-09 浙江大学 一种基于人体表观结构的人群数量预测方法
CN108133185B (zh) * 2017-12-21 2021-02-23 山东师范大学 基于轨迹数据判断行人关系的方法和系统
CN108133185A (zh) * 2017-12-21 2018-06-08 山东师范大学 基于轨迹数据判断行人关系的方法和系统
CN109460723A (zh) * 2018-10-26 2019-03-12 思百达物联网科技(北京)有限公司 对鼠情进行统计的方法、装置以及存储介质
CN111192604A (zh) * 2019-12-12 2020-05-22 秒针信息技术有限公司 录音设备控制方法和装置
CN111191524A (zh) * 2019-12-12 2020-05-22 华迪计算机集团有限公司 运动人群计数方法
CN113673392A (zh) * 2021-08-09 2021-11-19 南京理工大学 基于时空聚类的动目标检测方法
CN117911949A (zh) * 2024-01-19 2024-04-19 北京易融信达科技有限公司 社会稳定风险智能评估方法及系统
CN117911949B (zh) * 2024-01-19 2024-06-11 北京易融信达科技有限公司 一种风险智能评估方法及系统

Also Published As

Publication number Publication date
CN101976353B (zh) 2012-08-22

Similar Documents

Publication Publication Date Title
CN101976353B (zh) 低密度人群的统计方法及装置
CN102063613B (zh) 基于头部识别的人群计数方法及装置
WO2022126669A1 (zh) 一种基于视频行人识别的地铁人流网络融合方法及人流预测方法
CN103824070B (zh) 一种基于计算机视觉的快速行人检测方法
DE112013001858B4 (de) Mehrfachhinweis-Objekterkennung und -Analyse
CN101311947A (zh) 基于自然视频的实时智能监控方法
CN102842045A (zh) 一种基于组合特征的行人检测方法
CN104732236B (zh) 一种基于分层处理的人群异常行为智能检测方法
Badura et al. Intelligent traffic system: Cooperation of MANET and image processing
CN106919939B (zh) 一种交通标识牌跟踪识别方法及系统
CN107330390A (zh) 一种基于图像分析和深度学习的人数统计方法
CN111753651A (zh) 一种基于车站二维人群密度分析的地铁群体异常行为检测方法
Xu et al. Efficient CityCam-to-edge cooperative learning for vehicle counting in ITS
CN105321187A (zh) 基于人头检测的行人计数方法
CN107659754A (zh) 一种树叶扰动情况下监控视频的有效浓缩方法
CN103927875A (zh) 基于视频的交通溢流状态识别方法
Ma et al. Event detection using local binary pattern based dynamic textures
Vu et al. Traffic incident recognition using empirical deep convolutional neural networks model
CN113963373A (zh) 一种基于视频图像动态检测与跟踪算法系统及方法
CN105447463A (zh) 变电站基于人体特征识别的跨摄像机自动跟踪系统
CN113239718A (zh) 一种校园人流量统计方法
Han et al. Smart CCTV and the management of urban space
CN102034243B (zh) 一种获取视频图像中人群密度图的方法及装置
CN101877135A (zh) 一种基于背景重构的运动目标检测方法
Yang et al. Crowd density and counting estimation based on image textural feature

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
ASS Succession or assignment of patent right

Owner name: NETPOSA TECHNOLOGIES, LTD.

Free format text: FORMER OWNER: BEIJING ZANB SCIENCE + TECHNOLOGY CO., LTD.

Effective date: 20150806

C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20150806

Address after: 100102, Beijing, Chaoyang District, Tong Tung Street, No. 1, Wangjing SOHO tower, two, C, 26 floor

Patentee after: NETPOSA TECHNOLOGIES, Ltd.

Address before: 100048 Beijing city Haidian District Road No. 9, building 4, 5 layers of international subject

Patentee before: Beijing ZANB Technology Co.,Ltd.

PP01 Preservation of patent right

Effective date of registration: 20220726

Granted publication date: 20120822

PP01 Preservation of patent right