一种用于高超声速流场纹影显示中的激波定位方法技术

本发明专利技术公开了一种用于高超声速流场纹影显示中的激波定位方法，用于高超声速流场纹影显示技术；首先提取图像灰度值，采用归一化互相关函数算法，并基于设定的零点位置，计算发生了的像素偏差，最后将像素距离按比例换算成实际位移。本发明专利技术给出了灰度数组的相关性计算模型，能够精确描述激波在任一时刻发生的振荡位移偏差，实现了对大批量纹影试验数据的数字化定量处理，大幅提高纹影图像数据在该研究领域中的利用率。

A shock wave location method for schlieren display of hypersonic flow fields

The invention discloses a method for hypersonic shock wave position schlieren display, for hypersonic schlieren display technology; first extract the image gray value, the normalized cross-correlation function algorithm, and based on the zero position, calculate the pixel deviation occurs, the pixel distance according to the actual displacement scaling. The invention provides a calculation model of gray correlation array, which can accurately describe the displacement deviation of the oscillating shock occurred at a time, to realize the digital processing of large quantities of quantitative schlieren experimental data, greatly improve the utilization rate of the research in the field of schlieren image data.

【技术实现步骤摘要】
一种用于高超声速流场纹影显示中的激波定位方法
本专利技术涉及高超声速流场测量显示
，具体指代一种用于高超声速流场纹影显示中的激波定位方法。
技术介绍
高超声速飞行器是目前国内外的热点研究之一，随着美国X-43、X-51等高超声速飞行器相继进行的试飞试验，各航天大国争相发展高速临近空间飞行器。国内众多学者对高超声速飞行器也展开了大量研究。在高超声速
的研究过程中，流场结构的认识是必不可少的环节，激波之间的相互作用、复杂流场的相互干扰、激波串结构特性等均是目前的热点研究内容。在高超声速风洞试验研究中，流场监控和测量技术是认识流场结构的有效途径。其中纹影技术是较常用的流场显示技术之一。纹影测量属于非侵入式流场诊断技术，能够清晰、直观的显示流场特征，被广泛应用于高速流场的显示测量中。在高超声速进气道研究中，不起动现象是重要的学术问题之一。通过纹影及高速摄影技术，可以有效还原进气道从起动到不起动过程中流场结构的变化，对激波结构进行定性研究。然而通过对纹影图像的定量处理，如对激波位置的描述和定位，将可以从定量的角度分析激波振荡的位移特性及频率变化特征，更细节性地认识激波之间相互作用的机理及影响方式。
技术实现思路
针对于上述现有技术问题，本专利技术的目的在于提供一种用于高超声速流场纹影显示中的激波定位方法，将纹影图像中的激波位置，通过灰度提取并计算的方式给出精确的位移大小，实现大批量纹影数据的定量化数字处理，以提高对纹影图像数据的利用率。为达到上述目的，本专利技术的一种用于高超声速流场纹影显示中的激波定位方法，包括步骤：设定激波零点位置，沿一直线方向提取图像灰度值后采用归一化互相关函数算法，将计算得到的偏移像素值换算成实际发生的位移大小，实现对激波的位移定位；具体步骤如下：1)选取某一瞬时流场状态激波位置作为零点位置，设定图片中某一固定线段，保证该线段始终贯穿激波可能出现的任一位置；2)基于上述确定的线段位置及长度，提取所有纹影图片中该线段所扫过的灰度值数据，并生成二维数组；3)上述得到的二维数组按顺序将每一元素与定义零点位置的灰度元素进行计算任一偏差位移的归一化互相关函数，并挑选出使得归一化互相关性最大的偏差，将所有元素计算出的该偏差生成一组一维数组；4)将上述得到的一维数组乘以实际位移与像素位移的比率，得出任一时间激波发生的实际位移。本专利技术的有益效果：本专利技术基于图像灰度提取及归一化互相关函数算法的激波定位方式，能够精确定位任一时间刻激波偏移的位置，完成了纹影图像显示从定性分析到定量分析的转变，实现了海量纹影试验数据的数字化分析处理，大大提高纹影数据的研究利用率。附图说明图1为本实施例中纹影图片所定义的零点激波位置示意图；图2为激波偏移后在水平方向上产生的位移示意图；图3为实施例全部纹影图像进行激波定位的位移数据与同时采集的压力数据示意图；图4为图3数据在3s至3.5s之间的局部放大图。具体实施方式为了便于本领域技术人员的理解，下面结合实施例与附图对本专利技术作进一步的说明，实施方式提及的内容并非对本专利技术的限定。本专利技术的一种用于高超声速流场纹影显示中的激波定位方法，包括步骤：设定激波零点位置，沿一直线方向提取图像灰度值后采用归一化互相关函数算法，将计算得到的偏移像素值换算成实际发生的位移大小，实现对激波的位移定位；具体步骤如下：1)选取某一瞬时流场状态激波位置作为零点位置，设定图片中某一固定线段，保证该线段始终贯穿激波可能出现的任一位置；2)基于上述确定的线段位置及长度，提取所有纹影图片中该线段所扫过的灰度值数据，并生成二维数组；3)上述得到的二维数组按顺序将每一元素与定义零点位置的灰度元素进行计算任一偏差位移的归一化互相关函数，并挑选出使得归一化互相关性最大的偏差，将所有元素计算出的该偏差生成一组一维数组；4)将上述得到的一维数组乘以实际位移与像素位移的比率，得出任一时间激波发生的实际位移。将本专利技术方法应用于某一高超声速风洞试验采集得到的一组纹影图片中，图片数量大于3000张，具体实施步骤如下：(1)如图1，选取某一瞬时流场状态激波位置作为零点位置x0，设定图片中某一固定线段g0，保证该线段始终贯穿激波可能出现的任一位置；(2)基于上述确定的线段位置及长度，提取所有纹影图片中该线段所扫过的灰度值数据，并生成二维数组G[m][n]，其中m为纹影图像的总张数，n为每张图像所提取的灰度的数据量；(3)上述得到的二维数组G[m][n]按顺序将每一元素G[x][n]与定义零点位置的灰度元素G[0][n]进行计算任一偏差位移τ的归一化互相关函数：其中，E(Gx)、E(G0)是G[x][n]与G[0][n]的均值，σ(Gx)、σ(G0)是二者的标准差，挑选出使得归一化互相关性最大的偏差：将所有元素计算出的该偏差τmax生成一组一维数组SD[m]；(4)将上述得到的一维数组SD[m]乘以实际位移与像素位移的比率η，由此得出任一时间激波发生的实际位移SDx[t]。参照图2，本实施例中，x1是激波偏移后的位置，Δ1是激波发生的像素位移，某一纹影图像通过人为测量发现激波在水平方向上向前偏移55-59像素的距离，利用本专利技术的方法计算值为57。将全部纹影图像中激波发生的位移偏移量与同时采集的壁面压力做对比，如图3和4，在同一时间轴上，1是实施例中全部纹影图像定位出的激波位移数据，纵坐标代表像素位移大小，向前移动为负，向后移动为正；2是同时间段采集得到的壁面压力数据，纵坐标代表壁面绝对压力值大小；3是激波位移数据的局部展示，4是壁面压力数据的局部展示，3和4均表达了瞬态特征，随时间推移均发生周期性跳动；位移与压力有着一致的变化规律(压力数据为真实可靠数据，所发生的脉动代表真实的激波偏移，每个脉动都同时刻对应着计算得到的一个激波跳动位移)，本专利技术方法得出与压力一致的时域变化趋势，以上证实了本专利技术方法的准确性及实现大批量纹影图像处理的可行性。本专利技术具体应用途径很多，以上所述仅是本专利技术的优选实施方式，应当指出，对于本
的普通技术人员来说，在不脱离本专利技术原理的前提下，还可以作出若干改进，这些改进也应视为本专利技术的保护范围。本文档来自技高网...

【技术保护点】
一种用于高超声速流场纹影显示中的激波定位方法，其特征在于，包括步骤：设定激波零点位置，沿一直线方向提取图像灰度值后采用归一化互相关函数算法，将计算得到的偏移像素值换算成实际发生的位移大小，实现对激波的位移定位；具体步骤如下：1）选取某一瞬时流场状态激波位置作为零点位置，设定图片中某一固定线段，保证该线段始终贯穿激波可能出现的任一位置；2）基于上述确定的线段位置及长度，提取所有纹影图片中该线段所扫过的灰度值数据，并生成二维数组；3）上述得到的二维数组按顺序将每一元素与定义零点位置的灰度元素进行计算任一偏差位移的归一化互相关函数，并挑选出使得归一化互相关性最大的偏差，将所有元素计算出的该偏差生成一组一维数组；4）将上述得到的一维数组乘以实际位移与像素位移的比率，得出任一时间激波发生的实际位移。

【技术特征摘要】
1.一种用于高超声速流场纹影显示中的激波定位方法，其特征在于，包括步骤：设定激波零点位置，沿一直线方向提取图像灰度值后采用归一化互相关函数算法，将计算得到的偏移像素值换算成实际发生的位移大小，实现对激波的位移定位；具体步骤如下：1）选取某一瞬时流场状态激波位置作为零点位置，设定图片中某一固定线段，保证该线段始终贯穿激波可能出现的任一位置；2）基于上...

【专利技术属性】
技术研发人员：王成鹏，薛龙生，程克明，程川，焦运，赵俊淇，徐培，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏,32