CN104317352A - 一种自适应光学控制系统快速去倾斜分量处理方法 - Google Patents
一种自适应光学控制系统快速去倾斜分量处理方法 Download PDFInfo
- Publication number
- CN104317352A CN104317352A CN201410539953.5A CN201410539953A CN104317352A CN 104317352 A CN104317352 A CN 104317352A CN 201410539953 A CN201410539953 A CN 201410539953A CN 104317352 A CN104317352 A CN 104317352A
- Authority
- CN
- China
- Prior art keywords
- matrix
- slope
- restoration calculation
- row
- tilt component
- 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
Links
Abstract
本发明公开了一种自适应光学控制系统快速去倾斜分量处理方法,该方法将波前复原运算分解为复原运算1和复原运算2,然后将波前斜率计算与复原运算1安排为流水线结构,待斜率计算全部完成并得到倾斜量后,开始复原运算2,其间的矩阵相乘都将矩阵分割多个可以并行运算的模块单元,最后复原运算1的结果减去复原运算2结果即完成复原运算。本发明有效克服了过去自适应光学控制系统在斜率需要去倾斜分量时,斜率计算与复原运算不能实现流水线工作的困难,大大减小了复原运算的处理延时,节省了资源,降低了开发难度,有效提高了系统处理带宽。
Description
技术领域
本发明涉及一种自适应光学控制系统快速去倾斜分量处理方法,具体适用于自适应光学系统在斜率需要去倾斜分量的方式下减小波前复原的运算延时。
背景技术
在自适应光学系统中,为了有效校正大气湍流造成的动态波前误差,对波前处理机的计算能力和实时性要求非常高。通用的微型计算机已不能满足要求,必须根据自适应光学系统波前处理的特点研制专用的高速波前处理机。波前处理机的工作流程包括图像处理、斜率计算、复原运算、控制运算和D/A转换五个模块,随着子孔径数和驱动器单元数的不断增加,波前复原运算已经成为波前处理中运算量最大的部分。将子孔径波前斜率向量G乘波前复原矩阵D得到波前误差向量E,即:
E=DG (1)
设子孔径数为n,驱动器数为m,则D是m×2n的矩阵,G有2n个分量,为各子孔径的X和Y方向斜率的组合。
自适应光学系统在斜率需要去倾斜分量时,要先计算出倾斜分量即X方向斜率的平均值和Y方向斜率的平均值,G中各子孔径X和Y方向斜率要分别减去对应的和
目前的自适应光学系统在斜率不需要去倾斜分量时的主要处理方法是,将斜率计算与复原运算安排为流水线结构,当一部分子孔径计算完成以后即开始相应部分的复原运算,这样可以大大减小复原运算的处理延时。而在斜率需要去倾斜分量时的处理方法是,等待斜率全部计算完成以后,计算出倾斜分量和斜率减去和后再开始复原运算。这样复原运算无法采用流水线结构,运算耗时很长,资源消耗巨大。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提出一种用于自适应光学系统在斜率需要去倾斜分量的方式下减小波前复原运算延时的处理方法,该方法有效克服了过 去在斜率需要去倾斜分量的方式下斜率计算与复原运算不能实现流水线工作的困难,大大减小了复原运算的处理延时,节省了资源,降低了开发难度,有效提高了系统处理带宽。
本发明的技术方案如下:一种自适应光学控制系统快速去倾斜分量处理方法,将波前复原运算分解为复原运算1和复原运算2,然后将波前斜率计算与复原运算1安排为流水线结构,待斜率计算全部完成并得到倾斜量后,开始复原运算2,其间的矩阵相乘都将矩阵分割多个可以并行运算的模块单元,最后复原运算1的结果减去复原运算2结果即完成复原运算。
进一步的,所述的波前复原运算分解步骤如下:
①斜率矩阵为1列2N(N为子孔径数)行的矩阵,分解为1列2N行的原始斜率矩阵减去1列2N行的倾斜分量矩阵;
②复原运算即为M(驱动器单元数)行2N列的复原矩阵乘以斜率矩阵,可以分解为复原矩阵分别乘以原始斜率矩阵和倾斜分量矩阵的结果之差。
进一步的,所述的波前斜率计算与复原运算1的流水线结构为:以哈特曼传感器子孔径行为单位,当一个子孔径行的斜率计算完成后,随即开始该行斜率矩阵对应的复原运算;当下一行子孔径的斜率计算完成时,上一行子孔径斜率对应的复原运算已经完成。
进一步的,所述的复原运算2步骤如下:
①倾斜分量矩阵被分解为由两个倾斜量和重复组成的一列2N行的矩阵;
②复原矩阵乘以倾斜分量矩阵分解为复原矩阵每一行的奇数列乘以加上偶数列乘以
③由于复原矩阵可以提前得到的固定矩阵,所以将复原矩阵每一行的奇数列和偶数列分别累加得到这样复原矩阵简化为M行2列的矩阵;
④复原运算2即为M行2列的矩阵乘以2行1列的矩阵。
进一步的,所述的矩阵分割多个可以并行运算的模块单元是按行为单位,将复原矩阵分为多行并行进行乘累加。
本发明的原理在于:
一种自适应光学控制系统快速去倾斜分量处理方法:(1)波前复原运算分解,将复原矩阵与原始斜率减去倾斜分量后的斜率矩阵相乘的过程,分解为:复原矩阵与原始斜率矩阵的乘积(以下简称复原运算1)减去复原矩阵与倾斜分量矩阵的乘积(以下简 称复原运算2);(2)将波前斜率计算与复原运算1安排为以子孔径行图像为单位的流水线结构,一个子孔径行斜率计算完成后即开始进行对应部分的复原运算;(3)斜率计算全部完成,计算得到倾斜量后,开始复原运算2;(4)将复原运算2中的复原矩阵与倾斜分量矩阵相乘简化为:复原矩阵奇数列累加和、偶数列累加和组成的两列矩阵与倾斜量X、Y组成的两行矩阵相乘的过程;(5)矩阵分割,目的是将矩阵分割为多个可以并行运算的模块单元,进一步提高运算效率;(6)复原运算1的结果减去复原运算2的结果即得到整个复原运算的结果。
本发明与现有技术相比的优点在于:本发明通过一种巧妙的矩阵拆分和简化,有效克服了过去自适应光学系统在斜率需要去倾斜分量的方式下只能等待斜率计算全部完成以后才能开始复原运算所导致的运算量巨大,耗时很长的困难,尤其是解决了随着自适应光学系统子孔径数目和驱动器数目不断增加而导致复原矩阵运算量剧增,进一步导致波前处理的高实时性难于实现的难题。
附图说明
图1为本发明中的波前复原矩阵拆分示意图;
图2为本发明中的斜率计算与复原运算1流水线结构示意图;
图3为本发明中的复原运算2简化示意图。
具体实施方式
下面结合图1至图3具体说明本实施方式。本发明的一种自适应光学控制系统快速去倾斜分量处理方法,具体适用于自适应光学系统在斜率需要去倾斜分量的方式下减小波前复原运算延时,它通过以下步骤完成:
1、如图1所示,将波前复原运算分解为复原运算1和复原运算2,斜率矩阵为1列2N(N为子孔径数)行的矩阵,分解为1列2N行的原始斜率矩阵减去1列2N行的倾斜分量矩阵;复原运算即为M(驱动器单元数)行2N列的复原矩阵乘以斜率矩阵,可以分解为复原矩阵分别乘以原始斜率矩阵和倾斜分量矩阵的结果之差。
2、如图2所示,将波前斜率计算与复原运算1设计为流水线结构:以哈特曼传感器子孔径行为单位,当一个子孔径行的斜率计算完成后,随即开始该行斜率矩阵对应的复原运算;当下一行子孔径的斜率计算完成时,上一行子孔径斜率对应的复原运算已经完成。当一帧图像输出结束时,斜率计算延时仅为最后一个子孔径行的斜率计算延时ttitl,复原运算1的运算延时为最后一个子孔径行的复原运算延时tre1。
3、如图3所示,复原运算2通过以下步骤简化:
①倾斜分量矩阵被分解为由两个倾斜量和重复组成的一列2N行的矩阵;
②复原矩阵乘以倾斜分量矩阵分解为复原矩阵每一行的奇数列乘以加上偶数列乘以
③由于复原矩阵可以提前得到的固定矩阵,所以将复原矩阵每一行的奇数列和偶数列分别累加得到这样复原矩阵简化为M行2列的矩阵;
④复原运算2即为M行2列的矩阵乘以2行1列的矩阵,运算量大大减小,其运算延时为tre2。
4、矩阵分割多个可以并行运算的模块单元是按行为单位,将复原矩阵分为多行并行进行乘累加。
Claims (5)
1.一种自适应光学控制系统快速去倾斜分量处理方法,其特征在于:将波前复原运算分解为复原运算1和复原运算2,然后将波前斜率计算与复原运算1安排为流水线结构,待斜率计算全部完成并得到倾斜量后,开始复原运算2,其间的矩阵相乘都将矩阵分割多个可以并行运算的模块单元,最后复原运算1的结果减去复原运算2结果即完成复原运算。
2.根据权利要求1所述的一种自适应光学控制系统快速去倾斜分量处理方法,其特征在于:所述的波前复原运算1分解步骤如下:
①斜率矩阵为1列2N行的矩阵,分解为1列2N行的原始斜率矩阵减去1列2N行的倾斜分量矩阵,其中N为子孔径数;
②复原运算即为M行2N列的复原矩阵乘以斜率矩阵,其中M为驱动器单元数,可以分解为复原矩阵分别乘以原始斜率矩阵和倾斜分量矩阵的结果之差。
3.根据权利要求1所述的一种自适应光学控制系统快速去倾斜分量处理方法,其特征在于:所述的波前斜率计算与复原运算1的流水线结构为:以哈特曼传感器子孔径行为单位,当一个子孔径行的斜率计算完成后,随即开始该行斜率矩阵对应的复原运算;当下一行子孔径的斜率计算完成时,上一行子孔径斜率对应的复原运算已经完成。
4.根据权利要求1所述的一种自适应光学控制系统快速去倾斜分量处理方法,其特征在于:所述的复原运算2步骤如下:
①倾斜分量矩阵被分解为由两个倾斜量和重复组成的一列2N行的矩阵;
②复原矩阵乘以倾斜分量矩阵分解为复原矩阵每一行的奇数列乘以加上偶数列乘以
③由于复原矩阵可以提前得到的固定矩阵,所以将复原矩阵每一行的奇数列和偶数列分别累加得到这样复原矩阵简化为M行2列的矩阵;
④复原运算2即为M行2列的矩阵乘以2行1列的矩阵。
5.根据权利要求1所述的一种自适应光学控制系统快速去倾斜分量处理方法,其特征在于:所述的矩阵分割多个可以并行运算的模块单元是按行为单位,将复原矩阵分为多行并行进行乘累加。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410539953.5A CN104317352B (zh) | 2014-10-13 | 2014-10-13 | 一种自适应光学控制系统快速去倾斜分量处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410539953.5A CN104317352B (zh) | 2014-10-13 | 2014-10-13 | 一种自适应光学控制系统快速去倾斜分量处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104317352A true CN104317352A (zh) | 2015-01-28 |
CN104317352B CN104317352B (zh) | 2017-10-24 |
Family
ID=52372594
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410539953.5A Expired - Fee Related CN104317352B (zh) | 2014-10-13 | 2014-10-13 | 一种自适应光学控制系统快速去倾斜分量处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104317352B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106768395A (zh) * | 2016-12-14 | 2017-05-31 | 中国科学院光电技术研究所 | 一种自适应光学波前传感器与波前校正器对准误差精密测量方法 |
CN110231958A (zh) * | 2017-08-31 | 2019-09-13 | 北京中科寒武纪科技有限公司 | 一种矩阵乘向量运算方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102722217A (zh) * | 2012-06-05 | 2012-10-10 | 中国科学院光电技术研究所 | 一种千单元级自适应光学系统波前信号处理平台 |
CN103335950A (zh) * | 2013-06-24 | 2013-10-02 | 中国科学院光电技术研究所 | 一种测量大气湍流非等晕性波前误差及湍流特征参数的测量装置及方法 |
-
2014
- 2014-10-13 CN CN201410539953.5A patent/CN104317352B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102722217A (zh) * | 2012-06-05 | 2012-10-10 | 中国科学院光电技术研究所 | 一种千单元级自适应光学系统波前信号处理平台 |
CN103335950A (zh) * | 2013-06-24 | 2013-10-02 | 中国科学院光电技术研究所 | 一种测量大气湍流非等晕性波前误差及湍流特征参数的测量装置及方法 |
Non-Patent Citations (3)
Title |
---|
张浩田: "基于多核DSP的实时波前复原研究", 《中国优秀硕士学位论文全文数据库(电子期刊)基础科学辑》 * |
王春鸿: "61单元自适应光学系统实时波前处理技术研究", 《中国优秀硕士学位论文全文数据库(电子期刊)基础科学辑》 * |
贾建禄等: "基于FPGA的自适应光学系统波前处理机", 《光学精密工程》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106768395A (zh) * | 2016-12-14 | 2017-05-31 | 中国科学院光电技术研究所 | 一种自适应光学波前传感器与波前校正器对准误差精密测量方法 |
CN106768395B (zh) * | 2016-12-14 | 2019-10-08 | 中国科学院光电技术研究所 | 一种自适应光学波前传感器与波前校正器对准误差精密测量方法 |
CN110231958A (zh) * | 2017-08-31 | 2019-09-13 | 北京中科寒武纪科技有限公司 | 一种矩阵乘向量运算方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN104317352B (zh) | 2017-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105681628B (zh) | 一种卷积网络运算单元及可重构卷积神经网络处理器和实现图像去噪处理的方法 | |
CN111459877B (zh) | 基于FPGA加速的Winograd YOLOv2目标检测模型方法 | |
CN108205701B (zh) | 一种执行卷积计算的系统及方法 | |
CN103310443B (zh) | 高分辨率遥感影像快速处理方法及系统 | |
CN107085562B (zh) | 一种基于高效复用数据流的神经网络处理器及设计方法 | |
CN101763627B (zh) | 一种高斯模糊的实现方法和装置 | |
JP3228927B2 (ja) | プロセッサエレメント、プロセッシングユニット、プロセッサ、及びその演算処理方法 | |
CN107463990A (zh) | 一种卷积神经网络的fpga并行加速方法 | |
CN104035750A (zh) | 一种基于fpga的实时模板卷积实现方法 | |
CN106528490B (zh) | 一种fpga异构加速计算装置及系统 | |
CN108446253A (zh) | 一种针对神威体系架构的稀疏矩阵向量乘的并行计算方法 | |
CN104317352A (zh) | 一种自适应光学控制系统快速去倾斜分量处理方法 | |
CN103413273A (zh) | 一种基于gpu加速实现图像复原处理方法 | |
CN103455518A (zh) | 一种数据处理方法及装置 | |
CN113888524A (zh) | 缺陷检测模型训练方法、装置、设备及可读存储介质 | |
CN106228506A (zh) | 一种基于gpu的多帧图像并行处理的方法 | |
CN103902762A (zh) | 一种针对正定对称矩阵进行最小二乘方程求解的电路结构 | |
Takahashi et al. | Rank minimization approach to image inpainting using null space based alternating optimization | |
CN103279973A (zh) | 基于混合并行的立体图像匹配系统 | |
CN111260598A (zh) | 基于图像处理的Depthwise快速卷积系统及图像识别方法 | |
WO2023184754A1 (zh) | 可配置实时视差点云计算装置及方法 | |
CN110163793B (zh) | 卷积计算加速方法和装置 | |
CN110648313B (zh) | 基于fpga的激光条纹中心线拟合方法 | |
CN112580492A (zh) | 车辆检测方法及装置 | |
CN105373498A (zh) | 一种数据处理系统及方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171024 Termination date: 20211013 |