CN114757993A - 一种纹影图像的运动、参数识别方法及系统 - Google Patents

一种纹影图像的运动、参数识别方法及系统 Download PDF

Info

Publication number
CN114757993A
CN114757993A CN202210658880.6A CN202210658880A CN114757993A CN 114757993 A CN114757993 A CN 114757993A CN 202210658880 A CN202210658880 A CN 202210658880A CN 114757993 A CN114757993 A CN 114757993A
Authority
CN
China
Prior art keywords
light source
schlieren
motion
motion model
marking
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
CN202210658880.6A
Other languages
English (en)
Other versions
CN114757993B (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.)
Institute of Mechanics of CAS
Original Assignee
Institute of Mechanics of CAS
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 Institute of Mechanics of CAS filed Critical Institute of Mechanics of CAS
Priority to CN202210658880.6A priority Critical patent/CN114757993B/zh
Publication of CN114757993A publication Critical patent/CN114757993A/zh
Application granted granted Critical
Publication of CN114757993B publication Critical patent/CN114757993B/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/70Determining position or orientation of objects or cameras
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00Aerodynamic testing; Arrangements in or on wind tunnels
    • G01M9/06Measuring arrangements specially adapted for aerodynamic testing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/20Analysis of motion
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • G06T7/66Analysis of geometric attributes of image moments or centre of gravity
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30241Trajectory
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Geometry (AREA)
  • Multimedia (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Fluid Mechanics (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明提供了一种纹影图像的运动、参数识别方法及系统,在运动模型的壁面固定设有多个标记光源,每个标记光源的灯心安装方向与纹影光源的光路平行,且面向纹影相机用于标记运动模型的运动轨迹;风洞试验后,获得包含有运动模型和标记光源的原始纹影图像,基于标记光源在运动模型上相对于运动模型的空间位置坐标及运动模型的质心坐标分析原始纹影图像,获得风洞实验中运动模型的分离运动轨迹以及动力学参数。本发明在运动模型上安装标记光源,通过简单的坐标对应关系以及几何关系处理获得试验中运动模型的位移、俯仰角变化,解决了现有技术中测量难度大、多相机位受场地限制、标定复杂的问题。

Description

一种纹影图像的运动、参数识别方法及系统
技术领域
本发明涉及风洞特种试验领域,具体涉及一种纹影图像的运动、参数识别方法及系统,特别适用于超声速/高超声速风洞多体分离、自由飞试验时,能够获得模型分离轨迹和运动参数的识别方法。
背景技术
在航空航天基础研究领域,高超声速风洞动态试验技术为研究高超声速飞行器动态飞行特性、多体分离动态特性做出了很大的贡献。目前在超声速/高超风洞中进行的多体分离试验一般采用CTS轨迹模拟系统进行试验,为一种准静态的分离模拟试验,非真实飞行条件下的动态分离模拟试验。模型自由飞和动态分离试验虽然有CTS系统支撑,但是无法直接获得模型的运动参数以及受力状态等数据信息。而这些试验数据对分析飞行器的动态(分离)气动特性以及流场分析非常重要。
在风洞动态试验中,一般用于获得飞行器运动参数信息的方式有两种:一是摄影测量法,二是在模型内部安装运动加速度传感器和信号采集器。对于摄像测量法,需要采用多相机或者单相机拍摄,在不同的角度对同一目标点进行多次拍摄。标定好双相机之间几何关系或者单相机不同机位之间的几何关系,可以通过对目标物立体像对的分析获得目标的空间坐标。该方法在应用与超声速/高超声速风洞试验中,测量难度大、设备成本高、标定复杂且测量效率低。一方面需要精心设计多相机位,但是这往往受到风洞试验舱观察窗口以及场地的限制,另一方面需要在风洞试验中给运动模型进行打光,反射足够多的光用于相机镜头捕捉。最后还需要进行复杂的坐标转换以及图像处理方法才能得到目标点的三维空间坐标和随时间的运动轨迹。这种方法的测量精度与相机方位、相机台数以及处理方法密切相关。所以综合来说,这种方法对于风洞动态试验的模型运动参数计算和捕捉来说难度较大,成本过高,实施操作复杂。最后对于在运动模型中安装加速度传感器及其信号采集存储系统的方法来说,由于一般用于风洞动态分离或者自由飞试验的试验模型尺寸较小,无法全部放置其中测量所需的惯性加速度传感器以及信号采集存储系统,而且尽管模型做的足够大用于安装测量设备,但是会增加试验模型重量,影响或者偏离预期动态试验效果,会对试验产生重要影响。另一方面,这会增加试验成本。
以上这些情况限制下,就需要一种普适性强,实施简单方便、成本低、数据处理简便高效而且对风洞动态试验结果影响小的模型多体分离或自由飞轨迹和运动参数识别方法。针对以上问题,提出一种基于风洞试验纹影图像的多体分离、自由飞轨迹和运动参数识别方法。
发明内容
为了解决以上方案中的缺陷和不足,本发明提出了一种纹影图像的运动、参数识别方法及系统,该识别方法是针对在超声速/高超声速风洞中开展多体动态分离或者自由飞试验时,通过对纹影图像进行后期处理获得运动模型的运动参数和受力状态等信息,具有成本低、操作简单高效、对实验影响小等特点,解决了现有技术中测量难度大、多相机位受场地限制、标定复杂的问题。
在本发明的第一个方面,提供一种纹影图像的运动、参数识别方法,包括如下方法步骤:
步骤1、在运动模型的壁面固定设有多个标记光源,每个所述标记光源的灯心安装方向与纹影光源的光路平行,且面向纹影相机用于标记运动模型的运动轨迹;
步骤2、风洞试验后,获得包含有运动模型和标记光源的原始纹影图像,识别所述原始纹影图像中的标记光源的位置像素点,获得每张所述原始纹影图像中每个所述标记光源的像素坐标;
步骤3、利用步骤2处理后的纹影图像中某一时刻
Figure 269277DEST_PATH_IMAGE002
、上一时刻
Figure 739573DEST_PATH_IMAGE004
的每个所述标记光源的所述像素坐标和运动模型的质心坐标计算运动模型的俯仰角姿态变化,获得某一时刻相对于所述原始纹影图像的起始时刻的运动模型的姿态角的变化量;
步骤4、以步骤2处理后的纹影图片中的某一固定点作为原点,建立参考坐标系,将每个所述标记光源的所述像素坐标转化为物理空间坐标变化,结合步骤3的结果,获得风洞实验中运动模型的分离运动轨迹以及动力学参数。
在一个可能的实施方式中,步骤1中,每个所述标记光源的中轴线与所述纹影光源的光路平行实现每个所述标记光源的灯心安装方向与纹影光源的光路平行,且每个所述标记光源背向纹影光源安装在所述运动模型的壁面上用于标记运动模型的运动轨迹以及姿态。
在一个可能的实施方式中,步骤2中所述获得每张纹影图像中每个标记光源的像素坐标的具体内容为:
获得包含有运动模型和标记光源的原始纹影图像,对所有所述原始纹影图像进行批量灰度处理得到灰度图片,以判定出标记光源的位置像素点与背景像素点;继续将所述灰度图片转化为黑白二值图片,以得到每个标记光源的灯心的像素点坐标;获得每张纹影图片中每个标记光源的像素坐标。
在一个可能的实施方式中,步骤4中所述参考坐标系中,以风洞试验中试验气流的来流方向建立参考坐标系的横轴,垂直来流方向向上为参考坐标系的纵轴,得到每个所述标记光源在所述运动模型的运动过程中在所述参考坐标系中的坐标和所述运动模型的质心在参考坐标系中的坐标。
在一个可能的实施方式中,在所述参考坐标系中,结合所述运动模型在某一时刻
Figure 407315DEST_PATH_IMAGE005
相对于上一时刻的
Figure DEST_PATH_IMAGE006
的位移,得到所述运动模型的运动速度、加速度、受力和俯仰力矩。
在一个可能的实施方式中,步骤4中的所述动力学参数的精度随着标记光源的数量增多而增高。
在一个可能的实施方式中,步骤1中,所述标记光源的数量大于等于3个。
在本发明的第二个方面,提供一种应用了上述方法的纹影图像的运动、参数识别系统,包括:用于风洞试验的运动模型、纹影光源和纹影相机,在运动模型的壁面上固定设有3个及以上光线方向背对纹影光源且面向纹影相机镜头的标记光源,每个所述标记光源的中轴线与纹影光源的光路平行用于标记运动模型的运动轨迹,在运动模型的空腔内部固定有连接标记光源的供电电路;
风洞试验后,获得包含有运动模型和标记光源的原始纹影图像,基于标记光源在运动模型上相对于运动模型的空间位置坐标及运动模型的质心坐标分析所述原始纹影图像。
在一个可能的实施方式中,每个所述标记光源固定在运动模型壁面的一个贯穿孔内。
在一个可能的实施方式中,所述标记光源是LED灯、小电珠或小灯泡中的一种。
本发明和现有技术相比具有如下有益效果:本发明提供了一种纹影图像的运动、参数识别方法及系统,该方法系统在运动模型上安装了标记光源,通过简单的坐标对应关系以及几何关系处理获得试验中运动模型的位移、俯仰角变化,进而获得加速度以及受力状态,该方法操作简单,方便快捷、精度可靠,而且成本很低,对风洞试验效果几乎没有影响。
附图说明
为了更清楚地说明本发明的实施方式或现有技术中的技术方案,下面将对实施方式或现有技术描述中所需要使用的附图作简单地介绍。显而易见地,下面描述中的附图仅仅是示例性的,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图引申获得其它的实施附图。
图1为本发明实施例中风洞试验纹影图像的运动、参数识别方法的流程示意图;
图2为本发明中风洞试验纹影图像的运动、参数识别系统的结构示意图;
图中:1-标记光源,2-运动模型,3-纹影光源,4-纹影相机,5-观察窗。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种纹影图像的运动、参数识别方法,包括如下方法步骤:
步骤1、在运动模型的壁面固定设有多个标记光源,每个所述标记光源的灯心安装方向与纹影光源的光路平行,且面向纹影相机用于标记运动模型的运动轨迹。
在一个具体实施例中,标记光源数量一般大于等于3个,通过每个标记光源的中轴线与纹影光源的光路平行可以实现标记光源的灯心安装方向与纹影光源的光路平行的安装要求。标记光源的中轴线与纹影光路平行时,使得标记光源中间的最大光强亮度正朝向纹影相机拍摄,使得纹影相机捕捉到最大亮度的标记光源发光强度,方便后续轨迹识别处理。
此外,由图2可知,每个标记光源面向纹影相机即背向纹影光源,每个标记光源必须背向纹影光源安装在所述运动模型的壁面上,是为了在运动模型背向纹影光源的一面形成亮光标记点,可以标记运动模型在风洞试验中的运动轨迹。
在纹影拍摄系统中,由于运动模型会挡住纹影光源发出的光,所以模型在拍摄的纹影照片中是黑色的,而在运动模型背着纹影光源的一侧壁面上安装标记光源(优选白光LED灯),一方面可以清晰的标记出模型的运动轨迹以及姿态,一方面LED灯发出的白光与模型的黑色环境形成强烈对比,方便后续图像轨迹识别的处理。
本实施例中,标记光源优选为白光LED灯,以LED灯为实施例的流程图如图1所示,标记光源的数量视具体情况而定,数量越多,后处理运动的参数精度越高。
试验前,需要确定标记光源在运动模型上相对于运动模型的空间位置坐标
Figure DEST_PATH_IMAGE008
,以及运动模型的质心坐标
Figure DEST_PATH_IMAGE010
。在试验前需要确定纹影观察窗口的尺寸,比如观察窗的直径D或者长和宽。确保以上工作齐全后,开始风洞试验,获得包含运动模型和标记光源的纹影图像,纹影相机拍摄帧率为f。
步骤2、风洞试验后,获得包含有运动模型和标记光源的原始纹影图像,识别所述原始纹影图像中的标记光源的位置像素点,获得每张所述原始纹影图像中每个所述标记光源的像素坐标。
获得纹影图像的原始图片后,需要对试验过程中运动模型上安装的若干标记光源进行图像识别,由于标记光源发白光并且亮度明显区别于背景环境,所以该过程比较简单。首先对原始纹影图像进行批量处理,将原始彩色图片转变为灰度图片,然后对比得出图片中标记光源的像素点的灰度值,设定某一灰度阈值,能够判定出标记光源的位置像素点与背景像素点。根据设定的灰度阈值,将灰度图片转化为二值图片(即黑和白),黑色用来标记背景像素、白色用来标记运动模型上安装的标记光源位置。这样就能够很容易的得到各个标记光源白色像素点的坐标,然后对每个标记光源的像素点坐标做平均化处理,获得标记光源灯心(贯穿孔圆心)的像素点坐标,以此类推,获得每张纹影图片中每个标记光源的像素坐标。
步骤3、利用步骤2处理后的纹影图像中某一时刻
Figure 935248DEST_PATH_IMAGE002
、上一时刻
Figure 381142DEST_PATH_IMAGE004
的每个所述标记光源的所述像素坐标和运动模型的质心坐标计算运动模型的俯仰角姿态变化,获得某一时刻相对于所述原始纹影图像的起始时刻的运动模型的姿态角的变化量。
在获得模型运动过程中的每张图片中每个标记光源的像素坐标
Figure DEST_PATH_IMAGE012
之后,需要将标记光源的像素坐标转化为物理空间坐标变化(即位移和俯仰角变化)。由于试验的运动模型一般为面对称模型,在试验模拟中,运动模型一般只在俯仰方向存在明显的姿态角变化,在侧滑以及滚转方向上无明显位移或者姿态变化,所以通过纹影图像进行后处理获得轴向和法向位移以及俯仰角变化是足够的。
根据纹影相机的拍摄图片,获得纹影观察窗在纹影相机拍摄下的像素坐标范围,比如纹影观察窗直径D对应的像素坐标长度为D’,那么一个像素点的边长代表真实空间长度为D/D’。所以根据纹影图像中某一时刻
Figure DEST_PATH_IMAGE014
和上一时刻的
Figure DEST_PATH_IMAGE016
的每一个标记光源的像素坐标。以下为运动模型的位移、俯仰角、加速度以及受力状态计算过程及表达式。
定义两向量之间所成的夹角计算函数为
Figure DEST_PATH_IMAGE018
Figure DEST_PATH_IMAGE020
Figure DEST_PATH_IMAGE022
上述函数中,A、B为向量,x1、y1为向量A的分量,x2、y2为向量B的分量。该计算函数的定义为:
定义计算函数 f=getangle(x1,y1,x2,y2),函数名为getangle,输入的自变量为向量A的分量x1、y1和向量B的分量x2、y2。经过文中的计算公式可以获得计算输出两向量A、B之间的夹角 angle。dot(A,B)函数为求两向量的点积。
每个标记光源与运动模型的质心之间形成的直线与水平线所成夹角为
Figure DEST_PATH_IMAGE024
每个标记光源与运动模型的质心所成线段长度为
Figure DEST_PATH_IMAGE026
利用纹影图像中某一时刻
Figure 590931DEST_PATH_IMAGE005
和上一时刻的
Figure 796784DEST_PATH_IMAGE006
的每一个补光像素坐标计算运动模型的俯仰角姿态变化:采用最小二乘法确定时刻
Figure 38279DEST_PATH_IMAGE002
和上一时刻的
Figure 936964DEST_PATH_IMAGE004
的补光等拟合出的直线斜率
Figure DEST_PATH_IMAGE028
Figure DEST_PATH_IMAGE030
,然后通过斜率确定当前时刻
Figure 568803DEST_PATH_IMAGE002
与上一时刻
Figure 312768DEST_PATH_IMAGE031
的姿态角相对变化值:
Figure 425080DEST_PATH_IMAGE033
以此类推,获得每一时刻与前一时刻的相对姿态角变化量。然后对这些姿态角变化量进行累加取和,获得某一时刻相对于原始纹影图像的起始时刻的运动模型的姿态角的变化量为:
Figure 475426DEST_PATH_IMAGE035
步骤4、以步骤2处理后的纹影图片中的某一固定点作为原点,建立参考坐标系,将每个所述标记光源的所述像素坐标转化为物理空间坐标变化,结合步骤3的结果,获得风洞实验中运动模型的分离运动轨迹以及动力学参数。
取纹影图像中的某一固定点作为参考坐标系的原点,该原点在纹影图像中的像素坐标为
Figure 407610DEST_PATH_IMAGE037
。以风洞试验中气流的来流方向建立坐标系的横轴,垂直来流方向向上为坐标的纵轴,则每个标记光源在模型运动过程中在建立的参考坐标系的坐标为:
Figure 955266DEST_PATH_IMAGE039
则运动模型的质心在参考坐标系中的坐标为:
Figure 922085DEST_PATH_IMAGE041
(当质心位于该LED灯之前时,式中“±”取负号,反之取正号)
Figure 411841DEST_PATH_IMAGE043
(当质心位于该LED灯之下时,式中“±”取负号,反之取正号)
则运动模型的质心在时刻
Figure 831321DEST_PATH_IMAGE005
的坐标为:
Figure 182668DEST_PATH_IMAGE045
式中,n为LED灯个数。
则模型在某一时刻
Figure 987682DEST_PATH_IMAGE005
相对于上一时刻的
Figure 399072DEST_PATH_IMAGE006
的位移为:
Figure 305848DEST_PATH_IMAGE047
模型运动速度为:
Figure 929727DEST_PATH_IMAGE049
加速度为:
Figure 120406DEST_PATH_IMAGE051
受力为:
Figure 702697DEST_PATH_IMAGE053
俯仰力矩为:
Figure DEST_PATH_IMAGE055
其中m为模型质量,Izz为绕运动模型体轴系的Z轴的惯性矩。
经过以上步骤便可以基于原始纹影图像计算得到运动模型的分离轨迹以及运动参数,并且在JF12激波风洞两级入轨飞行器纵向级间分离轨道级模型运动参数识别中获得了成功的应用,成本低,方便快捷,且具有良好的识别精度。
如图2所示,本发明还提供一种应用了上述纹影图像的运动、参数识别方法的系统,包括:用于风洞试验的运动模型2、纹影光源3和纹影相机4,在运动模型2的壁面上固定设有3个及以上光线方向背对纹影光源3且面向纹影相机4镜头的标记光源1,每个所述标记光源1的中轴线与纹影光源3的光路平行用于标记运动模型2的运动轨迹,在运动模型2的空腔内部固定有连接标记光源1的供电电路;风洞试验后,获得包含有运动模型2和标记光源1的原始纹影图像,基于标记光源1在运动模型2上相对于运动模型2的空间位置坐标及运动模型2的质心坐标分析所述原始纹影图像,获得风洞实验中运动模型的分离运动轨迹以及动力学参数。
在一个具体的实施例中,在运动模型2的壁面一侧上打上直径为3~5 mm的若干贯穿孔,贯穿孔的中轴线与纹影光源3的光路平行,每个贯穿孔内安装固定有一个标记光源1,保证了标记光源1的中轴线与纹影光源3的光路平行。
标记光源可以选择体积小、质量小、电路简单、安装方便的光源,可以是LED灯、小电珠、小灯泡中的一种,一般情况下LED灯是最理想的选择。在不明显影响试验效果的情况下,也可以在运动模型2的壁面上安装其他发光体。
本实施例中,标记光源选择LED灯,试验前,准备好足够数量的白光LED灯(发光二极管),并且准备好串联电阻、电源(可以选择纽扣电池)、LED灯以及导线组成的串联电路。在关闭试验舱前,沿运动模型壁面的贯穿孔安装白光LED发光二极管,LED安装尽量保证灯的中轴线与孔的中轴线重合,然后用胶水粘合固定,电源电压需要保证LED发光亮度在试验高帧率的纹影相机拍摄下能够明显辨析(因为高速相机拍摄帧数越高,图片亮度越暗),即保证LED灯的足够发光亮度(可以选择串联多块纽扣电池)。然后在试验模型空腔内用胶带或者其他方式固定好LED等供电电路,防止试验时由于模型运动导致电路断开或者接触不良影响LED正常工作。
风洞试验后,获得包含有运动模型和标记光源的原始纹影图像,基于标记光源在运动模型上相对于运动模型的空间位置坐标及运动模型的质心坐标分析所述原始纹影图像,基于上述分析方法,可以获得风洞实验中运动模型的分离运动轨迹以及动力学参数。
本发明提供的识别方法及系统,通过在运动模型上安装标记光源,经过简单的坐标对应关系以及几何关系处理即可获得试验中运动模型的位移、俯仰角变化,进而获得加速度以及受力状态,该方法操作简单,方便快捷、精度可靠,而且成本很低,选择质量小的LED灯等标记光源,对风洞试验效果几乎没有影响。
以上实施例仅为本申请的示例性实施例,不用于限制本申请,本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内,对本申请做出各种修改或等同替换,这种修改或等同替换也应视为落在本申请的保护范围内。

Claims (10)

1.一种纹影图像的运动、参数识别方法,其特征在于,包括如下方法步骤:
步骤1、在运动模型的壁面固定设有多个标记光源,每个所述标记光源的灯心安装方向与纹影光源的光路平行,且面向纹影相机用于标记运动模型的运动轨迹;
步骤2、风洞试验后,获得包含有运动模型和标记光源的原始纹影图像,识别所述原始纹影图像中的标记光源的位置像素点,获得每张所述原始纹影图像中每个所述标记光源的像素坐标;
步骤3、利用步骤2处理后的纹影图像中某一时刻
Figure 60238DEST_PATH_IMAGE002
、上一时刻
Figure DEST_PATH_IMAGE003
的每个所述标记光源 的所述像素坐标和运动模型的质心坐标计算运动模型的俯仰角姿态变化,获得某一时刻相 对于所述原始纹影图像的起始时刻的运动模型的姿态角的变化量;
步骤4、以步骤2处理后的纹影图片中的某一固定点作为原点,建立参考坐标系,将每个所述标记光源的所述像素坐标转化为物理空间坐标变化,结合步骤3的结果,获得风洞实验中运动模型的分离运动轨迹以及动力学参数。
2.根据权利要求1所述的一种纹影图像的运动、参数识别方法,其特征在于,
步骤1中,每个所述标记光源的中轴线与所述纹影光源的光路平行实现每个所述标记光源的灯心安装方向与纹影光源的光路平行,且每个所述标记光源背向纹影光源安装在所述运动模型的壁面上用于标记运动模型的运动轨迹以及姿态。
3.根据权利要求1所述的一种纹影图像的运动、参数识别方法,其特征在于,
步骤2中所述获得每张纹影图像中每个标记光源的像素坐标的具体内容为:
获得包含有运动模型和标记光源的原始纹影图像,对所有所述原始纹影图像进行批量灰度处理得到灰度图片,以判定出标记光源的位置像素点与背景像素点;继续将所述灰度图片转化为黑白二值图片,以得到每个标记光源的灯心的像素点坐标;获得每张纹影图片中每个标记光源的像素坐标。
4.根据权利要求1所述的一种纹影图像的运动、参数识别方法,其特征在于,
步骤4中所述参考坐标系中,以风洞试验中试验气流的来流方向建立参考坐标系的横轴,垂直来流方向向上为参考坐标系的纵轴,得到每个所述标记光源在所述运动模型的运动过程中在所述参考坐标系中的坐标和所述运动模型的质心在参考坐标系中的坐标。
5.根据权利要求4所述的一种纹影图像的运动、参数识别方法,其特征在于,
在所述参考坐标系中,结合所述运动模型在某一时刻
Figure 61561DEST_PATH_IMAGE005
相对于上一时刻的
Figure DEST_PATH_IMAGE007
的位移,得到所述运动模型的运动速度、加速度、受力和俯仰力矩。
6.根据权利要求1所述的一种纹影图像的运动、参数识别方法,其特征在于,
步骤4中的所述动力学参数的精度随着标记光源的数量增多而增高。
7.根据权利要求1所述的一种纹影图像的运动、参数识别方法,其特征在于,
步骤1中,所述标记光源的数量大于等于3个。
8.一种纹影图像的运动、参数识别系统,其特征在于,包括:用于风洞试验的运动模型、纹影光源和纹影相机,在运动模型的壁面上固定设有3个及以上光线方向背对纹影光源且面向纹影相机镜头的标记光源,每个所述标记光源的中轴线与纹影光源的光路平行用于标记运动模型的运动轨迹,在运动模型的空腔内部固定有连接标记光源的供电电路;
风洞试验后,获得包含有运动模型和标记光源的原始纹影图像,基于标记光源在运动模型上相对于运动模型的空间位置坐标及运动模型的质心坐标分析所述原始纹影图像。
9.根据权利要求8所述的一种纹影图像的运动、参数识别系统,其特征在于,每个所述标记光源固定在运动模型壁面的一个贯穿孔内。
10.根据权利要求8-9任一项所述的一种纹影图像的运动、参数识别系统,其特征在于,所述标记光源是LED灯、小电珠或小灯泡中的一种。
CN202210658880.6A 2022-06-13 2022-06-13 一种纹影图像的运动、参数识别方法及系统 Active CN114757993B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210658880.6A CN114757993B (zh) 2022-06-13 2022-06-13 一种纹影图像的运动、参数识别方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210658880.6A CN114757993B (zh) 2022-06-13 2022-06-13 一种纹影图像的运动、参数识别方法及系统

Publications (2)

Publication Number Publication Date
CN114757993A true CN114757993A (zh) 2022-07-15
CN114757993B CN114757993B (zh) 2022-09-09

Family

ID=82336339

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210658880.6A Active CN114757993B (zh) 2022-06-13 2022-06-13 一种纹影图像的运动、参数识别方法及系统

Country Status (1)

Country Link
CN (1) CN114757993B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116929701A (zh) * 2023-09-15 2023-10-24 中国空气动力研究与发展中心低速空气动力研究所 一种翼型表面流动迹线测量方法及系统
CN118190338A (zh) * 2024-05-16 2024-06-14 中国航空工业集团公司沈阳空气动力研究所 一种适用于风洞的变体模型的动态变形测量方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030030820A1 (en) * 2001-06-29 2003-02-13 Samsung Corporation Marker used for measuring displacement of moving object and method therefor
CN103884486A (zh) * 2014-02-27 2014-06-25 中国科学院力学研究所 纹影测量成像系统及方法
CN106290256A (zh) * 2016-08-03 2017-01-04 中国空气动力研究与发展中心高速空气动力研究所 基于视频测量的定量背景纹影方法
CN107270875A (zh) * 2017-05-09 2017-10-20 南京航空航天大学 运动模糊效应下的视觉特征三维重建方法
CN112577708A (zh) * 2020-12-29 2021-03-30 中国航天空气动力技术研究院 一种基于背光成像的高焓激波风洞发动机推阻测量系统方法
CN112697388A (zh) * 2021-01-11 2021-04-23 中国空气动力研究与发展中心高速空气动力研究所 一种基于纹影图像的高超声速风洞模型姿态角测量方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030030820A1 (en) * 2001-06-29 2003-02-13 Samsung Corporation Marker used for measuring displacement of moving object and method therefor
CN103884486A (zh) * 2014-02-27 2014-06-25 中国科学院力学研究所 纹影测量成像系统及方法
CN106290256A (zh) * 2016-08-03 2017-01-04 中国空气动力研究与发展中心高速空气动力研究所 基于视频测量的定量背景纹影方法
CN107270875A (zh) * 2017-05-09 2017-10-20 南京航空航天大学 运动模糊效应下的视觉特征三维重建方法
CN112577708A (zh) * 2020-12-29 2021-03-30 中国航天空气动力技术研究院 一种基于背光成像的高焓激波风洞发动机推阻测量系统方法
CN112697388A (zh) * 2021-01-11 2021-04-23 中国空气动力研究与发展中心高速空气动力研究所 一种基于纹影图像的高超声速风洞模型姿态角测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YUNPENG WANG 等: ""Numerical Investigation of Supersonic Oscillatory Flow with Strong Interference over a Capsule-shaped Abort System"", 《TRANS. JAPAN SOC. AERO. SPACE SCI.》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116929701A (zh) * 2023-09-15 2023-10-24 中国空气动力研究与发展中心低速空气动力研究所 一种翼型表面流动迹线测量方法及系统
CN116929701B (zh) * 2023-09-15 2023-12-01 中国空气动力研究与发展中心低速空气动力研究所 一种翼型表面流动迹线测量方法及系统
CN118190338A (zh) * 2024-05-16 2024-06-14 中国航空工业集团公司沈阳空气动力研究所 一种适用于风洞的变体模型的动态变形测量方法
CN118190338B (zh) * 2024-05-16 2024-07-23 中国航空工业集团公司沈阳空气动力研究所 一种适用于风洞的变体模型的动态变形测量方法

Also Published As

Publication number Publication date
CN114757993B (zh) 2022-09-09

Similar Documents

Publication Publication Date Title
CN114757993B (zh) 一种纹影图像的运动、参数识别方法及系统
CN109272532B (zh) 基于双目视觉的模型位姿计算方法
CN109857144A (zh) 无人机、无人机控制系统及控制方法
CN111274959B (zh) 一种基于可变视场角的加油锥套位姿精确测量方法
CN106403827A (zh) 一种gis母线舱对地三维位移的测量装置和测量方法
CN110426178B (zh) 一种基于尾流示踪的风场测量方法及系统
CN108680165B (zh) 基于光学图像的目标飞行器姿态确定方法和装置
CN109781068B (zh) 面向空间应用的视觉测量系统地面模拟评估系统及方法
CN109063532A (zh) 一种基于无人机的野外失联人员搜寻方法
CN113237628A (zh) 一种低速风洞水平自由飞模型姿态测量方法
CN111551111A (zh) 一种基于标准球阵的零件特征机器人快速视觉定位方法
CN112197713A (zh) 直升机旋翼桨叶变形及表面压力同步测量装置及方法
CN118074435A (zh) 一种平面电机的动子位姿检测装置及位姿检测方法
CN106225668B (zh) 基于多折射模型的风洞投放物高速位姿测量方法
CN117629567A (zh) 一种高瞬时速度模型姿态角在线测量方法
CN101793905A (zh) 光纤式二维风速/风向测量装置及方法
CN212256370U (zh) 光学式运动捕捉系统
CN117972885A (zh) 一种基于仿真增强的空间智能感知数据生成方法
CN112859907A (zh) 少样本条件下基于三维特效仿真的火箭残骸高空检测方法
CN111862170A (zh) 光学式运动捕捉系统及方法
CN115359422A (zh) 一种高空抛物图像生成方法、装置及系统
Lindh Development and implementation of star tracker electronics
CN110006368B (zh) 一种基于视觉测量的绳索摆角测量方法及系统
CN109360241B (zh) 单相机测量三维起落架机轮中心位移方法
CN113358325A (zh) 一种风洞投放物位姿测量方法、装置、设备及存储介质

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