CN113139995B - 一种低成本的物体间光线遮挡检测和评估方法 - Google Patents

一种低成本的物体间光线遮挡检测和评估方法 Download PDF

Info

Publication number
CN113139995B
CN113139995B CN202110418561.3A CN202110418561A CN113139995B CN 113139995 B CN113139995 B CN 113139995B CN 202110418561 A CN202110418561 A CN 202110418561A CN 113139995 B CN113139995 B CN 113139995B
Authority
CN
China
Prior art keywords
objects
light
projection
evaluating
polyhedrons
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.)
Active
Application number
CN202110418561.3A
Other languages
English (en)
Other versions
CN113139995A (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.)
Lingyang Technology Hangzhou Co ltd
Original Assignee
Hangzhou Bozi Enterprise Management Partnership LP
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 Hangzhou Bozi Enterprise Management Partnership LP filed Critical Hangzhou Bozi Enterprise Management Partnership LP
Priority to CN202110418561.3A priority Critical patent/CN113139995B/zh
Publication of CN113139995A publication Critical patent/CN113139995A/zh
Application granted granted Critical
Publication of CN113139995B publication Critical patent/CN113139995B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/507Depth or shape recovery from shading
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/08Projecting images onto non-planar surfaces, e.g. geodetic screens
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Geometry (AREA)
  • Image Generation (AREA)

Abstract

本发明公开了一种低成本的物体间光线遮挡检测和评估方法,包括以下步骤:(1)使用多面体分别表示目标物体和遮挡物体,分别计算两个多面体的顶点,得到两个多面体的顶点集合;(2)确定光线方向和投影面;(3)计算两个多面体的顶点在投影面上的坐标;(4)判断投影面上的两个多边形是否相交,若相交,则物体间发生光线遮挡,否则没有发生光线遮挡;对于发生光线遮挡的,进一步评估阴影重叠程度。本发明能够快速、有效的检测给出的光线条件下物体的阴影是否重叠,并评估重叠程度,算法计算量小,便于在小型设备上实施。

Description

一种低成本的物体间光线遮挡检测和评估方法
技术领域
本发明属于图形处理领域,尤其是涉及一种低成本的物体间光线遮挡检测和评估方法。
背景技术
在图形处理领域,遮挡检测是影响制作的关键一步。遮挡检测算法根据原理不同,可以分为Z-buffer算法、基于角度的遮挡检测算法和基于高程的遮挡检测算法。
如公开号为CN108182675A的中国专利文献公开了一种声波照射随机起伏界面时的面元遮挡判断方法,包括将随机起伏界面投影到观测平面,并得到观测平面边长分别为Lx和Ly;设置观测平面像素分辨率Npix;利用Z-buffer算法以像素为单位将各个面元离散化,计算得到每个面元的遮挡率QM;设置面元遮挡判断阈值Q,根据遮挡率QM和阈值Q的大小关系进行面元遮挡判断。
基于角度的遮挡检测算法由Habib等提出,算法的基本原理是:在地底点和待检测点的水平连线上,每个地物点对应的投影光线与水平面都有一个夹角,根据夹角的变化来分析待检测点的可见性。在地物点和地底点的水平连线上,逐渐远离地底点的方向上,如果投影光线与水平面的夹角逐渐变小,则没有遮挡,如果在某一位置突然变大,而后又变小回复或小于原来的角度,则该段区域被遮挡。该方法可参考谢文寒等发表的论文“城市大比例尺真正摄影像阴影与遮挡问题的研究”。
基于高程的遮蔽检测算法的基本原理是:在对某一地物点进行可见性分析时,如果该地物点可见则该地物点与投影中心的连线在地形数据的上方。
上述三种现有的遮挡检测算法中,基于角度的遮挡检测算法简单明了,能适用于各种复杂的环境,不存在伪遮蔽、伪可见的问题,得到了广泛应用。然而,常规的基于角度的遮挡检测算法需要对物体的每个点都进行角度分析,频繁的计算是的算法在没有充分优化时计算量大,耗时长。
发明内容
为解决现有技术存在的上述问题,本发明提供了一种低成本的物体间光线遮挡检测和评估方法,能够快速、有效的检测给出的光线条件下,物体的阴影是否重叠,并评估重叠程度。
一种低成本的物体间光线遮挡检测和评估方法,包括以下步骤:
(1)使用多面体分别表示目标物体和遮挡物体,分别计算两个多面体的顶点,得到两个多面体的顶点集合;
(2)确定光线方向和投影面;
(3)计算两个多面体的顶点在投影面上的坐标;
(4)使用多边形表示坐标围绕的区域;
(5)判断投影面上的两个多边形是否相交,若相交,则物体间发生光线遮挡,否则没有发生光线遮挡;对于发生光线遮挡的,进一步评估阴影重叠程度。
在遮挡检测问题中,一般已知的信息包括,光线的方向信息,物体的三维信息以及用于显示物体阴影的投影面信息。本发明的方法,能够减少阴影区域和遮挡关系判断的计算量,有利于在低计算资源设备上部署实现阴影遮挡判断,在地图绘制、光伏发电、影像制作领域都可以使用。
步骤(1)中,所述的目标物体和遮挡物体用多面体表示的公式为:
Figure BDA0003026961510000021
其中A1,B1,C1,D1表示目标物体或遮挡物体第一个面的参数,一共有n个面,x,y,z表示多面体中任意一点的三维坐标。
多面体的顶点计算方法为:
(1-1)在公式(1)中任意挑选三个面,计算三个面的交点,如果该交点属于多面体,则记录到多面体的顶点集合中;
(1-2)不断重复上述步骤,直到找到多面体的所有顶点。
进一步的,步骤(2)中,光线方向向量为:
Figure BDA0003026961510000031
投影面为:
m1x+m2y+m3z=0
其中,m1,m2,m3表示投影面法向向量在三维坐标系下的值。
步骤(3)的具体过程为:
对于多面体顶点集合中的任意一点p,其在投影平面上的投影点p'之间满足以下关系:
Figure BDA0003026961510000032
其中,
Figure BDA0003026961510000033
表示从点p指向p'的向量,
Figure BDA0003026961510000034
表示光线方向向量,r表示比值;通过公式(2),分别求出两个多面体所有顶点在投影面上的投影点。
步骤(4)中,使用多边形表示坐标围绕的区域的具体形式为::
以顺时针的方向确定所有投影点组成的多边形,该多边形在投影平面内以二维形式表示:
Figure BDA0003026961510000035
其中G11,G12,g1表示多边形第一条边的参数,一共有N条边,x,y表示多边形中任意一点的二维坐标。
(步骤(5)中,判断投影面上的两个多边形是否相交的方法为:
假设目标物体和遮挡物体对应的多边形参数分别为Gα、gα和Gβ、gβ,求解以下优化问题:
min J
s.t.
Figure BDA0003026961510000041
Figure BDA0003026961510000042
λ≥0
μ≥0
其中,
Figure BDA0003026961510000043
表示两个物体的投影多边形的最小距离,当J>0时,表示两个物体没有光线遮挡,否侧存在光线遮挡;λ和μ是两个辅助变量。
步骤(5)中,评估阴影重叠程度的过程为:
当发生阴影重叠,假设目标物体的中心在投影面的坐标为(xα,yα),遮挡物体的中心在投影面的坐标为(xβ,yβ),阴影重叠度的判断指标为:
Figure BDA0003026961510000044
其中ε表示重叠指标因子,|J|表示物体将最小距离的绝对值,||(xα,yα)-(xβ,yβ)||表示两中心的距离。
与现有技术相比,本发明具有以下有益效果:
1、能够快速、有效的检测给出的光线条件下,物体的阴影是否重叠,并评估重叠程度。
2、相较于传统算法,本发明计算量小,有利于在小型设备上实施。
3、避免使用图形处理单元等高成本的计算设备,有效的降低了设备成本。
附图说明
图1为本发明方法的流程示意图;
图2为本发明实施例物体和光线的三维示意图;
图3为本发明实施例中投影面上阴影重叠情况。
具体实施方式
下面结合附图和实施例对本发明做进一步详细描述,需要指出的是,以下所述实施例旨在便于对本发明的理解,而对其不起任何限定作用。
本实施例以两立方体在光线下的投影为例,对一种低成本的物体间光线遮挡检测和评估方法做详细描述。
如图1所示,一种低成本的物体间光线遮挡检测和评估方法,包括以下步骤:
步骤1,使用多面体分别表示目标物体和遮挡物体,计算两个多面体的顶点。具体过程为:
(1-1)如图2所示,立方体1和立方体2可表示为
Figure BDA0003026961510000051
Figure BDA0003026961510000052
(1-2)在公式(3)中任意挑选三个面,计算三个面的交点,如果该交点属于多面体,则记录到多面体的顶点集合中;
(1-3)不断重复步骤(1-2),直到找到多面体的所有顶点,在本实施例中,立方体1和立方体2的顶点集合分别是
{(1,1,1),(2,1,1),(2,1.5,1),(1,1.5,1),(1,1,2.2),(2,1,2.2),(2,1.5,2.2),(1,1.5,2.2)}
{(1.5,2,1.5),(2.5,2,1.5),(2.5,2.5,1.5),(1.5,2.5,1.5),(1.5,2,2.7),(2.5,2,2.7),(2.5,2.5,2,7),(1.5,2.5,2.7)}
步骤2,确定光线方向和投影面。
在本实施例中,光线方向向量为
Figure BDA0003026961510000053
投影面为:
y=3.5
步骤3,计算两个多面体的顶点在投影面上的坐标。
具体过程为:
对于多面体顶点集合中的任意一点p,其在投影平面上的投影点p'之间满足以下关系:
Figure BDA0003026961510000061
其中,
Figure BDA0003026961510000062
表示从点p指向p'的向量,
Figure BDA0003026961510000063
表示光线方向向量,r表示比值。通过公式(4),求出多面体所有顶点在投影面上的投影点。
步骤4,判断投影面上的两个多边形是否相交,若相交,则物体间发生光线遮挡,评估阴影重叠程度,否则没有发生光线遮挡。
具体方法是:
(4-1)以顺时针的方向确定所有投影点组成的多边形,该多边形在投影平面内以二维形式表示:
Figure BDA0003026961510000064
Figure BDA0003026961510000065
其中G11,G12,g1表示多边形一条边的参数,一共有N条边,x,y表示多边形中任意一点的二维坐标;
(4-2)假设目标物体和遮挡物体对应的多边形参数分别为Gα、gα和Gβ、gβ,求解以下优化问题:
min J
s.t.
Figure BDA0003026961510000066
Figure BDA0003026961510000067
λ≥0
μ≥0
其中,
Figure BDA0003026961510000068
表示两个物体的投影多边形的最小距离,当J>0时,表示两个物体没有光线遮挡,否则,存在光线遮挡,在本实施例中,J=-0.86,即存在光线遮挡。
(4-3)本实施例中的两个立方体发生阴影重叠,如图3所示,立方体1中心在投影面的坐标为(1.5,1.6),遮挡物体的中心在投影面的坐标为(2,2.1),阴影重叠度的判断指标为:
Figure BDA0003026961510000071
以上所述的实施例对本发明的技术方案和有益效果进行了详细说明,应理解的是以上所述仅为本发明的具体实施例,并不用于限制本发明,凡在本发明的原则范围内所做的任何修改、补充和等同替换,均应包含在本发明的保护范围之内。

Claims (6)

1.一种低成本的物体间光线遮挡检测和评估方法,其特征在于,包括以下步骤:
(1)使用多面体分别表示目标物体和遮挡物体,分别计算两个多面体的顶点,得到两个多面体的顶点集合;
(2)确定光线方向和投影面;
(3)计算两个多面体的顶点在投影面上的坐标;
(4)使用多边形表示坐标围绕的区域;
(5)判断投影面上的两个多边形是否相交,若相交,则物体间发生光线遮挡,否则没有发生光线遮挡;对于发生光线遮挡的,进一步评估阴影重叠程度;
判断投影面上的两个多边形是否相交的方法为:
假设目标物体和遮挡物体对应的多边形参数分别为Gα、gα和Gβ、gβ,求解以下优化问题:
min J
s.t.
Figure FDA0003595705040000011
Figure FDA0003595705040000012
λ≥0
μ≥0
其中,
Figure FDA0003595705040000013
表示两个物体的投影多边形的最小距离,当J>0时,表示两个物体没有光线遮挡,否侧存在光线遮挡;λ和μ是两个辅助变量;
评估阴影重叠程度的过程为:
当发生阴影重叠,假设目标物体的中心在投影面的坐标为(xα,yα),遮挡物体的中心在投影面的坐标为(xβ,yβ),阴影重叠度的判断指标为:
Figure FDA0003595705040000014
其中ε表示重叠指标因子,|J|表示物体将最小距离的绝对值,||(xα,yα)-(xβ,yβ)||表示两中心的距离。
2.根据权利要求1所述的低成本的物体间光线遮挡检测和评估方法,其特征在于,步骤(1)中,所述的目标物体和遮挡物体用多面体表示的公式为:
Figure FDA0003595705040000021
其中A1,B1,C1,D1表示目标物体或遮挡物体第一个面的参数,一共有n个面,x,y,z表示多面体中任意一点的三维坐标。
3.根据权利要求2所述的低成本的物体间光线遮挡检测和评估方法,其特征在于,步骤(1)中,多面体的顶点计算方法为:
(1-1)在公式(1)中任意挑选三个面,计算三个面的交点,如果该交点属于多面体,则记录到多面体的顶点集合中;
(1-2)不断重复上述步骤,直到找到多面体的所有顶点。
4.根据权利要求1所述的低成本的物体间光线遮挡检测和评估方法,其特征在于,步骤(2)中,光线方向向量为:
Figure FDA0003595705040000022
投影面为:
m1x+m2y+m3z=0
其中,m1,m2,m3表示投影面法向向量在三维坐标系下的值。
5.根据权利要求1所述的低成本的物体间光线遮挡检测和评估方法,其特征在于,步骤(3)的具体过程为:
对于多面体顶点集合中的任意一点p,其在投影平面上的投影点p'之间满足以下关系:
Figure FDA0003595705040000023
其中,
Figure FDA0003595705040000024
表示从点p指向p'的向量,
Figure FDA0003595705040000025
表示光线方向向量,r表示比值;通过公式(2),分别求出两个多面体所有顶点在投影面上的投影点。
6.根据权利要求1所述的低成本的物体间光线遮挡检测和评估方法,其特征在于,步骤(4)中,使用多边形表示坐标围绕的区域的具体形式为:
以顺时针的方向确定所有投影点组成的多边形,该多边形在投影平面内以二维形式表示:
Figure FDA0003595705040000031
其中G11,G12,g1表示多边形第一条边的参数,一共有N条边,x,y表示多边形中任意一点的二维坐标。
CN202110418561.3A 2021-04-19 2021-04-19 一种低成本的物体间光线遮挡检测和评估方法 Active CN113139995B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110418561.3A CN113139995B (zh) 2021-04-19 2021-04-19 一种低成本的物体间光线遮挡检测和评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110418561.3A CN113139995B (zh) 2021-04-19 2021-04-19 一种低成本的物体间光线遮挡检测和评估方法

Publications (2)

Publication Number Publication Date
CN113139995A CN113139995A (zh) 2021-07-20
CN113139995B true CN113139995B (zh) 2022-06-21

Family

ID=76812698

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110418561.3A Active CN113139995B (zh) 2021-04-19 2021-04-19 一种低成本的物体间光线遮挡检测和评估方法

Country Status (1)

Country Link
CN (1) CN113139995B (zh)

Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2302589A1 (en) * 2009-09-01 2011-03-30 Fondazione Bruno Kessler Method for efficient target detection from images robust to occlusion
CN103576994A (zh) * 2013-06-27 2014-02-12 香港应用科技研究院有限公司 在输入触摸面板上使用传感模块确定真正触点的方法及系统
CN104599253A (zh) * 2015-01-30 2015-05-06 武汉大学 一种自然图像阴影消除方法
CN105335732A (zh) * 2015-11-17 2016-02-17 西安电子科技大学 基于分块及鉴别非负矩阵分解的有遮挡人脸识别方法
CN105912852A (zh) * 2016-04-08 2016-08-31 河海大学 一种基于距离势函数任意凸多边形块体离散单元法
CN106203429A (zh) * 2016-07-06 2016-12-07 西北工业大学 基于双目立体视觉复杂背景下的遮挡目标检测方法
CN107944094A (zh) * 2017-11-06 2018-04-20 中国航天空气动力技术研究院 一种复杂外形航天器投影面积的确定方法及系统
CN108629742A (zh) * 2018-03-29 2018-10-09 深圳飞马机器人科技有限公司 真正射影像阴影检测与补偿方法、装置及存储介质
CN109961043A (zh) * 2019-03-22 2019-07-02 广西北斗星测绘科技有限公司 一种基于无人机高分辨率影像的单木高度测量方法及系统
CN110288203A (zh) * 2019-05-31 2019-09-27 中国计量大学 一种障碍物阴影区域光伏阵列的辐照量损失率动态计算方法
CN110349246A (zh) * 2019-07-17 2019-10-18 广西师范大学 一种应用于光场绘制中降低视点的重构失真度的方法
WO2019201355A1 (en) * 2018-04-17 2019-10-24 Shanghaitech University Light field system occlusion removal
CN111597733A (zh) * 2020-06-03 2020-08-28 清华大学 一种计算空间结构中的各种构件的太阳辐射量的方法
CN111625955A (zh) * 2020-05-25 2020-09-04 浙江大学 塔式太阳能镜场阴影与遮挡效率的计算方法

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2302589A1 (en) * 2009-09-01 2011-03-30 Fondazione Bruno Kessler Method for efficient target detection from images robust to occlusion
CN103576994A (zh) * 2013-06-27 2014-02-12 香港应用科技研究院有限公司 在输入触摸面板上使用传感模块确定真正触点的方法及系统
CN104599253A (zh) * 2015-01-30 2015-05-06 武汉大学 一种自然图像阴影消除方法
CN105335732A (zh) * 2015-11-17 2016-02-17 西安电子科技大学 基于分块及鉴别非负矩阵分解的有遮挡人脸识别方法
CN105912852A (zh) * 2016-04-08 2016-08-31 河海大学 一种基于距离势函数任意凸多边形块体离散单元法
CN106203429A (zh) * 2016-07-06 2016-12-07 西北工业大学 基于双目立体视觉复杂背景下的遮挡目标检测方法
CN107944094A (zh) * 2017-11-06 2018-04-20 中国航天空气动力技术研究院 一种复杂外形航天器投影面积的确定方法及系统
CN108629742A (zh) * 2018-03-29 2018-10-09 深圳飞马机器人科技有限公司 真正射影像阴影检测与补偿方法、装置及存储介质
WO2019201355A1 (en) * 2018-04-17 2019-10-24 Shanghaitech University Light field system occlusion removal
CN109961043A (zh) * 2019-03-22 2019-07-02 广西北斗星测绘科技有限公司 一种基于无人机高分辨率影像的单木高度测量方法及系统
CN110288203A (zh) * 2019-05-31 2019-09-27 中国计量大学 一种障碍物阴影区域光伏阵列的辐照量损失率动态计算方法
CN110349246A (zh) * 2019-07-17 2019-10-18 广西师范大学 一种应用于光场绘制中降低视点的重构失真度的方法
CN111625955A (zh) * 2020-05-25 2020-09-04 浙江大学 塔式太阳能镜场阴影与遮挡效率的计算方法
CN111597733A (zh) * 2020-06-03 2020-08-28 清华大学 一种计算空间结构中的各种构件的太阳辐射量的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
融合 SGM 与 Z-Buffer 的真正射影像遮挡检测算法;化雪诚 等;《桂林电子科技大学学报》;20170430;第37卷(第2期);1-6 *
遮挡检测算法分析比较研究;王玉锋 等;《科技视界》;20151231(第26期);1-2 *

Also Published As

Publication number Publication date
CN113139995A (zh) 2021-07-20

Similar Documents

Publication Publication Date Title
Schauer et al. The peopleremover—removing dynamic objects from 3-d point cloud data by traversing a voxel occupancy grid
CN111046776B (zh) 基于深度相机的移动机器人行进路径障碍物检测的方法
Poulin et al. Lights from highlights and shadows
US8743114B2 (en) Methods and systems to determine conservative view cell occlusion
CN103473814B (zh) 一种基于gpu的三维几何图元拾取方法
CN109903372B (zh) 深度图超分辨率补全方法及高质量三维重建方法与系统
CN104246831A (zh) 地图显示装置
CN102496181A (zh) 面向规模化生产的真正射影像制作方法
CN112258649A (zh) 一种基于光线投射的三维城市太阳辐射计算方法
US6356264B1 (en) Image processing device and image processing method
CN107765257A (zh) 一种基于反射强度辅助外部校准的激光探测与测量方法
CN113139995B (zh) 一种低成本的物体间光线遮挡检测和评估方法
CN112215308B (zh) 一种吊装物体单阶检测方法、装置、电子设备及存储介质
CN101937577A (zh) 生成具有边界像素超采样效果的阴影的方法
Zhou et al. Segmentation algorithm for 3D LiDAR point cloud based on region clustering
Jang et al. View-Dependent Multiresolution Splatting of Non-Uniform Data.
Rau et al. Lod generation for 3d polyhedral building model
Yokoyama et al. 3D modeling for historical structure using terrestrial laser ranging data
CN113808243B (zh) 一种可形变雪地网格的绘制方法和装置
US20240046549A1 (en) Ray tracing method and apparatus based on attention for dynamic scenes
Mura et al. Reconstructing Complex Indoor Environments with Arbitrary Wall Orientations.
CN106845501A (zh) 一种针对输变电工程环境敏感区域识别的缓冲区分析方法
JP2001148028A (ja) 図形表示装置及びその方法
Kim et al. Automatic 3D city modeling using a digital map and panoramic images from a mobile mapping system
Yang et al. LiDAR data reduction assisted by optical image for 3D building reconstruction

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20230531

Address after: Room 302, Building 1, No. 13, Yingqian Road, Meixian Section, Meixian Village, Guali Town, Xiaoshan District, Hangzhou City, Zhejiang Province, 311243

Patentee after: Lingyang Technology (Hangzhou) Co.,Ltd.

Address before: 310053 room 31232, 3rd floor, building 1, 88 Puyan street, Binjiang District, Hangzhou City, Zhejiang Province

Patentee before: Hangzhou Bozi enterprise management partnership (L.P.)

TR01 Transfer of patent right