太阳城集团

  • / 9
  • 下载费用:30 金币  

一种基于人头检测的计数方法.pdf

摘要
申请专利号:

CN201611235768.2

申请日:

2016.12.28

公开号:

太阳城集团CN106779043A

公开日:

2017.05.31

当前法律状态:

实审

有效性:

审中

法律详情: 实质审查的生效IPC(主分类):G06M 11/00申请日:20161228|||公开
IPC分类号: G06M11/00; G06K9/00; G06K9/32; G06K9/46 主分类号: G06M11/00
申请人: 南京沃顿物联网科技有限公司
发明人: 孙成建; 杜辉
地址: 210038 江苏省南京市栖霞区燕子矶街道和燕路371号
优先权:
专利代理机构: 南京瑞弘专利商标事务所(普通合伙) 32249 代理人: 吴旭
PDF完整版下载: PDF下载
法律状态
申请(专利)号:

CN201611235768.2

授权太阳城集团号:

|||

法律状态太阳城集团日:

2017.06.23|||2017.05.31

法律状态类型:

实质审查的生效|||公开

摘要

本发明公开了一种基于人头检测的计数方法,主要功能是能够自动实现公交车、办公室和商场等进、出口场所行人的人数统计,便于相关人员进行高效的统计分析。本方法能够根据人头部的颜色特征和类似椭圆特征,采用三帧差分法能够比较准确的检测出每一帧图像中的行人的头部轮廓,并与给定的最小轮廓面积S进行比较,筛选出符合要求的轮廓;在进、出人数统计环节,运用两个计数器分别统计进、出总人数,并且以人头计数的方式避免了行人之间由于相互遮挡而造成的统计误差。根据实验确定,本方法能够实现3到5个人同时进、出情况下的人数统计,且准确率最高可以达到90%左右。

权利要求书

1.一种基于人头检测的计数方法,其特征在于,包括以下步骤:
步骤(1),通过摄像机采集实时图像F,并设定帧计数器C记录视频的帧数;
步骤(2),如果C=1,则对图像F进行初始化;否则,通过金字塔放缩算法使得每一帧图
像F的高度和宽度减半,得到缩放后的每一帧图像SF;
步骤(3),对缩放后的每一帧图像SF进行图像灰度转换;
步骤(4),对经步骤(3)操作后的每一帧图像SF,根据人头的颜色特性与类似椭圆特性,
采用三帧差分法进行人头检测;
步骤(5),对人头检测后的每一帧图像SF依次进行膨胀、腐蚀操作,初步提取出对应的头
部轮廓;
步骤(6),采用混合高斯背景建模进行像素匹配,如果匹配成功,则被判断为背景像素,
否则被判断为前景像素,即头部轮廓;像素匹配的过程中同时进行背景更新;
步骤(7),对每一帧图像SF中检测到的头部轮廓进行相应的面积筛选,对于大于给定最
小轮廓面积S的头部轮廓绘制最小外接矩形;
步骤(8),定义每一帧图像SF的最小外接矩形内部轮廓的状态:o表示到当前为止图像SF
中原有目标轮廓的状态,n表示目标轮廓新进入SF中的状态,g表示目标轮廓离开SF的状态,m
表示目标轮廓之间的合并,s表示目标轮廓之间的分离;
步骤(9),设定行人进、出方向标志f以及进、出人数计数器I与O,以所有最小外接矩形
内的人头数目之和为计数标准,实现行人的进、出统计。
2.根据权利要求1所述的一种基于人头检测的计数方法,其特征在于,所述步骤(1)中,
摄像机采集实时图像F为原视频序列的奇数帧图像。
3.根据权利要求1所述的一种基于人头检测的计数方法,其特征在于,所述步骤(3)具
体包括以下步骤:
(31),采用颜色空间转换函数CvtColor(.)将缩放过后的图像SF分别转换为前景灰度图
像Fr与背景灰度图像Bk;
(32),采用灰度矩阵转换函数CvConvert(.)将前景灰度图像Fr转换为前景灰度矩阵FM,
背景灰度图像Bk转换为背景灰度矩阵BM;
(33),采用阈值函数CvThreshod(.)将前景灰度图像Fr转换为二值前景灰度矩阵FMB。
4.根据权利要求1所述的一种基于人头检测的计数方法,其特征在于,所述步骤(4)中
三帧差分法的阈值DT确定方法为:DT=h*w*p;其中h,w,p分别为图像SF的高度、宽度以及百
分比;定义st表示选取连续三帧图像的顺序标志,st=0表示选取第一帧图像,st=1表示选
取第二帧图像,st=2表示选取第三帧图像;三帧差分法每一帧的选取方式包括以下具体步
骤:
(41),如果二值前景灰度矩阵FMB的非零像素个数大于三帧差分法的阈值DT,同时st=0,
表示第一帧的异物大于阈值,则需要将当前帧SF存储为第一帧F1,继续采集第二帧;
(42),否则,如果二值前景灰度矩阵FMB的非零像素个数小于三帧差分法的阈值DT,同时
st=1,表示第一帧的异物大于阈值,而第二帧没有;
(43),否则,如果二值前景灰度矩阵FMB的非零像素个数大于三帧差分法的阈值DT,同时
st=1,表示第一帧和第二帧的异物都大于阈值,且第一帧和第二帧连续,则继续采集第3
帧;
(44),否则,如果前景灰度二值矩阵FMB的非零像素个数小于三帧差分法的阈值DT,同时
st=2,表示第一帧和第二帧的异物都大于阈值,而第三帧没有;
(45),否则,如果前景灰度二值矩阵FMB的非零像素个数大于三帧差分法的阈值DT,同时
st=2,且第二帧和第三帧连续,则三帧差分法选取成功。
5.根据权利要求1所述的一种基于人头检测的计数方法,其特征在于,所述步骤(7)对
头部轮廓进行面积筛选包括以下具体步骤:
(71),获取每一帧图像SF中每个头部轮廓的面积A,如果面积A大于给定的最小轮廓面积
S,则为符合要求的头部轮廓;
(72),根据SF的高度h和宽度w确定每一帧图像SF的感兴趣区域,对于分辨率为h*w尺寸
的SF,所述感兴趣区域的尺寸设置为0.5~0.8h*0.5~0.8w;
(73),采用轮廓绘制函数cvBoundingRect(.)对符合要求的头部轮廓寻找最小外接拟
合矩形,cX=cvRound(r.w/2+r.x),cY=cvRound(r.h/2+r.y);其中r.x与r.y分别表示头部
轮廓的最小外接矩形的原点坐标,r.w与r.h分别表示最小外接矩形的宽度和高度,cvRound
(.)表示四舍五入算法,cX,cY分别代表外接拟合矩形中心的横坐标和纵坐标;
(74),如果头部轮廓的最小外接拟合矩形的中心坐标cX,cY同时位于所述感兴趣区域
内,则绘制该头部轮廓的最小外接拟合矩形。
6.根据权利要求5所述的一种基于人头检测的计数方法,其特征在于,所述步骤(9)中,
只有当行人进入或离开SF时,并且头部轮廓的最小外接拟合矩形的中心坐标的数值大于或
小于感兴趣区域的中心坐标时才开始计数,包括以下具体步骤:
(91),感兴趣区域的面积分别由w/2-DX/2,w/2+DX/2,h/2-DY/2,h/2+DY/2四条边围成的
矩形区域,其中DX,DY分别表示感兴趣区域的横坐标和纵坐标;
(92),如果头部轮廓的最小外接拟合矩形中心坐标cX的取值先大于cvRound(w/2-DX/
2),再大于cvRound(w/2),则表示进,设置进出标志f=1;
(93),如果头部轮廓的最小外接拟合矩形中心坐标cX的取值先小于cvRound(w/2+DX/
2),再小于cvRound(w/2),则表示出,设置进出标志f=2;
(94),当某个目标的状态为n,同时行人进出标志f=1,则表示行人为进,因此I=I1+i;
(95),当某个目标的状态为g,同时行人进出标志f=2,则表示行人为出,因此O=O1+i;
其中i为某个轮廓内的人数,I为到当前为止,进入感兴趣区域内的总人数,I1表示前一
刻的进入感兴趣区域内目标总人数,O为到当前为止,离开感兴趣区域内的总人数,O1表示
前一刻的离开感兴趣区域内目标总人数。

说明书

一种基于人头检测的计数方法

技术领域

本发明涉及一种基于人头检测的计数方法及其装置,该装置能够自动实现公交车
和办公场所行人的进、出人数统计,便于相关人员进行高效的统计分析。

背景技术

随着社会的进步,人们对生活的需求越来越高,为了节省人力,具有自动化性能的
设备显得越来越重要,人数统计装置就是生活中一个很重要的方面。一个性能良好的人数
统计装置给人们的出行带来了很大的方便。用户可以在掌握监控区域实时动态太阳城集团的同
时,还能构及时得到现场准确的人数和人群流量太阳城集团,统计比较不同时期客流量,可以评估
营销、促销策略的合理性;特殊时段、特殊区域可起到警戒作用,避免造成的财产损失;同
时,人数统计还是大型商场、购物中心、连锁店、机场、车站、博物馆、展览馆等公共场所在管
理和决策方面不可缺少的数据。对于零售业而言,人流量更是非常基础的指标。因其与销售
量直接的正比关系。因此,具有良好客流量统计性能的系统在生活中具有很重要的作用。

发明内容

发明目的:针对上述现有技术,提出一种基于人头检测的计数方法,能够准确的统
计出公共场所行人的进、出人数。

技术方案:一种基于人头检测的计数方法,包括以下步骤:

步骤(1),通过摄像机采集实时图像F,并设定帧计数器C记录视频的帧数;

步骤(2),如果C=1,则对图像F进行初始化;否则,通过金字塔放缩算法使得每一
帧图像F的高度和宽度减半,得到缩放后的每一帧图像SF;

步骤(3),对缩放后的每一帧图像SF进行图像灰度转换;

步骤(4),对经步骤(3)操作后的每一帧图像SF,根据人头的颜色特性与类似椭圆
特性,采用三帧差分法进行人头检测;

步骤(5),对人头检测后的每一帧图像SF依次进行膨胀、腐蚀操作,初步提取出对
应的头部轮廓;

步骤(6),采用混合高斯背景建模进行像素匹配,如果匹配成功,则被判断为背景
像素,否则被判断为前景像素,即头部轮廓;像素匹配的过程中同时进行背景更新;

步骤(7),对每一帧图像SF中检测到的头部轮廓进行相应的面积筛选,对于大于给
定最小轮廓面积S的头部轮廓绘制最小外接矩形;

步骤(8),定义每一帧图像SF的最小外接矩形内部轮廓的状态:o表示到当前为止
图像SF中原有目标轮廓的状态,n表示目标轮廓新进入SF中的状态,g表示目标轮廓离开SF的
状态,m表示目标轮廓之间的合并,s表示目标轮廓之间的分离;

步骤(9),设定行人进、出方向标志f以及进、出人数计数器I与O,以所有最小外接
矩形内的人头数目之和为计数标准,实现行人的进、出统计。

进一步的,所述步骤(1)中,摄像机采集实时图像F为原视频序列的奇数帧图像。

进一步的,所述步骤(3)具体包括以下步骤:

(31),采用颜色空间转换函数CvtColor(.)将缩放过后的图像SF分别转换为前景
灰度图像Fr与背景灰度图像Bk;

(32),采用灰度矩阵转换函数CvConvert(.)将前景灰度图像Fr转换为前景灰度矩
阵FM,背景灰度图像Bk转换为背景灰度矩阵BM;

(33),采用阈值函数CvThreshod(.)将前景灰度图像Fr转换为二值前景灰度矩阵
FMB。

进一步的,所述步骤(4)中三帧差分法的阈值DT确定方法为:DT=h*w*p;其中h,w,p
分别为图像SF的高度、宽度以及百分比;定义st表示选取连续三帧图像的顺序标志,st=0表
示选取第一帧图像,st=1表示选取第二帧图像,st=2表示选取第三帧图像;三帧差分法每
一帧的选取方式包括以下具体步骤:

(41),如果二值前景灰度矩阵FMB的非零像素个数大于三帧差分法的阈值DT,同时
st=0,表示第一帧的异物大于阈值,则需要将当前帧SF存储为第一帧F1,继续采集第二帧;

(42),否则,如果二值前景灰度矩阵FMB的非零像素个数小于三帧差分法的阈值DT,
同时st=1,表示第一帧的异物大于阈值,而第二帧没有;

(43),否则,如果二值前景灰度矩阵FMB的非零像素个数大于三帧差分法的阈值DT,
同时st=1,表示第一帧和第二帧的异物都大于阈值,且第一帧和第二帧连续,则继续采集
第3帧;

(44),否则,如果前景灰度二值矩阵FMB的非零像素个数小于三帧差分法的阈值DT,
同时st=2,表示第一帧和第二帧的异物都大于阈值,而第三帧没有;

(45),否则,如果前景灰度二值矩阵FMB的非零像素个数大于三帧差分法的阈值DT,
同时st=2,且第二帧和第三帧连续,则三帧差分法选取成功。

进一步的,所述步骤(7)对头部轮廓进行面积筛选包括以下具体步骤:

(71),获取每一帧图像SF中每个头部轮廓的面积A,如果面积A大于给定的最小轮
廓面积S,则为符合要求的头部轮廓;

(72),根据SF的高度h和宽度w确定每一帧图像SF的感兴趣区域,对于分辨率为h*w
尺寸的SF,所述感兴趣区域的尺寸设置为0.5~0.8h*0.5~0.8w;

(73),采用轮廓绘制函数cvBoundingRect(.)对符合要求的头部轮廓寻找最小外
接拟合矩形,cX=cvRound(r.w/2+r.x),cY=cvRound(r.h/2+r.y);其中r.x与r.y分别表示
头部轮廓的最小外接矩形的原点坐标,r.w与r.h分别表示最小外接矩形的宽度和高度,
cvRound(.)表示四舍五入算法,cX,cY分别代表外接拟合矩形中心的横坐标和纵坐标;

(74),如果头部轮廓的最小外接拟合矩形的中心坐标cX,cY同时位于所述感兴趣区
域内,则绘制该头部轮廓的最小外接拟合矩形。

进一步的,所述步骤(9)中,只有当行人进入或离开SF时,并且头部轮廓的最小外
接拟合矩形的中心坐标的数值大于或小于感兴趣区域的中心坐标时才开始计数,包括以下
具体步骤:

(91),感兴趣区域的面积分别由w/2-DX/2,w/2+DX/2,h/2-DY/2,h/2+DY/2四条边围
成的矩形区域,其中DX,DY分别表示感兴趣区域的横坐标和纵坐标;

(92),如果头部轮廓的最小外接拟合矩形中心坐标cX的取值先大于cvRound(w/2-
DX/2),再大于cvRound(w/2),则表示进,设置进出标志f=1;

(93),如果头部轮廓的最小外接拟合矩形中心坐标cX的取值先小于cvRound(w/2+
DX/2),再小于cvRound(w/2),则表示出,设置进出标志f=2;

(94),当某个目标的状态为n,同时行人进出标志f=1,则表示行人为进,因此I=
I1+i;

(95),当某个目标的状态为g,同时行人进出标志f=2,则表示行人为出,因此O=
O1+i;

其中i为某个轮廓内的人数,I为到当前为止,进入感兴趣区域内的总人数,I1表示
前一刻的进入感兴趣区域内目标总人数,O为到当前为止,离开感兴趣区域内的总人数,O1
表示前一刻的离开感兴趣区域内目标总人数。

有益效果:本发明主要实现一种基于人头检测的计数方法,主要功能是能够自动
实现公交车、办公室和商场等进、出口场所行人的人数统计,便于相关人员进行高效的统计
分析。本发明的有益效果主要体现在以下几个方面:第一,本发明所用的视频是通过摄像头
对准固定的场景俯视拍摄的,主要是把摄像头垂直向下固定在门上,因此背景不会发生太
大变化;第二,根据人头部的颜色特征和类似椭圆特征,采用三帧差分法能够比较准确的检
测出每一帧图像中的行人的头部轮廓,并与给定的最小轮廓面积S进行比较,筛选出符合要
求的轮廓。第三,为了保证检测的准确性,同时对检测后的每一帧图像再次进行混合高斯背
景建模,即便是在背景相对固定的情况下,仍然进行一定的背景更新,这在一定程度上防止
了人头的漏检与误检,提高了检测准确度。第四,在进、出人数统计环节,运用两个计数器分
别统计进、出总人数,并且以人头计数的方式避免了行人之间由于相互遮挡而造成的统计
误差。经实验确定,本发明能够实现3到5个人同时进、出情况下的人数统计,且准确率最高
可以达到90%左右。

附图说明

图1为本发明方法的流程图。

具体实施方式

下面结合附图对本发明做更进一步的解释。

如图1所示,一种基于人头检测的计数方法,包括以下步骤:

步骤(1),通过摄像机采集实时图像F,并设定帧计数器C记录视频的帧数。

步骤(2),如果C=1,则对图像F进行初始化,即对处理过程中所涉及到的中间环节
进行申请内存;否则,通过金字塔放缩算法使得每一帧图像F的高度和宽度减半,得到缩放
后的每一帧图像SF,这样做的目的是减小检测区域的大小,提高检测效率。

步骤(3),对缩放后的每一帧图像SF进行图像灰度转换;具体为:

(31),采用颜色空间转换函数CvtColor(.)将缩放过后的图像SF分别转换为前景
灰度图像Fr与背景灰度图像Bk,以便进行人头识别。

(32),采用灰度矩阵转换函数CvConvert(.)将前景灰度图像Fr转换为前景灰度矩
阵FM,背景灰度图像Bk转换为背景灰度矩阵BM。

(33),采用阈值函数CvThreshod(.)将前景灰度图像Fr转换为二值前景灰度矩阵
FMB,实现轮廓提取。

步骤(4),对经步骤(3)操作后的每一帧图像SF,根据人头的颜色特性与类似椭圆
特性,采用三帧差分法进行人头检测。三帧差分法的阈值DT确定方法为:DT=h*w*p;其中h,
w,p分别为图像SF的高度、宽度以及百分比;定义st表示选取连续三帧图像的顺序标志,st=
0表示选取第一帧图像,st=1表示选取第二帧图像,st=2表示选取第三帧图像;三帧差分法
每一帧的选取方式包括以下具体步骤:

(41),如果二值前景灰度矩阵FMB的非零像素个数大于三帧差分法的阈值DT,同时
st=0,表示第一帧的异物大于阈值,则需要将当前帧SF存储为第一帧F1,继续采集第二帧。

(42),否则,如果二值前景灰度矩阵FMB的非零像素个数小于三帧差分法的阈值DT,
同时st=1,表示第一帧的异物大于阈值,而第二帧没有。

(43),否则,如果二值前景灰度矩阵FMB的非零像素个数大于三帧差分法的阈值DT,
同时st=1,表示第一帧和第二帧的异物都大于阈值,且第一帧和第二帧连续,则继续采集
第3帧。

(44),否则,如果前景灰度二值矩阵FMB的非零像素个数小于三帧差分法的阈值DT,
同时st=2,表示第一帧和第二帧的异物都大于阈值,而第三帧没有。

(45),否则,如果前景灰度二值矩阵FMB的非零像素个数大于三帧差分法的阈值DT,
同时st=2,且第二帧和第三帧连续,则三帧差分法选取成功。

步骤(5),对人头检测后的每一帧图像SF依次进行膨胀、腐蚀操作,初步提取出对
应的头部轮廓。

步骤(6),为了防止图像SF中人头部轮廓的漏检与误检,进一步采用混合高斯背景
建模进行像素匹配,如果匹配成功,则被判断为背景像素,否则被判断为前景像素,即头部
轮廓;像素匹配的过程中同时进行背景更新。

步骤(7),对每一帧图像SF中检测到的头部轮廓进行相应的面积筛选,对于大于给
定最小轮廓面积S的头部轮廓绘制最小外接矩形。对头部轮廓进行面积筛选包括以下具体
步骤:

(71),获取每一帧图像SF中每个头部轮廓的面积A,如果面积A大于给定的最小轮
廓面积S,则为符合要求的头部轮廓。

(72),根据SF的高度h和宽度w确定每一帧图像SF的感兴趣区域;对于分辨率为h*w
尺寸的SF,所述感兴趣区域的尺寸设置为0.5~0.8h*0.5~0.8w,例如对于分辨率为
600x800尺寸的SF图像,感兴趣区域一般可以设置为400x600。

(73),采用轮廓绘制函数cvBoundingRect(.)对符合要求的头部轮廓寻找最小外
接拟合矩形,cX=cvRound(r.w/2+r.x),cY=cvRound(r.h/2+r.y);其中r.x与r.y分别表示
头部轮廓的最小外接矩形的原点坐标,r.w与r.h分别表示最小外接矩形的宽度和高度,
cvRound(.)表示四舍五入算法,cX,cY分别代表外接拟合矩形中心的横坐标和纵坐标。

(74),如果头部轮廓的最小外接拟合矩形的中心坐标cX,cY同时位于所述感兴趣区
域内,则绘制该头部轮廓的最小外接拟合矩形,即:

cX>cvRound(w/2-DX/2),

且cX<cvRound(w/2+DX/2),

cY>cvRound(h/2-DY/2),

且cY<cvRound(h/2+DY/2),

其中,DX,DY分别表示感兴趣区域的横坐标和纵坐标。

步骤(8),定义每一帧图像SF的最小外接矩形内部轮廓的状态:o表示到当前为止
图像SF中原有目标轮廓的状态,n表示目标轮廓新进入SF中的状态,g表示目标轮廓离开SF的
状态,m表示目标轮廓之间的合并,s表示目标轮廓之间的分离。其中,如果某个头部轮廓的
最小外接矩形的状态为o,则该目标为SF中原有的目标,此时不更新人数太阳城集团;如果某个头
部轮廓的最小外接矩形的状态为n,则该目标为SF中新进入的目标,此时更新人数太阳城集团;如
果某个头部轮廓的最小外接矩形的状态为g,则该目标为SF中离开的目标,此时更新人数信
息;如果某个头部轮廓的最小外接矩形的状态为m,则表示SF中发生目标融合,即若干个最
小外接矩形轮廓融合为一个,此时不更新人数太阳城集团,同时最小外接矩形轮廓内的总人数为
被融合的轮廓的总人数之和;如果某个头部轮廓的最小外接矩形的状态为s,则表示SF中发
生目标分离,即一个最小外接轮廓矩形分离成若干个;此时不更新人数太阳城集团,同时分离后最
小外接矩形轮廓内的总人数仍为未分离之前的总人数之和。

步骤(9),设定行人进、出方向标志f以及进、出人数计数器I与O,以所有最小外接
矩形内的人头数目之和为计数标准,实现行人的进、出统计。具体为,只有当行人进入或离
开SF时,并且头部轮廓的最小外接拟合矩形的中心坐标的数值大于或小于感兴趣区域的中
心坐标时才开始计数,包括以下具体步骤:

(91),感兴趣区域的面积分别由w/2-DX/2,w/2+DX/2,h/2-DY/2,h/2+DY/2四条边围
成的矩形区域,其中DX,DY分别表示感兴趣区域的横坐标和纵坐标;

(92),如果头部轮廓的最小外接拟合矩形中心坐标cX的取值先大于cvRound(w/2-
DX/2),再大于cvRound(w/2),则表示进,设置进出标志f=1;

(93),如果头部轮廓的最小外接拟合矩形中心坐标cX的取值先小于cvRound(w/2+
DX/2),再小于cvRound(w/2),则表示出,设置进出标志f=2;

(94),当某个目标的状态为n,同时行人进出标志f=1,则表示行人为进,因此I=
I1+i;

(95),当某个目标的状态为g,同时行人进出标志f=2,则表示行人为出,因此O=
O1+i;

其中i为某个轮廓内的人数,I为到当前为止,进入感兴趣区域内的总人数,I1表示
前一刻的进入感兴趣区域内目标总人数,O为到当前为止,离开感兴趣区域内的总人数,O1
表示前一刻的离开感兴趣区域内目标总人数。

优选的,步骤(1)中,摄像机采集实时图像F为原视频序列的奇数帧图像。

以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人
员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应
视为本发明的保护范围。

关 键 词:
一种 基于 人头 检测 计数 方法
  专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
太阳城集团本文
本文标题:一种基于人头检测的计数方法.pdf
链接地址:http://zh228.com/p-6019790.html
太阳城集团我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们

copyright@ 2017-2018 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备17046363号-1 
 


收起
展开
葡京赌场|welcome document.write ('');