CN105976317A - 一种图像空间退化模拟方法及系统 - Google Patents
一种图像空间退化模拟方法及系统 Download PDFInfo
- Publication number
- CN105976317A CN105976317A CN201610274396.8A CN201610274396A CN105976317A CN 105976317 A CN105976317 A CN 105976317A CN 201610274396 A CN201610274396 A CN 201610274396A CN 105976317 A CN105976317 A CN 105976317A
- Authority
- CN
- China
- Prior art keywords
- image
- function
- mtf
- matrix
- spatial resolution
- 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 53
- 238000004088 simulation Methods 0.000 title claims abstract description 53
- 230000015556 catabolic process Effects 0.000 title claims abstract description 46
- 238000006731 degradation reaction Methods 0.000 title claims abstract description 46
- 239000011159 matrix material Substances 0.000 claims description 102
- 238000012546 transfer Methods 0.000 claims description 48
- 230000000007 visual effect Effects 0.000 claims description 41
- 230000000694 effects Effects 0.000 claims description 10
- 230000007850 degeneration Effects 0.000 claims description 7
- 230000009466 transformation Effects 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 abstract description 7
- 238000011160 research Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 14
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4092—Image resolution transcoding, e.g. by using client-server architectures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- 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/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- 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/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种图像空间退化模拟方法,该方法是依据待模拟图像奈奎斯特频率下的MTF值,以及基础图像和待模拟图像的空间分辨率值,针对基础图像进行空间退化模拟,获得空间退化后的待模拟图像。本发明通过考虑MTF或PSF对于图像空间实际分辨能力的影响,解决了现有技术中由于采用像元求平均值来实现空间退化模拟的方法得到的图像分辨率无法满足新型遥感器图像数据模拟精度或无法为模型算法研究提供需要的图像数据的问题。
Description
技术领域
本发明涉及遥感传感器成像技术领域,更具体的说,是涉及一种图像空间退化模拟方法及系统。
背景技术
随着遥感技术和计算机技术的快速发展,遥感成像产品在各行各业中起着越来越重要的作用,因此遥感成像分辨率越来越高,数据量也越来越大。当需要模拟特定空间分辨率的图像时,一般是以较高空间分辨率的图像数据为基础图像,对其进行计算处理获得对应的空间分辨率图像,这个过程为图像的空间退化模拟,它一般可用于新研制遥感器的图像数据模拟或为模型算法研究提供需要的图像数据。
现有技术中,图像空间退化模拟根据较高的空间分辨率图像与待模拟图像的空间分辨率差异,针对较高的空间分辨率图像数据,确定其像元与对应待模拟图像的像元关系,通过像元的空间退化计算,进而通过像元求平均,来实现空间退化模拟,例如:若待模拟图像空间分辨率为10m,基础图像空间分辨率为5m,则将基础图像上4个像元求平均作为待模拟图像上的一个像元来实现图像空间退化模拟。但是这种像元求平均值来实现空间退化模拟的方法得到的图像分辨率无法满足新型遥感器的图像数据模拟精度或无法为模型算法研究提供需要的图像数据。
发明内容
有鉴于此,本发明提供了一种图像空间退化模拟方法及系统,以解决现有技术中由于采用像元求平均值来实现空间退化模拟的方法得到的图像分辨率无法满足新型遥感器的图像数据模拟精度或无法为模型算法研究提供需要的图像数据的问题。
为实现上述目的,本发明提供如下技术方案:
本发明公开了一种图像空间退化模拟方法,包括:
获取基础图像数据、基础图像空间分辨率、待模拟图像奈奎斯特频率下的调制传递函数的参数值以及待模拟图像空间分辨率;
根据所述待模拟图像奈奎斯特频率下的调制传递函数的参数值,求得有效瞬时视场角和调制传递函数的傅里叶变换函数,并对所述傅里叶变换函数的参数值求解;
根据所述基础图像空间分辨率、所述待模拟图像空间分辨率以及所述傅里叶变换函数求解的参数值构建点扩展函数矩阵;
将所述点扩展函数矩阵对应的曲线作为权重,与所述基础图像数据进行空间卷积计算,得到图像空间退化后的模拟图像。
优选的,所述根据所述待模拟图像奈奎斯特频率下的调制传递函数的参数值,求得有效瞬时视场角和调制传递函数的傅里叶变换函数,并对所述傅里叶变换函数求解的参数值具体为:
获取所述待模拟图像奈奎斯特频率下的调制传递函数的参数值对应参数向量值;
预设所述调制传递函数为高斯型函数,根据所述调制传递函数的参数值对应参数向量值对每个波段依次进行所述傅里叶变换函数的参数系数求解;
根据所述调制传递函数和所述参数系数求解对应的所述有效瞬时视场角EIFOV;
根据所述有效瞬时视场角EIFOV求解所述傅里叶变换函数的参数值。
优选的,所述根据所述基础图像空间分辨率、所述待模拟图像空间分辨率以及所述傅里叶变换函数求解的参数值构建点扩展函数矩阵具体为:
根据所述基础图像空间分辨率和所述待模拟图像空间分辨率得到空间分辨率关系因子;
根据所述空间分辨率关系因子和所述有效瞬时视场角求得所述点扩展函数矩阵的大小;
根据所述点扩展函数矩阵的大小和所述有效瞬时视场角求得所述点扩展函数矩阵中的自变量对应的输入矩阵;
根据所述输入矩阵和所述调制传递函数的参数值求得所述点扩展函数矩阵,并对所述点扩展函数矩阵进行归一化,得到归一化的点扩展函数矩阵。
优选的,所述将所述点扩展函数矩阵对应的曲线作为权重,与所述基础图像数据进行空间卷积计算,得到图像空间退化后的模拟图像具体为:
根据所述空间分辨率关系因子和所述点扩展函数矩阵,求解输出图像行列数;
根据所述有效瞬时视场角和所述输出图像行列数确定输出图像对应输入图像的像元增量;
根据所述像元增量和所述归一化的点扩展函数矩阵进行图像空间退化模拟,得到图像空间退化后的模拟图像。
在上述公开的方法的基础上,本发明还公开了一种图像空间退化模拟装置,包括:
获取模块,用于获取基础图像数据、基础图像空间分辨率、待模拟图像奈奎斯特频率下的调制传递函数的参数值以及待模拟图像空间分辨率;
求解模块,用于根据所述待模拟图像奈奎斯特频率下的调制传递函数的参数值,求得有效瞬时视场角和调制传递函数的傅里叶变换函数,并对所述傅里叶变换函数的参数值求解;
构建模块,用于根据所述基础图像空间分辨率、所述待模拟图像空间分辨率以及所述傅里叶变换函数求解的参数值构建点扩展函数矩阵;
模型退化模块,用于将所述点扩展函数矩阵对应的曲线作为权重,与所述基础图像数据进行空间卷积计算,得到图像空间退化后的模拟图像。
优选的,所述求解模块包括:
第一获取模块,用于获取所述待模拟图像奈奎斯特频率下的调制传递函数的参数值对应参数向量值;
第一求解模块,用于预设所述调制传递函数为高斯型函数,根据所述调制传递函数的参数值对应参数向量值对每个波段依次进行所述傅里叶变换函数的参数系数求解;
第二求解模块,用于根据所述调制传递函数和所述参数系数求解对应的所述有效瞬时视场角EIFOV;
第三求解模块,用于根据所述有效瞬时视场角EIFOV求解所述傅里叶变换函数的参数值。
优选的,所述构建模块包括:
第四求解模块,用于根据所述基础图像空间分辨率和所述待模拟图像空间分辨率得到空间分辨率关系因子;
第五求解模块,用于根据所述空间分辨率关系因子和所述有效瞬时视场角求得所述点扩展函数矩阵的大小;
第六求解模块,用于根据所述点扩展函数矩阵的大小和所述有效瞬时视场角求得所述点扩展函数矩阵中的自变量对应的输入矩阵;
归一化模块,用于根据所述输入矩阵和所述调制传递函数的参数值求得所述点扩展函数矩阵,并对所述点扩展函数矩阵进行归一化,得到归一化的点扩展函数矩阵。
优选的,所述模型退化模块包括:
第七求解模块,用于根据所述空间分辨率关系因子和所述点扩展函数矩阵,求解输出图像行列数;
第八求解模块,用于根据所述有效瞬时视场角和所述输出图像行列数确定输出图像对应输入图像的像元增量;
退化模块,用于根据所述像元增量和所述归一化的点扩展函数矩阵进行图像空间退化模拟,得到图像空间退化后的模拟图像。
经由上述的技术方案可知,与现有技术相比,本发明公开了一种图像空间退化模拟方法,该方法是依据待模拟图像奈奎斯特(Nyquist)频率下的MTF值,以及基础图像和待模拟图像的空间分辨率值,针对基础图像进行空间退化模拟,获得空间退化后的待模拟图像。处理步骤为:首先,基于待模拟图像奈奎斯特频率下的MTF值,求得其有效瞬时视场角,进而求得MTF函数的傅里叶变换函数(即PSF函数)的参数;其次,根据基础图像的空间分辨率、待模拟图像的空间分辨率,计到空间比值,利用PSF函数参数,计算得到PSF二维矩阵;最后,将PSF二维矩阵作为权重,与输入图像进行卷积计算,进行逐像元的空间退化计算,最终获得考虑MTF效应空间退化后的待模拟图像,上述方法通过考虑MTF或PSF对于图像空间实际分辨能力的影响,解决了现有技术中由于采用像元求平均值来实现空间退化模拟的方法得到的图像分辨率无法满足新型遥感器的图像数据模拟精度或无法为模型算法研究提供需要的图像数据的问题。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1为本发明实施例公开的一种图像空间退化模拟方法流程图;
图2为本发明实施例公开的图1中步骤S102的方法流程图;
图3为本发明实施例公开的图1中步骤S103的方法流程图;
图4为本发明实施例公开的图1中步骤S104的方法流程图;
图5为本发明实施例公开的一种图像空间退化模拟装置的结构示意图;
图6为本发明实施例公开的图5中的求解模块的具体结构示意图;
图7为本发明实施例公开的图5中的构建模块的具体结构示意图;
图8为本发明实施例公开的图5中的模型退化模块的具体结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
传感器的PSF(point spread function,点扩展函数)或MTF(modulationtransfer function,调制传递函数)是传感器成像性能的重要指标,其中PSF是信号扩散与削弱在空间域的描述,而MTF是PSF在频率域的描述。MTF或PSF是用一个函数去表达,反映的是系统在对目标成像过程中的信号扩散与削弱程度,是决定图像质量及其应用效能的重要参量。因此,在对图像进行空间退化模拟,应考虑MTF或PSF对于图像空间实际分辨能力的影响。
具体的,本发明采用以下实施例手段进行实现:
获取图像模拟对应的基础图像imgin(nl,ns,nb),其中nl表示图像行数,ns表示图像列数,nb为图像波段数;图像模拟对应的基础图像的空间分辨率SRin;待模拟图像,即空间退化后图像的空间分辨率SRout;待模拟图像的奈奎斯特频率下的MTF参数值mtf_n(nb)。
预先设定MTF函数为高斯型,即函数根据给定的奈奎斯特频率下的MTF参数值,求得系数c,第k波段对应的c(k)为:
通过计算得到有效瞬时视场角EIFOV,作为图像实际的分辨率:
其中,x0.5(k)为上述高斯函数值等于0.5时候对应的x值即:
根据EIFOV,求得图像的二维PSF函数,该函数为高斯概率密度函数:
该函数公式中对应到波段k的σ(k)为:
具体的,根据基础图像、待模拟图像的空间分辨率,以及每个波段的有效瞬时视场角,构建二维PSF函数矩阵。对于波段k有:
求得空间分辨率比例关系:
二维PSF函数矩阵为方阵,其中一维单元数nk(k)为:
temp(k)=ceil(2*EIFOV(k)*scale+1),
其中,ceil表示向上取整。如果temp(k)为奇数,则nk(k)=temp(k);如果temp(k)为偶数,则nk(k)=temp(k)+1。
求得函数矩阵中自变量x(k),y(k)对应的输入矩阵X(k),Y(k),为nk(k)*nk(k)的方阵:
x1(k)=-(nk(k)-1)/2:1:(nk(k)-1)/2
X(1,:)=X(2,:)=...=X(nk(k),:)=x1(k)*2*EIFOV(k)/nk(k)
Y(k)=X(k)
根据自变量x(k),y(k)对应的输入矩阵X(k),Y(k),求得二维PSF概率密度函数矩阵为:
其中,Z(k)为nk(k)*nk(k)的方阵。
对PSF概率密度函数矩阵Z(k)求和:
factor(k)=sum(Z(k))
对PSF概率密度函数矩阵Z(k)进行归一化:
Z1(k)=Z(k)/factor(k)
优选的,确认输出图像imgout的大小,列数nsout,行数nlout为:
nsout=floor(ns/scale)
nlout=floor(nl/scale)
其中,floor表示向下取整的意思,如floor(10.5)为10。
最后,针对输出图像的每个点基于PSF概率密度函数矩阵Z(k)进行卷积求解,实现基于MTF参数的图像空间退化模拟。
主要的,针对波段k有:
确定像元增量:
step(k)=2*EIFOV(k)*scale/nk(k)
Z1(k)作为核函数对图像进行卷积计算,对于第k波段,针对输出图像的第m,n像元,通过以下过程进行计算得到空间退化后的结果:
求得该像元对应基础图像上的空间位置:
xx=(scale+1)/2+(m-1)*scale;
yy=(scale+1)/2+(n-1)*scale;
以xx,yy为中心位置,求得二维PSF函数每个点在基础图像上的对应的空间位置,对于该矩阵中第ki,kj点的像元位置(nx(ki),ny(kj))有:
var1(k)=(nk(k)-1)/2
var(k)=(nk(k)-1)/2
nx(ki)=xx+(ki-1-var1(k))*step(k),ki=1:nk(k)
ny(kj)=yy+(kj-1-var1(k))*step(k),kj=1:nk(k)
求得矩阵中像元位置(nx(ki),ny(kj))在基础图像对应的像元值,对于第k波段imgin(nx(ki),ny(kj),k)为:
x1=floor(nx(ki));
x2=ceil(nx(ki));
y1=floor(ny(kj));
y2=ceil(ny(kj));
dx=x2-nx(ki);
dy=y2-ny(kj);
tmp=dx*imgin(x1,y1,k)+(1-dx)*imgin(x2,y1,k);
具体的,imgin(nx(ki),ny(kj),k)=dy*tmp+(1-dy)*(dx*imgin(x1,y2,k)+(1-dx)*imgin(x2,y2,k))。
以归一化后的二维PSF函数矩阵Z1(nk(k),nk(k))为核函数,通过卷积运算求解空间退化后的imgout(m,n,k)的像元值:
请参阅附图1,图1为本发明公开的一种图像空间退化模拟方法流程图。如图1所示,本发明公开了一种图像空间退化模拟方法,该方法具体包括如下步骤:
S101、获取基础图像数据、基础图像空间分辨率、待模拟图像奈奎斯特频率下的调制传递函数的参数值以及待模拟图像空间分辨率。
其中,图像模拟对应的基础图像imgin,三维矩阵,nl*ns*nb,其中nl表示图像行数,ns表示图像列数,nb为图像波段数;图像模拟对应的基础图像的空间分辨率变量SRin;待模拟图像,即空间退化后图像的空间分辨率变量SRout;待模拟图像的奈奎斯特频率下的MTF参数值向量mtf_n,长度为波段数nb。
S102、根据待模拟图像奈奎斯特频率下的调制传递函数的参数值,求得有效瞬时视场角和调制传递函数的傅里叶变换函数,并对傅里叶变换函数的参数值求解。
基于奈奎斯特频率下的MTF值,求得有效瞬时视场角,进而求得MTF函数的傅里叶变换函数,即PSF函数。
S103、根据基础图像空间分辨率待模拟图像空间分辨率以及傅里叶变换函数求解的参数值构建点扩展函数矩阵。
S104、将点扩展函数矩阵对应的曲线作为权重,与基础图像数据进行空间卷积计算,得到图像空间退化后的模拟图像。
根据基础图像的空间分辨率、待模拟图像的空间分辨率比值,将PSF函数曲线作为权重,与输入图像进行卷积计算,进行逐像元的空间退化计算,最终获得空间退化后的模拟图像。
具体的,请参阅附图2,图2为图1中步骤S102的方法流程图。具体的,步骤S102为如下步骤:
S201、获取待模拟图像奈奎斯特频率下的调制传递函数的参数值对应参数向量值。
S202、预设调制传递函数为高斯型函数,根据调制传递函数的参数值对应参数向量值对每个波段依次进行傅里叶变换函数的参数系数求解。
S203、根据调制传递函数和参数系数求解对应的有效瞬时视场角EIFOV。
S204、根据有效瞬时视场角EIFOV求解傅里叶变换函数的参数值。
上述步骤可以描述为以下内容:
输入待模拟图像的奈奎斯特频率下的MTF参数值向量mtf_n(nb)。
设定MTF函数为高斯型,即函数根据给定的奈奎斯特频率下的MTF参数值,对每个波段依次进行PSF函数参数系数求解,第k波段对应的c(k)为:
计算得到有效瞬时视场角EIFOV,作为图像实际的分辨率:
其中,x0.5(k)为上述高斯函数值等于0.5时候对应的x值即:
根据EIFOV,求得图像的二维PSF函数,该函数为高斯概率密度函数该函数公式中对应到波段k的σ(k)为:
具体的,请参阅附图3,图3为图1中步骤S103的方法流程图。具体的,步骤S103为如下步骤:
S301、根据基础图像空间分辨率和待模拟图像空间分辨率得到空间分辨率关系因子。
S302、根据空间分辨率关系因子和有效瞬时视场角求得点扩展函数矩阵的大小。
S303、根据点扩展函数矩阵的大小和有效瞬时视场角求得点扩展函数矩阵中的自变量对应的输入矩阵。
S304、根据输入矩阵和调制传递函数的参数值求得点扩展函数矩阵,并对点扩展函数矩阵进行归一化,得到归一化的点扩展函数矩阵。
上述步骤可以描述为以下内容:
输入基础图像的空间分辨率变量SRin;输入待模拟图像,即空间退化后图像的空间分辨率变量SRout;输入有效瞬时视场角EIFOV(nb);输入二维PSF函数参数σ(nb)。
求得空间分辨率关系因子:
求得二维PSF函数矩阵大小,该矩阵方阵,对于第k波段,其某一维单元数nk(k)为:
temp(k)=ceil(2*EIFOV(k)*scale+1),
其中ceil表示向上取整;如果temp(k)为奇数,则nk(k)=temp(k);如果temp(k)为偶数,则nk(k)=temp(k)+1。
求得二维PSF函数矩阵中自变量x,y对应的输入矩阵X、Y,为nk(k)*nk(k)的方阵,对于第k波段X(k)、Y(k)为:
x1(k)=-(nk(k)-1)/2:1:(nk(k)-1)/2
X(1,:)=X(2,:)=...=X(nk(k),:)=x1(k)*2*EIFOV(k)/nk(k)
Y(k)=X(k)
求得二维PSF函数矩阵,对于第k波段,根据自变量x(k),y(k)对应的输入矩阵X(k),Y(k),有:
其中,Z(k)为nk(k)*nk(k)的方阵。
对PSF函数矩阵进行归一化,具体为:
对PSF函数矩阵求和,对于第k波段有:
factor(k)=sum(Z(k))
对PSF概率密度函数矩阵进行归一化,对于第k波段有:
Z1(k)=Z(k)/factor(k)
最后得到归一化的点扩展函数矩阵。
具体的,请参阅附图4,图4为图1中步骤S104的方法流程图。具体的,步骤S104为如下步骤:
S401、根据空间分辨率关系因子和点扩展函数矩阵,求解输出图像行列数。
S402、根据有效瞬时视场角和输出图像行列数确定输出图像对应输入图像的像元增量。
S403、根据像元增量和归一化的点扩展函数矩阵进行图像空间退化模拟,得到图像空间退化后的模拟图像。
上述步骤可以描述为以下内容:
根据空间分辨率关系因子scale,输入基础图像矩阵imgin(nl,ns,nb)、每个波段的有效瞬时视场角EIFOV(nb)、每个波段的归一化后的二维PSF函数矩阵,对于第k波段矩阵为Z1(nk(k),nk(k))。
求解输出图像行列数,数nlout分别为:
nlout=floor(nl/scale)
其中,floor表示向下取整的意思,如floor(10.5)为10。
确定输出图像单个像元对应输入图像的像元增量:
step(k)=2*EIFOV(k)*scale/nk(k)
逐像元的求得待模拟图像像元值,进行基于MTF参数的图像空间退化模拟。对于第k波段,针对输出图像的第m,n像元,通过以下过程进行计算得到空间退化后的结果:
求得该像元对应基础图像上的空间位置:
xx=(scale+1)/2+(m-1)*scale;
yy=(scale+1)/2+(n-1)*scale;
以xx,yy为中心位置,求得二维PSF函数每个点在基础图像上的对应的空间位置,对于该矩阵中第ki,kj点的像元位置(nx(ki),ny(kj))有:
var1(k)=(nk(k)-1)/2
nx(ki)=xx+(ki-1-var1(k))*step(k),ki=1:nk(k)
ny(kj)=yy+(kj-1-var1(k))*step(k),kj=1:nk(k)
求得矩阵中像元位置(nx(ki),ny(kj))在基础图像对应的像元值,对于第k波段imgin(nx(ki),ny(kj),k)为:
x1=floor(nx(ki));
x2=ceil(nx(ki));
y1=floor(ny(kj));
y2=ceil(ny(kj));
dx=x2-nx(ki);
dy=y2-ny(kj);
tmp=dx*imgin(x1,y1,k)+(1-dx)*imgin(x2,y1,k);
imgin(nx(ki),
ny(kj),k)=dy*tmp+(1-dy)*(dx*imgin(x1,y2,k)+(1-dx)*imgin(x2,y2,k));
以归一化后的二维PSF函数矩阵Z1(nk(k),nk(k))为核函数,通过卷积运算求解空间退化后的imgout(m,n,k)的像元值:
本发明公开了一种图像空间退化模拟方法,该方法是依据待模拟图像奈奎斯特(Nyquist)频率下的MTF值,以及基础图像和待模拟图像的空间分辨率值,针对基础图像进行空间退化模拟,获得空间退化后的待模拟图像。本发明通过考虑MTF或PSF对于图像空间实际分辨能力的影响,解决了现有技术中由于采用像元求平均值来实现空间退化模拟的方法得到的图像分辨率无法满足新型遥感器的图像数据模拟精度或无法为模型算法研究提供需要的图像数据的问题。
在上述公开的方法的基础上,请参阅附图5,图5为本发明公开的一种图像空间退化模拟装置的结构示意图。本发明还公开了一种图像空间退化模拟装置,具体的该装置包括:
获取模块501,用于获取基础图像数据、基础图像空间分辨率、待模拟图像奈奎斯特频率下的调制传递函数的参数值以及待模拟图像空间分辨率。
求解模块502,用于根据待模拟图像奈奎斯特频率下的调制传递函数的参数值,求得有效瞬时视场角和调制传递函数的傅里叶变换函数,并对傅里叶变换函数的参数值求解。
构建模块503,用于根据基础图像空间分辨率、待模拟图像空间分辨率以及傅里叶变换函数求解的参数值构建点扩展函数矩阵。
模型退化模块504,用于将点扩展函数矩阵对应的曲线作为权重,与基础图像数据进行空间卷积计算,得到图像空间退化后的模拟图像。
本实施例的装置,可以用于执行图1所示方法实施例的技术方案,其实现原理和技术效果类似,此处不再赘述。
具体的,请参阅附图6,图6为图5中的求解模块的具体结构示意图,所述求解模块502包括:
第一获取模块601,用于获取待模拟图像奈奎斯特频率下的调制传递函数的参数值对应参数向量值。
第一求解模块602,用于预设调制传递函数为高斯型函数,根据调制传递函数的参数值对应参数向量值对每个波段依次进行傅里叶变换函数的参数系数求解。
第二求解模块603,用于根据调制传递函数和参数系数求解对应的有效瞬时视场角EIFOV。
第三求解模块604,用于根据有效瞬时视场角EIFOV求解傅里叶变换函数的参数值。
本实施例的装置,可以用于执行图2所示方法实施例的技术方案,其实现原理和技术效果类似,此处不再赘述。
具体的,请参阅附图7,图7为图5中的构建模块的具体结构示意图,所述构建模块503包括:
第四求解模块701,用于根据基础图像空间分辨率和待模拟图像空间分辨率得到空间分辨率关系因子。
第五求解模块702,用于根据空间分辨率关系因子和有效瞬时视场角求得点扩展函数矩阵的大小。
第六求解模块703,用于根据点扩展函数矩阵的大小和有效瞬时视场角求得点扩展函数矩阵中的自变量对应的输入矩阵。
归一化模块704,用于根据输入矩阵和调制传递函数的参数值求得点扩展函数矩阵,并对点扩展函数矩阵进行归一化,得到归一化的点扩展函数矩阵。
本实施例的装置,可以用于执行图3所示方法实施例的技术方案,其实现原理和技术效果类似,此处不再赘述。
具体的,请参阅附图8,图8为图5中的模型退化模块的具体结构示意图,所述模型退化模块504包括:
第七求解模块801,用于根据空间分辨率关系因子和点扩展函数矩阵,求解输出图像行列数。
第八求解模块802,用于根据有效瞬时视场角和输出图像行列数确定输出图像对应输入图像的像元增量。
退化模块803,用于根据像元增量和归一化的点扩展函数矩阵进行图像空间退化模拟,得到图像空间退化后的模拟图像。
本实施例的装置,可以用于执行图4所示方法实施例的技术方案,其实现原理和技术效果类似,此处不再赘述。
综上所述,本发明公开了一种图像空间退化模拟方法,该方法是依据待模拟图像奈奎斯特(Nyquist)频率下的MTF值,以及基础图像和待模拟图像的空间分辨率值,针对基础图像进行空间退化模拟,获得空间退化后的待模拟图像。处理步骤为:首先,基于待模拟图像奈奎斯特频率下的MTF值,求得其有效瞬时视场角,进而求得MTF函数的傅里叶变换函数(即PSF函数)的参数;其次,根据基础图像的空间分辨率、待模拟图像的空间分辨率,计到空间比值,利用PSF函数参数,计算得到PSF二维矩阵;最后,将PSF二维矩阵作为权重,与输入图像进行卷积计算,进行逐像元的空间退化计算,最终获得考虑MTF效应空间退化后的待模拟图像,上述方法通过考虑MTF或PSF对于图像空间实际分辨能力的影响,解决了现有技术中由于采用像元求平均值来实现空间退化模拟的方法得到的图像分辨率无法满足新型遥感器的图像数据模拟精度或无法为模型算法研究提供需要的图像数据的问题。另外,本发明还公开了一种图像空间退化模拟装置。
需要说明的是,本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (8)
1.一种图像空间退化模拟方法,其特征在于,包括:
获取基础图像数据、基础图像空间分辨率、待模拟图像奈奎斯特频率下的调制传递函数的参数值以及待模拟图像空间分辨率;
根据所述待模拟图像奈奎斯特频率下的调制传递函数的参数值,求得有效瞬时视场角和调制传递函数的傅里叶变换函数,并对所述傅里叶变换函数的参数值求解;
根据所述基础图像空间分辨率、所述待模拟图像空间分辨率以及所述傅里叶变换函数求解的参数值构建点扩展函数矩阵;
将所述点扩展函数矩阵对应的曲线作为权重,与所述基础图像数据进行空间卷积计算,得到图像空间退化后的模拟图像。
2.根据权利要求1所述的方法,其特征在于,所述根据所述待模拟图像奈奎斯特频率下的调制传递函数的参数值,求得有效瞬时视场角和调制传递函数的傅里叶变换函数,并对所述傅里叶变换函数求解的参数值具体为:
获取所述待模拟图像奈奎斯特频率下的调制传递函数的参数值对应参数向量值;
预设所述调制传递函数为高斯型函数,根据所述调制传递函数的参数值对应参数向量值对每个波段依次进行所述傅里叶变换函数的参数系数求解;
根据所述调制传递函数和所述参数系数求解对应的所述有效瞬时视场角EIFOV;
根据所述有效瞬时视场角EIFOV求解所述傅里叶变换函数的参数值。
3.根据权利要求1所述的方法,其特征在于,所述根据所述基础图像空间分辨率、所述待模拟图像空间分辨率以及所述傅里叶变换函数求解的参数值构建点扩展函数矩阵具体为:
根据所述基础图像空间分辨率和所述待模拟图像空间分辨率得到空间分辨率关系因子;
根据所述空间分辨率关系因子和所述有效瞬时视场角求得所述点扩展函数矩阵的大小;
根据所述点扩展函数矩阵的大小和所述有效瞬时视场角求得所述点扩展函数矩阵中的自变量对应的输入矩阵;
根据所述输入矩阵和所述调制传递函数的参数值求得所述点扩展函数矩阵,并对所述点扩展函数矩阵进行归一化,得到归一化的点扩展函数矩阵。
4.根据权利要求3所述方法,其特征在于,所述将所述点扩展函数矩阵对应的曲线作为权重,与所述基础图像数据进行空间卷积计算,得到图像空间退化后的模拟图像具体为:
根据所述空间分辨率关系因子和所述点扩展函数矩阵,求解输出图像行列数;
根据所述有效瞬时视场角和所述输出图像行列数确定输出图像对应输入图像的像元增量;
根据所述像元增量和所述归一化的点扩展函数矩阵进行图像空间退化模拟,得到图像空间退化后的模拟图像。
5.一种图像空间退化模拟装置,其特征在于,包括:
获取模块,用于获取基础图像数据、基础图像空间分辨率、待模拟图像奈奎斯特频率下的调制传递函数的参数值以及待模拟图像空间分辨率;
求解模块,用于根据所述待模拟图像奈奎斯特频率下的调制传递函数的参数值,求得有效瞬时视场角和调制传递函数的傅里叶变换函数,并对所述傅里叶变换函数的参数值求解;
构建模块,用于根据所述基础图像空间分辨率、所述待模拟图像空间分辨率以及所述傅里叶变换函数求解的参数值构建点扩展函数矩阵;
模型退化模块,用于将所述点扩展函数矩阵对应的曲线作为权重,与所述基础图像数据进行空间卷积计算,得到图像空间退化后的模拟图像。
6.根据权利要求5所述的装置,其特征在于,所述求解模块包括:
第一获取模块,用于获取所述待模拟图像奈奎斯特频率下的调制传递函数的参数值对应参数向量值;
第一求解模块,用于预设所述调制传递函数为高斯型函数,根据所述调制传递函数的参数值对应参数向量值对每个波段依次进行所述傅里叶变换函数的参数系数求解;
第二求解模块,用于根据所述调制传递函数和所述参数系数求解对应的所述有效瞬时视场角EIFOV;
第三求解模块,用于根据所述有效瞬时视场角EIFOV求解所述傅里叶变换函数的参数值。
7.根据权利要求5所述的装置,其特征在于,所述构建模块包括:
第四求解模块,用于根据所述基础图像空间分辨率和所述待模拟图像空间分辨率得到空间分辨率关系因子;
第五求解模块,用于根据所述空间分辨率关系因子和所述有效瞬时视场角求得所述点扩展函数矩阵的大小;
第六求解模块,用于根据所述点扩展函数矩阵的大小和所述有效瞬时视场角求得所述点扩展函数矩阵中的自变量对应的输入矩阵;
归一化模块,用于根据所述输入矩阵和所述调制传递函数的参数值求得所述点扩展函数矩阵,并对所述点扩展函数矩阵进行归一化,得到归一化的点扩展函数矩阵。
8.根据权利要求7所述装置,其特征在于,所述模型退化模块包括:
第七求解模块,用于根据所述空间分辨率关系因子和所述点扩展函数矩阵,求解输出图像行列数;
第八求解模块,用于根据所述有效瞬时视场角和所述输出图像行列数确定输出图像对应输入图像的像元增量;
退化模块,用于根据所述像元增量和所述归一化的点扩展函数矩阵进行图像空间退化模拟,得到图像空间退化后的模拟图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610274396.8A CN105976317B (zh) | 2016-04-28 | 2016-04-28 | 一种图像空间退化模拟方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610274396.8A CN105976317B (zh) | 2016-04-28 | 2016-04-28 | 一种图像空间退化模拟方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105976317A true CN105976317A (zh) | 2016-09-28 |
CN105976317B CN105976317B (zh) | 2019-03-19 |
Family
ID=56993797
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610274396.8A Active CN105976317B (zh) | 2016-04-28 | 2016-04-28 | 一种图像空间退化模拟方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105976317B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110460765A (zh) * | 2018-05-02 | 2019-11-15 | 展讯通信(上海)有限公司 | 空间频率响应曲线的获取方法、装置及电子设备 |
CN110954303A (zh) * | 2019-11-20 | 2020-04-03 | 东南大学 | 一种基于高分辨率遥感图像基准的mtf自动测算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101441764A (zh) * | 2008-12-31 | 2009-05-27 | 中国资源卫星应用中心 | 一种mtfc遥感图像复原方法 |
CN104298844A (zh) * | 2014-05-23 | 2015-01-21 | 中国科学院光电研究院 | 获取点阵法测量光学遥感载荷在轨mtf测量精度的方法 |
WO2015195417A1 (en) * | 2014-06-20 | 2015-12-23 | Rambus Inc. | Systems and methods for lensed and lensless optical sensing |
CN105427267A (zh) * | 2015-10-28 | 2016-03-23 | 中国科学院光电研究院 | 获取刃边法测量光学遥感载荷在轨mtf精度的方法 |
-
2016
- 2016-04-28 CN CN201610274396.8A patent/CN105976317B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101441764A (zh) * | 2008-12-31 | 2009-05-27 | 中国资源卫星应用中心 | 一种mtfc遥感图像复原方法 |
CN104298844A (zh) * | 2014-05-23 | 2015-01-21 | 中国科学院光电研究院 | 获取点阵法测量光学遥感载荷在轨mtf测量精度的方法 |
WO2015195417A1 (en) * | 2014-06-20 | 2015-12-23 | Rambus Inc. | Systems and methods for lensed and lensless optical sensing |
CN105427267A (zh) * | 2015-10-28 | 2016-03-23 | 中国科学院光电研究院 | 获取刃边法测量光学遥感载荷在轨mtf精度的方法 |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110460765A (zh) * | 2018-05-02 | 2019-11-15 | 展讯通信(上海)有限公司 | 空间频率响应曲线的获取方法、装置及电子设备 |
CN110460765B (zh) * | 2018-05-02 | 2020-12-18 | 展讯通信(上海)有限公司 | 空间频率响应曲线的获取方法、装置及电子设备 |
CN110954303A (zh) * | 2019-11-20 | 2020-04-03 | 东南大学 | 一种基于高分辨率遥感图像基准的mtf自动测算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105976317B (zh) | 2019-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3557484B1 (en) | Neural network convolution operation device and method | |
US10839483B2 (en) | Method of converting low-resolution image to high-resolution image and image conversion device performing method | |
US7991240B2 (en) | Methods, systems and apparatuses for modeling optical images | |
US20090257680A1 (en) | Method and Device for Video Stitching | |
CN104052979B (zh) | 用于图像处理的装置和技术 | |
CN107016210B (zh) | 一种菲涅尔衍射光学系统的成像质量仿真方法 | |
CN103985085A (zh) | 图像超分辨率放大的方法和装置 | |
CN102682424B (zh) | 基于边缘方向差值的图像放大处理方法 | |
CN110838085B (zh) | 一种图像的超分辨率重建方法、装置及电子设备 | |
Oppermann et al. | Reconstructing signals from noisy data with unknown signal and noise covariance | |
CN101790069A (zh) | 一种基于图像边缘方向的尺度变换方法 | |
Pratley et al. | A Fast and Exact w-stacking and w-projection Hybrid Algorithm for Wide-field Interferometric Imaging | |
CN105869141A (zh) | 图像配准方法及装置 | |
US9251575B2 (en) | Image processing apparatus, image pickup apparatus, image processing method, and image processing program | |
CN105976317A (zh) | 一种图像空间退化模拟方法及系统 | |
US20200334791A1 (en) | Image down-scaling with pixel sets selected via blue noise sampling | |
US10521918B2 (en) | Method and device for filtering texture, using patch shift | |
CN106470293B (zh) | 图像处理装置及图像处理方法 | |
US20150254815A1 (en) | Image downsampling apparatus and method | |
CN104574315B (zh) | 基于光强传输矩阵的光学系统成像复原方法 | |
Krumm et al. | Sampled-grating and crossed-grating models of moire patterns from digital imaging | |
KR101582578B1 (ko) | 그래픽 처리 장치 및 방법 | |
US20150213578A1 (en) | Method for electronic zoom with sub-pixel offset | |
CN112907450B (zh) | 三维时序图像处理方法、装置、计算机设备和存储介质 | |
CN113126944B (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 |