采用仿射变换的桨叶变形摄像测量Newton-Raphson算法.pdfVIP

第33卷第3期 振动、测试与诊断 V01．33No．3 2013年6月 ournalof Jun．2013 J Vibration．MeasurementDiagnosis 采用仿射变换的桨叶变形摄像测量 Newton—Raphson算法。 高亚东， 刘荣明， 王华明 (南京航空航天大学直升机旋翼动力学国防科技重点实验室南京，210016) 摘要 变形前后的图像中选取匹配特征点，计算仿射变换参数；然后，根据仿射变换位移模式计算目标子区中心点的位 移，将其用作N—R算法的初值。为降低计算量，采用全场插值，插值运算次数降为常规插值的1％左右。仿真测试 表明，基于仿射变换的N—R算法在精度保持不变的情况下显著减少计算时间，克服了常规N—R算法计算量大、效 率低、只适用于无刚体转动目标等缺陷，为开发直升机桨叶变形摄像测量系统提供了算法基础。 关键词直升机；摄像测量；特征点匹配；仿射变换；Newton—Raphson算法 中图分类号P234．1；TP391 引 言 1 Newton—Raphson算法 桨叶变形的准确测量对直升机动力学精细化设 N—R算法考虑模板窗口的变形，可同时获得位 计具有重要价值。常用的电阻应变片法对被测物 移和应变，摒弃对目标子区零应变的假设。图1给 体一表面质量、温度、电磁干扰等因素有较高要求，且 出了图像子区在变形前后的映射关系，其中点P(x。， 难以在旋转桨叶上实现。探索利用数字图像相关技 Y。)为样本子区的中心，P7(z7。，Y7。)为P点在变形图 术实现非接触式桨叶变形测量具有重要的理论意义 像子区中的对应点。点Q(x，y)为图像子区内任一 与实用价值。在数字图像相关技术中，寻求图像间 点，它在变形后图像子区的对应点为d(z7，Y7)，对 最佳匹配的相关运算是核心。其中，如何提高计算 (z，y)变形后的坐标(z7，Y7)[21可表示为 精度降低计算量，是图像相关算法的研究重点。目 前，相关算法主要有双参数法、粗细搜索法、十字搜索 法、(Newton-Raphson，N—R)法等。由于双参数法和 粗细搜索法的计算量较大，应用较少。十字搜索法搜 索速度很快，但对亚像素搜索精度不高[1]。N-R算法 可同时获得位移和应变信息，需要较准确的初值估 计。通常用整像素搜索结果作为其初值，但整像素搜 索方法计算量大，对于有刚体转动的测量目标适用性 差。针对上述变形摄像测量中的难题，笔者尝试采用 人机交互的方法，在变形前后图像上选取三对匹配 的特征点，求得仿射变换的变换系数，根据仿射变换 的位移模式计算得到目标子区中心点的位移作为 N—R算法的初值，以克服整像素搜索方法的不足， 为桨叶变形摄像测量系统提供关键算法。 图1 变形前后子区示意图 * (2010ZD52044)；江苏高校优势学科建设工程资助项目 收稿日期：2011-11-11；修改稿收到日期：2013—02—20 万方数据 第3期 为N—R算法的初值。 r：一z+‰+如△z+‰△y (1) lY’=Y+uo+让△z+royAy 仿射变换包含了刚体平移、旋转、剪切和正应变 目标子区的位移和变形可用其中心点的位移， 等力学量，可用一阶形函数表示。平面仿射变换有 “。，可。表示图像子区中心点在X，Y方向的位移；Ax，6个自由度，对应有6个矩阵元素[4]，可表示为 △y分别为Q点和P点在．717，Y方向的距离；“：，“，， ㈤ u：，口。表示子区的位移梯度，综合反映了子区的伸