一种图像空间退化模拟方法及系统的制作方法
【专利摘要】本发明公开了一种图像空间退化模拟方法，该方法是依据待模拟图像奈奎斯特频率下的MTF值，以及基础图像和待模拟图像的空间分辨率值，针对基础图像进行空间退化模拟，获得空间退化后的待模拟图像。本发明通过考虑MTF或PSF对于图像空间实际分辨能力的影响，解决了现有技术中由于采用像元求平均值来实现空间退化模拟的方法得到的图像分辨率无法满足新型遥感器图像数据模拟精度或无法为模型算法研究提供需要的图像数据的问题。
【专利说明】
一种图像空间退化模拟方法及系统
技术领域
[0001] 本发明涉及遥感传感器成像技术领域，更具体的说，是涉及一种图像空间退化模 拟方法及系统。
【背景技术】
[0002] 随着遥感技术和计算机技术的快速发展，遥感成像产品在各行各业中起着越来越 重要的作用，因此遥感成像分辨率越来越高，数据量也越来越大。当需要模拟特定空间分辨 率的图像时，一般是以较高空间分辨率的图像数据为基础图像，对其进行计算处理获得对 应的空间分辨率图像，这个过程为图像的空间退化模拟，它一般可用于新研制遥感器的图 像数据模拟或为模型算法研究提供需要的图像数据。
[0003] 现有技术中，图像空间退化模拟根据较高的空间分辨率图像与待模拟图像的空间 分辨率差异，针对较高的空间分辨率图像数据，确定其像元与对应待模拟图像的像元关系， 通过像元的空间退化计算，进而通过像元求平均，来实现空间退化模拟，例如：若待模拟图 像空间分辨率为l〇m，基础图像空间分辨率为5m，则将基础图像上4个像元求平均作为待模 拟图像上的一个像元来实现图像空间退化模拟。但是这种像元求平均值来实现空间退化模 拟的方法得到的图像分辨率无法满足新型遥感器的图像数据模拟精度或无法为模型算法 研究提供需要的图像数据。
【发明内容】
[0004] 有鉴于此，本发明提供了一种图像空间退化模拟方法及系统，以解决现有技术中 由于采用像元求平均值来实现空间退化模拟的方法得到的图像分辨率无法满足新型遥感 器的图像数据模拟精度或无法为模型算法研究提供需要的图像数据的问题。
[0005] 为实现上述目的，本发明提供如下技术方案：
[0006] 本发明公开了 一种图像空间退化模拟方法，包括：
[0007] 获取基础图像数据、基础图像空间分辨率、待模拟图像奈奎斯特频率下的调制传 递函数的参数值以及待模拟图像空间分辨率；
[0008] 根据所述待模拟图像奈奎斯特频率下的调制传递函数的参数值，求得有效瞬时视 场角和调制传递函数的傅里叶变换函数，并对所述傅里叶变换函数的参数值求解；
[0009] 根据所述基础图像空间分辨率、所述待模拟图像空间分辨率以及所述傅里叶变换 函数求解的参数值构建点扩展函数矩阵；
[0010] 将所述点扩展函数矩阵对应的曲线作为权重，与所述基础图像数据进行空间卷积 计算，得到图像空间退化后的模拟图像。
[0011]优选的，所述根据所述待模拟图像奈奎斯特频率下的调制传递函数的参数值，求 得有效瞬时视场角和调制传递函数的傅里叶变换函数，并对所述傅里叶变换函数求解的参 数值具体为：
[0012]获取所述待模拟图像奈奎斯特频率下的调制传递函数的参数值对应参数向量值；
[0013] 预设所述调制传递函数为高斯型函数，根据所述调制传递函数的参数值对应参数 向量值对每个波段依次进行所述傅里叶变换函数的参数系数求解；
[0014] 根据所述调制传递函数和所述参数系数求解对应的所述有效瞬时视场角EIF0V; [0015]根据所述有效瞬时视场角EIF0V求解所述傅里叶变换函数的参数值。
[0016]优选的，所述根据所述基础图像空间分辨率、所述待模拟图像空间分辨率以及所 述傅里叶变换函数求解的参数值构建点扩展函数矩阵具体为：
[0017] 根据所述基础图像空间分辨率和所述待模拟图像空间分辨率得到空间分辨率关 系因子；
[0018] 根据所述空间分辨率关系因子和所述有效瞬时视场角求得所述点扩展函数矩阵 的大小；
[0019] 根据所述点扩展函数矩阵的大小和所述有效瞬时视场角求得所述点扩展函数矩 阵中的自变量对应的输入矩阵；
[0020] 根据所述输入矩阵和所述调制传递函数的参数值求得所述点扩展函数矩阵，并对 所述点扩展函数矩阵进行归一化，得到归一化的点扩展函数矩阵。
[0021] 优选的，所述将所述点扩展函数矩阵对应的曲线作为权重，与所述基础图像数据 进行空间卷积计算，得到图像空间退化后的模拟图像具体为：
[0022] 根据所述空间分辨率关系因子和所述点扩展函数矩阵，求解输出图像行列数； [0023]根据所述有效瞬时视场角和所述输出图像行列数确定输出图像对应输入图像的 像元增量；
[0024] 根据所述像元增量和所述归一化的点扩展函数矩阵进行图像空间退化模拟，得到 图像空间退化后的模拟图像。
[0025] 在上述公开的方法的基础上，本发明还公开了一种图像空间退化模拟装置，包括：
[0026] 获取模块，用于获取基础图像数据、基础图像空间分辨率、待模拟图像奈奎斯特频 率下的调制传递函数的参数值以及待模拟图像空间分辨率；
[0027] 求解模块，用于根据所述待模拟图像奈奎斯特频率下的调制传递函数的参数值， 求得有效瞬时视场角和调制传递函数的傅里叶变换函数，并对所述傅里叶变换函数的参数 值求解；
[0028] 构建模块，用于根据所述基础图像空间分辨率、所述待模拟图像空间分辨率以及 所述傅里叶变换函数求解的参数值构建点扩展函数矩阵；
[0029] 模型退化模块，用于将所述点扩展函数矩阵对应的曲线作为权重，与所述基础图 像数据进行空间卷积计算，得到图像空间退化后的模拟图像。
[0030] 优选的，所述求解模块包括：
[0031] 第一获取模块，用于获取所述待模拟图像奈奎斯特频率下的调制传递函数的参数 值对应参数向量值；
[0032] 第一求解模块，用于预设所述调制传递函数为高斯型函数，根据所述调制传递函 数的参数值对应参数向量值对每个波段依次进行所述傅里叶变换函数的参数系数求解； [0033]第二求解模块，用于根据所述调制传递函数和所述参数系数求解对应的所述有效 瞬时视场角EIF0V;
[0034]第三求解模块，用于根据所述有效瞬时视场角EIF0V求解所述傅里叶变换函数的 参数值。
[0035] 优选的，所述构建模块包括：
[0036] 第四求解模块，用于根据所述基础图像空间分辨率和所述待模拟图像空间分辨率 得到空间分辨率关系因子；
[0037] 第五求解模块，用于根据所述空间分辨率关系因子和所述有效瞬时视场角求得所 述点扩展函数矩阵的大小；
[0038]第六求解模块，用于根据所述点扩展函数矩阵的大小和所述有效瞬时视场角求得 所述点扩展函数矩阵中的自变量对应的输入矩阵；
[0039]归一化模块，用于根据所述输入矩阵和所述调制传递函数的参数值求得所述点扩 展函数矩阵，并对所述点扩展函数矩阵进行归一化，得到归一化的点扩展函数矩阵。
[0040] 优选的，所述模型退化模块包括：
[0041] 第七求解模块，用于根据所述空间分辨率关系因子和所述点扩展函数矩阵，求解 输出图像行列数；
[0042] 第八求解模块，用于根据所述有效瞬时视场角和所述输出图像行列数确定输出图 像对应输入图像的像元增量；
[0043] 退化模块，用于根据所述像元增量和所述归一化的点扩展函数矩阵进行图像空间 退化模拟，得到图像空间退化后的模拟图像。
[0044] 经由上述的技术方案可知，与现有技术相比，本发明公开了一种图像空间退化模 拟方法，该方法是依据待模拟图像奈奎斯特(Ny qu i s t)频率下的MTF值，以及基础图像和待 模拟图像的空间分辨率值，针对基础图像进行空间退化模拟，获得空间退化后的待模拟图 像。处理步骤为:首先，基于待模拟图像奈奎斯特频率下的MTF值，求得其有效瞬时视场角， 进而求得MTF函数的傅里叶变换函数（即PSF函数）的参数;其次，根据基础图像的空间分辨 率、待模拟图像的空间分辨率，计到空间比值，利用PSF函数参数，计算得到PSF二维矩阵;最 后，将PSF二维矩阵作为权重，与输入图像进行卷积计算，进行逐像元的空间退化计算，最终 获得考虑MTF效应空间退化后的待模拟图像，上述方法通过考虑MTF或PSF对于图像空间实 际分辨能力的影响，解决了现有技术中由于采用像元求平均值来实现空间退化模拟的方法 得到的图像分辨率无法满足新型遥感器的图像数据模拟精度或无法为模型算法研究提供 需要的图像数据的问题。
【附图说明】
[0045] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本 发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据 提供的附图获得其他的附图。
[0046] 图1为本发明实施例公开的一种图像空间退化模拟方法流程图；
[0047]图2为本发明实施例公开的图1中步骤S102的方法流程图；
[0048]图3为本发明实施例公开的图1中步骤S103的方法流程图；
[0049] 图4为本发明实施例公开的图1中步骤S104的方法流程图；
[0050] 图5为本发明实施例公开的一种图像空间退化模拟装置的结构示意图；
[0051] 图6为本发明实施例公开的图5中的求解模块的具体结构示意图；
[0052] 图7为本发明实施例公开的图5中的构建模块的具体结构示意图；
[0053] 图8为本发明实施例公开的图5中的模型退化模块的具体结构示意图。
【具体实施方式】
[0054]下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0055] 传感器的PSF(po int spread function，点扩展函数）或MTF(modulation transfer function，调制传递函数)是传感器成像性能的重要指标，其中PSF是信号扩散与 肖IJ弱在空间域的描述，而MTF是PSF在频率域的描述。MTF或PSF是用一个函数去表达，反映的 是系统在对目标成像过程中的信号扩散与削弱程度，是决定图像质量及其应用效能的重要 参量。因此，在对图像进行空间退化模拟，应考虑MTF或PSF对于图像空间实际分辨能力的影 响。
[0056] 具体的，本发明采用以下实施例手段进行实现：
[0057]获取图像模拟对应的基础图像imgin(nl，ns，nb)，其中nl表示图像行数，ns表示图 像列数，nb为图像波段数；图像模拟对应的基础图像的空间分辨率SRin;待模拟图像，即空 间退化后图像的空间分辨率SRout;待模拟图像的奈奎斯特频率下的MTF参数值mtf_n(nb)。
[0058] 预先设定MTF函数为高斯型，即函数/(x：> = eW ^艮据给定的奈奎斯特频率下的 MTF参数值，求得系数c，第k波段对应的c (k)为：
[0059]
[0060]通过计算得到有效瞬时视场角EIF0V，作为图像实际的分辨率：
[0061]
[0062] 其中，XQ.5(k)为上述高斯函数值等于0.5时候对应的X值即：
[0063]
[0064] 根据EIF0V，求得图像的二维PSF函数，该函数为高斯概率密度函数：
[0065]
[0066]该函数公式中对应到波段k的〇 (k)为：
[0067]
[0068] 具体的，根据基础图像、待模拟图像的空间分辨率，以及每个波段的有效瞬时视场 角，构建二维PSF函数矩阵。对于波段k有：
[0069]求得空间分辨率比例关系：
[0070]
[0071] 二维PSF函数矩阵为方阵，其中一维单元数nk(k)为：
[0072] temp(k)=ceil(2*EIF0V(k)*scale+l)，
[0073] 其中，ceil表示向上取整。如果temp(k)为奇数，贝ljnk(k) = temp(k);如果temp(k) 为偶数，贝1Jnk(k)=temp(k)+1。
[0074] 求得函数矩阵中自变量x(k)，y(k)对应的输入矩阵X(k)，Y(k)，为nk(k)*nk(k)的 方阵：
[0075] xl(k)=-(nk(k)-l)/2:l:(nk(k)-l)/2
[0076] X(l, :)=X(2，：）= ··· =X(nk(k)，：）=xl(k)*2*EIF0V(k)/nk(k)
[0077] Y(k)=X(k)
[0078] 根据自变量x(k)，y(k)对应的输入矩阵X(k)，Y(k)，求得二维PSF概率密度函数矩 阵为：
[0079]
[0080]其中，Z(k)为 nk(k)*nk(k)的方阵。
[0081 ] 对PSF概率密度函数矩阵Z(k)求和：
[0082] f actor (k) =sum(Z(k))
[0083] 对PSF概率密度函数矩阵Z (k)进行归一化：
[0084] Zl(k)=Z(k)/factor(k)
[0085]优选的，确认输出图像imgout的大小，列数nsout，行数nlout为：
[0086] nsout = floor (ns/scale)
[0087] nlout = floor(nl/scale)
[0088] 其中，floor表示向下取整的意思，如floor(10.5)为10。
[0089] 最后，针对输出图像的每个点基于PSF概率密度函数矩阵Z(k)进行卷积求解，实现 基于MTF参数的图像空间退化模拟。
[0090] 主要的，针对波段k有：
[0091] 确定像元增量：
[0092] step(k) =2*EIF0V(k)*scale/nk(k)
[0093] Zl(k)作为核函数对图像进行卷积计算，对于第k波段，针对输出图像的第m，n像 元，通过以下过程进行计算得到空间退化后的结果：
[0094] 求得该像元对应基础图像上的空间位置：
[0095] xx= (scale+1 )/2+(m_l )*scale;
[0096] yy = (scale+1 )/2+( n-1 )*scale ；
[0097] 以xx，yy为中心位置，求得二维PSF函数每个点在基础图像上的对应的空间位置， 对于该矩阵中第ki，kj点的像元位置(nx(ki)，ny(kj))有：
[0098] varl(k)=(nk(k)-l)/2
[0099] var(k) = (nk(k)-l)/2
[0100] nx(ki) =xx+(ki-1-varl(k) )*step(k)，ki = l :nk(k)
[0101 ] ny(kj) =yy+(kj-1-varl(k) )*step(k) ,kj = l ：nk(k)
[0102]求得矩阵中像元位置(nx (ki)，ny (kj))在基础图像对应的像元值，对于第k波段 imgin(nx(ki)，ny(kj)，k)为：
[0103] xl =f loor(nx(ki))；
[0104] x2 = ceil(nx(ki))；
[0105] yl = f loor(ny(kj))；
[0106] y2 = ceil(ny(kj))；
[0107] dx = x2-nx(ki)；
[0108] dy = y2-ny(kj)；
[0109] tmp = dx*imgin(xl，yl，k) + ( l_dx)*imgin(x2，yl，k);
[0110] 具体的，imgin(nx(ki)，ny(kj)，k) =dy*tmp+( l_dy)*(dx*imgin(xl，y2，k) + (l_ dx)*imgin(x2，y2，k))〇
[0111] 以归一化后的二维PSF函数矩阵Z1 (nk(k)，nk(k))为核函数，通过卷积运算求解空 间退化后的imgout (m，η，k)的像元值：
[0112]
[0113]请参阅附图1，图1为本发明公开的一种图像空间退化模拟方法流程图。如图1所 示，本发明公开了 一种图像空间退化模拟方法，该方法具体包括如下步骤：
[0114] S101、获取基础图像数据、基础图像空间分辨率、待模拟图像奈奎斯特频率下的调 制传递函数的参数值以及待模拟图像空间分辨率。
[0??5] 其中，图像模拟对应的基础图像imgin，三维矩阵，nl*ns*nb，其中nl表示图像行 数，ns表示图像列数，nb为图像波段数；图像模拟对应的基础图像的空间分辨率变量SRin; 待模拟图像，即空间退化后图像的空间分辨率变量SRout;待模拟图像的奈奎斯特频率下的 MTF参数值向量mtf _n，长度为波段数nb。
[0116] S102、根据待模拟图像奈奎斯特频率下的调制传递函数的参数值，求得有效瞬时 视场角和调制传递函数的傅里叶变换函数，并对傅里叶变换函数的参数值求解。
[0117] 基于奈奎斯特频率下的MTF值，求得有效瞬时视场角，进而求得MTF函数的傅里叶 变换函数，即PSF函数。
[0118] S103、根据基础图像空间分辨率待模拟图像空间分辨率以及傅里叶变换函数求解 的参数值构建点扩展函数矩阵。
[0119] S104、将点扩展函数矩阵对应的曲线作为权重，与基础图像数据进行空间卷积计 算，得到图像空间退化后的模拟图像。
[0120] 根据基础图像的空间分辨率、待模拟图像的空间分辨率比值，将PSF函数曲线作为 权重，与输入图像进行卷积计算，进行逐像元的空间退化计算，最终获得空间退化后的模拟 图像。
[0121] 具体的，请参阅附图2,图2为图1中步骤S102的方法流程图。具体的，步骤S102为如 下步骤：
[0122] S201、获取待模拟图像奈奎斯特频率下的调制传递函数的参数值对应参数向量 值。
[0123] S202、预设调制传递函数为高斯型函数，根据调制传递函数的参数值对应参数向 量值对每个波段依次进行傅里叶变换函数的参数系数求解。
[0124] S203、根据调制传递函数和参数系数求解对应的有效瞬时视场角EIFOV。
[0125] S204、根据有效瞬时视场角EIF0V求解傅里叶变换函数的参数值。
[0126] 上述步骤可以描述为以下内容：
[0127] 输入待模拟图像的奈奎斯特频率下的MTF参数值向量mtf_n(nb)。
[0128] 设定MTF函数为高斯型，即函数/(x) = e^/e2，根据给定的奈奎斯特频率下的MTF参 数值，对每个波段依次进行PSF函数参数系数求解，第k波段对应的c (k)为：
[0129]
[0130]计算得到有效瞬时视场角EIF0V，作为图像实际的分辨率：
[0131]
[0132] 其中，XQ.5(k)为上述高斯函数值等于0.5时候对应的X值即：
[0133]
[0134] 根据EIF0V，求得图像的二维PSF函数，该函数为高斯概率密度函数2 = (?_$^该 .， 函数公式中对应到波段沾勺〇 (k)为：
[0135]
[0136] 具体的，请参阅附图3,图3为图1中步骤S103的方法流程图。具体的，步骤S103为如 下步骤：
[0137] S301、根据基础图像空间分辨率和待模拟图像空间分辨率得到空间分辨率关系因 子。
[0138] S302、根据空间分辨率关系因子和有效瞬时视场角求得点扩展函数矩阵的大小。
[0139] S303、根据点扩展函数矩阵的大小和有效瞬时视场角求得点扩展函数矩阵中的自 变量对应的输入矩阵。
[0140] S304、根据输入矩阵和调制传递函数的参数值求得点扩展函数矩阵，并对点扩展 函数矩阵进行归一化，得到归一化的点扩展函数矩阵。
[0141] 上述步骤可以描述为以下内容：
[0142] 输入基础图像的空间分辨率变量SRin;输入待模拟图像，即空间退化后图像的空 间分辨率变量SRout;输入有效瞬时视场角EIFOV(nb);输入二维PSF函数参数〇(nb)。
[0143] 求得空间分辨率关系因子：
[0144]
[0145] 求得二维PSF函数矩阵大小，该矩阵方阵，对于第k波段，其某一维单元数nk(k)为：
[0146] temp(k) =ceil(2*EIF0V(k)*scale+l)，
[0147] 其中06；11表示向上取整;如果七611^1(1〇为奇数，贝1|111^(1〇 =七611^1(1〇;如果七611^1(1〇为 偶数，贝1Jnk(k) =temp(k)+l。
[0148] 求得二维PSF函数矩阵中自变量x，y对应的输入矩阵X、Y，为nk(k)*nk(k)的方阵， 对于第k波段X(k)、Y(k)为：
[0149] xl(k)=-(nk(k)-l)/2:l:(nk(k)-l)/2
[0150] X(l，：）=X(2，：）= ...=X(nk(k)，：）=xl (k)*2*EIF0V(k)/nk(k)
[0151] Y(k)=X(k)
[0152] 求得二维PSF函数矩阵，对于第k波段，根据自变量x(k)，y(k)对应的输入矩阵X (k)，Y(k)，有：
[0153]
其中，Z(k)为 nk(k)*nk(k)的方阵。
[0154] 对PSF函数矩阵进行归一化，具体为：
[0155] 对PSF函数矩阵求和，对于第k波段有：
[0156] f actor (k) =sum(Z(k))
[0157] 对PSF概率密度函数矩阵进行归一化，对于第k波段有：
[0158] Zl(k)=Z(k)/factor(k)
[0159] 最后得到归一化的点扩展函数矩阵。
[0160]具体的，请参阅附图4,图4为图1中步骤S104的方法流程图。具体的，步骤S104为如 下步骤：
[0161] S401、根据空间分辨率关系因子和点扩展函数矩阵，求解输出图像行列数。
[0162] S402、根据有效瞬时视场角和输出图像行列数确定输出图像对应输入图像的像元 增量。
[0163] S403、根据像元增量和归一化的点扩展函数矩阵进行图像空间退化模拟，得到图 像空间退化后的模拟图像。
[0164] 上述步骤可以描述为以下内容：
[0165] 根据空间分辨率关系因子scale，输入基础图像矩阵imgin(nl，ns，nb)、每个波段 的有效瞬时视场角EIFOV(nb)、每个波段的归一化后的二维PSF函数矩阵，对于第k波段矩阵 为Zl(nk(k)，nk(k))〇
[οι 66]求解输出图像行列数，imgojg澈f&胃織亍数nlout分别为：
[0167] nlout = f loor(nl/scale)
[0168] 其中，floor表示向下取整的意思，如floor (10.5)为10。
[0169] 确定输出图像单个像元对应输入图像的像元增量：
[0170] step(k) =2*EIF0V(k)*scale/nk(k)
[0171] 逐像元的求得待模拟图像像元值，进行基于MTF参数的图像空间退化模拟。对于第 k波段，针对输出图像的第m，n像元，通过以下过程进行计算得到空间退化后的结果：
[0172] 求得该像元对应基础图像上的空间位置：
[0173] xx= (scale+1 )/2+(m_l )*scale;
[0174] yy = (scale+1 )/2+(n_l )*scale ;
[0175] 以xx，yy为中心位置，求得二维PSF函数每个点在基础图像上的对应的空间位置， 对于该矩阵中第ki，kj点的像元位置(nx(ki)，ny(kj))有：
[0176] varl (k) = (nk(k)-1)/2
[0177] nx(ki) =xx+(ki-l_varl(k) )*step(k)，ki = l :nk(k)
[0178] ny(kj) = yy+(kj-l_varl(k))*step(k)，kj = 1 :nk(k)
[0179] 求得矩阵中像元位置(nx (ki)，ny (kj))在基础图像对应的像元值，对于第k波段 imgin(nx(ki)，ny(kj)，k)为：
[0180] xl =f loor(nx(ki))；
[0181] x2 = ceil(nx(ki))；
[0182] yl = f loor(ny(kj))；
[0183] y2 = ceil(ny(kj))；
[0184] dx = x2_nx(ki);
[0185] dy = y2-ny(kj)；
[0186] tmp = dx*imgin(xl，yl，k) + ( l_dx)*imgin(x2，yl，k);
[0187] imgin(nx(ki),
[0188] ny(kj)，k)=dy*tmp+(l_dy)*(dx*imgin(xl，y2，k) + (l_dx)*imgin(x2，y2，k));
[0189] 以归一化后的二维PSF函数矩阵Z1 (nk(k)，nk(k))为核函数，通过卷积运算求解空 间退化后的imgout (m，η，k)的像元值：
[0190]
[0191] 本发明公开了一种图像空间退化模拟方法，该方法是依据待模拟图像奈奎斯特 (Nyqui st)频率下的MTF值，以及基础图像和待模拟图像的空间分辨率值，针对基础图像进 行空间退化模拟，获得空间退化后的待模拟图像。本发明通过考虑MTF或PSF对于图像空间 实际分辨能力的影响，解决了现有技术中由于采用像元求平均值来实现空间退化模拟的方 法得到的图像分辨率无法满足新型遥感器的图像数据模拟精度或无法为模型算法研究提 供需要的图像数据的问题。
[0192] 在上述公开的方法的基础上，请参阅附图5,图5为本发明公开的一种图像空间退 化模拟装置的结构示意图。本发明还公开了一种图像空间退化模拟装置，具体的该装置包 括：
[0193] 获取模块501，用于获取基础图像数据、基础图像空间分辨率、待模拟图像奈奎斯 特频率下的调制传递函数的参数值以及待模拟图像空间分辨率。
[0194] 求解模块502,用于根据待模拟图像奈奎斯特频率下的调制传递函数的参数值，求 得有效瞬时视场角和调制传递函数的傅里叶变换函数，并对傅里叶变换函数的参数值求 解。
[0195] 构建模块503,用于根据基础图像空间分辨率、待模拟图像空间分辨率以及傅里叶 变换函数求解的参数值构建点扩展函数矩阵。
[0196] 模型退化模块504,用于将点扩展函数矩阵对应的曲线作为权重，与基础图像数据 进行空间卷积计算，得到图像空间退化后的模拟图像。
[0197] 本实施例的装置，可以用于执行图1所示方法实施例的技术方案，其实现原理和技 术效果类似，此处不再赘述。
[0198] 具体的，请参阅附图6,图6为图5中的求解模块的具体结构示意图，所述求解模块 502包括：
[0199] 第一获取模块601，用于获取待模拟图像奈奎斯特频率下的调制传递函数的参数 值对应参数向量值。
[0200] 第一求解模块602,用于预设调制传递函数为高斯型函数，根据调制传递函数的参 数值对应参数向量值对每个波段依次进行傅里叶变换函数的参数系数求解。
[0201] 第二求解模块603,用于根据调制传递函数和参数系数求解对应的有效瞬时视场 角EIF0V。
[0202]第三求解模块604,用于根据有效瞬时视场角EIF0V求解傅里叶变换函数的参数 值。
[0203] 本实施例的装置，可以用于执行图2所示方法实施例的技术方案，其实现原理和技 术效果类似，此处不再赘述。
[0204] 具体的，请参阅附图7,图7为图5中的构建模块的具体结构示意图，所述构建模块 503包括：
[0205]第四求解模块701，用于根据基础图像空间分辨率和待模拟图像空间分辨率得到 空间分辨率关系因子。
[0206]第五求解模块702,用于根据空间分辨率关系因子和有效瞬时视场角求得点扩展 函数矩阵的大小。
[0207]第六求解模块703,用于根据点扩展函数矩阵的大小和有效瞬时视场角求得点扩 展函数矩阵中的自变量对应的输入矩阵。
[0208]归一化模块704,用于根据输入矩阵和调制传递函数的参数值求得点扩展函数矩 阵，并对点扩展函数矩阵进行归一化，得到归一化的点扩展函数矩阵。
[0209] 本实施例的装置，可以用于执行图3所示方法实施例的技术方案，其实现原理和技 术效果类似，此处不再赘述。
[0210] 具体的，请参阅附图8,图8为图5中的模型退化模块的具体结构示意图，所述模型 退化模块504包括：
[0211] 第七求解模块801，用于根据空间分辨率关系因子和点扩展函数矩阵，求解输出图 像行列数。
[0212] 第八求解模块802,用于根据有效瞬时视场角和输出图像行列数确定输出图像对 应输入图像的像元增量。
[0213] 退化模块803,用于根据像元增量和归一化的点扩展函数矩阵进行图像空间退化 模拟，得到图像空间退化后的模拟图像。
[0214] 本实施例的装置，可以用于执行图4所示方法实施例的技术方案，其实现原理和技 术效果类似，此处不再赘述。
[0215] 综上所述，本发明公开了一种图像空间退化模拟方法，该方法是依据待模拟图像 奈奎斯特(Nyqui st)频率下的MTF值，以及基础图像和待模拟图像的空间分辨率值，针对基 础图像进行空间退化模拟，获得空间退化后的待模拟图像。处理步骤为:首先，基于待模拟 图像奈奎斯特频率下的MTF值，求得其有效瞬时视场角，进而求得MTF函数的傅里叶变换函 数(即PSF函数)的参数;其次，根据基础图像的空间分辨率、待模拟图像的空间分辨率，计到 空间比值，利用PSF函数参数，计算得到PSF二维矩阵;最后，将PSF二维矩阵作为权重，与输 入图像进行卷积计算，进行逐像元的空间退化计算，最终获得考虑MTF效应空间退化后的待 模拟图像，上述方法通过考虑MTF或PSF对于图像空间实际分辨能力的影响，解决了现有技 术中由于采用像元求平均值来实现空间退化模拟的方法得到的图像分辨率无法满足新型 遥感器的图像数据模拟精度或无法为模型算法研究提供需要的图像数据的问题。另外，本 发明还公开了一种图像空间退化模拟装置。
[0216] 需要说明的是，本说明书中的各个实施例均采用递进的方式描述，每个实施例重 点说明的都是与其他实施例的不同之处，各个实施例之间相同相似的部分互相参见即可。
[0217] 对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。 对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的 一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明 将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一 致的最宽的范围。
【主权项】
1. 一种图像空间退化模拟方法，其特征在于，包括： 获取基础图像数据、基础图像空间分辨率、待模拟图像奈奎斯特频率下的调制传递函 数的参数值以及待模拟图像空间分辨率； 根据所述待模拟图像奈奎斯特频率下的调制传递函数的参数值，求得有效瞬时视场角 和调制传递函数的傅里叶变换函数，并对所述傅里叶变换函数的参数值求解； 根据所述基础图像空间分辨率、所述待模拟图像空间分辨率以及所述傅里叶变换函数 求解的参数值构建点扩展函数矩阵； 将所述点扩展函数矩阵对应的曲线作为权重，与所述基础图像数据进行空间卷积计 算，得到图像空间退化后的模拟图像。2. 根据权利要求1所述的方法，其特征在于，所述根据所述待模拟图像奈奎斯特频率下 的调制传递函数的参数值，求得有效瞬时视场角和调制传递函数的傅里叶变换函数，并对 所述傅里叶变换函数求解的参数值具体为： 获取所述待模拟图像奈奎斯特频率下的调制传递函数的参数值对应参数向量值； 预设所述调制传递函数为高斯型函数，根据所述调制传递函数的参数值对应参数向量 值对每个波段依次进行所述傅里叶变换函数的参数系数求解； 根据所述调制传递函数和所述参数系数求解对应的所述有效瞬时视场角EIFOV; 根据所述有效瞬时视场角EIFOV求解所述傅里叶变换函数的参数值。3. 根据权利要求1所述的方法，其特征在于，所述根据所述基础图像空间分辨率、所述 待模拟图像空间分辨率以及所述傅里叶变换函数求解的参数值构建点扩展函数矩阵具体 为： 根据所述基础图像空间分辨率和所述待模拟图像空间分辨率得到空间分辨率关系因 子； 根据所述空间分辨率关系因子和所述有效瞬时视场角求得所述点扩展函数矩阵的大 小； 根据所述点扩展函数矩阵的大小和所述有效瞬时视场角求得所述点扩展函数矩阵中 的自变量对应的输入矩阵； 根据所述输入矩阵和所述调制传递函数的参数值求得所述点扩展函数矩阵，并对所述 点扩展函数矩阵进行归一化，得到归一化的点扩展函数矩阵。4. 根据权利要求3所述方法，其特征在于，所述将所述点扩展函数矩阵对应的曲线作为 权重，与所述基础图像数据进行空间卷积计算，得到图像空间退化后的模拟图像具体为： 根据所述空间分辨率关系因子和所述点扩展函数矩阵，求解输出图像行列数； 根据所述有效瞬时视场角和所述输出图像行列数确定输出图像对应输入图像的像元 增量； 根据所述像元增量和所述归一化的点扩展函数矩阵进行图像空间退化模拟，得到图像 空间退化后的模拟图像。5. -种图像空间退化模拟装置，其特征在于，包括： 获取模块，用于获取基础图像数据、基础图像空间分辨率、待模拟图像奈奎斯特频率下 的调制传递函数的参数值以及待模拟图像空间分辨率； 求解模块，用于根据所述待模拟图像奈奎斯特频率下的调制传递函数的参数值，求得 有效瞬时视场角和调制传递函数的傅里叶变换函数，并对所述傅里叶变换函数的参数值求 解； 构建模块，用于根据所述基础图像空间分辨率、所述待模拟图像空间分辨率以及所述 傅里叶变换函数求解的参数值构建点扩展函数矩阵； 模型退化模块，用于将所述点扩展函数矩阵对应的曲线作为权重，与所述基础图像数 据进行空间卷积计算，得到图像空间退化后的模拟图像。6. 根据权利要求5所述的装置，其特征在于，所述求解模块包括： 第一获取模块，用于获取所述待模拟图像奈奎斯特频率下的调制传递函数的参数值对 应参数向量值； 第一求解模块，用于预设所述调制传递函数为高斯型函数，根据所述调制传递函数的 参数值对应参数向量值对每个波段依次进行所述傅里叶变换函数的参数系数求解； 第二求解模块，用于根据所述调制传递函数和所述参数系数求解对应的所述有效瞬时 视场角EIFOV; 第三求解模块，用于根据所述有效瞬时视场角EIFOV求解所述傅里叶变换函数的参数 值。7. 根据权利要求5所述的装置，其特征在于，所述构建模块包括： 第四求解模块，用于根据所述基础图像空间分辨率和所述待模拟图像空间分辨率得到 空间分辨率关系因子； 第五求解模块，用于根据所述空间分辨率关系因子和所述有效瞬时视场角求得所述点 扩展函数矩阵的大小； 第六求解模块，用于根据所述点扩展函数矩阵的大小和所述有效瞬时视场角求得所述 点扩展函数矩阵中的自变量对应的输入矩阵； 归一化模块，用于根据所述输入矩阵和所述调制传递函数的参数值求得所述点扩展函 数矩阵，并对所述点扩展函数矩阵进行归一化，得到归一化的点扩展函数矩阵。8. 根据权利要求7所述装置，其特征在于，所述模型退化模块包括： 第七求解模块，用于根据所述空间分辨率关系因子和所述点扩展函数矩阵，求解输出 图像行列数； 第八求解模块，用于根据所述有效瞬时视场角和所述输出图像行列数确定输出图像对 应输入图像的像元增量； 退化模块，用于根据所述像元增量和所述归一化的点扩展函数矩阵进行图像空间退化 模拟，得到图像空间退化后的模拟图像。
【文档编号】G06T5/50GK105976317SQ201610274396
【公开日】2016年9月28日
【申请日】2016年4月28日
【发明人】张文娟, 张兵, 高连如, 陈正超, 高建威
【申请人】中国科学院遥感与数字地球研究所