• 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
    • / 13
    • 下载费用:30 金币  

    重庆时时彩预测推荐: 一种基于双目立体视觉的排爆机器人空间定位方法.pdf

    关 键 词:
    一种 基于 双目 立体 视觉 机器人 空间 定位 方法
      专利查询网所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
    摘要
    申请专利号:

    CN201510076626.5

    申请日:

    2015.02.11

    公开号:

    CN104680528A

    公开日:

    2015.06.03

    当前法律状态:

    驳回

    有效性:

    无权

    法律详情: 发明专利申请公布后的驳回IPC(主分类):G06T 7/00申请公布日:20150603|||实质审查的生效IPC(主分类):G06T 7/00申请日:20150211|||公开
    IPC分类号: G06T7/00 主分类号: G06T7/00
    申请人: 广州霞光技研有限公司; 中国人民解放军63908部队
    发明人: 蒋梁中; 马振书; 张勇; 马东玺; 杜峰坡; 罗磊; 徐学锦
    地址: 510000广东省广州市花都区炭步镇安泰路2号002
    优先权:
    专利代理机构: 北京市盈科律师事务所11344 代理人: 马丽丽
    PDF完整版下载: PDF下载
    法律状态
    申请(专利)号:

    CN201510076626.5

    授权公告号:

    ||||||

    法律状态公告日:

    2018.07.20|||2015.07.01|||2015.06.03

    法律状态类型:

    发明专利申请公布后的驳回|||实质审查的生效|||公开

    摘要

    本发明涉及一种基于双目立体视觉的排爆机器人空间定位方法,其包括:(1)根据实际的图像像素坐标结合摄像机内参数及畸变参数求出理想投影点归一化坐标;(2)根据理想投影点的归一化坐标和摄像机外参数获取左右投影点分别与摄像机光心形成的两射线的方程;(3)计算所述两射线的公垂线;(4)计算公垂线的长度并将公垂线的中点作为目标点的三维坐标。本发明提出的基于双目立体视觉的排爆机器人空间定位方法,能够提高排爆机器人对目标点三维坐标的求取速度和定位精度,有效地解决了目前排爆机器人对目标点三维坐标的求取速度和定位精度较低的技术问题。

    权利要求书

    权利要求书
    1.  一种基于双目立体视觉的排爆机器人空间定位方法,其特征在于,包括以下步骤(1) 至步骤(4):
    (1)根据实际的图像像素坐标结合摄像机内参数及畸变参数求出理想投影点归一化坐 标;
    (2)根据理想投影点的归一化坐标和摄像机外参数获取左右投影点分别与摄像机光心形 成的两射线的方程;
    (3)计算所述两射线的公垂线;
    (4)计算公垂线的长度并将公垂线的中点作为目标点的三维坐标。

    2.  根据权利要求1所述的基于双目立体视觉的排爆机器人空间定位方法,其特征在于, 所述步骤(1)包括步骤(11)和步骤(12):
    (11)令为理想投影点归一化坐标矢量:
    x n → = X c / Z c Y c / Z c = x n y n - - - ( 1 ) ]]>
    其中:Xc、Yc、Zc为摄像机光心的三维坐标,xn、yn为理想投影点归一化坐标;
    令r2=xn2+yn2,根据摄像机的非线性模型,考虑径向畸变和切向畸变,实际投影点的归 一化坐标矢量可表示如下:
    x d → = x d y d = ( 1 + k c 1 r 2 + k c 2 r 4 ) x n → + dx → - - - ( 2 ) ]]>
    dx → = 2 k c 3 x n y n + k c 4 ( r 2 + 2 x n 2 ) k c 3 ( r 2 + 2 y n 2 ) + 2 k c 4 x n y n - - - ( 3 ) ]]>
    其中:r为光心坐标在图像平面投影到原点的距离,d代表畸变后的坐标,xd和yd为图 像畸变后的二维坐标;为切向畸变矢量,kc1,kc2,kc3,kc4为镜头畸变系数;
    把式(2)和式(3)代入以下式(4)可得投影点的实际图像像素坐标(ud,vd)
    u = 1 dx ( x + α c y ) + u 0 ]]>   (4)
    v = y dy + v 0 ]]>
    u d = f x x d + f x α c y d + u 0 v d = f y y d + v 0 - - - ( 5 ) ]]>
    其中:u和v为目标点图像像素坐标,dx为一个像素在x方向上的物理尺寸,dy为一个 像素在y方向上的物理尺寸,(u0,v0)为光心的像素坐标;fx和fy为摄像机的焦距,αc为变 形因子,其与u,v像素轴的角度有关;
    (12)的求解过程:
    根据式(1)至式(5),根据标定所得出的镜头畸变系数,利用迭代法,可以由(ud,vd) 求得具体步骤如下:
    首先求得
    x d = ( u d - u 0 - α c ( v d - v 0 ) / f y ) / f x y d = ( v d - v 0 ) / f y - - - ( 6 ) ]]>
    由(xd,yd)经过如下迭代函数式可得(xn yn)T
    x n ( k + 1 ) = ( x d - 2 * k c 3 * x n ( k ) * y n ( k ) - k c 4 * ( r 2 ( k ) + 2 * x n ( k ) * x n ( k ) ) ) / ( 1 + k c 1 * r 2 ( k ) + k c 2 * r 2 ( k ) * r 2 ( k ) ) y n ( k + 1 ) = ( y d - k c 3 * ( r 2 ( k ) + 2 * y n ( k ) * y n ( k ) ) - 2 * k c 4 * x n ( k ) * y n ( k ) ) / ( 1 + k c 1 * r 2 ( k ) + k c 2 * r 2 ( k ) * r 2 ( k ) ) - - - ( 7 ) ]]>
    其中: r n ( k ) = ( x n ( k ) ) 2 + ( y n ( k ) ) 2 . ]]>

    3.  根据权利要求2所述的基于双目立体视觉的排爆机器人空间定位方法,其特征在于, 所述步骤(2)包括步骤(21)和步骤(22):
    (21)令世界坐标系与右摄像机坐标系重合无旋转,则两摄像机的外参数分别为:
    [R T]right=[I 0],其中,I为单位矩阵,R为右摄像机坐标系相对于左摄像坐标系的 旋转矩阵;T为右摄像机坐标系相对于左摄像坐标系的平移矩阵;
    [R T]left=[R T],其中[R T]为右摄像机坐标系相对于左摄像坐标系的旋转和平移矩 阵;
    (22)设左、右摄像机的光心分别为Ol和Or,场景中P点在左摄像机坐标系中的像点为 pl,空间点为Pl,空间点Pl坐标为(Xl,Yl,Zl),在右摄像机坐标系中的像点为pr,空间点为Pr, 坐标为(Xr,Yr,Zr),可根据Ol、pl、Or、pr四点的空间位置求所述两射线和
    令在世界坐标系下,则各矢量在右摄像机 坐系中的定义如下:



    4.  根据权利要求3所述的基于双目立体视觉的排爆机器人空间定位方法,其特征在于, 所述步骤(3)计算所述两射线的公垂线具体如下:

    其中,为所述两射线与的公垂线。

    5.  根据权利要求4所述的基于双目立体视觉的排爆机器人空间定位方法,其特征在于, 所述步骤(4)具体计算步骤如下:
    因为为射线与的公垂线,则可得:

    其中:⊙代表两向量的点乘;
    由于R,T可通过标定求得,(Xl/Zl Yl/Zl 1)T、(Xr/Zr Yr/Zr 1)T可根据式(7) 求得,因此在式(11)中Zl、Zr为仅有的两个未知量,式(11)有且存在唯一解,至此 Pr(Xr,Yr,Zr)Pl(Xl,Yl,Zl)完全确定,取公垂线的中点最为最终点P的坐标:
    P = P r + ( R * P l + T ) 2 - - - ( 12 ) . ]]>

    说明书

    说明书一种基于双目立体视觉的排爆机器人空间定位方法
    技术领域
    本发明涉及机器视觉领域,具体涉及一种基于双目立体视觉的排爆机器人空间定位方法。
    背景技术
    目前单目视觉系统排爆机器人采用一台摄像机对目标进行测量,其单目视觉系统对深度 信息的恢复能力较弱,多目视觉系统排爆机器人采用多台摄像机对目标进行测量,其多目测 量精度较高但是图像传输速度慢,计算量大,测量耗时多。同时由于排爆机器人工作环境的 特殊性,图像与数据信息的有线传输不仅增加了排爆的负担,有碍其灵活工作还给其自身安 全性带来隐患。且排爆机器人空间定位时,由于摄像机模型的和投影点的近似误差的存在, 左右投影点与摄像机光心的射线方程不交于一点,给目标点的三维求距带来了麻烦。因此, 现有技术中排爆机器人对目标点三维坐标的求取速度和定位精度都较低。
    发明内容
    本发明针对现有技术中排爆机器人对目标点三维坐标的求取速度和定位精度较低的技术 问题,提出一种基于双目立体视觉的排爆机器人空间定位方法。
    本发明提出的一种基于双目立体视觉的排爆机器人空间定位方法,其包括以下步骤(1) 至步骤(4):
    (1)根据实际的图像像素坐标结合摄像机内参数及畸变参数求出理想投影点归一化坐 标;
    (2)根据理想投影点的归一化坐标和摄像机外参数获取左右投影点分别与摄像机光心形 成的两射线的方程;
    (3)计算所述两射线的公垂线;
    (4)计算公垂线的长度并将公垂线的中点作为目标点的三维坐标。
    进一步,所述步骤(1)包括步骤(11)和步骤(12):
    (11)令为理想投影点归一化坐标矢量:
    x n → = X c / Z c Y c / Z c = x n y n - - - ( 1 ) ]]>
    其中:Xc、Yc、Zc为摄像机光心的三维坐标,xn、yn为理想投影点归一化坐标;
    令r2=xn2+yn2,根据摄像机的非线性模型,考虑径向畸变和切向畸变,实际投影点的归 一化坐标矢量可表示如下:
    x d → = x d y d = ( 1 + k c 1 r 2 + k c 2 r 4 ) x n → + d x → - - - ( 2 ) ]]>
    dx → = 2 k c 3 x n y n + k c 4 ( r 2 + 2 x n 2 ) k c 3 ( r 2 + 2 y n 2 ) + 2 k c 4 x n y n - - - ( 3 ) ]]>
    其中:r为光心坐标在图像平面投影到原点的距离,d代表畸变后的坐标,xd和yd为图 像畸变后的二维坐标;为切向畸变矢量,kc1,kc2,kc3,kc4为镜头畸变系数;
    把式(2)和式(3)代入以下式(4)可得投影点的实际图像像素坐标(ud,vd)
    u = 1 dx ( x + α c y ) + u 0 ]]>
    v = y dy + v 0 - - - ( 4 ) ]]>
    u d = f x x d + f x α c y d + u 0 v d = f y y d + v 0 - - - ( 5 ) ]]>
    其中:u和v为目标点图像像素坐标,dx为一个像素在x方向上的物理尺寸,dy为一个 像素在y方向上的物理尺寸,(u0,v0)为光心的像素坐标;fx和fy为摄像机的焦距,αc为变 形因子,其与u,v像素轴的角度有关;
    (12)的求解过程:
    根据式(1)至式(5),根据标定所得出的镜头畸变系数,利用迭代法,可以由(ud,vd) 求得具体步骤如下:
    首先求得
    x d = ( u d - u 0 - α c ( v d - v 0 ) / f y ) / f x y d = ( v d - v 0 ) / f y - - - ( 6 ) ]]>
    由(xd,yd)经过如下迭代函数式可得
    x n ( k + 1 ) = ( x d - 2 * k c 3 * x n ( k ) * y n ( k ) - k c 4 * ( r 2 ( k ) + 2 * x n ( k ) * x n ( k ) ) ) / ( 1 + k c 1 * r 2 ( k ) + k c 2 * r 2 ( k ) * r 2 ( k ) ) y n ( k + 1 ) = ( y d - k c 3 * ( r 2 ( k ) + 2 * y n ( k ) * y n ( k ) ) - 2 * k c 4 * x n ( k ) * y n ( k ) ) / ( 1 + k c 1 * r 2 ( k ) + k c 2 * r 2 ( k ) * r 2 ( k ) ) - - - ( 7 ) ]]>
    其中: r n ( k ) = ( x n ( k ) ) 2 + ( y n ( k ) ) 2 . ]]>
    进一步,所述步骤(2)包括步骤(21)和步骤(22):
    (21)令世界坐标系与右摄像机坐标系重合无旋转,则两摄像机的外参数分别为:
    [R T]right=[I 0],其中,I为单位矩阵,R为右摄像机坐标系相对于左摄像坐标系的 旋转矩阵;T为右摄像机坐标系相对于左摄像坐标系的平移矩阵;
    [R T]left=[R T],其中[R T]为右摄像机坐标系相对于左摄像坐标系的旋转和平移矩 阵;
    (22)设左、右摄像机的光心分别为Ol和Or,场景中P点在左摄像机坐标系中的像点为 pl,空间点为Pl,空间点Pl坐标为(Xl,Yl,Zl),在右摄像机坐标系中的像点为pr,空间点为Pr, 坐标为(Xr,Yr,Zr),可根据Ol、pl、Or、pr四点的空间位置求所述两射线和
    令在世界坐标系下,则各矢量在右摄像机 坐系中的定义如下:


    进一步,所述步骤(3)计算所述两射线的公垂线具体如下:

    其中,为所述两射线与的公垂线。
    进一步,所述步骤(4)具体计算步骤如下:
    因为为射线与的公垂线,则可得:

    其中:⊙代表两向量的点乘;
    由于R,T可通过标定求得,可根据式(7) 求得,因此在式(11)中Zl、Zr为仅有的两个未知量,式(11)有且存在唯一解,至此 Pr(Xr,Yr,Zr)Pl(Xl,Yl,Zl)完全确定,取公垂线的中点最为最终点P的坐标:
    P = P r + ( R * P l + T ) 2 - - - ( 12 ) ]]>
    有益效果:相对于现有技术中的单目或多目视觉系统排爆机器人,本发明提出的基于双 目立体视觉的排爆机器人空间定位方法,能够提高排爆机器人对目标点三维坐标的求取速度 和定位精度,有效地解决了目前排爆机器人对目标点三维坐标的求取速度和定位精度较低的 技术问题。
    附图说明
    图1是用于辅助理解本发明实施例方案的一种基于双目立体视觉的排爆机器人图像采集 系统的工作原理示意图。
    图2是本发明提出的基于双目立体视觉的排爆机器人空间定位方法实施例的计算原理示 意图。
    图3是本发明提出的基于双目立体视觉的排爆机器人空间定位方法实施例整体流程图。
    图4是发明提出的基于双目立体视觉的排爆机器人空间定位方法实施例详细流程图。
    具体实施方式
    为了便于本领域技术人员理解,下面将结合附图以及实施例对本发明进行进一步描述。
    在介绍本发明提出的基于双目立体视觉的排爆机器人空间定位方法实施例之前,作为辅 助理解,以下先介绍一种基于双目立体视觉的排爆机器人图像采集系统。
    请参阅图1,所述基于双目立体视觉的排爆机器人图像采集系统,包括左右两个摄像机、 无线图像收发器、图像采集卡、服务器、数传电台等。两个摄像机固定安装在排爆机器人机 械臂的小臂上作为机器人计算物体目标的“双眼”用来拍摄空间场景图像;摄像机采集的图 像通过无线图像传送设备传送到后台服务器,后台接收到图像后,通过安装在计算机中的图 像采集卡,便可以将图像送入计算机进行处理,处理结果通过数传电台再传送给下位机 PC/104。
    所述的无线图像收发器包括图像发送器和图像接收器,图像发送器与两个摄像机相连安 装在排爆机器人上,图像发送器将摄像机传来的图像通过无线传输至图像接收器,
    所述的图像采集卡安装在服务器内,是将图像接受器传来的模拟图像转换为计算机能识 别的数字图像,然后在对应的程序窗口中进行显示。
    所述基于双目立体视觉的排爆机器人图像采集系统具有以下优点:1、较单目视觉排爆机 器人恢复三维信息能力强;2、较多目排爆机器人处理速度快,定位精度高;3、由于采用了 无线传输方式,减轻了排爆机器人的负担,增加了其工作的灵活性,保障了其危险场合工作 的安全性。
    以下介绍本发明提出的基于双目立体视觉的排爆机器人空间定位方法实施例。
    请参阅图2至图4,本发明提出的基于双目立体视觉的排爆机器人空间定位方法实施例 包括以下步骤(1)至步骤(4):
    (1)根据实际的图像像素坐标结合摄像机内参数及畸变参数求出理想投影点归一化坐 标;
    (2)根据理想投影点的归一化坐标和摄像机外参数获取左右投影点分别与摄像机光心形 成的两射线的方程;
    (3)计算所述两射线的公垂线;
    (4)计算公垂线的长度并将公垂线的中点作为目标点的三维坐标。
    为了更好地实现本实施例的目的,本实施例对以上步骤(1)至步骤(4)做了具体的算 法设计。
    具体的,所述步骤(1)包括步骤(11)和步骤(12):
    (11)令为理想投影点归一化坐标矢量:
    x n → = X c / Z c Y c / Z c = x n y n - - - ( 1 ) ]]>
    其中:Xc、Yc、Zc为摄像机光心的三维坐标,xn、yn为理想投影点归一化坐标;
    令r2=xn2+yn2,根据摄像机的非线性模型,考虑径向畸变(本实施例中,由于径向畸变 较小,因此可以忽略)和切向(即切线方向)畸变,实际投影点的归一化坐标矢量可表示 如下:
    x d → = x d y d = ( 1 + k c 1 r 2 + k c 2 r 4 ) x n → + d x → - - - ( 2 ) ]]>
    dx → = 2 k c 3 x n y n + k c 4 ( r 2 + 2 x n 2 ) k c 3 ( r 2 + 2 y n 2 ) + 2 k c 4 x n y n - - - ( 3 ) ]]>
    其中:r为光心坐标在图像平面投影到原点的距离,d代表畸变后的坐标,xd和yd为图 像畸变后的二维坐标;为切向畸变矢量,kc1,kc2,kc3,kc4为镜头畸变系数;
    把式(2)和式(3)代入以下式(4)可得投影点的实际图像像素坐标(ud,vd)
    u = 1 dx ( x + α c y ) + u 0 ]]>
    v = y dy + v 0 - - - ( 4 ) ]]>
    u d = f x x d + f x α c y d + u 0 v d = f y y d + v 0 - - - ( 5 ) ]]>
    其中:u和v为目标点图像像素坐标,dx为一个像素在x方向上的物理尺寸,dy为一个 像素在y方向上的物理尺寸,(u0,v0)为光心的像素坐标;fx和fy为摄像机的焦距,αc为变 形因子,其与u,v像素轴的角度有关;
    (12)的求解过程:
    在三维计算的实际操作中,获得的投影点的图像像素坐标是经过畸变后的实际图像像素 坐标,我们需要根据实际图像像素坐标得到理想投影点归一化的图像平面坐标矢量,以便于 后面的几何求距;
    根据式(1)至式(5),根据标定所得出的镜头畸变系数,利用迭代法,可以由(ud,vd) 求得具体步骤如下:
    首先求得
    x d = ( u d - u 0 - α c ( v d - v 0 ) / f y ) / f x y d = ( v d - v 0 ) / f y - - - ( 6 ) ]]>
    由(xd,yd)经过如下迭代函数式可得
    x n ( k + 1 ) = ( x d - 2 * k c 3 * x n ( k ) * y n ( k ) - k c 4 * ( r 2 ( k ) + 2 * x n ( k ) * x n ( k ) ) ) / ( 1 + k c 1 * r 2 ( k ) + k c 2 * r 2 ( k ) * r 2 ( k ) ) y n ( k + 1 ) = ( y d - k c 3 * ( r 2 ( k ) + 2 * y n ( k ) * y n ( k ) ) - 2 * k c 4 * x n ( k ) * y n ( k ) ) / ( 1 + k c 1 * r 2 ( k ) + k c 2 * r 2 ( k ) * r 2 ( k ) ) - - - ( 7 ) ]]>
    其中: r n ( k ) = ( x n ( k ) ) 2 + ( y n ( k ) ) 2 ; ]]>
    在进行迭代时,可以设置迭代初值为迭代次数为20次。
    具体的,所述步骤(2)包括步骤(21)和步骤(22):
    (21)令世界坐标系与右摄像机坐标系重合无旋转,则两摄像机的外参数分别为:
    [R T]right=[I 0],其中,I为单位矩阵,R为右摄像机坐标系相对于左摄像坐标系的 旋转矩阵;T为右摄像机坐标系相对于左摄像坐标系的平移矩阵;
    [R T]left=[R T],其中[R T]为右摄像机坐标系相对于左摄像坐标系的旋转和平移矩 阵;
    (22)设左、右摄像机的光心分别为Ol和Or,场景中P点在左摄像机坐标系中的像点为 pl,空间点为Pl,空间点Pl坐标为(Xl,Yl,Zl),在右摄像机坐标系中的像点为pr,空间点为Pr, 坐标为(Xr,Yr,Zr),可根据Ol、pl、Or、pr四点的空间位置求所述两射线和
    令在世界坐标系(即右摄像机坐标系)下, 则各矢量在右摄像机坐系中的定义如下:


    具体的,所述步骤(3)计算所述两射线的公垂线具体如下:

    其中,为所述两射线与的公垂线。
    具体的,所述步骤(4)计算步骤如下:
    因为为射线与的公垂线,则可得:

    其中:⊙代表两向量的点乘;
    由于R,T可通过标定求得,可根据式(7) 求得,因此在式(11)中Zl、Zr为仅有的两个未知量,式(11)有且存在唯一解,至此 Pr(Xr,Yr,Zr)Pl(Xl,Yl,Zl)完全确定,取公垂线的中点最为最终点P的坐标:
    P = P r + ( R * P l + T ) 2 - - - ( 12 ) . ]]>
    本实施例提出的基于双目立体视觉的排爆机器人空间定位方法,能够提高排爆机器人对 目标点三维坐标的求取速度和定位精度,有效地解决了目前排爆机器人对目标点三维坐标的 求取速度和定位精度较低的技术问题。
    以上所述实施例仅表达了本发明的几种实施方式,其描述较为具体和详细,但并不能因 此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在 不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的?;し段?。 因此,本发明专利的?;し段вσ运饺ɡ笪?。

    关于本文
    本文标题:一种基于双目立体视觉的排爆机器人空间定位方法.pdf
    链接地址://www.4mum.com.cn/p-5886949.html
    关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

    [email protected] 2017-2018 www.4mum.com.cn网站版权所有
    经营许可证编号:粤ICP备17046363号-1 
     


    收起
    展开
  • 四川郎酒股份有限公司获第十二届人民企业社会责任奖年度环保奖 2019-05-13
  • 银保监会新规剑指大企业多头融资和过度融资 2019-05-12
  • 韩国再提4国联合申办世界杯 中国网友无视:我们自己来 2019-05-11
  • 中国人为什么一定要买房? 2019-05-11
  • 十九大精神进校园:风正扬帆当有为 勇做时代弄潮儿 2019-05-10
  • 粽叶飘香幸福邻里——廊坊市举办“我们的节日·端午”主题活动 2019-05-09
  • 太原设禁鸣路段 设备在测试中 2019-05-09
  • 拜耳医药保健有限公司获第十二届人民企业社会责任奖年度企业奖 2019-05-08
  • “港独”没出路!“梁天琦们”该醒醒了 2019-05-07
  • 陈卫平:中国文化内涵包含三方面 文化复兴表现在其中 2019-05-06
  • 人民日报客户端辟谣:“合成军装照”产品请放心使用 2019-05-05
  • 【十九大·理论新视野】为什么要“建设现代化经济体系”?   2019-05-04
  • 聚焦2017年乌鲁木齐市老城区改造提升工程 2019-05-04
  • 【专家谈】上合组织——构建区域命运共同体的有力实践者 2019-05-03
  • 【华商侃车NO.192】 亲!楼市火爆,别忘了买车位啊! 2019-05-03
  • kk盘锦棋牌外挂 昨天体彩浙江11选5开奖 浙江11选5选号技巧 中国福利彩票36选7开奖结果 双色球+7码复式破解 河北十一选五手机走势 貔喜游戏脉动棋牌官网 巨丰投资股票推荐 体彩河北十一选五开奖结果 新浪棋牌下载 黑马股票推荐11月 博远棋牌最新版 铁龙物流股票行情 搜河北十一选五走势图 福彩3d组三组六怎么看 代买彩票