一种可减弱大气湍流影响的成像方法及装置

摘要

本发明公开了一种可减弱大气湍流影响的成像方法及装置，方法步骤包括：采集成像目标的四维光场分布及退化图像，解算全视场各视角方向相位波前畸变，估计求解成像系统的系统成像点扩散函数，在频域对退化图像进行解卷积运算，求解得到成像目标理想成像的傅里叶变换，然后将成像目标理想成像的傅里叶变换进行傅里叶逆变换，得到成像目标接近于衍射极限的理想成像；装置包括主透镜、分光器件、成像信息采集高帧频感光CCD、微透镜阵列和光场信息采集低噪声高帧频感光CCD。本发明具有探测视场大、成本低、动态范围大、图像分辨率高、适用范围广的优点。

说明书

技术领域

本发明涉及高空目标成像的自适应光学和图像处理领域，具体涉及一种可减弱大气湍流 影响的成像方法及装置。

背景技术

大气是一种不稳定的随机介质，在成像系统中，由于“大气湍流”这种动态干扰的影响， 光学望远镜的实际分辨率远远达不到理论上所预期的光学衍射极限。一个口径达数米的大型 无补偿光学望远镜，在地面通过大气层对天空及天文目标进行观察时，由于大气湍流的影响， 实际成像分辨率不会超过口径仅为0.1～0.2m的小望远镜，并且观察到的图像存在一定的模 糊和抖动。

自适应光学技术具有克服光学系统各种误差因素影响的能力，使其适用于望远镜和光学 遥感器的高分辨率成像，自适应光学技术具有克服光学系统各种误差因素影响、保持系统始 终工作在最佳状态的能力，使其适用于望远镜和光学遥感器的高分辨率成像，在激光大气传 输、激光加工、眼科科学、航空航天及通信领域展现了一定的应用前景。尽管自适应光学技 术能够实时校正光束的波前畸变，但是由于自身设计、计算机处理能力、闭环伺服带宽、波 前观测数据误差以及噪声等因素影响，自适应光学对大气湍流的补偿仅仅是部分的、不充分 的，目标的高频信息仍然受到严重的抑制和衰减。另外，自适应光学系统还存在光能利用率 低(哈特曼分光，暗弱目标使用受限)、结构复杂(特别是针对扩展目标成像)、成本昂贵等 特点，并且其有限的工作带宽，使其不适用于气动光学效应的校正。图像后处理技术是利用 望远镜和数字图像采集系统记录目标大气湍流退化图像，记录图像再经计算机复原，获得目 标的清晰图像，主要包括散斑成像、幸运成像和盲解卷积等。图像后处理技术利用望远镜和 数字图像采集系统记录目标大气湍流退化图像，记录图像再经计算机复原，获得目标的清晰 图像，由于需要大量目标的先验性知识且易受成像噪声影响，在实际应用中受到了很大的限 制。散斑成像和幸运成像对图像的拍摄条件有一定的要求，并且需要大量的样本信息；盲解 卷积需要知道图像大量的先验信息。另外，图像后处理技术需要进行大量的数据处理，计算 复杂很难做到实时或者准实时，在对高速运动目标成像时，受到了很大的限制。

发明内容

本发明要解决的技术问题：针对目前自适应光学技术和传统图像后处理技术存在的问题， 提供一种探测视场大、成本低、动态范围大、图像分辨率高、适用范围广的可减弱大气湍流 影响的成像方法及装置。

为了解决上述技术问题，本发明采用的技术方案为：

本发明提供一种可减弱大气湍流影响的成像方法，步骤包括：

1)采集成像目标的四维光场分布I(u,v,x,y)、受大气湍流影响的退化图像g(x,y)；

2)利用成像目标的四维光场分布I(u,v,x,y)解算大气湍流导致的全视场各视角方向相位 波前畸变

3)根据大气湍流导致的全视场各视角方向相位波前畸变估计成像系统的相干成 像点扩散函数h(x_i，y_i)，利用估计得到的相干成像点扩散函数h(x_i，y_i)求解成像系统的系统成 像点扩散函数h_I(x,y)；

4)根据退化图像g(x,y)的傅里叶变换G(f_x,f_y)、系统成像点扩散函数h_I(x,y)的傅里叶 变换H_I(f_x,f_y)在频域对退化图像g(x,y)进行解卷积运算，求解得到成像目标理想成像的傅里 叶变换I(f_x,f_y)，然后将成像目标理想成像的傅里叶变换I(f_x,f_y)进行傅里叶逆变换，得到 成像目标接近于衍射极限的理想成像i(x,y)。

优选地，所述步骤2)的详细步骤包括：

2.1)针对成像目标的四维光场分布I(u,v,x,y)进行积分运算，得到位于主望远镜焦面上 的原始参考图像L(u,v)；

2.2)在所述原始参考图像L(u,v)上选择若干关心的视角，选择一个视角作为当前视角θ_i， 跳转执行步骤2.3)；

2.3)以当前视角θ_i的中心位置(u_i,v_i)为中心，设置窗口 u∈(u_i-M,u_i+M),v∈(v_i-N,v_i+N)对成像目标的四维光场分布I(u,v,x,y)和原始参考图像 L(u,v)进行裁减，得到当前视角(u_i,v_i)的光场四维矩阵数据和基准参考图像 L_θi(u,v)，其中整数M为窗口的最大长度、整数N为窗口的最大宽度，u表示设置窗口的长 度，v表示设置窗口的宽度；

2.4)对窗口数据进行重组，固定坐标(x,y)，得到子光瞳对成像目标的一幅子图像I_xy(u,v)， 重复获取各子光瞳对成像目标的一幅子图像I_xy(u,v)，生成成像目标在当前视角θ_i的中心位置 (u_i,v_i)上的成像阵列；

2.5)针对每一幅子图像I_xy(u,v)，将当前视角θ_i的中心位置(u_i,v_i)的光场四维矩阵数据 获取各个视角上x和y方向上的波前斜率信息，最终得到由各个子图像 I_xy(u,v)的x和y两个方向上的波前斜率信息构成的两组波前斜率矩阵S_x和S_y；

2.6)判断是否所有视角均已经处理完毕，如果已经处理完毕，则跳转执行步骤2.7)；否 则选择下一个视角作为当前视角θ_i，跳转执行步骤2.3)；

2.7)得到所有视角上的两组波前斜率矩阵S_x和S_y，结合现阶段较为成熟的将波前复原 算法将两组波前斜率矩阵S_x和S_y代入波前复原算法，复原出全视场不同视角上的相位畸变的 分布

优选地，所述步骤2.1)中进行积分运算的函数表达式如式(1)所示；

L(u,v)＝∫∫I(u,v,x,y)dxdy(1)

式(1)中，L(u,v)为原始参考图像，I(u,v,x,y)为成像目标的四维光场分布。

优选地，所述步骤2.5)中运算获取各个视角上x和y方向上的波前斜率信息的函 数表达式如式(2)所示；

C(u,v)＝F^-1{[F{L_xy(u,v)}]*[F{L_θi(u,v)}]}(2)

式(2)中，C(u,v)为当前视角θ_i的子图像I_xy(u,v)与基准参考图像L_θi(u,v)的互相关值， C(x^*,y^*)为C(u,v)的峰值，(x^*,y^*)为子图像I_xy(u,v)与基准参考图像L_θi(u,v)之间的偏移量， (x^*,y^*)为对应于x和y方向上的波前斜率信息，F{Lx_y(u,v)}为当前视角θ_i的子图像I_xy(u,v)的 傅里叶变换，F{L_θi(u,v)}为基准参考图像L_θi(u,v)的傅里叶变换，F^-1为傅里叶逆变换， Ix_y(u,v)为当前视角θ_i的子图像，L_θi(u,v)为基准参考图像。

优选地，所述步骤3)中估计得到的系统相干成像点扩散函数如式(3)所示；

h(x_i,y_i)＝F{t(x,y)}(3)

式(3)中，h(x_i,y_i)为相干成像点扩散函数，(x_i,y_i)为成像面坐标，t(x,y)为成像系统 的复振幅透过率函数，(x,y)为成像系统入曈面坐标，成像系统的复振幅透过率函数t(x,y)的 计算函数表达式如式(4)所示；

式(4)中，t(x,y)为成像系统的复振幅透过率函数，(x,y)为成像系统入曈面坐标，P(x,y) 为孔径函数，为相位透过率函数。

优选地，所述步骤3)中求解系统成像点扩散函数的函数表达式如式(5)所示；

h_I(x_i,y_i)＝|h(x_i,y_i)|²(5)

式(5)中，h_I(x_i,y_i)为系统成像点扩散函数，h(x_i,y_i)为相干成像点扩散函数。

优选地，所述步骤4)中进行解卷积运算的函数表达式如式(6)所示；

$\left(\begin{array}{c}I(f_x,f_y)=\frac{G(f_x,f_y).{H^{*}}_I(f_x,f_y)}{{\left|H_I(f_x,f_y)\right|}^2}\\ G(f_x,f_y)=H_I(f_x,f_y).I(f_x,f_y)+N(f_x,f_y)\\ g(x,y)=i(x,y)*h_I(x,y)+n(x,y)\end{array}\right)---\left(6\right)$

式(6)中，I(f_x,f_y)为成像目标的理想成像i(x,y)的傅里叶变换，G(f_x,f_y)为退化图 像g(x,y)的傅里叶变换，H_I(f_x,f_y)为系统成像点扩散函数h_I(x,y)的傅里叶变换，H^*_I(f_x,f_y) 为H_I(f_x,f_y)的复共轭，I(f_x,f_y)为成像目标的理想成像i(x,y)的傅里叶变换，N(f_x,f_y)为 退化图像g(x,y)中的噪声成分n(x,y)的傅里叶变换，g(x,y)为退化图像，i(x,y)为成像目标 的理想成像，h_I(x_i,y_i)为系统成像点扩散函数，n(x,y)为噪声成分。

优选地，所述步骤4)中进行解卷积运算之前，还包括下述步骤：根据噪声与真实图像 频谱的分布特点对退化图像g(x,y)进行预处理以对高频噪声进行区分化规整并对其进行抑 制，在保证噪声抑制的前提下保持图像的低频部分、抑制图像解卷积中的高频噪声影响。

本发明还提供一种可减弱大气湍流影响的成像装置，包括主透镜、分光器件、成像信息 采集高帧频感光CCD、微透镜阵列和光场信息采集低噪声高帧频感光CCD，所述主透镜、微 透镜阵列、光场信息采集低噪声高帧频感光CCD三者依次平行布置，所述分光器件倾斜布置 于主透镜、微透镜阵列之间，所述微透镜阵列位于主透镜的焦平面上，所述光场信息采集低 噪声高帧频感光CCD布置于主透镜的共轭位置上，所述分光器件分别将来自主透镜的光分成 两路，一路通过微透镜阵列在光场信息采集低噪声高帧频感光CCD上成像得到成像目标的四 维光场分布、另一路在光场信息采集低噪声高帧频感光CCD上成像得到成像目标的受大气湍 流影响的退化图像。

优选地，所述主透镜和微透镜阵列具有相同的相对孔径。

本发明可减弱大气湍流影响的成像方法具有下述优点：

1、本发明充分利用了光场相机结构在波前探测方面的优势，在波前探测方面，不受信标 扩展性限制，同时具有大视场波前探测能力，通过一套光学系统、一次曝光来获取目标的光 场信息与成像，一次性获得多视角对应的波前畸变，具有结构简单、视场大的优点，同时绕 开了多传感器采样时同步性和均匀性的问题。

2、本发明不受传统自适应光学成像系统矫正器件如变形镜的机械结构限制，具有更大的 动态矫正范围，在图像复原方面，采取数字恢复方案，处理速度可在并行运算环境下达到准 实时，在实际应用中结构简单且具有更宽广的适用范围。

3、本发明中采取的光学成像系统分辨率不受微透镜尺寸限制，而取决于感光器件像素尺 寸，成像分辨率远高于已有的光场成像系统。

本发明可减弱大气湍流影响的成像装置具有下述优点：本发明可减弱大气湍流影响的成 像装置包括主透镜、分光器件、成像信息采集高帧频感光CCD、微透镜阵列和光场信息采集 低噪声高帧频感光CCD，通过分光器件分别将来自主透镜的光分成两路，一路通过微透镜阵 列在光场信息采集低噪声高帧频感光CCD上成像得到成像目标的四维光场分布、另一路在光 场信息采集低噪声高帧频感光CCD上成像得到成像目标的受大气湍流影响的退化图像，从而 能够为本发明可减弱大气湍流影响的成像方法提供基础数据，具有结构简单、成本低、适用 范围广的优点。

附图说明

图1为本发明实施例方法的基本流程示意图。

图2为本发明实施例装置的结构示意图。

图3为本发明实施例装置的成像应用原理示意图。

具体实施方式

如图1所示，本实施例可减弱大气湍流影响的成像方法的步骤包括：

1)采集成像目标的四维光场分布I(u,v,x,y)、受大气湍流影响的退化图像g(x,y)；

2)利用成像目标的四维光场分布I(u,v,x,y)解算大气湍流导致的全视场各视角方向相位 波前畸变

3)根据大气湍流导致的全视场各视角方向相位波前畸变估计成像系统的相干成 像点扩散函数h(x_i，y_i)，利用估计得到的相干成像点扩散函数h(x_i，y_i)求解成像系统的系统成 像点扩散函数h_I(x,y)；

4)根据退化图像g(x,y)的傅里叶变换G(f_x,f_y)、系统成像点扩散函数h_I(x,y)的傅里叶 变换H_I(f_x,f_y)在频域对退化图像g(x,y)进行解卷积运算，求解得到成像目标理想成像的傅里 叶变换I(f_x，f_y)，然后将成像目标理想成像的傅里叶变换I(f_x，f_y)进行傅里叶逆变换，得到 成像目标接近于衍射极限的理想成像i(x,y)。

本实施例中，步骤2)的详细步骤包括：

2.1)针对成像目标的四维光场分布I(u,v,x,y)进行积分运算，得到位于主望远镜焦面上 的原始参考图像L(u,v)；

2.2)在所述原始参考图像L(u,v)上选择若干关心的视角，选择一个视角作为当前视角θ_i， 跳转执行步骤2.3)；

2.3)以当前视角θ_i的中心位置(u_i,v_i)为中心，设置窗口 u∈(u_i-M,u_i+M),v∈(v_i-N,v_i+N)对成像目标的四维光场分布I(u,v,x,y)和原始参考图像 L(u,v)进行裁减，得到当前视角(u_i,v_i)的光场四维矩阵数据和基准参考图像 L_θi(u,v)，其中整数M为窗口的最大长度、整数N为窗口的最大宽度，u表示设置窗口的长 度，v表示设置窗口的宽度；

2.4)对窗口数据进行重组，固定坐标(x,y)，得到子光瞳对成像目标的一幅子图像I_xy(u,v)， 重复获取各子光瞳对成像目标的一幅子图像I_xy(u,v)，生成成像目标在当前视角θ_i的中心位置 (u_i,v_i)上的成像阵列；

2.5)针对每一幅子图像I_xy(u,v)，将当前视角θ_i的中心位置(u_i,v_i)的光场四维矩阵数据 获取各个视角上x和y方向上的波前斜率信息，最终得到由各个子图像 I_xy(u,v)的x和y两个方向上的波前斜率信息构成的两组波前斜率矩阵S_x和S_y；

2.6)判断是否所有视角均已经处理完毕，如果已经处理完毕，则跳转执行步骤2.7)；否 则选择下一个视角作为当前视角θ_i，跳转执行步骤2.3)；

2.7)得到所有视角上的两组波前斜率矩阵S_x和S_y，结合现阶段较为成熟的将波前复原 算法将两组波前斜率矩阵S_x和S_y代入波前复原算法，复原出全视场不同视角上的相位畸变的 分布由于大气湍流导致的非等晕性，各视角方向上的波前信息不同，使用窗口裁减 成像目标的光场数据，窗口大小与大气等晕区域对应，每个窗口对应一个视角，每组窗口数 据解算出两个波前斜率矩阵S_x,S_y，将波前斜率矩阵代入波前复原算法，复原出全视场不同视 角上的相位畸变的分布

本实施例中，步骤2.1)中进行积分运算的函数表达式如式(1)所示；

L(u,v)＝∫∫I(u,v,x,y)dxdy(1)

式(1)中，L(u,v)为原始参考图像，I(u,v,x,y)为成像目标的四维光场分布。

本实施例中，步骤2.5)中运算获取各个视角上x和y方向上的波前斜率信息的函 数表达式如式(2)所示；

C(u,v)＝F^-1{[F{L_xy(u,v)}]*[F{L_θi(u,v)}]}(2)

式(2)中，C(u,v)为当前视角θ_i的子图像I_xy(u,v)与基准参考图像L_θi(u,v)的互相关值， C(x^*,y^*)为C(u,v)的峰值，(x^*,y^*)为子图像I_xy(u,v)与基准参考图像L_θi(u,v)之间的偏移量， (x^*,y^*)为对应于x和y方向上的波前斜率信息，F{L_xy(u,v)}为当前视角θ_i的子图像I_xy(u,v)的 傅里叶变换，F{L_θi(u,v)}为基准参考图像L_θi(u,v)的傅里叶变换，F^-1为傅里叶逆变换， I_xy(u,v)为当前视角θ_i的子图像，L_θi(u,v)为基准参考图像。

由于大气湍流导致的相位畸变的非等晕性，各视角方向上的波前信息不同。单一视角的 相位分布可以用于估计这一视场角度上系统相干成像点扩散函数h(x_i,y_i)；本实施例 中，步骤3)中估计得到的系统相干成像点扩散函数如式(3)所示；

h(x_i,y_i)＝F{t(x,y)}(3)

式(3)中，h(x_i,y_i)为相干成像点扩散函数，(x_i,y_i)为成像面坐标，t(x,y)为成像系统 的复振幅透过率函数，(x,y)为成像系统入曈面坐标，成像系统的复振幅透过率函数t(x,y)的 计算函数表达式如式(4)所示；

式(4)中，t(x,y)为成像系统的复振幅透过率函数，(x,y)为成像系统入曈面坐标，P(x,y) 为孔径函数，为相位透过率函数。参见式(4)可知，成像系统的复振幅透过率函数 t(x,y)由孔径函数P(x,y)与相位透过率函数两部分组成。

本实施例中，步骤3)中求解系统成像点扩散函数的函数表达式如式(5)所示；

h_I(x_i,y_i)＝|h(x_i,y_i)|²(5)

式(5)中，h_I(x_i,y_i)为系统成像点扩散函数，h(x_i,y_i)为相干成像点扩散函数。对所有 的视角θ_i都进行上述估计，就获得了全视场扩展目标成像的点扩散函数。对卫星图像的不同 区域，如卫星电池板、卫星翼等分别进行处理，区域的大小与等晕区大小对应，利用获取的 全视场各个视角方向上波前信息，对处于对应等晕区域的图像的成像点扩散函数进行估计。

本实施例中，步骤4)中还包括根据噪声与真实图像频谱的分布特点，对退化图像进行 预处理，对高频噪声进行区分化规整并对其进行抑制，在保证噪声抑制的前提下保持图像的 低频部分，抑制图像解卷积中的高频噪声影响。然后利用估计的不同等晕区域的成像点扩散 函数，对处于同一等晕区域的图像进行解卷积操作，针对退化图像g(x,y)恢复达到接近理论 衍射极限分辨率的图像i(x,y)。步骤4)中进行解卷积运算的函数表达式如式(6)所示；

$\left(\begin{array}{c}I(f_x,f_y)=\frac{G(f_x,f_y).{H^{*}}_I(f_x,f_y)}{{\left|H_I(f_x,f_y)\right|}^2}\\ G(f_x,f_y)=H_I(f_x,f_y).I(f_x,f_y)+N(f_x,f_y)\\ g(x,y)=i(x,y)*h_I(x,y)+n(x,y)\end{array}\right)---\left(6\right)$

式(6)中，I(f_x,f_y)为成像目标的理想成像i(x,y)的傅里叶变换，G(f_x,f_y)为退化图 像g(x,y)的傅里叶变换，H_I(f_x,f_y)为系统成像点扩散函数h_I(x,y)的傅里叶变换，H^*_I(f_x,f_y) 为H_I(f_x,f_y)的复共轭，I(f_x,f_y)为成像目标的理想成像i(x,y)的傅里叶变换，N(f_x,f_y)为 退化图像g(x,y)中的噪声成分n(x,y)的傅里叶变换，g(x,y)为退化图像，i(x,y)为成像目标 的理想成像，h_I(x_i,y_i)为系统成像点扩散函数，n(x,y)为噪声成分。本实施例根据噪声信息 与成像信息在频域中分布的差异性，在频域中对噪声进行规整，抑制噪声对解卷积操作的影 响，带入获取的目标成像傅立叶变换G(f_x,f_y)和成像点扩散函数傅立叶变换H_I(f_x,f_y)，利用 式(6)在频域对图像进行解卷积运算，求解得到成像目标理想成像的傅里叶变换I(f_x,f_y)， 然后将成像目标理想成像的傅里叶变换I(f_x,f_y)进行傅里叶逆变换，得到成像目标接近于衍 射极限的理想成像i(x,y)。

综上所述，本实施例的实现方法是利用光场采集结构获取成像目标的四维光场分布 I(u,v,x,y)，以此解算大气湍流导致的波前畸变通过解算的波前数据对系统成像点 扩散函数h_I(x,y)进行估计，结合成像结构采集的退化图像g(x,y)，利用有效的图像处理技术， 快速的实现图像高清晰复原，有效的减弱大气湍流对大口径望远镜成像质量造成的退化影响， 同时具有大视场波前探测能力和较大像差矫正动态范围，适用范围广，不需要成像目标先验 知识，能够准实时进行图像清晰化，不受波前探测信标扩展性限制而直接使用成像目标本身 作为波前探测信标，可通过一套光学系统、一次曝光来获得全视场的由大气湍流带来的相位 畸变信息，并对图像进行清晰化复原，探测视场大、成本低、动态范围大、图像分辨率高、 适用范围广的优点。本实施例综合了光场相机大视场波前测量的优势和传统成像的优势，可 实现宽视场目标的高清晰成像，在大口径望远镜清晰成像方面具有巨大的应用前景。

如图2所示，本实施例可减弱大气湍流影响的成像装置包括主透镜1、分光器件2、成像 信息采集高帧频感光CCD3、微透镜阵列4和光场信息采集低噪声高帧频感光CCD5，主透镜 1、微透镜阵列4、光场信息采集低噪声高帧频感光CCD5三者依次平行布置，分光器件2倾 斜布置于主透镜1、微透镜阵列4之间，微透镜阵列4位于主透镜1的焦平面上，光场信息 采集低噪声高帧频感光CCD5布置于主透镜1的共轭位置上，分光器件2分别将来自主透镜 1的光分成两路，一路通过微透镜阵列4在光场信息采集低噪声高帧频感光CCD5上成像得 到成像目标的四维光场分布、另一路在光场信息采集低噪声高帧频感光CCD5上成像得到成 像目标的受大气湍流影响的退化图像。本实施例中，主透镜1和微透镜阵列4具有相同的相 对孔径，以确保成像器件的像素得到最大限度的利用。

如图3所示，以卫星6的清晰成像为例，利用大口径望远镜对成像目标进行成像时，由 于大气湍流导致波前畸变7的影响，成像质量严重退化，必须同步获取所有视角的波前畸变 7的信息并对其进行补偿矫正，才能获得目标的高分辨率成像。在获取成像目标的四维光场 分布与成像图像时，以卫星6本身作为波前探测信标，将其拉入主透镜1(探测望远镜)的 视场内，将主透镜1对焦至飞机(成像目标)，再将微透镜阵4列置于主透镜1的焦面上(二 者具有相同的相对孔径以确保成像器件的像素得到最大限度的利用)，将微透镜阵4的微透镜 焦面成像到光场信息采集低噪声高帧频感光CCD5，构成大口径光场采集系统，在光场信息 采集低噪声高帧频感光CCD5上，单个像素的强度可以表示为I(u,v,x,y)(其中u,v为微透镜 单元的坐标，x,y为瞳面像内部的坐标，与瞳面坐标一一对应)。这种结构下，在光场信息采 集低噪声高帧频感光CCD5采集目标四维光场分布信息同时，成像信息采集高帧频感光CCD3 将获取目标的模糊成像信息。在采集成像目标的四维光场分布I(u,v,x,y)、受大气湍流影响 的退化图像g(x,y)的基础上，即可应用本实施例前述可减弱大气湍流影响的成像方法获得减 弱大气湍流影响的、成像目标接近于衍射极限的理想成像i(x,y)。

以上所述仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡 属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通 技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，这些改进和润饰也应视为本 发明的保护范围。