一种高分辨率光学遥感卫星相机最佳焦面确定方法
【技术领域】
[0001] 本发明涉及一种在轨高分辨率光学遥感卫星相机最佳焦面确定方法，属于航天遥 感图像处理技术领域。
【背景技术】
[0002] 高分辨率卫星显著的特点是长焦距，分辨率优于1米的遥感卫星等效焦距一般接 近10米，例如美国研制的IK0N0S、QuickBird卫星、中国研制的高分二号卫星等效焦距分别 达到10米、8.8米及7.8米。卫星完成在轨测试交付用户时，卫星焦面处于最佳位置，获取的 图像质量达到最优水平。由于长焦相机焦深浅，卫星在轨运行期间，更容易受到温度、大气 压力以及应力释放等因素的影响，引起相机成像面偏离焦面，即造成离焦。高分辨率卫星相 机焦距越长，焦深越浅，越容易离焦。离焦图像数字值(Digital Number，DN)趋于单一灰度， 纹理模糊，信息含量减少。为了避免高分辨率卫星相机离焦工作，卫星在轨运行期间，需要 定期监测图像质量，通过分析图像清晰度判读相机是否离焦。一旦离焦即启动相机调焦程 序，调焦的目的是通过调整焦平面的位置，保证相机工作在最佳焦面上，获取清晰度最高的 图像。
[0003] 可以在地面采用图像判读、分析的技术分析卫星在轨图像清晰度，以判断相机是 否离焦并进一步确定最佳焦面。目前分析卫星在轨图像清晰度的算法没有针对性，一般采 用图像标准差、峰值信噪比等算法，上述算法均针对整幅图像，当包含2-3种反差大的均匀 地物时，纹理信息缺乏，图像标准差或峰值信噪比很大，不能精确反映地物的真实纹理特 征，降低了参考图像和待检测图像对比结果的可信度，另外，现有技术采用相对独立的指标 对最佳焦面进行确定，没有综合利用各指标综合评价，降低了最佳焦面确定的准确性。
【发明内容】
[0004] 本发明解决的技术问题是：克服现有技术的不足，提供一种高分辨率光学遥感卫 星相机最佳焦面确定方法，保证了待检测图像和参考图像对比结果的可信度，提高了最佳 焦面确定的准确性。
[0005] 本发明的技术方案是:一种高分辨率光学遥感卫星相机最佳焦面确定方法，包括 如下步骤：
[0006] (1)在光学遥感卫星相机处于最佳焦面时，从某区域地面站接收的图像中选择纹 理清晰的经过正射校正的正视图像作为参考图像；
[0007] (2)在光学遥感卫星相机在轨运行期间过程中，从该区域地面站接收的图像中选 择与参考图像同样位置的图像作为待检测图像，所述待检测图像是经过正射校正的正视图 像或侧摆角度小于5度的图像，且能够与参考图像几何配准；
[0008] (3)判断待检测图像和参考图像亮度是否一致，如果一致，进入步骤(4);如果不一 致，依据参考图像对待检测图像进行亮度归一化处理，使待检测图像平均亮度与参考图像 一致，进入步骤(4);
[0009] (4)计算参考图像的分块标准差及梯度；
[0010] (5)计算待检测图像的分块标准差及梯度；
[0011] (6)确定待检测图像和参考图像的分块标准差比值〇scale以及梯度比值g scale，利用
[0012] (7)当r大于0.9，即认为光学遥感卫星相机处于最佳焦面；当0.75 < r < 0.9时，需 要继续选择待检测图像进行确定；当r小于0.75时，光学遥感卫星相机处于离焦状态，启动 调焦工作，进入步骤(8);
[0013] (8)在调焦经历的每个焦面下，从该区域地面站接收的图像中选择与参考图像同 样位置的图像作为待检测图像，执行步骤(4) 一 (6)，计算每个焦面下待检测图像和参考图 像的图像指标比例因子，直到图像指标比例因子大于0.9时，认为该图像指标比例因子对应 的焦面为最佳焦面。
[0014] 所述步骤(3)中依据参考图像对待检测图像进行亮度归一化处理的方法为：
[0015] (2.1)计算参考图像的亮度均值/:
[0017] f (i，j)为参考图像坐标(i，j)处的图像壳度，M*N为参考图像整幅图像尺寸；
[0018] (2.2)计算待检测图像的亮度均值
[0020] g (i，j)为待检测图像坐标(i，j)处的图像亮度；
[0021] (2.3)利用参考图像的亮度均值/和待检测图像的亮度均值I，对待检测图像中的 每个点进行亮度归一化：
[0022] g(U) = ：g(i,jj*f/g
[0023] 8'（1，」)为待检测图像坐标(1，」)处归一化后的图像亮度。
[0024] 所述步骤(4)或步骤(5)中计算分块标准差的方法为：
[0025] (3.1)利用以下公式计算分块后第p块图像的亮度均值μΡ
[0027] m*n为分块后每一子块的尺寸，fp(i，j)为第ρ块图像中坐标(i，j)处的图像亮度；
[0028] (3.2)利用以下公式计算第ρ块图像的标准差％
[0030] (3.3)计算整幅图像各子块标准差的均值孓,J即为整幅图像的分块标准差
[0032] k为子块总数，即k= [M/m]*[N/n]，其中，[M/m]为M/m四舍五入后的整数值，[N/n] 为N/n四舍五入后的整数值，M*N为整幅图像尺寸。
[0033]所述步骤(4)或步骤(5)中利用如下公式计算图像梯度g:
[0035] 其中，f(i，j)为坐标(i，j)处的图像亮度，f(i，j-l)为坐标处的图像亮度， f (i-1，j)为坐标(i-1，j)处的图像亮度。
[0036] 本发明与现有技术相比的优点如下：
[0037] (1)本发明方法采用分块标准差及图像梯度精确获取图像的纹理细节，综合分块 标准差与图像梯度两个指标判断最佳焦面，相较于现有技术中针对整幅图像的图像标准 差、峰值信噪比等算法以及各指标独立评价的方法，本发明方法科学可靠，针对性强，保证 了待检测图像和参考图像对比结果的可信度，提高了最佳焦面确定的准确性。
[0038] (2)利用本发明方法既能够在卫星在轨运行过程中确定相机是否处于最佳焦面或 离焦，又可以在离焦后调焦过程中确定最佳焦面，方法简单可靠，能够有效保证光学遥感卫 星相机工作在最佳焦面上，获取清晰度高的图像。
[0039] (3)本发明在卫星运行过程中确定相机是否处于最佳焦面或离焦，以及在离焦后 调焦过程中确定最佳焦面时，利用待检测图像和参考图像的亮度均值，对待检测图像中的 每个点进行亮度归一化，从而抑制甚至消除不同季节太阳光照的影响，使待检测图像与参 考图像壳度一致，提尚最佳焦面研判的精准度。
【附图说明】
[0040] 图1为本发明方法流程图；
[0041] 图2为高分一号卫星获取的拉萨图像，其中（a)为参考图像，（b)为待检测图像，（c) 为归一化后的待检测图像；
[0042]图3为高分二号卫星获取的敦煌图像，其中(a)为参考图像，（b)为待检测图像。
【具体实施方式】
[0043]下面结合附图和具体实施例对本发明作进一步详细的描述：
[0044] 本发明针对现有技术判断相机是否离焦以及如何确定最佳焦面方面存在的不足， 提出了一种最佳焦面确定方法，实现两个目的：1)卫星在轨运行过程中图像分析，确定相机 是否处于最佳焦面或离焦;2)卫星相机调焦过程中图像分析、确定最佳焦面。该方法基本思 路为:首先采集受天气影响较小、地面景物相对丰富且地面站可以直接接收到的典型区域 (例如敦煌、喀什、银川）图像，起到参考图像的目的。卫星在轨运行期间，一旦获取该区域的 有效图像，即可作为待检测图像，针对纹理、细节等反映离焦程度的参数，与参考图像对比、 分析，判读是否离焦及估计离焦的大致程度;一旦离焦启动调焦程序，对目标区成像，计算 目标区图像的关键指标，与参考图像的指标对比，直到数值接近即可认定为最佳焦面。
[0045] 如图1所示，本发明提出一种高分辨率光学遥感卫星相机最佳焦面确定方法，包括 以下步骤：
[0046] (1)在光学遥感卫星相机处于最佳焦面时（如在轨测试结束时），从某区域地面站 接收的图像中选择纹理清晰的经过正射校正的正视图像作为参考图像。
[0047]选取的参考图像应该能够反映卫星相机正常工作状态，具有典型性，选取时应满 足以下条件：
[0048] 1)该区域受天气影响较小，晴天多阴天少，有利于获取该区域有效图像；
[0049] 2)该区域地面站可以直接接收，容易获取该区域图像；
[0050] 3)在不同经度尽量均匀分布；
[0051] 4)图像不少于20幅，有条件可在50幅以上，满足时序化监测要求；
[0052] 5)景物目标相对固定，变化较大需要重新采集；
[0053] 6)图像尽量反映不同时相；
[0054] 7)反映相机波谱特性，全色+多光谱同步采集；
[0055] 8)地面景物相对丰富，一般为城市，例如敦煌、喀什、银川及我国西部近邻区域；
[0056] 9)正视成像(无侧摆），图像质量得到保证；
[0057] 10)获取的图像成像条件大气良好，图像纹理清晰；
[0058] 11 )Geotif f格式，图像具有几何地理信息；
[0059] 12)图像经过正射校正，有利于几何匹配；
[0060] 13)图像尺寸1024像元* 1024像元-2048像元*2048像元。
[0061] (2)在光学遥感卫星相机在轨运行期间过程中，从该区域地面站接收的图像中选 择与参考图像同样位置的图像作为待检测图像。
[0062]待检测图像所在区域及成像满足以下条件：
[0063] 1)成像区域包含参考图像；
[0064] 2)正视成像或放宽条件至侧摆角度小于5度；
[0065] 3)Geotiff格式，能够与参考图像几何配准；
[0066] 4)图像经过正射校正；
[0067] 5)图像尺寸与参考图像一致。
[0068]待检测图像的幅数与参考图像一致，然后从这多幅图像中，选择一组参考图像和 相对应的待检测图像纹理都清晰的图像对进行最佳焦面确定。
[0069] (3)判断待检测图像和参考图像亮度是否一致，如果一致，进入步骤(4);如果不一 致，依据参考图像对待检测图像进行亮度归一化处理，使待检测图像平均亮度与参考图像 一致，进入步骤(4);
[0070] 亮度归一化处理的方法为：
[0071] (3.1)计算参考图像的亮度均值
[0073