CN116847067A - 成像系统评估方法和装置 - Google Patents
成像系统评估方法和装置 Download PDFInfo
- Publication number
- CN116847067A CN116847067A CN202311099702.5A CN202311099702A CN116847067A CN 116847067 A CN116847067 A CN 116847067A CN 202311099702 A CN202311099702 A CN 202311099702A CN 116847067 A CN116847067 A CN 116847067A
- Authority
- CN
- China
- Prior art keywords
- value
- transfer function
- modulation transfer
- parameter
- preset
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 116
- 238000011156 evaluation Methods 0.000 title claims abstract description 71
- 238000012546 transfer Methods 0.000 claims abstract description 327
- 238000012360 testing method Methods 0.000 claims abstract description 162
- 238000001228 spectrum Methods 0.000 claims abstract description 151
- 238000000034 method Methods 0.000 claims abstract description 125
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 29
- 230000006870 function Effects 0.000 claims description 331
- 241001237160 Kallima inachus Species 0.000 claims description 37
- 238000004364 calculation method Methods 0.000 claims description 33
- 230000003595 spectral effect Effects 0.000 claims description 23
- 230000008859 change Effects 0.000 claims description 21
- 238000004590 computer program Methods 0.000 claims description 10
- 238000004458 analytical method Methods 0.000 claims description 9
- 230000001131 transforming effect Effects 0.000 claims description 5
- 230000009286 beneficial effect Effects 0.000 abstract description 9
- 230000008569 process Effects 0.000 description 33
- 238000012545 processing Methods 0.000 description 16
- 239000011159 matrix material Substances 0.000 description 14
- 238000010586 diagram Methods 0.000 description 11
- 230000009466 transformation Effects 0.000 description 5
- 238000011161 development Methods 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000012854 evaluation process Methods 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 238000011430 maximum method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 241000283715 Damaliscus lunatus Species 0.000 description 1
- 101100481875 Drosophila melanogaster topi gene Proteins 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000003631 expected effect Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N17/00—Diagnosis, testing or measuring for television systems or their details
- H04N17/002—Diagnosis, testing or measuring for television systems or their details for television cameras
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- General Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)
Abstract
本申请提供一种成像系统评估方法和装置,可以使得用于分析成像系统成像能力的第一调制传递函数曲线,更接近基于测试图像得到的多个第一调制传递函数取值的变化趋势,有利于减小误差、提高评估结果的可信度。具体包括:获取来自第二电子设备的测试图像;根据测试图像的第一频谱数据和测试图像的第一标准频谱数据,得到测试图像在频率域上的多个第一调制传递函数取值;利用第一算法对多个第一调制传递函数取值进行拟合,得到第一调制传递函数曲线,在第一调制传递函数曲线中,频率为零对应的调制传递函数取值为第一调制传递函数曲线的实际取值;对第一调制传递函数曲线进行分析,得到待评估成像系统的评估结果。
Description
技术领域
本申请涉及电子设备技术领域,尤其涉及一种成像系统评估方法和装置。
背景技术
随着电子设备的更迭发展,人们已经习惯于通过电子设备拍照、录制视频来记录生活,电子设备成像系统的能力可以通过图像分辨率、对比度等参数来体现,对用户观感而言可以表现为图像对所拍摄实物的还原程度、图像表达的纹理细节程度等。
可见,电子设备成像系统的能力会对用户的使用体验造成影响,故而厂商会在研发过程中或电子设备出厂之前对该电子设备的成像系统进行评估,以尽量保证出厂的电子设备的成像系统能力满足用户需求。
然而,现有的成像系统评估方法误差较大,得出的评估结果可信度低。
发明内容
本申请提供了一种成像系统评估方法和装置,可以使得用于分析成像系统成像能力的第一调制传递函数曲线,更接近基于测试图像得到的多个第一调制传递函数取值的变化趋势,有利于减小误差、提高评估结果的可信度。
第一方面,本申请提供一种成像系统评估方法,应用于第一电子设备,包括:获取来自第二电子设备的测试图像,所述测试图像是通过所述第二电子设备中的待评估成像系统拍摄预设图像得到的;根据所述测试图像的第一频谱数据和所述测试图像的第一标准频谱数据,得到所述测试图像在频率域上的多个第一调制传递函数取值;利用第一算法对所述多个第一调制传递函数取值进行拟合,得到第一调制传递函数曲线,所述第一调制传递函数曲线为调制传递函数取值与频率之间的关系曲线,在所述第一调制传递函数曲线中,频率为零对应的调制传递函数取值为所述第一调制传递函数曲线的实际取值;对所述第一调制传递函数曲线进行分析,得到所述待评估成像系统的评估结果。
在本申请实施例中,利用第一算法对多个第一调制函数取值进行拟合所得到的第一调制传递函数曲线,其在频率为零时对应的调制传递函数取值为该第一调制传递函数曲线的实际取值,该第一调制传递函数曲线,与基于测试图像得到的离散的多个第一调制传递函数取值在频率域上的变化趋势、以及各个频率所对应的调制传递函数取值均接近,基于该第一调制传递函数曲线进行分析得到的分析结果有利于反映测试图像的真实图像质量,即有利于反映待评估成像系统的真实情况,减小了对待评估成像系统的评估误差,使得评估结果可信度更高。
结合第一方面,在第一方面的某些实现方式中,所述第一标准频谱数据是基于第一参数的第一取值、第二参数的第一取值、所述测试图像的频谱坐标和所述测试图像的像素大小计算得到的;所述方法还包括:根据所述预设图像的第二频谱数据和所述预设图像的第二标准频谱数据,得到所述预设图像在频率域上的多个第二调制传递函数取值,所述第二标准频谱数据是基于所述第一参数的第一取值、所述第二参数的第一取值、所述预设图像的频谱坐标和所述预设图像的像素大小计算得到的;利用所述第一算法对所述多个第二调制传递函数取值进行拟合,得到第二调制传递函数曲线,所述第二调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于预设阈值。
应理解,理想条件下的第一调制传递函数曲线应当是调制传递函数取值为1的直线,在具体实现时并非如此的原因,除了待评估成像系统本身的性能因素之外,还可能是因为第一标准频谱数据是基于经验公式得到的,用该第一标准频谱数据表示测试图像真实的标准频谱数据存在误差。基于经验公式计算第一标准频谱数据虽然可以通过直接设置参数以使测试图像与预设图像的像素大小对等,以免去预设图像与测试图像像素不相等时复杂的像素配准过程,但是基于公式中的第一参数和第二参数的含义,其取值理应随待评估成像系统所参照的预设图像的不同而变化,故而第一参数的取值和第二参数取值相对于当前所采用的预设图像的准确性,可能会影响第一标准频谱数据计算结果的准确性,进而可能会对得到的第一调制传递函数曲线造成影响。
故而,在本申请实施例中,基于已知的预设图像,计算预设图像对应的第二调制传递函数曲线,并将第二调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于预设阈值时对应的、计算预设图像的第二标准频谱数据时用到的第一参数的第一取值和第二参数的第一取值,作为计算第一标准频谱数据时的第一参数的取值和第二参数的取值,以减小第一参数的取值和第二参数的取值选择不当造成的误差、提高后续成像系统的评估精度。同时,本申请实施例提供的方法,总能使得测试图像的第一标准频谱数据是基于当前所使用的预设图像的得到的,基于当前使用的预设图像校准第一参数和第二参数的取值,有利于避免因使用不同的预设图像而对同一个待评估成像系统造成误差,有利于提高成像系统评估结果的可信度。
结合第一方面,在第一方面的某些实现方式中,在根据所述测试图像的第一频谱数据和所述测试图像的第一标准频谱数据,得到所述测试图像在频率域上的多个第一调制传递函数取值之前,所述方法还包括:根据所述预设图像的第二频谱数据和所述预设图像的第三标准频谱数据,得到所述预设图像在频率域上的多个第三调制传递函数取值,所述第三标准频谱数据是基于所述第一参数的第二取值、所述第二参数的第二取值、所述预设图像的频谱坐标和所述预设图像的像素大小计算得到的;利用所述第一算法对所述多个第三调制传递函数取值进行拟合,得到第三调制传递函数曲线;在所述第三调制传递函数曲线在所述预设频率范围内的调制传递函数取值与1的差值的绝对值存在大于所述预设阈值的情况下,调整所述第一参数的取值,得到所述第一参数的第一取值;和/或,调整所述第二参数的取值,得到所述第二参数的第一取值。
在本申请实施例中,在计算测试图像的第一标准频谱数据之前,若预设图像的第三调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值存在大于预设阈值的情况,则需要对第一参数和/或第二参数的取值进行调整。本申请实施例提供的方法,可以在初始计算的预设图像的第三调制传递函数曲线不满足上述条件的情况下,通过调整第一参数和第二参数使其满足该条件后,再将此时的第一参数的第一取值和第二参数的第二取值用于第一标准频谱数据的计算,有利于减小因参数选取不当而造成的误差,有利于提高成像系统评估的精度。
结合第一方面,在第一方面的某些实现方式中,在所述第三调制传递函数曲线中,频率为零对应的调制传递函数取值为所述第三调制传递函数曲线的实际取值。
结合第一方面,在第一方面的某些实现方式中,判断所述第三调制传递函数曲线在预设频率范围内的变化率是否为零;若否,则调整所述第一参数的取值,以更新所述第三调制传递函数曲线;确定更新的第三调制传递函数曲线在所述预设频率范围内的变化率为零时所述第一参数的目标取值,并将所述第一参数的目标取值确定为所述第一参数的第一取值;判断所述更新的第三调制传递函数曲线中,预设频率的调制传递函数取值与1的差值的绝对值是否小于或等于所述预设阈值,所述预设频率在所述预设频率范围内;若否,则调整所述第二参数的取值,以再次更新所述第三调制传递函数曲线;确定再次更新的第三调制传递函数曲线中所述预设频率的调制传递函数取值与1的差值的绝对值小于或等于所述预设阈值时所述第二参数的目标取值,并将所述第二参数的目标取值确定为所述第二参数的第一取值。
结合第一方面,在第一方面的某些实现方式中,所述方法还包括:若所述更新的第三调制传递函数曲线中,所述预设频率的调制传递函数取值与1的差值的绝对值小于或等于所述预设阈值,则将所述第二参数的第二取值确定为所述第二参数的第一取值。
结合第一方面,在第一方面的某些实现方式中,所述方法还包括:若所述第三调制传递函数曲线在所述预设频率范围内的变化率为零,则将所述第一参数的第二取值确定为所述第一参数的第一取值;调整所述第二参数的取值,以更新所述第三调制传递函数曲线;确定更新后的第三调制传递函数曲线中所述预设频率的调制传递函数取值与1的差值的绝对值小于或等于所述预设阈值时所述第二参数的目标取值,并将所述第二参数的目标取值确定为所述第二参数的第一取值。
应理解,在第三调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值存在大于预设阈值的情况下,上述“在预设频率范围内的变化率为零”和“预设频率对应的调制传递函数取值与1的差值的绝对值小于或等于预设阈值”的两级条件,该第三调制传递函数曲线至少不满足其一。
在本申请实施例中,将“在预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于预设阈值”的判断条件拆解为“在预设频率范围内的变化率为零”和“预设频率对应的调制传递函数取值与1的差值的绝对值小于或等于预设阈值”的两级条件,有利于快速定位可调整的参数,节约设备算力,减小参数调整过程所消耗的时间,有利于提高整个成像系统评估过程的效率。
结合第一方面,在第一方面的某些实现方式中,预设所述第三调制传递函数曲线中频率为零对应的调制传递函数取值为1。
结合第一方面,在第一方面的某些实现方式中,调整所述第一参数的取值和所述第二参数的取值,以更新所述第三调制传递函数曲线;确定更新的第三调制传递函数曲线在所述预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于所述预设阈值时所述第一参数的目标取值和所述第二参数的目标取值,并将所述第一参数的目标取值确定为所述第一参数的第一取值、所述第二参数的目标取值确定为所述第二参数的第一取值。
本申请提供的方法,在第三调制传递函数曲线中频率为零对应的调制传递函数取值为1时,也能实现对第一参数和第二参数的调整,使得更新的第三调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于预设阈值,有提高获得测试图像的第一标准频谱数据的准确性。
结合第一方面,在第一方面的某些实现方式中,所述方法还包括:对所述测试图像进行傅里叶变换,得到所述测试图像的第一频谱数据。
应理解,通过傅里叶变换,可以将图像从基于像素坐标和灰度值描述的空间域,转换至基于空间频率、频谱坐标和频谱强度描述的频率域,以便于后续对图像质量的分析、对成像系统的评估。
可选地,本申请实施例所使用的傅里叶变换可以包括快速傅里叶变换、二维傅里叶变换等任意可实现本申请的傅里叶变换的变形或拓展,本申请所采用的傅里叶变换的形式不作具体限定。
结合第一方面,在第一方面的某些实现方式中,在所述测试图像为彩色图像的情况下,在对所述测试图像进行傅里叶变换之前,所述方法还包括:将所述测试图像转换为灰度图像;所述对所述测试图像进行傅里叶变换,得到所述测试图像的第一频谱数据,包括:对所述灰度图像进行傅里叶变换,得到所述测试图像的第一频谱数据。
可选地,可以采用加权法、均值法和最大值法中的任一个,将所述测试图像转换为灰度图像。
在本申请实施例中,将彩色图像转换为灰度图像之后,可以将彩色的测试图像中每个像素的R、G、B三个亮度值转换为一个灰度值,可以减小后续傅里叶变换的次数,简化计算过程,提高成像系统评估效率。
结合第一方面,在第一方面的某些实现方式中,所述对所述第一调制传递函数曲线进行分析,得到所述待评估成像系统的评估结果,包括:采用纹理锐度值分析法、十分之一幅度法和半功率法中的任一个,对所述第一调制传递函数曲线进行分析,得到所述待评估成像系统的评估结果。
结合第一方面,在第一方面的某些实现方式中,所述预设图像为枯叶图。
应理解,任意具有丰富纹理细节的图像均可以作为预设图像来对成像系统进行评估,本申请对预设图像的具体图像类型和内容不做具体限定。
结合第一方面,在第一方面的某些实现方式中,所述第一算法为最小二乘法、多项式拟合法、样条差值法和非参数拟合法中的任一个。
第二方面,本申请提供一种成像系统评估方法,应用于第一电子设备,包括:根据预设图像的第二频谱数据和所述预设图像的第二标准频谱数据,得到所述预设图像在频率域上的多个第二调制传递函数取值,所述第二标准频谱数据是基于第一参数的第一取值、第二参数的第一取值、所述预设图像的频谱坐标和所述预设图像的像素大小计算得到的;利用第一算法对所述多个第二调制传递函数取值进行拟合,得到第二调制传递函数曲线,所述第二调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值小于或等于预设阈值;获取来自第二电子设备的测试图像,所述测试图像是通过所述第二电子设备中的待评估成像系统拍摄所述预设图像得到的;根据所述测试图像的第一频谱数据和所述测试图像的第一标准频谱数据,得到所述测试图像在频率域上的多个第一调制传递函数取值,所述第一标准频谱数据是基于第一参数的第一取值、第二参数的第一取值、所述测试图像的频谱坐标和所述测试图像的像素大小计算得到;利用所述第一算法对所述多个第一调制传递函数取值进行拟合,得到第一调制传递函数曲线,所述第一调制传递函数曲线为调制传递函数取值与频率之间的关系曲线;对所述第一调制传递函数曲线进行分析,得到所述待评估成像系统的评估结果。
在本申请实施例中,在对测试图像进行计算之前,先基于预设图像得到第二调制传递函数曲线计算时所用到的第一参数的第一取值、第二参数的第二取值,且该第二调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值小于或等于预设阈值。本申请提供的方法,可以减小因为第一参数和第二参数的选取不当所造成的第一标准频谱数据的误差,也有利于保证在使用不同的预设图像对同一待评估成像系统进行评估时第一参数和第二参数取值的合理性,有利于保障成像系统评估结果的准确性。
结合第二方面,在第二方面的某些实现方式中,在根据所述测试图像的第一频谱数据和所述测试图像的第一标准频谱数据,得到所述测试图像在频率域上的多个第一调制传递函数取值之前,所述方法还包括:根据所述预设图像的第二频谱数据和所述预设图像的第三标准频谱数据,得到所述预设图像在频率域上的多个第三调制传递函数取值,所述第三标准频谱数据是基于所述第一参数的第二取值、所述第二参数的第二取值、所述预设图像的频谱坐标和所述预设图像的像素大小计算得到的;利用所述第一算法对所述多个第三调制传递函数取值进行拟合,得到第三调制传递函数曲线;在所述第三调制传递函数曲线在所述预设频率范围内的调制传递函数取值与1的差值的绝对值存在大于所述预设阈值的情况下,调整所述第一参数的取值,得到所述第一参数的第一取值;和/或,调整所述第二参数的取值,得到所述第二参数的第一取值。
应理解,本申请实施例中,所述“调整所述第一参数的取值,得到所述第一参数的第一取值;和/或,调整所述第二参数的取值,得到所述第二参数的第一取值”的方式及有益效果与上述第一方面相同,在此不再赘述。
第三方面,本申请提供一种成像系统评估装置,包括用于实现上述第一方面或第二方面的方法的模块。
第四方面,本申请提供了一种电子设备,包括:处理器和存储器,该处理器用于读取存储器中存储的指令,以执行上述第一方面或第二方面中的任一种可能实现方式中的方法。
可选地,处理器为一个或多个,存储器为一个或多个。
可选地,存储器可以与处理器集成在一起,或者存储器与处理器分离设置。
在具体实现过程中,存储器可以为非瞬时性(non-transitory)存储器,例如只读存储器(read only memory,ROM),其可以与处理器集成在同一块芯片上,也可以分别设置在不同的芯片上,本申请实施例对存储器的类型以及存储器与处理器的设置方式不做限定。
上述第四方面中的电子设备可以是一个芯片,该处理器可以通过硬件来实现也可以通过软件来实现,当通过硬件实现时,该处理器可以是逻辑电路、集成电路等;当通过软件来实现时,该处理器可以是一个通用处理器,通过读取存储器中存储的软件代码来实现,该存储器可以集成在处理器中,可以位于该处理器之外,独立存在。
第五方面,本申请提供了一种计算机可读存储介质,该计算机可读存储介质存储有计算机程序(也可以称为代码,或指令)当其在计算机上运行时,使得计算机执行上述第一方面或第二方面中的任一种可能实现方式中的方法。
第六方面,本申请提供了一种计算机程序产品,计算机程序产品包括:计算机程序(也可以称为代码,或指令),当计算机程序被运行时,使得计算机执行上述第一方面或第二方面中的任一种可能实现方式中的方法。
附图说明
图1为基于现有技术得到的测试图像的调制传递函数曲线示意图;
图2为本申请所适用的场景图;
图3为本申请实施例提供的一种成像系统评估方法的示意性流程图;
图4为本申请实施例提供的测试图像的第一调制传递函数曲线示意图;
图5为本申请实施例提供的另一种成像系统评估方法的示意性流程图;
图6为本申请实施例提供的一种成像系统评估方法的第一阶段的示意性流程图;
图7为本申请实施例提供的一种成像系统评估方法的第二阶段的示意性流程图;
图8为本申请实施例提供的预设图像的调制传递函数曲线在参数调整过程中的变化示意图;
图9为本申请实施例提供的一种成像系统评估方法的第三阶段的示意性流程图;
图10为本申请实施例提供的测试图像的调制传递函数曲线与其真实调制传递函数曲线的示意性对比图;
图11为本申请实施例提供的另一种成像系统评估方法的示意性流程图;
图12为本申请实施例提供的一种成像系统评估装置的示意性框图;
图13为本申请实施例提供的另一种成像系统评估装置的示意性框图。
具体实施方式
下面将结合附图,对本申请中的技术方案进行描述。
为了便于清楚描述本申请实施例的技术方案,在本申请的实施例中,采用了“第一”、“第二”等字样对功能和作用基本相同的相同项或相似项进行区分。本领域技术人员可以理解“第一”、“第二”等字样并不对数量和执行次序进行限定,并且“第一”、“第二”等字样也并不限定一定不同。
需要说明的是,本申请实施例中,“示例性地”或者“例如”等词用于表示作例子、例证或说明。本申请中被描述为“示例性地”或者“例如”的任何实施例或设计方案不应被解释为比其他实施例或设计方案更优选或更具优势。确切而言,使用“示例性地”或者“例如”等词旨在以具体方式呈现相关概念。
本申请实施例中,“至少一个”是指一个或者多个,“多个”是指两个或两个以上。“和/或”,描述关联对象的关联关系,表示可以存在三种关系,例如,A和/或B,可以表示:单独存在A,同时存在A和B,单独存在B的情况,其中A,B可以是单数或者复数。字符“/”一般表示前后关联对象是一种“或”的关系。“以下至少一项(个)”或其类似表达,是指的这些项中的任意组合,包括单项(个)或复数项(个)的任意组合。例如,a,b,或c中的至少一项(个),可以表示:a,b,c,a-b,a--c,b-c,或a-b-c,其中a,b,c可以是单个,也可以是多个。
随着电子设备的更迭发展,人们已经习惯于通过电子设备拍照、录制视频来记录生活,电子设备成像系统的能力可以通过图像分辨率、对比度等参数来体现,对用户观感而言可以表现为图像对所拍摄实物的还原程度、图像表达的纹理细节程度等。
可见,电子设备成像系统的能力会对用户的使用体验造成影响,故而厂商会在研发过程中或电子设备出厂之前对该电子设备的成像系统进行评估,以尽量保证出厂的电子设备的成像系统能力满足用户需求。
调制传递函数(Modulation Transfer Function,MTF)可以用于描述成像系统对不同频率的图像细节的传递能力,故而经常被用于评估成像系统的成像能力。在通过调制传递函数对成像系统进行评估时,需要基于该成像系统所拍摄的测试图像计算不同频率对应的多个调制传递函数取值,并基于该多个调制传递函数取值拟合出随频率变化的调制传递函数曲线,进而对调制传递函数曲线进行分析以得到评估结果。但是现有技术中,在拟合测试图像的调制传递函数曲线时,预设了该曲线中频率为零时的调制传递函数的取值为1,得到的调制传递函数曲线如图1所示。参照图1可以看出,拟合出的测试图像的调制传递函数曲线与测试图像的多个调制传递函数取值的原始数据的变化趋势在低频阶段差异较大,这就可能导致基于该测试图像的调制传递函数曲线所进行的、关于成像系统能力的分析误差较大,进而导致得到的评估结果可信度较低。
有鉴于此,本申请提供一种成像系统评估方法,采用第一算法对待评估成像系统所拍摄的测试图像的多个第一调制传递函数取值进行拟合,拟合得到的第一调制传递函数曲线中,频率为零对应的调制传递函数取值为该第一调制传递函数曲线的实际取值,使得拟合得到的第一调制传递函数曲线更接近原始的多个第一调制传递函数取值的变化趋势,有利于减小得到该第一调制传递函数曲线的过程对成像系统能力评估造成的误差,有利于提高评估结果的可信度。
图2为本申请所适用的应用场景200的示意图,至少包括第一电子设备201和第二电子设备202。其中,第二电子设备202包括待评估的成像系统,第一电子设备201中包括用于实现本申请所提供的技术方案的计算机程序代码。
在一种可能的实现方式中,第二电子设备202可以实现测试图像的自动拍摄,并将所拍摄的测试图像传输至第一电子设备201,第一电子设备201可以运行所存储的计算机程序代码,以通过第二电子设备所拍摄的测试图像对其成像系统进行评估。
在一种可能的实现方式中,上述场景中还可以包括灯光控制设备(图2中未示出),该灯光控制设备包括自动控制灯光的程序代码,可以根据预设的灯光参数自动调节灯光,以使该第二电子设备202所拍摄的测试图像曝光度达到预期效果。
可选地,第一电子设备201可以是笔记本电脑、工作站、服务器等任意运算能力可满足本申请所提供的技术方案的程序代码运行要求的设备,本申请对此不作限定。
可选地,第二电子设备202可以是手机、平板电脑、照相机、摄像机等任意具有成像功能的设备,本申请对此也不作限定。
下面结合图3至图10对本申请提供的一种成像系统评估方法进行详细说明。
图3为本申请实施例提供的一种成像系统评估方法300的示意性流程图。该方法300由第一电子设备执行,第一电子设备可以是笔记本电脑、工作站、服务器等,但本申请对此不做限定。方法300具体包括如下步骤:
S301、获取来自第二电子设备的测试图像,测试图像是通过第二电子设备中的待评估成像系统拍摄预设图像得到的。
S302、根据测试图像的第一频谱数据和测试图像的第一标准频谱数据,得到测试图像在频率域上的多个第一调制传递函数取值。
S303、利用第一算法对多个第一调制传递函数取值进行拟合,得到第一调制传递函数曲线,第一调制传递函数曲线为调制传递函数取值与频率之间的关系曲线,在第一调制传递函数曲线中,频率为零对应的调制传递函数取值为第一调制传递函数曲线的实际取值。
S304、对第一调制传递函数曲线进行分析,得到待评估成像系统的评估结果。
应理解,第二电子设备可以是手机、平板电脑等包括待评估成像系统的设备,也可以是待评估成像系统本身,本申请对此不作限定。
在本申请实施例中,利用第一算法对多个第一调制函数取值进行拟合所得到的第一调制传递函数曲线,其在频率为零时对应的调制传递函数取值为该第一调制传递函数曲线的实际取值,如图4所示,图4中所呈现的第一调制传递函数曲线,与拟合前的多个离散的第一调制传递函数取值在整个频率域上的均有较高的重合度。通过本申请提供的方法得到的第一调制传递函数曲线,与基于测试图像得到的离散的多个第一调制传递函数取值在频率域上的变化趋势、以及各个频率所对应的调制传递函数取值均更接近,基于该第一调制传递函数曲线进行分析得到的分析结果更有利于反映测试图像的真实图像质量,即有利于反映待评估成像系统的真实情况,减小了对待评估成像系统的评估误差,使得评估结果可信度更高。
作为一个可选的实施例,上述S302中的第一频谱数据是对测试图像进行傅里叶变换得到的。
在一种可能的实现方式中,若测试图像为彩色图像,则将该测试图像的R、G、B三个通道的图像数据均进行傅里叶变换后,以得到测试图像的第一频谱数据。
在另一种可能的实现方式中,若测试图像为彩色图像,则在对测试图像进行傅里叶变换之前,则采用加权法、均值法和最大值法中的任一个,将测试图像转换为灰度图像,再对灰度图像进行傅里叶变换,得到第一频谱数据。
将彩色图像转换为灰度图像之后,可以将彩色的测试图像中每个像素的R、G、B三个亮度值转换为一个灰度值,可以减小后续傅里叶变换的次数,简化计算过程,提高成像系统评估效率。
应理解,通过傅里叶变换,可以将图像从基于像素坐标和灰度值描述的空间域,转换至基于空间频率、频谱坐标和频谱强度描述的频率域,以便于后续对图像质量的分析、对成像系统的评估。经过傅里叶变换所得到的第一频谱数据,是与表示测试图像灰度值的矩阵同型的矩阵形式。
还应理解,本申请实施例中所涉及的频率指的是图像的频率,也可称为空间频率,反映图像的像素对应的灰度值在不同像素位置的变化情况。
可选地,本申请实施例所使用的傅里叶变换可以包括快速傅里叶变换、二维傅里叶变换等任意可实现本申请的傅里叶变换的变形或拓展,本申请所采用的傅里叶变换的形式不作具体限定。
在一种可能的实现方式中,为了节约设备算力,简化计算流程,提高整个成像系统评估过程的效率,上述S302中的第一标准频谱数据是基于第一参数的第一取值、第二参数的第一取值、测试图像的频谱坐标和测试图像的像素大小,根据如下公式(1)和公式(2)计算得到。
(1)
(2)
其中,第一参数为经验参数,第二参数/>为描述图像中各个像素的灰度值大小差异的参数,m和n分别表示图像的频谱坐标的横轴的取值和纵轴的取值,N为表示图像的像素大小的参数。
应理解,在计算测试图像的第一标准频谱数据时,上述公式(2)中N的取值为测试图像的像素大小,m的取值和n的取值分别对应测试图像的频谱坐标的取值范围,而由于测试图像是待评估成像系统拍摄预设图像得到的,其第一标准频谱数据应与预设图像相关,在计算第一标准频谱数据时,第一参数和第二参数/>是与预设图像相关的参数。计算测试图像的第一标准频谱数据需要遍历该测试图像的每一个频谱坐标,得到的第一标准频谱数据为与上述第一频谱数据同型的矩阵。
还应理解,公式(1)和公式(2)为计算图像频谱数据的经验公式,根据公式(1)和公式(2)计算得到的频谱数据可以认为是图像的标准频谱数据。基于公式(1)和公式(2)计算测试图像的第一标准频谱数据的方式,相较于将预设图像与测试图像进行像素配准后再进行傅里叶变换以得到测试图像的第一标准频谱数据的传统方式,其计算过程简单,能够节约设备算力同时提高成像系统评估的效率。
在一种可能的实现方式中,上述S302~S303中描述的得到第一调制传递函数曲线的过程,包括:基于傅里叶变换得到的第一频谱数据矩阵中的每一个元素取模后,与基于上述公式(1)和公式(2)得到的第一标准频谱数据矩阵中的每个元素取模后的比值,为测试图像的二维调制传递函数取值,其数据形式为与第一频谱数据和第一标准频谱数据均同型的矩阵。然后,将该二维调制传递函数取值中,频率相同的多个值加和求平均,得到测试图像在频率域上的多个第一调制传递函数取值,多个第一调制传递函数取值的数据形式为一个一维矩阵。再利用第一算法对多个第一调至传递函数取值进行拟合,即可得到第一调制传递函数曲线。
通过上述得到第一调制传递函数曲线的过程可知,理想条件下的第一调制传递函数曲线应当是调制传递函数取值为1的直线,但事实并非如此的原因,除了待评估成像系统本身的性能因素之外,还可能是因为第一标准频谱数据是基于经验公式得到的,用该第一标准频谱数据表示测试图像真实的标准频谱数据存在误差。基于经验公式计算第一标准频谱数据虽然可以通过直接设置N为测试图像的像素大小,以免去预设图像与测试图像像素不相等时复杂的像素配准过程,但是基于公式中的第一参数和第二参数的含义,其取值理应随待评估成像系统所参照的预设图像的不同而变化,故而第一参数的取值和第二参数取值相对于当前所采用的预设图像的准确性,可能会影响第一标准频谱数据计算结果的准确性,进而可能会对得到的第一调制传递函数曲线造成影响。
又因为,测试图像是成像系统拍摄预设图像得到的,且已知第一参数和第二参数若取值合理,预设图像对应的第二调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值应均小于或等于预设阈值。其中,第二调制传递函数曲线的获取方式包括:基于第一参数的第一取值、第二参数的第一取值、预设图像的频谱坐标和预设图像的像素大小计算得到第二标准频谱数据;根据预设图像的第二频谱数据和预设图像的第二标准频谱数据,得到预设图像在频率域上的多个第二调制传递函数取值;利用第一算法对多个第二调制传递函数取值进行拟合,得到第二调制传递函数曲线。
故而,作为一个可选的实施例,将计算预设图像的第二标准频谱数据时用到的第一参数的第一取值和第二参数的第一取值,作为计算第一标准频谱数据时的第一参数的取值和第二参数的取值,以减小第一参数的取值和第二参数的取值选择不当造成的误差、提高后续成像系统的评估精度,该第一参数的第一取值和第二参数的第一取值应当使得预设图像对应的第二调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于预设阈值。
可选地,预设频率范围可以是指低频范围,具体可以是0~0.1周期每像素(cycles/pixel);预设阈值可以是0.01,但本申请对此不作具体限定。
作为一个可选的实施例,在上述S302之前,若在计算预设图像对应的调制传递函数曲线时,预设图像的第三标准频谱数据是基于第一参数的第二取值、第二参数的第二取值、预设图像的频谱坐标和预设图像的像素大小计算得到的,则会导致根据预设图像的第二频谱数据和预设图像的第三标准频谱数据,得到预设图像在频率域上的多个第三调制传递函数取值后,利用第一算法对多个第三调制传递函数取值进行拟合所得到第三调制传递函数曲线,其在预设频率范围内的调制传递函数取值与1的差值的绝对值存在大于预设阈值的情况。
故而在此种情况下,需要调整第一参数的取值,得到第一参数的第一取值;和/或,调整第二参数的取值,得到第二参数的第一取值,以校准计算第一标准频谱数据时第一参数和第二参数的取值。具体包括如下两种调整方式:
方式一是在第三调制传递函数曲线中频率为零对应的调制传递函数取值为第三调制传递函数曲线的实际取值时所采用的方式,具体包括如下步骤:
步骤一:判断第三调制传递函数曲线在预设频率范围内的变化率是否为零;若否,则执行步骤二,若是则执行步骤五。
步骤二:调整第一参数的取值,以更新第三调制传递函数曲线,确定更新的第三调制传递函数曲线在预设频率范围内的变化率为零时第一参数的目标取值,并将第一参数的目标取值确定为第一参数的第一取值。
步骤三:判断更新的第三调制传递函数曲线中,预设频率的调制传递函数取值与1的差值的绝对值是否小于或等于预设阈值,预设频率在预设频率范围内;若否,则执行步骤四,若是,则将第二参数的第二取值确定为第二参数的第一取值。
步骤四:调整第二参数的取值,以再次更新第三调制传递函数曲线,确定再次更新的第三调制传递函数曲线中预设频率的调制传递函数取值与1的差值的绝对值小于或等于预设阈值时第二参数的目标取值,并将第二参数的目标取值确定为第二参数的第一取值。
步骤五:将第一参数的第二取值确定为第一参数的第一取值,并调整第二参数的取值,以更新第三调制传递函数曲线;确定此次更新的第三调制传递函数曲线中预设频率的调制传递函数取值与1的差值的绝对值小于或等于预设阈值时第二参数的目标取值,并将第二参数的目标取值确定为第二参数的第一取值。
应理解,上述步骤中的步骤一和步骤三的判断至少有一个判断结果为否。
还应理解,本实施例中,将第三调制传递函数曲线应满足的条件,拆解为步骤一和步骤三的两级条件的目的是简化参数调整的过程,通过逐级的调整可以缩短参数调整所用的时间,以提高本申请所提供的成像系统评估方法的效率。
方式二是在第三调制传递函数曲线中频率为零对应的调制传递函数取值是预设值1时所采用的方式,具体包括:调整第一参数的取值和第二参数的取值,以更新第三调制传递函数曲线;确定更新的第三调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于预设阈值时第一参数的目标取值和第二参数的目标取值,并将第一参数的目标取值确定为第一参数的第一取值、第二参数的目标取值确定为第二参数的第一取值。
在一种可能的实现方式中,方式二对应的调整方式可以是先调整第二参数的取值,以使调整第二参数的取值后的第三调制传递函数曲线在预设频率范围的取值在1附近波动,然后在此基础上调整第一参数的取值,以缩小调制函数取值波动的幅度,通过循环多次的调整第一参数的取值和第二参数的取值,直至调整后的第三调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于预设阈值。
值得注意的是,不论是方式一还是方式二,在第一参数和/或第二参数得到一个新的取值之后,均需要重新计算新的第三标准频谱数据、在频率域上的新的多个第三调制传递函数取值、以及基于新的多个第三调制传递函数取值拟合出新的第三调制传递函数曲线,直至其能满足在预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于预设阈值的条件,其中,每次得到新的第三调制传递函数曲线的方式与上述得到第一调制传递函数曲线的方式类似,不再赘述。
作为一个可选的实施例,上述303中用于对多个第一调制传递函数取值个拟合的第一算法为最小二乘法、多项式拟合法、样条差值法和非参数拟合法中的任一个。
作为一个可选的实施例,上述304的一种实现方式包括:采用纹理锐度值分析法、十分之一幅度法和半功率法中的任一个或多个,对第一调制传递函数曲线进行分析,得到待评估成像系统的评估结果。
示例性地,通过如下所示的纹理锐度(Texture Acutance)值计算公式(3)对第一调制传递函数曲线进行纹理锐度值计算,得到的纹理锐度值若小于第一阈值,则可以认为测试图像不够清晰,该待评估成像系统的对比度不足;若得到的纹理锐度值大于第二阈值,则可以认为测试图像清晰,该待评估成像系统的对比度满足要求。可选地,第一阈值可以是0.7,第二阈值可以是0.9,但本申请对此不作具体限定。
(3)
其中,即为第一调制传递函数曲线对应的函数,u为图像的频率,取值范围为0~0.5cycles/pixel,/>为对比敏感度函数(Contrast Sensitivity Function,CSF),其公式为/>,其中,a、b、c为常数。示例性地,a、b、c的取值分别可以设为75、0.2、0.8。
作为一个可选的实施例,预设图像可以是枯叶图。枯叶图包含丰富的纹理细节且与自然图像相近,可用于对成像系统的评估。但应理解,任意具有丰富纹理细节的图像均可以作为预设图像来对成像系统进行评估,本申请对预设图像的具体图像类型和内容不做具体限定。
下面以预设图像为枯叶图、第一电子设备为笔记本电脑、第二电子设备为手机为例,对本申请的具体实现进行详细说明。在下文中,调制传递函数曲线简称为MTF曲线,调制传递函数取值简称为MTF值。
如图5所示,本申请提供的成像系统评估方法500包括三个阶段。
第一阶段:计算基于枯叶图的MTF曲线1。
第二阶段:判断MTF曲线1在预设在预设频率范围内的MTF值与1的差值的绝对值是否均小于或等于预设阈值。若是,则进入第三阶段;若否,则对第一参数和第二参数进行调整,直至满足该条件后,进入第三阶段。
第三阶段:计算基于测试图像的MTF曲线2,并对MTF曲线2进行分析得到待评估成像系统的评估结果。
其中,第一阶段的具体步骤如图6所示,具体包括:
S601、获取枯叶图,采用傅里叶变换,计算得到枯叶图的频谱数据。
应理解,该枯叶图可以预存于笔记本电脑中,也可以由笔记本电脑从标准公布网站查询、获取。
可选地,枯叶图的频谱数据可以记为。/>是一个二维矩阵,矩阵中的元素/>表示的是该枯叶图中不同频率对应的频谱强度,i、j的取值范围均为从1到P的正整数,P表示枯叶图的像素大小。
S602、获取枯叶图的像素大小和频谱坐标,并预设的第一参数取值为/>、第二参数/>取值为/>,计算得到该枯叶图的标准频谱数据。
可选地,该枯叶图的标准频谱数据可以记为,/>与/>为同型矩阵,二者的不同在于,/>是基于傅里叶变换得到的,而/>是基于经验公式得到的,两个矩阵中元素的含义相同,元素取值因计算方式的不同而存在差异。
在一种可能的实现方式中,中的每一个元素均可以根据上述公式(1)和公式(2)的标准频谱公式计算得到。
示例性地,若该枯叶图的像素大小为300*300,图像的频谱坐标原点位于图像的中心,在计算该枯叶图的标准频谱数据时,N即为300,m和n的取值为[-150,150]内的整数,第一参数取值为/>、第二参数/>取值为/>。图像的频谱坐标(-150,150)对应的标准频谱数据即为矩阵/>中的元素/>,图像的频谱坐标(150,150)对应的标准频谱数据即为矩阵/>中的元素/>。/>
将枯叶图中每一个频谱坐标均带入上述公式(1)和公式(2)进行计算,即可得到该枯叶图的标准频谱数据。
S603、利用枯叶图的频谱数据和枯叶图的标准频谱数据,计算得到枯叶图的二维MTF数据。
在一种可能的实现方式中,枯叶图的二维MTF数据记为,/>与/>、/>均为同型矩阵。
具体地,根据如下公式(4)计算得到二维MTF数据,即中的元素取模值后,与/>中对应位置元素的模值的比值,即得到矩阵/>对应位置的元素。
(4)
S604、对枯叶图的二维MTF数据中,频率相同的数据加和求平均,得到的多个值为枯叶图的一维MTF数据。
应理解,枯叶图的二维MTF数据中存在多个不同频率对应的数据,每个频率对应的数据个数不唯一,频率相同的数据也即是距离频谱坐标原点的欧氏距离相同的数据。枯叶图的一维MTF数据为多个频率各自对应的唯一值形成的一维矩阵。
S605、采用第一算法对枯叶图的一维MTF数据进行拟合,得到基于枯叶图的MTF曲线1,MTF曲线1为MTF值随频率变化的曲线。
第二阶段的具体步骤如图7所示,包括:
S701、判断MTF曲线1在预设频率范围内的MTF值与1的差值的绝对值是否均小于或等于预设阈值。
若是,则执行S702;若否,则对第一参数,和/或,第二参数/>进行调整。
S702、确定第一参数取值为/>、第二参数/>取值为/>,执行后续第三阶段对应的步骤。
若在拟合MTF曲线1时,频率为零的MTF值为MTF曲线1的实际值,其参数调整过程执行S703~S704~S705~S706的流程,S703~S707~S708的流程,或者,S703~S704~S705~S709的流程。
若在拟合MTF曲线1时,预设频率为零的MTF值等于1,则可以循环调整第一参数和第二参数/>的取值,直至调整后的MTF曲线1在预设频率范围内的MTF值与1的差值的绝对值均小于或等于预设阈值,获取此时第一参数/>的取值和第二参数/>的取值。
应理解,不论执行了上述哪一种流程,满足S701所描述的条件时的基于同一幅枯叶图得到的MTF曲线是相同的。
S703、判断MTF曲线1在该预设频率范围内的变化率是否为零。
若否,则执行S704;若是,则执行S707。
在一种可能的实现方式中,在一些评估精度要求不高的场景下,该MTF曲线1在预设频率范围内的变化率的绝对值在一个预设范围内时,也可认为其满足条件。
S704、调整第一参数的取值以更新MTF曲线1,直至更新的MTF曲线1在预设频率范围内的变化率为零时,获取此时第一参数/>的取值/>。
S705、判断该更新的MTF曲线1中,预设频率对应的MTF值与1的差值的绝对值是否小于或等于预设阈值。
应理解,预设频率为预设频率范围内的任一频率,在更新的MTF曲线1在预设频率范围内的变化率为零时,选择预设范围内的任意一个频率对应的MTF值,判断其是否小于或等于预设阈值,可以节约服务器算力,同时提高成像系统的评估效率。
若否,则执行S706;若是,则执行S709。
S706、调整第二参数的取值,以更新该更新的MTF曲线1,直至再次更新的MTF曲线1中,预设频率对应的MTF值与1的差值的绝对值均小于或等于预设阈值时,获取此时第二参数/>的取值/>。
S707、确定第一参数的取值为/>。
S708、调整第二参数的取值,以更新MTF曲线1,直至此次更新的MTF曲线1中,预设频率对应的MTF值与1的差值的绝对值小于或等于预设阈值时,获取此时第二参数/>的取值。
S709、获取此时第二参数的取值/>。
容易理解的是,枯叶图对应的MTF曲线同时满足“在预设频率范围内的变化率为零”且“预设频率对应的MTF值与1的差值的绝对值小于或等于预设阈值”,即满足S701中描述的“在预设频率范围内的MTF值与1的差值的绝对值均小于或等于预设阈值”。在S701中判断出MTF曲线1不满足判断条件时,后续设置的“在预设频率范围内的变化率是否为零”、“预设频率对应的MTF值与1的差值的绝对值是否小于或等于预设阈值”的判断并未改变S6的核心目的,仅是为了简化参数调整的过程。
图8示出了在执行S703~S704~S705~S706的流程过程中枯叶图的MTF曲线的变化过程。图8中a示出了MTF曲线1的形态以及MTF曲线1的原始数据,该原始数据即多个一维MTF值。在图8中的a所示出的MTF曲线1不满足预设频率范围内的MTF值与1的差值的绝对值小于或等于预设阈值,其初始第一参数的取值为/>,第二参数/>的取值为/>。然后,调整第一参数/>的取值为/>得到图8中的b,继续调整第二参数/>的取值为/>得到图8中的c。图8中的c满足预设频率范围内的MTF值与1的差值的绝对值小于或等于预设阈值。
应理解,上述各步骤的序号的大小并不意味着执行顺序的先后,各步骤的执行顺序应以其功能和内在逻辑确定。
第三阶段的具体步骤如图9所示,包括:
S901、获取来自手机的测试图像,该测试图像是由手机对枯叶图进行拍摄得到的。
S902、采用傅里叶变换,计算得到测试图像的频谱数据。
S903、获取测试图像的像素大小和频谱坐标,设置的第一参数和第二参数/>的取值为第二阶段所得到的取值,计算该测试图像的标准频谱数据。
在一种可能的实现方式中,S903在S702之后执行,在计算测试图像的标准频谱数据时,第一参数取值为/>、第二参数/>取值为/>。
在一种可能的实现方式中,基于枯叶图的MTF曲线1不能直接满足S701中描述的判断条件,随后执行了S703~S704~S705~S706的参数调整流程,则S903在S706之后执行,在计算测试图像的标准频谱数据时,确定第一参数取值为/>、第二参数/>取值为/>。
在一种可能的实现方式中,基于枯叶图的MTF曲线1不能直接满足S701中描述的判断条件,随后执行了S703~S707~S708的参数调整流程,则S903在S708之后执行,在计算测试图像的标准频谱数据时,确定第一参数取值为/>、第二参数/>取值为/>。
在一种可能的实现方式中,基于枯叶图的MTF曲线1不能直接满足S701中描述的判断条件,随后执行了S703~S704~S705~S709的参数调整流程,则S903在S709之后执行,在计算测试图像的标准频谱数据时,确定第一参数取值为/>、第二参数/>取值为/>。
值得注意的是,在运用上述公式(1)和公式(2)计算测试图像的标准频谱数据时,参数N应代入测试图像的像素大小,对应的,频谱坐标m和n也应为测试图像的频谱坐标。
S904、利用测试图像的频谱数据和测试图像的标准频谱数据,计算得到测试图像的二维MTF数据。
S905、对测试图像的二维MTF数据中,频率相同的数据加和求平均,得到的多个值为测试图像的一维MTF数据。
S906、对测试图像的一维MTF数据进行拟合,得到基于测试图像的MTF曲线2,MTF曲线2为MTF值随频率变化的曲线,MTF曲线2中频率为零对应的MTF值为MTF曲线的实际取值。
应理解,可以基于两次傅里叶变换的MTF曲线计算方法来对本申请所提供的成像系统进行验证。预设图像经过像素配准后再经过傅里叶变换得到的测试图像的标准频谱数据,在根据傅里叶变换得到测试图像的频谱数据,两次傅里叶变换所得到的测试图像的频谱数据和标准频谱数据所MTF曲线计算方法复杂但具有较高的可信度,图10中曲线K为根据两次傅里叶变换得到测试图像的MTF曲线;本申请中,通过对测试图像作一次傅里叶变换得到图像的频谱数据,再基于经验公式得到测试图像的标准频谱数据,以此数据得到的测试图像的MTF曲线如图10中的L所示;可以看出图10中曲线K和曲线L具有很高的相似度,可见,本申请提供的方法在减少计算量的同时,保持了计算精度。
S907、对该MTF曲线2进行分析,得到对手机成像系统的评估结果。
应理解,上述各步骤的序号的大小并不意味着执行顺序的先后,各步骤的执行顺序应以其功能和内在逻辑确定。
本申请还提供了另一种成像系统评估方法1100,将预设频率范围内的调制传递函数取值与1的差值的绝对值小于或等于预设阈值的、预设图像的第二调制传递函数曲线计算过程中所用到的第一参数的第一取值和第二参数的第二取值,用于计算得到测试图像的第一调制传递函数曲线,以使测试图像在计算第一标准频谱数据时的第一参数和第二参数的取值合理,减小由于参数选择不准确而造成的误差,提高成像系统评估的精度。
方法1100具体包括如下步骤:
S1101、根据预设图像的第二频谱数据和预设图像的第二标准频谱数据,得到预设图像在频率域上的多个第二调制传递函数取值,第二标准频谱数据是基于第一参数的第一取值、第二参数的第一取值、预设图像的频谱坐标和预设图像的像素大小计算得到的。
S1102、利用第一算法对多个第二调制传递函数取值进行拟合,得到第二调制传递函数曲线,第二调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值小于或等于预设阈值。
S1103、获取来自第二电子设备的测试图像,测试图像是通过第二电子设备中的待评估成像系统拍摄预设图像得到的。
S1104、根据测试图像的第一频谱数据和测试图像的第一标准频谱数据,得到测试图像在频率域上的多个第一调制传递函数取值,第一标准频谱数据是基于第一参数的第一取值、第二参数的第一取值、测试图像的频谱坐标和测试图像的像素大小计算得到。
S1105、利用第一算法对多个第一调制传递函数取值进行拟合,得到第一调制传递函数曲线,第一调制传递函数曲线为调制传递函数取值与频率之间的关系曲线;对第一调制传递函数曲线进行分析,得到待评估成像系统的评估结果。
应理解,第一标准频谱数据和第二标准频谱数据的计算公式均为上述公式(1)和公式(2),由于测试图像是待评估成像系统拍摄预设图像得到的,将预设频率范围内的调制传递函数取值与1的差值的绝对值小于或等于预设阈值的预设图像的第二调制传递函数曲线计算过程中所使用的第一参数的第一取值和第二参数的第二取值,作为计算测试图像的第一调制传递函数曲线时的参数取值,有利于减小第一调制传递函数曲线计算过程中参数取值与当前所使用的预设图像相关性差而产生的误差,有利于提高成像系统的评估精度。
还应理解,在基于不同的预设图像对同一待评估成像系统进行评估时,在计算预设图像各自对应的测试图像的调制传递函数曲线时,可以分别基于各自对应的预设图像对参数进行校准后再计算,以减小采用的预设图像不同而对同一待评估成像系统造成的误差。示例性地,待评估成像系统拍摄预设图像1得到测试图像1-1,则将预设频率范围内的调制传递函数取值与1的差值的绝对值小于或等于预设阈值的预设图像1的调制传递函数曲线1计算过程中所使用的第一参数的取值和第二参数的取值,作为计算测试图像1-1的调制传递函数曲线时的参数取值;待评估成像系统拍摄预设图像2得到测试图像2-1,则将预设频率范围内的调制传递函数取值与1的差值的绝对值小于或等于预设阈值的预设图像2的制传递函数曲线2计算过程中所使用的第一参数的取值和第二参数的取值,作为计算测试图像2-1的调制传递函数曲线时的参数取值。
值得注意的是,上述第一调制传递函数曲线和第二调制传递函数曲线在频率为零时的调制传递函数取值可以是第一调制传递函数曲线和第二调制传递函数曲线对应的实际取值,也可以是预设值1,本申请对此不作限定。
作为一个可选的实施例,若在上述S1104之前,计算预设图像的第二标准频谱数据时没有直接选取到第一参数的第一取值、第二参数的第一取值,则计算得到的预设图像的调制传递函数曲线,是预设频率范围内的调制传递函数取值与1的差值的绝对值存在大于预设阈值的第三调制传递函数曲线,其中,第三调制传递函数曲线是利用第一算法对多个第三调制传递函数取值进行拟合得到的,多个第三调制传递函数取值是根据预设图像的第二频谱数据和预设图像的第三标准频谱数据得到的,第三标准频谱数据是基于第一参数的第二取值、第二参数的第二取值、预设图像的频谱坐标和预设图像的像素大小计算得到的。
在此种情况下,需要调整第一参数的取值,得到第一参数的第一取值;和/或,调整第二参数的取值,得到第二参数的第一取值,其调整过程,与上述图7所描述的调整过程一致,此处不再赘述。
上文结合图2至图11,详细描述了本申请实施例的成像系统评估方法,下面结合图12和图13,详细描述本申请实施例的成像系统评估装置。
图12为本申请提供的一种成像系统评估装置1200,包括获取模块1201和处理模块1202。其中,获取模块1201用于:获取来自第二电子设备的测试图像,测试图像是通过第二电子设备中的待评估成像系统拍摄预设图像得到的;处理模块1202用于:根据测试图像的第一频谱数据和测试图像的第一标准频谱数据,得到测试图像在频率域上的多个第一调制传递函数取值;利用第一算法对多个第一调制传递函数取值进行拟合,得到第一调制传递函数曲线,第一调制传递函数曲线为调制传递函数取值与频率之间的关系曲线,在第一调制传递函数曲线中,频率为零对应的调制传递函数取值为第一调制传递函数曲线的实际取值;以及,对第一调制传递函数曲线进行分析,得到待评估成像系统的评估结果。
可选地,处理模块1202还用于:基于第一参数的第一取值、第二参数的第一取值、测试图像的频谱坐标和测试图像的像素大小计算得到第一标准频谱数据;根据预设图像的第二频谱数据和预设图像的第二标准频谱数据,得到预设图像在频率域上的多个第二调制传递函数取值,第二标准频谱数据是基于第一参数的第一取值、第二参数的第一取值、预设图像的频谱坐标和预设图像的像素大小计算得到的;利用第一算法对多个第二调制传递函数取值进行拟合,得到第二调制传递函数曲线,第二调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于预设阈值。
可选地,处理模块1202还用于:根据预设图像的第二频谱数据和预设图像的第三标准频谱数据,得到预设图像在频率域上的多个第三调制传递函数取值,第三标准频谱数据是基于第一参数的第二取值、第二参数的第二取值、预设图像的频谱坐标和预设图像的像素大小计算得到的;利用第一算法对多个第三调制传递函数取值进行拟合,得到第三调制传递函数曲线;在第三调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值存在大于预设阈值的情况下,调整第一参数的取值,得到第一参数的第一取值;和/或,调整第二参数的取值,得到第二参数的第一取值。
可选地,在第三调制传递函数曲线中,频率为零对应的调制传递函数取值为第三调制传递函数曲线的实际取值。
可选地,处理模块1202具体用于:判断第三调制传递函数曲线在预设频率范围内的变化率是否为零;若否,则调整第一参数的取值,以更新第三调制传递函数曲线;确定更新的第三调制传递函数曲线在预设频率范围内的变化率为零时第一参数的目标取值,并将第一参数的目标取值确定为第一参数的第一取值;判断更新的第三调制传递函数曲线中,预设频率的调制传递函数取值与1的差值的绝对值是否小于或等于预设阈值,预设频率在预设频率范围内;若否,则调整第二参数的取值,以再次更新第三调制传递函数曲线;确定再次更新的第三调制传递函数曲线中预设频率的调制传递函数取值与1的差值的绝对值小于或等于预设阈值时第二参数的目标取值,并将第二参数的目标取值确定为第二参数的第一取值。
可选地,处理模块1202还具体用于:若更新的第三调制传递函数曲线中,预设频率的调制传递函数取值与1的差值的绝对值小于或等于预设阈值,则将第二参数的第二取值确定为第二参数的第一取值。
可选地,处理模块1202还具体用于:若第三调制传递函数曲线在预设频率范围内的变化率为零,则将第一参数的第二取值确定为第一参数的第一取值;调整第二参数的取值,以更新第三调制传递函数曲线;确定更新后的第三调制传递函数曲线中预设频率的调制传递函数取值与1的差值的绝对值小于或等于预设阈值时第二参数的目标取值,并将第二参数的目标取值确定为第二参数的第一取值。
可选地,预设第三调制传递函数曲线中频率为零对应的调制传递函数取值为1。
可选地,处理模块1202还具体用于:调整第一参数的取值和第二参数的取值,以更新第三调制传递函数曲线;确定更新的第三调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于预设阈值时第一参数的目标取值和第二参数的目标取值,并将第一参数的目标取值确定为第一参数的第一取值、第二参数的目标取值确定为第二参数的第一取值。
可选地,处理模块1202还具体用于:对测试图像进行傅里叶变换,得到测试图像的第一频谱数据。
可选地,处理模块1202具体用于:将测试图像转换为灰度图像;对灰度图像进行傅里叶变换,得到测试图像的第一频谱数据。
可选地,处理模块1202具体用于:采用纹理锐度值分析法、十分之一幅度法和半功率法中的任一个,对第一调制传递函数曲线进行分析,得到待评估成像系统的评估结果。
可选地,预设图像为枯叶图。
可选地,第一算法为最小二乘法、多项式拟合法、样条差值法和非参数拟合法中的任一个。
应理解,这里的装置1200以功能模块的形式体现。这里的术语“模块”可以指应用特有集成电路(application specific integrated circuit,ASIC)、电子电路、用于执行一个或多个软件或固件程序的处理器(例如共享处理器、专有处理器或组处理器等)和存储器、合并逻辑电路和/或其它支持所描述的功能的合适组件。在一个可选的例子中,本领域技术人员可以理解,成像系统评估装置1200可以用于执行上述方法实施例中与第一电子设备对应的各个流程和/或步骤,为避免重复,在此不再赘述。
上述成像系统评估装置1200具有实现上述方法中第一电子设备执行的相应步骤的功能;上述功能可以通过硬件实现,也可以通过硬件执行相应的软件实现。该硬件或软件包括一个或多个与上述功能相对应的模块。
在本申请的实施例中,图12中的成像系统评估装置1200也可以是芯片,例如:SOC。对应地,处理模块1202可以是该芯片的处理电路,在此不做限定。
图13是本申请实施例提供的另一种成像系统评估装置1300。该成像系统评估装置1300包括处理器1301、收发器1302和存储器1303。其中,处理器1301、收发器1302和存储器1303通过内部连接通路互相通信,该存储器1303用于存储指令,该处理器1301用于执行该存储器1303存储的指令,以控制该收发器1302发送信号和/或接收信号。
应理解,成像系统评估装置1300可以具体为上述实施例中的第一电子设备,并且可以用于执行上述方法实施例中与第一电子设备对应的各个步骤和/或流程。可选地,该存储器1303可以包括只读存储器和随机存取存储器,并向处理器提供指令和数据。存储器的一部分还可以包括非易失性随机存取存储器。例如,存储器还可以存储设备类型的信息。该处理器1301可以用于执行存储器中存储的指令,并且当该处理器1301执行存储器中存储的指令时,该处理器1301用于执行上述与该第一电子设备对应的方法实施例的各个步骤和/或流程。该收发器1302可以包括发射器和接收器,该发射器可以用于实现上述收发器对应的用于执行发送动作的各个步骤和/或流程,该接收器可以用于实现上述收发器对应的用于执行接收动作的各个步骤和/或流程。
应理解,在本申请实施例中,该处理器可以是中央处理单元(central processingunit,CPU),该处理器还可以是其他通用处理器、数字信号处理器(DSP)、专用集成电路(ASIC)、现场可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。
在实现过程中,上述方法的各步骤可以通过处理器中的硬件的集成逻辑电路或者软件形式的指令完成。结合本申请实施例所公开的方法的步骤可以直接体现为硬件处理器执行完成,或者用处理器中的硬件及软件模块组合执行完成。软件模块可以位于随机存储器,闪存、只读存储器,可编程只读存储器或者电可擦写可编程存储器、寄存器等本领域成熟的存储介质中。该存储介质位于存储器,处理器执行存储器中的指令,结合其硬件完成上述方法的步骤。为避免重复,这里不再详细描述。
本申请还提供了一种计算机可读存储介质,该计算机可读存储介质用于存储计算机程序,该计算机程序用于实现上述实施例中与第一电子设备对应的方法。
本申请还提供了一种计算机程序产品,该计算机程序产品包括计算机程序(也可以称为代码,或指令),当该计算机程序在计算机上运行时,该计算机可以执行上述实施例所示的第一电子设备对应的方法。
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的模块及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本申请的范围。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统、装置和模块的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
在本申请所提供的几个实施例中,应该理解到,所揭露的系统、装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,模块的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个模块或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或模块的间接耦合或通信连接,可以是电性,机械或其它的形式。
作为分离部件说明的模块可以是或者也可以不是物理上分开的,作为模块显示的部件可以是或者也可以不是物理模块,即可以位于一个地方,或者也可以分布到多个网络模块上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。
另外,在本申请各个实施例中的各功能模块可以集成在一个处理模块中,也可以是各个模块单独物理存在,也可以两个或两个以上模块集成在一个模块中。
功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本申请各个实施例所描述的方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(read-only memory,ROM)、随机存取存储器(random access memory,RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上,仅为本申请的具体实施方式,但本申请实施例的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请实施例揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本申请实施例的保护范围之内。因此,本申请实施例的保护范围应以权利要求的保护范围为准。
Claims (17)
1.一种成像系统评估方法,其特征在于,应用于第一电子设备,包括:
获取来自第二电子设备的测试图像,所述测试图像是通过所述第二电子设备中的待评估成像系统拍摄预设图像得到的;
根据所述测试图像的第一频谱数据和所述测试图像的第一标准频谱数据,得到所述测试图像在频率域上的多个第一调制传递函数取值;
利用第一算法对所述多个第一调制传递函数取值进行拟合,得到第一调制传递函数曲线,所述第一调制传递函数曲线为调制传递函数取值与频率之间的关系曲线,在所述第一调制传递函数曲线中,频率为零对应的调制传递函数取值为所述第一调制传递函数曲线的实际取值;
对所述第一调制传递函数曲线进行分析,得到所述待评估成像系统的评估结果。
2.根据权利要求1所述的方法,其特征在于,所述第一标准频谱数据是基于第一参数的第一取值、第二参数的第一取值、所述测试图像的频谱坐标和所述测试图像的像素大小计算得到的;
所述方法还包括:
根据所述预设图像的第二频谱数据和所述预设图像的第二标准频谱数据,得到所述预设图像在频率域上的多个第二调制传递函数取值,所述第二标准频谱数据是基于所述第一参数的第一取值、所述第二参数的第一取值、所述预设图像的频谱坐标和所述预设图像的像素大小计算得到的;
利用所述第一算法对所述多个第二调制传递函数取值进行拟合,得到第二调制传递函数曲线,所述第二调制传递函数曲线在预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于预设阈值。
3.根据权利要求2所述的方法,其特征在于,在根据所述测试图像的第一频谱数据和所述测试图像的第一标准频谱数据,得到所述测试图像在频率域上的多个第一调制传递函数取值之前,所述方法还包括:
根据所述预设图像的第二频谱数据和所述预设图像的第三标准频谱数据,得到所述预设图像在频率域上的多个第三调制传递函数取值,所述第三标准频谱数据是基于所述第一参数的第二取值、所述第二参数的第二取值、所述预设图像的频谱坐标和所述预设图像的像素大小计算得到的;
利用所述第一算法对所述多个第三调制传递函数取值进行拟合,得到第三调制传递函数曲线;
在所述第三调制传递函数曲线在所述预设频率范围内的调制传递函数取值与1的差值的绝对值存在大于所述预设阈值的情况下,调整所述第一参数的取值,得到所述第一参数的第一取值,和/或,调整所述第二参数的取值,得到所述第二参数的第一取值。
4.根据权利要求3所述的方法,其特征在于,在所述第三调制传递函数曲线中,频率为零对应的调制传递函数取值为所述第三调制传递函数曲线的实际取值。
5.根据权利要求4所述的方法,其特征在于,所述调整所述第一参数的取值,得到所述第一参数的第一取值,和/或,调整所述第二参数的取值,得到所述第二参数的第一取值,包括:
判断所述第三调制传递函数曲线在预设频率范围内的变化率是否为零;
若否,则调整所述第一参数的取值,以更新所述第三调制传递函数曲线;
确定更新的第三调制传递函数曲线在所述预设频率范围内的变化率为零时所述第一参数的目标取值,并将所述第一参数的目标取值确定为所述第一参数的第一取值;
判断所述更新的第三调制传递函数曲线中,预设频率的调制传递函数取值与1的差值的绝对值是否小于或等于所述预设阈值,所述预设频率在所述预设频率范围内;
若否,则调整所述第二参数的取值,以再次更新所述第三调制传递函数曲线;
确定再次更新的第三调制传递函数曲线中所述预设频率的调制传递函数取值与1的差值的绝对值小于或等于所述预设阈值时所述第二参数的目标取值,并将所述第二参数的目标取值确定为所述第二参数的第一取值。
6.根据权利要求5所述的方法,其特征在于,所述方法还包括:
若所述更新的第三调制传递函数曲线中,所述预设频率的调制传递函数取值与1的差值的绝对值小于或等于所述预设阈值,则将所述第二参数的第二取值确定为所述第二参数的第一取值。
7.根据权利要求5所述的方法,其特征在于,所述方法还包括:
若所述第三调制传递函数曲线在所述预设频率范围内的变化率为零,则将所述第一参数的第二取值确定为所述第一参数的第一取值;
调整所述第二参数的取值,以更新所述第三调制传递函数曲线;
确定更新后的第三调制传递函数曲线中所述预设频率的调制传递函数取值与1的差值的绝对值小于或等于所述预设阈值时所述第二参数的目标取值,并将所述第二参数的目标取值确定为所述第二参数的第一取值。
8.根据权利要求3所述的方法,其特征在于,预设所述第三调制传递函数曲线中频率为零对应的调制传递函数取值为1。
9.根据权利要求8所述的方法,其特征在于,所述调整所述第一参数的取值,得到所述第一参数的第一取值,和/或,调整所述第二参数的取值,得到所述第二参数的第一取值,包括:
调整所述第一参数的取值和所述第二参数的取值,以更新所述第三调制传递函数曲线;
确定更新的第三调制传递函数曲线在所述预设频率范围内的调制传递函数取值与1的差值的绝对值均小于或等于所述预设阈值时所述第一参数的目标取值和所述第二参数的目标取值,并将所述第一参数的目标取值确定为所述第一参数的第一取值、所述第二参数的目标取值确定为所述第二参数的第一取值。
10.根据权利要求1所述的方法,其特征在于,所述方法还包括:
对所述测试图像进行傅里叶变换,得到所述测试图像的第一频谱数据。
11.根据权利要求10所述的方法,其特征在于,在所述测试图像为彩色图像的情况下,在对所述测试图像进行傅里叶变换之前,所述方法还包括:
将所述测试图像转换为灰度图像;
所述对所述测试图像进行傅里叶变换,得到所述测试图像的第一频谱数据,包括:
对所述灰度图像进行傅里叶变换,得到所述测试图像的第一频谱数据。
12.根据权利要求1所述的方法,其特征在于,所述对所述第一调制传递函数曲线进行分析,得到所述待评估成像系统的评估结果,包括:
采用纹理锐度值分析法、十分之一幅度法和半功率法中的任一个,对所述第一调制传递函数曲线进行分析,得到所述待评估成像系统的评估结果。
13.根据权利要求1所述的方法,其特征在于,所述预设图像为枯叶图。
14.根据权利要求1所述的方法,其特征在于,所述第一算法为最小二乘法、多项式拟合法、样条差值法和非参数拟合法中的任一个。
15.一种成像系统评估装置,其特征在于,包括:用于实现如权利要求1至14中任一项所述的方法的模块。
16.一种电子设备,其特征在于,包括:处理器,所述处理器与存储器耦合,所述存储器用于存储计算机程序,当所述处理器调用所述计算机程序时,使得所述电子设备执行如权利要求1至14中任一项所述的方法。
17.一种计算机可读存储介质,其特征在于,用于存储计算机程序,所述计算机程序包括用于实现如权利要求1至14中任一项所述的方法的指令。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311099702.5A CN116847067B (zh) | 2023-08-30 | 2023-08-30 | 成像系统评估方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311099702.5A CN116847067B (zh) | 2023-08-30 | 2023-08-30 | 成像系统评估方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116847067A true CN116847067A (zh) | 2023-10-03 |
CN116847067B CN116847067B (zh) | 2023-11-24 |
Family
ID=88163792
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311099702.5A Active CN116847067B (zh) | 2023-08-30 | 2023-08-30 | 成像系统评估方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116847067B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004134861A (ja) * | 2002-10-08 | 2004-04-30 | Nikon Corp | 解像度評価方法、解像度評価プログラム、および光学機器 |
CN106840617A (zh) * | 2017-03-02 | 2017-06-13 | 西安工业大学 | 变频栅条靶标及其光电成像系统动态传递函数测量方法 |
CN109587475A (zh) * | 2019-01-03 | 2019-04-05 | 北京理工大学 | 数字相机动态分辨率测量技术与装置 |
CN114577444A (zh) * | 2022-03-04 | 2022-06-03 | 合肥视涯技术有限公司 | 一种调制传递函数的检测方法及检测装置 |
-
2023
- 2023-08-30 CN CN202311099702.5A patent/CN116847067B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004134861A (ja) * | 2002-10-08 | 2004-04-30 | Nikon Corp | 解像度評価方法、解像度評価プログラム、および光学機器 |
CN106840617A (zh) * | 2017-03-02 | 2017-06-13 | 西安工业大学 | 变频栅条靶标及其光电成像系统动态传递函数测量方法 |
CN109587475A (zh) * | 2019-01-03 | 2019-04-05 | 北京理工大学 | 数字相机动态分辨率测量技术与装置 |
CN114577444A (zh) * | 2022-03-04 | 2022-06-03 | 合肥视涯技术有限公司 | 一种调制传递函数的检测方法及检测装置 |
Also Published As
Publication number | Publication date |
---|---|
CN116847067B (zh) | 2023-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20230419437A1 (en) | Systems and methods for fusing images | |
US10091479B2 (en) | Hardware-based convolutional color correction in digital images | |
US10949958B2 (en) | Fast fourier color constancy | |
US9672604B2 (en) | Convolutional color correction | |
US9749551B2 (en) | Noise models for image processing | |
US9979942B2 (en) | Per pixel color correction filtering | |
US10013764B2 (en) | Local adaptive histogram equalization | |
EP2582125B1 (en) | Use of noise-optimized selection criteria to calculate scene white points | |
US20160300337A1 (en) | Image fusion method and image processing apparatus | |
JP2016525840A (ja) | 画像アーティファクトを補正するシステム及び方法 | |
CN109672871B (zh) | 白平衡信息同步方法、装置及计算机可读介质 | |
US10554914B1 (en) | Adjusting confidence values for correcting pixel defects | |
US20190082154A1 (en) | Image processing method and device for auto white balance | |
CN116158087A (zh) | 一种多相机色彩一致性校正方法和装置 | |
CN112132925A (zh) | 用于重建水下图像颜色的方法和装置 | |
CN116847067B (zh) | 成像系统评估方法和装置 | |
JP6745317B2 (ja) | 投影システムにおける測色を補償する方法、その方法を実行する装置、コンピュータ可読記憶媒体、および、投影システム | |
CN115426487A (zh) | 色彩校正矩阵调整方法、装置、电子设备及可读存储介质 | |
US11232289B2 (en) | Face identification method and terminal device using the same |
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 |