CN107085843A - 用于估计光学系统中的调制传递函数的系统和方法 - Google Patents
用于估计光学系统中的调制传递函数的系统和方法 Download PDFInfo
- Publication number
- CN107085843A CN107085843A CN201611207952.6A CN201611207952A CN107085843A CN 107085843 A CN107085843 A CN 107085843A CN 201611207952 A CN201611207952 A CN 201611207952A CN 107085843 A CN107085843 A CN 107085843A
- Authority
- CN
- China
- Prior art keywords
- multinomial
- edge
- mtf
- image
- optical system
- 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
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000003287 optical effect Effects 0.000 title claims abstract description 27
- 238000012546 transfer Methods 0.000 title claims abstract description 9
- 238000001514 detection method Methods 0.000 claims abstract description 10
- 238000004088 simulation Methods 0.000 claims abstract description 6
- 230000008859 change Effects 0.000 claims abstract description 5
- 230000006870 function Effects 0.000 claims description 16
- 230000015654 memory Effects 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 description 28
- 238000004590 computer program Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 7
- 239000000306 component Substances 0.000 description 5
- 238000005070 sampling Methods 0.000 description 5
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000011514 reflex Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 230000000712 assembly Effects 0.000 description 2
- 238000000429 assembly Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000012533 medium component Substances 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
- G01M11/0292—Testing optical properties of objectives by measuring the optical modulation transfer function
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0004—Industrial image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30168—Image quality inspection
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Quality & Reliability (AREA)
- Studio Devices (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了用于估计光学系统中的调制传递函数的系统和方法。一种针对图像确定调制传递函数(MTF)的方法包括:接收通过光学系统捕获的图像;对图像中的列或行执行边缘检测以计算多个边缘点;计算拟合于所计算的边缘点的多个多项式,多个多项式中的每个在次数上是变化的;从多个多项式中选择用于表示所检测的边缘的多项式;以及基于所选择的多项式估计所述MTF。
Description
背景技术
在光学系统的技术领域中,光学设计师通常基于称为系统的调制传递函数(MTF)的测量结果而在系统之间比较性能。MTF是适用于量化系统在分辨率和对比度方面的总体成像性能的最佳工具之一,并且被广泛使用于评估像球面单透镜那样简单的光学组件以及像多元远心成像透镜组那样复杂的光学组件。ISO12233标准是用于确定光学系统的MTF的最常用方法。然而,ISO12233无法在各种不同条件下进行准确的MTF估计,例如,当被应用于高度失真的图像时。
发明内容
根据所公开主题的实施例,一种针对图像确定调制传递函数(MTF)的新方法,包括:对图像中的行进行边缘检测以计算多个边缘点;计算拟合于所计算的边缘点的多个多项式,多个多项式中的每个在次数上有所变化;从多个多项式中选择多项式以表示检测到的边缘;以及基于选择的多项式来估计MTF。
根据所公开主题的另一个实施例,一种系统包括存储器,其存储指令;处理器,其被配置成基于指令而对由光学系统所捕捉的数字图像执行过程,过程包括以下的操作:对图像中的列或者行进行边缘检测以计算多个边缘点;计算拟合所计算的边缘点的多个多项式,多个多项式中的每个在多项式次数上有所变化;从多个多项式中选择表示所检测到的边缘的多项式;以及基于所选择的多项式来所述MTF。
根据所公开主题的实施例,一种用于针对图像确定调制传递函数(MTF)的方法,包括:对图像中的行进行边缘检测以计算多个边缘点;计算拟合所计算的边缘点的多个多项式,多个多项式中的每个在次数上有所变化;从多个多项式中选择表示所检测到的边缘的多项式;以及基于所选择的多项式来估计MTF。
所公开主题的附加特征、优点以及实施例可以考虑以下详细描述、附图以及权利要求来加以阐明。此外,将理解的是,前述发明内容以及以下详细说明是说明性的并且旨在提供进一步解释,而并不限制权利要求的范围。
附图说明
被包括以提供对所公开主题的进一步理解的附图,被结合并构成本说明书的一部分。附图还图示出所公开主题的实施例并且与详细描述一起用以阐释所公开主题的实施例的原理。并未意欲示出比所公开主题的基本理解以及可以实践其的各种方式所必需的更详细的结构细节。
图1示出合成地生成的圆形边缘。
图2示出根据所公开主题的实施例的对图1的突出部分的MTF计算和根据ISO12233算法的其他MTF计算以及校正MTF曲线(SINC),其是已知的,因为图1中的边缘是完美的边缘。
图3示出从数字单镜头反光(DSLR)相机所采集的图像。
图4示出根据所公开主题的实施例的对图3的突出部分的MTF计算和根据ISO12233标准的MTF计算。
图5示出从数字单镜头反光(DSLR)相机所采集的图像。
图6示出根据所公开主题的实施例的对图5的突出部分的MTF计算和根据ISO12233标准的MTF计算。
图7示出从数字单镜头反光(DSLR)相机所采集的图像。
图8示出根据所公开主题的实施例的对图7的突出部分的MTF计算和根据ISO12233标准的MTF计算。
图9示出从数字单镜头反光(DSLR)相机所采集的图像。
图10示出根据所公开主题的实施例的对图9的突出部分的MTF计算。
图11示出根据所公开主题的实施例的计算设备。
图12示出根据所公开主题的实施例的用于确定MTF的流程图过程。
具体实施方式
将参照附图对本公开的各个方面或特征予以描述,其中通篇中使用相似的附图标记来表示相似的元素。在本说明书中,阐明诸多细节,以便提供对本公开的全面理解。然而,应理解的是,所公开主题的某些方面可以在没有这些具体细节的情况下或者采用其他方法、组件、材料等来实践。在其他实例中,以框图的形式示出公知的结构和设备,以便于描述本主题公开。
调制传递函数(MTF)是用于测量图像质量的最重要参数之一。光学设计师和工程师经常引用MTF数据,特别是在成败视对特定对象进行成像的准确性而定的应用中。ISO12233标准在几种重要的成像条件下不能为光学系统提供准确的MTF结果,并且在一些实例中,完全不能产生任何结果。例如,在基于高度失真的图像来测量MTF中,该标准并不准确。例如,通过宽视场(WFOV)相机,时常捕捉到这种类型的图像。在WFOV相机中,由镜头所致的桶形失真使物面中的直边缘在像面中变得弯曲。失真的图像对主要依赖于像面中的直边缘来估计MTF的ISO12233标准是有问题的。如下所示,图像失真可能导致并发的问题严重到无法使用ISO12233标准。
有可能直接从检查和分析光学系统来测量MTF,然而,在许多实例中,光学系统不可用于检查,仅可访问由光学系统所摄取的图像。当图像过于失真而不能使用ISO12233时,一种常规的技术是使用计算机通过数字方式使图像去失真,并且将ISO12233应用于改变的图像。然而,数字方式的去失真引入赝像并且导致测量出加工的图像而非原始图像的MTF,这降低了计算的准确性。另一种常规技术是尝试捕捉某物的图像,使其在图像的桶形失真中所呈现的自然曲线相反的方向上弯曲,以便获得直边缘。然而,该方法耗费时间的并且需要反复试错。
所公开的发明通过提供一种基于先前常规方法使用非常有问题的图像来确定光学系统的高度准确的MTF的方法,提供用于光学系统技术领域的显著改进。本文中所公开的实施例提供了一种改进的新方法来确定MTF,甚至从失真的图像来确定MTF。创建每个多项式具有不同的多项式次数的多项式的集合来估计在图像中所检测到的边缘。从该集合当中,选择多项式并且将其使用于估计MTF,产生比常规技术更加准确的图形。
本文中所公开的技术将被称为自由形态(FF)型MTF算法。图12示出FF型MTF算法的示例流程图1200。在操作100中,对图像的每个行进行边缘检测,例如,使用诸如Canny边缘检测的技术或者修改的复杂度更低的Canny边缘检测技术。在修改的Canny边缘检测中,独立地对每个图像行运行边缘检测。对每个行求微分,x'(i)=x(i+1)-x(i)。结果数据是沿图像行的梯度。再用(Gaussian)低通滤波器使微分数据平滑。再用自适应噪声电平对数据设定阈值(将小于某个值的任何值均设定成0)。再将该设定阈值的数据归一化到0至1。针对计算MTF的感兴趣区域中的所有行,计算该Canny型边缘轮廓(在以上步骤中)。能够对图像列计算边缘,例如,以确定垂直的MTF。
在操作200,使用在操作100所估计的边缘点以生成一个或多个多项式,这些多项式被使用于提高边缘位置的准确性。在操作300,算法从一个或多个多项式的集合中选择多项式。这能够通过分析每个多项式拟合中的残差来完成。随着拟合边缘的多项式次数增加,残差减少。当将残差绘制为多项式次数的函数时,其迅速下降,然后趋于平缓,最终降至0。通过尝试选择没有过拟合的次数最高的多项式,选出用于估计边缘的多项式。这通过经由测量相继的多项式拟合的残差之差而估计何时发生过拟合来完成。当该值小于阈值量时,所述算法确定已经发生过拟合。所选择的多项式能够具有至少高于二次的次数。
应指出的是,针对第N行(以及隐含第N个点),N-1次多项式会完美地拟合全部点,残差为零。然而,完美拟合因数据的过拟合而不利于计算准确的MTF。为了避免选择过拟合的多项式,能够针对每个多项式计算MTF曲线,以确定过拟合的不利影响的程度。
在操作400,将每个行移位由多项式所确定的量,以创建超采样的边缘扩展函数。该超采样的边缘扩展函数能够与ISO12233标准兼容。在操作500,基于边缘的位置,合并(bin)并平均边缘。能够合并数据,使得超采样的线扩展函数具有高于图像中的原始行数据的分辨率。例如,分辨率能够是原始行数据的分辨率的至少四倍。在操作600,对超采样的线扩展函数求微分,通过Hamming滤波器对其进行加窗,最后对其进行Fourier变换和标准化,以确定MTF。在操作700,能够基于MTF来调整光学系统的配置。
下面在算法1中示出所公开的算法的伪代码的示例实施例。所公开的算法的计算复杂性是O(num columns*(num rows)2,其中,num columns是图像中的列号,并且num rows是图像中的行号:
算法1
图1示出合成地生成的圆形边缘,并且图2示出所公开的FF算法的结果、其与ISO12233相比如何以及其与理论上预期的MTF值sinc(xfx)=sin(πxfx)/πxfx相比如何。应指出,fx是空间频率并且x是像素采样距离。如图可见,ISO12233标准基于圆形边缘产生不准确的结果,而所公开的FF算法对ISO12233标准提供显著改进,其测量结果非常接近地反映理论值。
图3示出从DSLR相机所采集的图像。该图像中的灰色框示出用于计算图4中所示的MTF的倾斜的直边缘。在此情形下,ISO12233标准因图像中存在直边缘而可行。ISO12233标准以及所公开的FF算法产生相同的测量结果。
图5示出从相同的DSLR所采集的另一个图像。该图像中的灰色框示出表示失真直边缘的凸形边缘。在图6中示出从该边缘算出的MTF。在此,边缘的略微失真导致ISO12233标准对光学系统的MTF计算与所公开的FF算法的计算之间的差异。
图7示出从相同的佳能(Canon)DSLR所采集的另一个图像。该图像中的灰色框示出表示失真直边缘的凹形边缘。在图8中示出从该边缘算出的MTF。如图可见,ISO12233标准不适于基于凹形边缘而做出MTF判定,而与之相反,所公开的FF算法产生一致并且更加准确的结果。
图9示出从相同的DSLR相机所采集的另一个图像。该图像中的灰色框示出表示失真直边缘的多项式边缘。在图10中示出由FF算法所算出的该边缘的MTF。应指出,所公开的FF算法提供了一种用于评估光学系统的解决方案,尽管常规的ISO12233标准完全没有产生该边缘的结果。
在本文中,将“MTF50”灵敏度定义为在MTF等于0.5情况下的空间频率,并且将灵敏度定义为标准偏差除以正在测量的统计数值的均值。在图1中,计算的对边缘选择的MTF50灵敏度,即如果将虚线框上移或下移几十行,针对所公开的FF算法以及ISO12233标准均约为5%。然而,ISO12233标准低估MTF50约6倍的数量级。
在图3、5、7和9中,即用使用高端鱼眼镜头的DSLR相机所捕捉的图像的表现,针对所公开的FF算法,对边缘选择的灵敏度小于5%,而针对常规的ISO12233标准,变化则大于51%。此外,ISO12233标准低估MTF50约1.6倍的数量级。针对许多边缘,ISO12233没有产生结果。通过所公开的FF算法的改进克服ISO12233标准对于在这种类型的图像中产生一致的准确的结果的有问题的失败。
所公开的FF算法能够被使用于估计系统的MTF,即使针对包括任意平滑边缘的复杂失真的图像亦然。这允许跨视场表征光学系统的MTF。另外,所公开的算法跨直线且失真的样本图像是一致的,获得与ISO12233标准对直边缘的相同的结果。
与ISO12233标准相比,所公开的FF算法对于低亮度的,即低信噪比(SNR)的图像更具鲁棒性。与ISO12233相比,所公开的FF算法对图像中的噪点也更具鲁棒性。通过使用Canny型边缘估计算法或者降低图像中噪点数量并且减少图像数据集的其他边缘检测算法能够在处理噪点时改进该算法。
因为MTF是适用于量化光学系统的成像性能的最佳工具之一,了解系统内的每个镜头和传感器的MTF允许设计师例如在针对特定分辨率进行优化时做出适当的选择。宽视场(WFOV)系统的设计师现在能够在各种配置中使用所公开的FF算法来确定MTF并且准确地选择提供所需结果的组合。
本发明所公开主题的实施例可以通过并且使用多种组件和网络架构来实施。图11是适用于实施本发明所公开主题的实施例的示例计算设备20。设备20可以是例如桌面型计算机或膝上型计算机或者诸如智能电话、平板型计算机等移动计算设备。设备20可以包括总线21,其使计算机20的主要组件互连,诸如中央处理器24,诸如随机存取存储器(RAM)、只读存储器(ROM)、闪存RAM等的存储器27,诸如显示屏的用户显示器22,可以包括一个或多个控制器以及诸如键盘、鼠标、触摸屏等的相关用户输入设备的用户输入接口26,诸如硬盘驱动器、闪存等的固定存储器23,可操作成控制和接收光盘、闪存驱动器等的可移动介质组件25以及可操作成经由合适的网络连接与一个或多个远程设备通信的网络接口29。
总线21允许中央处理器24与一个或多个存储器组件之间的通信,如前所述,一个或多个存储器组件可以包括RAM、ROM以及其他存储器。典型地,RAM是向其中加载指令、操作系统以及应用程序的主存储器。ROM或闪存存储器组件除代码之外还能够包含基本输入输出系统(BIOS),其控制诸如与外围组件的交互的基本的硬件操作。与计算机20同驻的应用通常被存储在诸如硬盘驱动器(例如固定储存器23)、光驱、软盘或者其他存储介质的计算机可读介质上或者经由该计算机可读介质来访问。中央处理器24能够执行存储在一个或多个存储器组件上的指令。指令能够包括上述关于所公开的FF型MTF算法的操作。
固定储存器23可以与计算机20集成为一体或者可以是单独的并且通过其他接口来接入。网络接口29可以经由有线或无线连接而提供到远程服务器的直接连接。网络接口29可以使用本领域技术人员容易理解的任何合适的技术和协议来提供这样的连接,包括数字蜂窝电话、WiFi、蓝牙(R)、近场等。例如,网络接口29可以允许计算机经由一个或多个局域、广域或其他通信网络与其他计算机通信,如下更加详细地描述。
许多其他设备或组件(未示出)可以通过类似的方式来连接(例如,文档扫描仪、数码相机等等)。反之,在图11中所示的全部组件对于实施本公开而言不是必须存在的。组件能够通过不同于所示的方式互连。本领域技术人员容易理解诸如图11中所示的计算机的操作,在本申请中不再赘述。用于实施本公开的代码能够被存储在诸如存储器27、固定储存器23、可移动介质25中的一个或多个的计算机可读存储介质中,或者存储在远程储存位置上。
更一般地,本发明所公开主题的各种实施例可以包括或实现为计算机实施的过程的形式以及用于实践那些过程的装置。这些实施例还可以通过计算机程序产品的形式来实现,该计算机程序产品具有包含在诸如软盘、CD-ROM、硬盘驱动器、USB(通用串行总线)驱动器或者任何其他的机器可读存储介质的非暂时性和/或有形介质中实现的指令的计算机程序代码,所述非暂时性和/或有形介质,使得当将该计算机程序代码加载至计算机中并由该计算机执行时,该计算机成为用于实践所公开主题的实施例的装置。这些实施例也可以通过计算机程序代码的形式来实现,计算机程序代码例如或者被存储在存储介质中、加载至计算机中和/或由该计算机执行,或者通过光纤或者经由电磁辐射经诸如电线或电缆的一些传送介质来传送,使得当计算机程序代码被加载至计算机中并且由该计算机执行时,计算机成为用于实践所公开主题的实施例的装置。当在通用微处理器上实施时,计算机程序代码段配置微处理器以创建特定的逻辑电路。
在一些配置中,存储在计算机可读存储介质上的计算机可读的指令集可以由通用处理器来实施,该指令集可以将通用处理器或者包含通用处理器的设备转换成配置用于实施或实行指令的专用设备。这些实施例可以使用硬件来实施,该硬件可以包括处理器,诸如通用微处理器和/或专用集成电路(ASIC),其在硬件和/或固件中实现根据所公开主题的实施例的全部或部分技术。处理器可以被耦合至诸如RAM、ROM、快闪存储器、硬盘或者能够存储电子信息的任何其他设备的存储器。存储器可以存储适于由处理器执行以进行根据所公开主题的实施例的技术的指令。
为了解释目的,参照具体实施例来描述前述内容。然而,以上说明性讨论并非旨在穷举或者将所公开主题的实施例限定成所公开的确定形式。鉴于以上教导,可能有许多修改和变型。已对多个实施例加以选择和描述,以便阐释所公开主题的实施例的原理及其实际应用,由此使得本领域技术人员能够充分利用那些实施例以及可能适于特定用途的具有各种修改的各种实施例。
Claims (18)
1.一种基于通过光学系统捕获的图像来确定所述光学系统的调制传递函数(MTF)的方法,所述方法包括:
接收通过所述光学系统捕获的图像;
对所述图像中的列或行执行边缘检测以计算多个边缘点;
计算拟合于所计算的多个边缘点的多个多项式,所述多个多项式中的每一个在多项式次数上是变化的;
从所述多个多项式中选择用于表示所检测的边缘的多项式;以及
基于所选择的多项式来估计所述MTF。
2.根据权利要求1所述的方法,进一步包括基于所估计的MTF来配置光学系统配置。
3.根据权利要求1所述的方法,其中,选择所述多项式包括:
针对所述多个多项式中的每一个,确定残差值;以及
确定在所述多个多项式中的每一个之间的残差值的差;
确定所述多个多项式当中的、导致所述残差值的降低小于阈值量的最高次数的多项式。
4.根据权利要求2所述的方法,其中,所选择的多项式具有高于二次的次数。
5.根据权利要求1所述的方法,其中,估计所述MTF包括:
基于所选择的多项式来创建线扩展函数;以及
合并所述线扩展函数的边缘点数据。
6.根据权利要求5所述的方法,进一步包括:合并所述边缘点数据,使得所述线扩展函数具有是来自所述图像的原始数据至少四倍高的分辨率。
7.根据权利要求1所述的方法,其中,所述边缘是凹形边缘。
8.根据权利要求1所述的方法,其中,所述边缘是圆形边缘。
9.根据权利要求1所述的方法,其中,所述边缘是任意平滑边缘。
10.一种系统,包括:
存储器,所述存储器存储指令;
处理器,所述处理器被配置成基于所述指令来对由光学系统捕获的数字图像执行过程,所述过程包括操作:
接收通过所述光学系统捕获的图像;
对所述图像中的列或行执行边缘检测以计算多个边缘点;
计算拟合于所计算的多个边缘点的多个多项式,所述多个多项式中的每一个在多项式次数上是变化的;
从所述多个多项式中选择用于表示所检测的边缘的多项式;以及
基于所选择的多项式估计MTF。
11.根据权利要求10所述的系统,其中,选择所述多项式的操作包括:
针对所述多个多项式中的每一个,确定残差值;以及
确定在所述多个多项式中的每一个之间的残差值的差;
确定所述多个多项式当中的、导致所述残差值的降低小于阈值量的最高次数的多项式。
12.根据权利要求11所述的系统,其中,所选择的多项式具有高于二次的次数。
13.根据权利要求10所述的系统,其中,估计所述MTF的操作包括:
基于所选择的多项式来创建线扩展函数;以及
合并所述线扩展函数的边缘点数据。
14.根据权利要求13所述的系统,所述过程进一步包括合并所述边缘点数据,使得所述线扩展函数具有是来自所述图像的原始数据的至少高四倍高的分辨率。
15.根据权利要求10所述的系统,其中,所述边缘是凹形边缘。
16.根据权利要求10所述的系统,其中,所述边缘是凸形边缘。
17.根据权利要求10所述的系统,其中,所述边缘是圆形边缘。
18.根据权利要求10所述的系统,其中,所述边缘是任意平滑边缘。
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201662294786P | 2016-02-12 | 2016-02-12 | |
US62/294,786 | 2016-02-12 | ||
US15/290,313 | 2016-10-11 | ||
US15/290,313 US10067029B2 (en) | 2016-02-12 | 2016-10-11 | Systems and methods for estimating modulation transfer function in an optical system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107085843A true CN107085843A (zh) | 2017-08-22 |
CN107085843B CN107085843B (zh) | 2020-10-13 |
Family
ID=58284395
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611207952.6A Active CN107085843B (zh) | 2016-02-12 | 2016-12-23 | 用于估计光学系统中的调制传递函数的系统和方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US10067029B2 (zh) |
CN (1) | CN107085843B (zh) |
DE (2) | DE202016107422U1 (zh) |
GB (1) | GB2547318B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107024339B (zh) * | 2017-04-21 | 2023-10-20 | 小艾帮帮(杭州)科技有限公司 | 一种头戴显示设备的测试装置及方法 |
CN110097563A (zh) * | 2019-03-18 | 2019-08-06 | 深圳蓝韵医学影像有限公司 | 数字x线成像系统中调制传递函数的测算方法与装置 |
US11228699B2 (en) * | 2020-04-15 | 2022-01-18 | Labsphere, Inc. | Systems and methods for determining a minimum resolvable distance of an imaging system |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050276513A1 (en) * | 2004-05-25 | 2005-12-15 | Nokia Corporation | Method for characterizing a digital imaging system |
CN101635050A (zh) * | 2009-06-26 | 2010-01-27 | 武汉大学 | 一种图像复原方法 |
CN102779333A (zh) * | 2012-07-10 | 2012-11-14 | 武汉大学 | 一种基于卡尔曼滤波的光学影像复原方法 |
US20140205140A1 (en) * | 2013-01-24 | 2014-07-24 | Kineticor, Inc. | Systems, devices, and methods for tracking moving targets |
US20150093043A1 (en) * | 2013-10-02 | 2015-04-02 | National Applied Research Laboratories | Method of Evaluating Image Correlation with Speckle Patter |
CN104537646A (zh) * | 2014-12-12 | 2015-04-22 | 南京理工大学 | 遥感图像的多角度自动mtf估计方法 |
JP2015094701A (ja) * | 2013-11-13 | 2015-05-18 | 日本放送協会 | Mtf測定装置およびmtf測定プログラム |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
NL9001265A (nl) | 1990-06-05 | 1992-01-02 | Philips Nv | Werkwijze voor bepaling van een modulatie-overdrachtfunctie van een digitaal beeldvormend systeem. |
WO2010044939A1 (en) * | 2008-10-17 | 2010-04-22 | Honda Motor Co., Ltd. | Structure and motion with stereo using lines |
DE112013000938B4 (de) * | 2012-02-10 | 2021-04-22 | Zachary T. Bonefas | System und verfahren zum materialtransport mit einem oder mehreren bildgebungsgeräten |
EP3074949A2 (en) * | 2013-11-27 | 2016-10-05 | Universidad Politecnica De Madrid | Method and system for determining the prognosis of a patient suffering from pulmonary embolism |
-
2016
- 2016-10-11 US US15/290,313 patent/US10067029B2/en active Active
- 2016-12-19 GB GB1621558.4A patent/GB2547318B/en active Active
- 2016-12-23 CN CN201611207952.6A patent/CN107085843B/zh active Active
- 2016-12-28 DE DE202016107422.3U patent/DE202016107422U1/de active Active
- 2016-12-28 DE DE102016125758.3A patent/DE102016125758B4/de active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050276513A1 (en) * | 2004-05-25 | 2005-12-15 | Nokia Corporation | Method for characterizing a digital imaging system |
CN101635050A (zh) * | 2009-06-26 | 2010-01-27 | 武汉大学 | 一种图像复原方法 |
CN102779333A (zh) * | 2012-07-10 | 2012-11-14 | 武汉大学 | 一种基于卡尔曼滤波的光学影像复原方法 |
US20140205140A1 (en) * | 2013-01-24 | 2014-07-24 | Kineticor, Inc. | Systems, devices, and methods for tracking moving targets |
US20150093043A1 (en) * | 2013-10-02 | 2015-04-02 | National Applied Research Laboratories | Method of Evaluating Image Correlation with Speckle Patter |
JP2015094701A (ja) * | 2013-11-13 | 2015-05-18 | 日本放送協会 | Mtf測定装置およびmtf測定プログラム |
CN104537646A (zh) * | 2014-12-12 | 2015-04-22 | 南京理工大学 | 遥感图像的多角度自动mtf估计方法 |
Non-Patent Citations (2)
Title |
---|
JUNQING CHEN等: "Imaging sensor Modulation Transfer Function estimation", 《2010 IEEE INTERNATIONAL CONFERENCE ON IMAGE PROCESSING》 * |
陈起行等: "刃边法边缘扩散函数最优提取方法", 《激光与光电子学进展》 * |
Also Published As
Publication number | Publication date |
---|---|
DE102016125758B4 (de) | 2022-09-29 |
DE202016107422U1 (de) | 2017-05-24 |
GB2547318B (en) | 2020-08-26 |
US20170234765A1 (en) | 2017-08-17 |
GB201621558D0 (en) | 2017-02-01 |
CN107085843B (zh) | 2020-10-13 |
DE102016125758A1 (de) | 2017-08-17 |
US10067029B2 (en) | 2018-09-04 |
GB2547318A (en) | 2017-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pereyra | Proximal markov chain monte carlo algorithms | |
US9412030B2 (en) | Visualization of defects in a frame of image data | |
US8948513B2 (en) | Blurring based content recognizer | |
US10963676B2 (en) | Image processing method and apparatus | |
US9811885B2 (en) | Detecting glare in a frame of image data | |
CN108961184B (zh) | 一种深度图像的校正方法、装置及设备 | |
US9235063B2 (en) | Lens modeling | |
CN104637046B (zh) | 图像检测方法和装置 | |
JP2017130929A (ja) | 撮像装置により取得された文書画像の補正方法及び補正装置 | |
CN107424160A (zh) | 通过视觉系统查找图像中线的系统和方法 | |
CN105225222B (zh) | 对不同图像集的感知视觉质量的自动评估 | |
CN108764247B (zh) | 基于稠密连接的深度学习物体检测方法及装置 | |
US10852208B2 (en) | Quantifying gas leak rates using frame images acquired by a camera | |
CN111275660B (zh) | 一种平板显示器缺陷检测方法及装置 | |
US11164306B2 (en) | Visualization of inspection results | |
US10332244B2 (en) | Methods and apparatuses for estimating an ambiguity of an image | |
CN110738235A (zh) | 肺结核判定方法、装置、计算机设备及存储介质 | |
CN107085843A (zh) | 用于估计光学系统中的调制传递函数的系统和方法 | |
US8873839B2 (en) | Apparatus of learning recognition dictionary, and method of learning recognition dictionary | |
CN111369523A (zh) | 显微图像中细胞堆叠的检测方法、系统、设备及介质 | |
CN109035167A (zh) | 对图像中的多个人脸进行处理的方法、装置、设备和介质 | |
Halder et al. | Model-free prediction of atmospheric warp based on artificial neural network | |
CN115760653A (zh) | 图像校正方法、装置、设备及可读存储介质 | |
CN114663681A (zh) | 用于对指针式仪表进行读数的方法及其相关产品 | |
EP4174794A1 (en) | Anomaly detection in images |
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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: American California Applicant after: Google limited liability company Address before: American California Applicant before: Google Inc. |
|
GR01 | Patent grant | ||
GR01 | Patent grant |