CN104732531B - 一种高分辨率遥感图像信噪比曲线自适应获取方法 - Google Patents
一种高分辨率遥感图像信噪比曲线自适应获取方法 Download PDFInfo
- Publication number
- CN104732531B CN104732531B CN201510106719.8A CN201510106719A CN104732531B CN 104732531 B CN104732531 B CN 104732531B CN 201510106719 A CN201510106719 A CN 201510106719A CN 104732531 B CN104732531 B CN 104732531B
- Authority
- CN
- China
- Prior art keywords
- noise
- signal
- image
- mask
- point
- 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.)
- Expired - Fee Related
Links
Landscapes
- Image Processing (AREA)
Abstract
本发明一种高分辨率遥感图像信噪比曲线自适应获取方法,包括以下步骤:对整幅图像进行傅里叶变换,通过不同频率区间的掩膜,求得功率谱密度趋于稳定的频率区间;掩膜后的噪声频谱傅里叶反变换到空间域,对信号图像和噪声图像进行非重叠的小块划分,分别求得每个小块的平均值和方差;去除信号平均值和噪声方差的野值,逐步迭代拟合得到噪声方差与信号平均值的线性关系,求得直线的斜率和截距;整个灰度值范围内的信号除以噪声标准偏差得到信噪比曲线。本发明针对单幅遥感图像本身,高精度获取整个灰度值范围内的信噪比曲线。解决了光学遥感卫星信噪比在轨监测需要的均匀区域数量非常多,耗费时间很长,难以实现业务化监测的难题。
Description
技术领域
本发明涉及一种高分辨率遥感图像信噪比曲线自适应获取方法,属于高分辨率遥感卫星成像质量监测与评价技术领域,特别适用于高分辨率遥感卫星在轨成像质量业务化监测。
背景技术
卫星在轨图像的信噪比(SNR)作为光学遥感卫星的关键性能指标之一,特别对于高分辨率遥感卫星来说,信噪比直接影响了卫星的图像质量,因此,能够准确获得卫星图像的信噪比具有重要意义。目前,国内外主要通过铺设靶标法和典型目标法定期测量光学遥感卫星的信噪比。
(1)铺设靶标法:在地面铺设不同反射率的均匀靶标,在不同的太阳高度角成像条件下,卫星在轨对靶标多次成像,经过多次测量插值拟合之后得到信噪比曲线。但是,该方法不可能做到铺设所有反射率的靶标,因此,得不到整个量化灰度范围内的信噪比曲线;同时,受卫星飞行轨道和天气等条件的限制,很难保证获取靶标图像的有效性,导致了信噪比的测量周期长。
(2)典型目标法:对具有典型特征的地物目标(比如马路、河流、桥梁、圆屋顶等)进行成像,充分利用这些刃边和点目标信息,类似于布设的靶标,通过特征提取分析成像***的信噪比。该方法对目标有一定条件的限制,由于地物的非绝对均匀性,因此,比较难找到满足测量要求的典型目标。
总之,铺设靶标与典型目标法只能通过有限数量的反射率靶标和多个典型目标计算分析得到几个特定灰度值的信噪比值,难以得到整个量化灰度范围内的信噪比曲线。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供一种高分辨率遥感图像信噪比曲线自适应获取方法,针对任意高分辨率光学遥感卫星图像,不仅能够得到某个灰度值的信噪比,同时可以准确获取整个灰度值范围内的信噪比曲线。解决了目前光学遥感卫星信噪比在轨监测依赖铺设靶标或寻找典型目标,而导致的信噪比监测周期长、难以实现业务化监测以及难以获得信噪比曲线的问题。
本发明的技术方案是:一种高分辨率遥感图像信噪比曲线自适应获取方法,步骤如下:
1)在高分辨率遥感图像g(x,y)上建立二维空间直角坐标系;以g(x,y)的左上角为原点、X轴正方向向下、Y轴正方向向右建立二维空间直角坐标系,x=0,1,…,M-1,y=0,1,…,N-1,g(x,y)的大小为M×N个像素;M、N为正整数;
2)对高分辨率遥感图像g(x,y)进行傅里叶变换,在得到的傅里叶变换上建立二维频率直角坐标系;以傅里叶变换的左上角为原点、U轴正方向向下、V轴正方向向右建立二维频率直角坐标系,将g(x,y)的傅里叶变换的零频中心变换到二维频率直角坐标系上的(M/2,N/2),得到零频中心平移后的傅里叶变换,表示为G(u,v),u,v为频率变量,u=0,1,…,M-1,v=0,1,…,N-1;
3)在二维频率直角坐标系中构造频谱掩膜核函数一K1(u,v)和频谱掩膜核函数二K2(u,v);其中,K1(u,v)表示为:
该函数在区间[R1,Rmax]范围内的值为1,称为非零区间一,在其它区间该函数的值为0;
K2(u,v)表示为:
该函数在区间[R2,Rmax]范围内的值为1,称为非零区间二,在其它区间该函数的值为0;
非零区间一[R1,Rmax]和非零区间二[R2,Rmax]中的R1、R2、Rmax的关系为:
4)对上述零频中心平移后的傅里叶变换G(u,v)分别与K1(u,v)和K2(u,v)进行点乘运算,得到掩膜频谱一和掩膜频谱二,分别将掩膜频谱一和掩膜频谱二的高频中心平移到二维频率直角坐标系的(M/2,N/2),得到高频中心平移后的掩膜频谱一N1(u,v)和掩膜频谱二N2(u,v),分别对N1(u,v)和N2(u,v)进行傅里叶反变换,将得到的二维空间直角坐标系上的反变换结果分别作为掩膜后的噪声图像一n1(x,y)和掩膜后的噪声图像二n2(x,y);
5)对高分辨率遥感图像g(x,y)进行分块,块大小为m×m个像素,共得到个图像块,计算所有图像块的灰度平均值mg(x,y)(i)和灰度标准方差 表示上取整;
分别对掩膜后的噪声图像一n1(x,y)和掩膜后的噪声图像二n2(x,y)进行分块,块大小为m×m个像素,都得到个噪声图像块,计算所有噪声图像块的噪声均值和噪声方差
将噪声方差与分别乘以系数λ1、λ2,得到转化的噪声方差和表示为:
其中,λ1=(M×N)/Γ1,λ2=(M×N)/Γ2,Γ1、Γ2分别表示掩膜频谱一和掩膜频谱二中的非零个数;其中m取值为8或16;
根据转化得到的噪声方差和则图像块的噪声方差σn(x,y) 2(i)计算公式如下:
6)根据上述所得到的所有图像块的灰度平均值mg(x,y)(i)与灰度标准方差σg(x,y)(i)剔除不满足条件的图像块,通过逐步迭代拟合得到信号与噪声方差的直线方程,进而得到整个灰度范围内的信噪比曲线,具体步骤如下:
(6.1)遍历所有图像块,判断每个图像块的灰度平均值mg(x,y)(i)与标准方差σg(x,y)(i)是否满足关系式如果不满足则剔除该图像块,否则保留该图像块;
(6.2)将保留下来的每个图像块的灰度平均值与其对应的噪声图像块的噪声方差作为一个点,所有图像块的灰度平均值与其对应的噪声图像块的噪声方差构成点集(mg(x,y)(j),σn(x,y) 2(j)),其中,点集中的点个数小于划分的图像块数即
(6.3)利用最小二乘法对点集(mg(x,y)(j),σn(x,y) 2(j))进行直线拟合,并计算点集中的各个点到直线的距离,如果距离小于设定的阈值,则保留点集中的该点,否则剔除点集中的该点,该步骤中剔除的点集中的点个数记为η;
(6.4)重复步骤(6.2)、步骤(6.3),直到剔除的点集中的点个数η=0,利用点集中的剩余的点拟合直线方程,p、q为拟合系数,则得到信号S与噪声方差σ2的关系式可表示为σ2=p×S+q,信噪比即信号S与噪声标准方差σ的比值,其计算公式如下:
7)对于量化级数为L的遥感图像,其对应的灰度值范围为[0,2L-1]之间的整数,将每一个灰度值作为信号S,通过公式计算每个灰度值的信噪比,最终得到整个灰度值范围内的信噪比曲线。
所述R1、R2、Rmax的取值范围一般为0.52≤R1≤0.56,0.60≤R2≤0.64,0.66≤Rmax≤0.7。
本发明与现有技术相比有益效果为:
(1)本发明提供一种高分辨率遥感图像信噪比曲线自适应获取方法,较传统的遥感卫星依赖铺设靶标或寻找典型目标进行信噪比在轨监测的方法,该方法仅需要单幅遥感图像就可以获得某个灰度值的信噪比,同时可以高精度获取图像整个灰度值范围内的信噪比曲线,实现遥感卫星的信噪比业务化监测。
(2)本发明提供一种基于高频掩膜的遥感图像信噪比曲线自适应获取方法,通过自适应高频掩膜,除了获取高斯分布的白噪声之外,还能获取与信号相关的泊松分布的散弹噪声。
附图说明
图1是本发明高分辨率遥感图像信噪比曲线自适应获取方法流程图;
图2是本发明信号能量与噪声示意图;
图3是本发明构造的频谱掩膜核函数一和频谱掩膜核函数二的示意图;
图4a是本发明实施例采用的遥感图像,图4b是对遥感图像进行模糊的MTF卷积核,图4c是对图4a进行MTF卷积模糊后的图像;
图5a是本发明频谱掩膜区间一对应的掩膜核图像,图5b是掩膜后的噪声图像,图5c是频谱掩膜区间二对应的掩膜核图像,图5d是掩膜后的噪声图像;
图6是本发明所有图像块的信号均值与噪声方差关系图;
图7是本发明基于最小二乘法去野值后的信号均值与噪声方差,以及信号均值与噪声方差拟合得到的直线方程;
图8a是本发明实施例中加入随机噪声的100次试验,每次试验的测量值与理论值的2阶范数误差图,图8b与图8c给出了第一次、第二次试验的测量值与理论值曲线图。
具体实施方式
本发明一种高分辨率遥感图像信噪比曲线自适应获取方法,具体步骤如图1所示,该方法由以下步骤实现:
1、在高分辨率遥感图像g(x,y)上建立二维空间直角坐标系;以g(x,y)的左上角为原点、X轴正方向向下、Y轴正方向向右建立二维空间直角坐标系,x=0,1,…,M-1,y=0,1,…,N-1,g(x,y)的大小为M×N个像素;M、N为正整数。
2、对高分辨率遥感图像g(x,y)进行傅里叶变换,在得到的傅里叶变换上建立二维频率直角坐标系;以傅里叶变换的左上角为原点、U轴正方向向下、V轴正方向向右建立二维频率直角坐标系,将g(x,y)的傅里叶变换的零频中心变换到二维频率直角坐标系上的(M/2,N/2),得到零频中心平移后的傅里叶变换,表示为G(u,v),u,v为频率变量,u=0,1,…,M-1,v=0,1,…,N-1。将得到的傅里叶坐标变换进行零频坐标平移后,为了确保移动后的坐标为整数,要求M和N为偶数。
根据遥感图像的退化原理,在轨获得的高分辨率遥感图像g(x,y)由清晰图像f(x,y)与模糊核h(x,y)进行模糊卷积,再加上噪声n(x,y)所得,即:
g(x,y)=h(x,y)*f(x,y)+n(x,y)
因此,上述退化原理在频域的傅里叶变换关系表示为:
G(u,v)=H(u,v)F(u,v)+N(u,v)
对于通常的卫星图像,尤其是高分辨率遥感卫星图像,其模糊核H(u,v)的高频部分是非常小的,遥感图像通过模糊核作用过后的H(u,v)F(u,v)在高频部分的能量非常弱,可以忽略不计。如图2所示,在高频R>R0处,噪声部分的能量N(u,v)占主导作用。根据这一点可以设计噪声方差σ2的快速获取算法,从g(x,y)中获取噪声n(x,y)的方差σ2。
3、在二维频率直角坐标系中构造频谱掩膜核函数一K1(u,v)和频谱掩膜核函数二K2(u,v);其中,K1(u,v)表示为:
该函数在区间[R1,Rmax]范围内的值为1,称为非零区间一,在其它区间该函数的值为0。
K2(u,v)表示为:
该函数在区间[R2,Rmax]范围内的值为1,称为非零区间二,在其它区间该函数的值为0。
所述频谱掩膜核函数一K1(u,v)的非零区间[R1,Rmax]和频谱掩膜核函数二K2(u,v)的非零区间[R2,Rmax]中的R1、R2、Rmax的关系约束为:
根据经验值,R1、R2、Rmax的取值范围一般为0.52≤R1≤0.56,0.60≤R2≤0.64,0.66≤Rmax≤0.7。
图3是本发明构造的频谱掩膜核函数一和频谱掩膜核函数二的示意图,选取两个高频区间的频谱掩膜核函数,其主要目的是防止选取的不够准确的高频下限值R1或R2会使得G(u,v)·K1(u,v)或G(u,v)·K2(u,v)在频率高频处衰减得不够厉害,噪声部分的能量N(u,v)占主导作用的假设不正确,从而导致获取出来的噪声方差往往要比实际的σ2值要偏大。因此,为了解决这一问题,采取了选择两个不同下界的策略[R1,Rmax]和[R2,Rmax];而设定上界Rmax的主要目的是为了避免边界效应带来的影响。
4、对上述零频中心平移后的傅里叶变换G(u,v)分别与K1(u,v)和K2(u,v)进行点乘运算,得到掩膜频谱一和掩膜频谱二,分别将掩膜频谱一和掩膜频谱二的高频中心平移到二维频率直角坐标系的(M/2,N/2),得到高频中心平移后的掩膜频谱一N1(u,v)和掩膜频谱二N2(u,v),分别对N1(u,v)和N2(u,v)进行傅里叶反变换,将得到的二维空间直角坐标系上的反变换结果分别作为掩膜后的噪声图像一n1(x,y)和掩膜后的噪声图像二n2(x,y),噪声图像的大小都等于M×N个像素。
对G(u,v)分别利用频谱掩膜核函数一K1(u,v)和频谱掩膜核函数二K2(u,v)进行点乘运算,即对应元素相乘,表示为G(u,v)·K1(u,v)=[H(u,v)F(u,v)+N(u,v)]·K1(u,v)和G(u,v)·K2(u,v)=[H(u,v)F(u,v)+N(u,v)]·K2(u,v),由于在高频部分H(u,v)F(u,v)对应的能量非常弱,本发明中频域中的圆环频谱掩膜函数半径的区间范围[R1,Rmax]和[R2,Rmax]都是高频区间,所以[H(u,v)F(u,v)]·K1(u,v)≈0和[H(u,v)F(u,v)]·K2(u,v)≈0。本部分也证明了遥感图像在高频部分占主导作用的主要是噪声部分。
5、对高分辨率遥感图像g(x,y)进行分块,块大小为m×m个像素,共得到个图像块,计算所有图像块的灰度平均值mg(x,y)(i)和标准方差 表示上取整;m一般取值8或16。
分别对掩膜后的噪声图像一n1(x,y)和掩膜后的噪声图像二n2(x,y)进行分块,块大小为m×m个像素,都得到个噪声图像块,计算所有噪声图像块的噪声均值和噪声方差
由于噪声方差仅在频谱掩膜区间得到的,在频谱掩膜函数为零的区间被舍弃,因此,需要转化到整个区间,即将噪声方差与分别乘以系数λ1、λ2,得到转化的噪声方差和表示为
其中,λ1=(M×N)/Γ1,λ2=(M×N)/Γ2,Γ1、Γ2分别表示掩膜频谱一和掩膜频谱二中的非零个数;其中m取值为8或16。
根据转化得到的噪声方差和则图像块的噪声方差σn(x,y) 2(i)计算公式如下:
为了防止向外插值引起的不稳定现象,对上述计算得到的σn(x,y) 2(i)的大小加以限制,要求即最终得到的σn(x,y) 2(i)被限制在的范围内。
6、根据上述所得到的所有图像块的灰度平均值mg(x,y)(i)与灰度标准方差σg(x,y)(i)剔除不满足条件的图像块,通过逐步迭代拟合得到信号与噪声方差的直线方程,进而得到整个灰度范围内的信噪比曲线,具体步骤如下:
(1)遍历所有图像块,判断每个图像块的灰度平均值mg(x,y)(i)与标准方差σg(x,y)(i)是否满足关系式如果不满足则剔除该图像块,否则保留该图像块;
(2)将保留下来的每个图像块的灰度平均值与其对应的噪声图像块的噪声方差作为一个点,所有图像块的灰度平均值与其对应的噪声图像块的噪声方差构成点集(mg(x,y)(j),σn(x,y) 2(j)),其中,点集中的点个数小于划分的图像块数即
(3)利用最小二乘法对点集(mg(x,y)(j),σn(x,y) 2(j))进行直线拟合,并计算点集中的各个点到直线的距离,如果距离小于设定的阈值δ,则保留点集中的该点,否则剔除点集中的该点,该步骤中剔除的点集中的点个数记为η,阈值δ取值范围为0.2≤δ≤0.5;
(4)重复步骤(2)、步骤(3),直到剔除的点集中的点个数η=0,利用点集中的剩余的点拟合直线方程,p、q为拟合系数,则得到信号S与噪声方差σ2的关系式可表示为σ2=p×S+q,信噪比即信号S与噪声标准方差σ的比值,其计算公式如下:
遥感图像的噪声部分主要包括两部分,一部分为与信号大小相关的满足泊松分布的散弹噪声(系数为p,该系数与信号大小成正比例关系),另外一部分为与信号大小无关的满足高斯分布的白噪声q。总的噪声方差满足:
σ2=p×S+q
对光学遥感卫星获取的任意图像,在不需要典型均匀场景的情况下,通过对设定的高频区域进行非零掩膜,而其他频率区间进行零掩膜,利用多个非重叠的分块计算得到与信号无关的高斯分布白噪声和与信号相关的泊松分布散弹噪声,线性拟合分析得到任意一幅图像的信噪比曲线。
7、对于量化级数为L的遥感图像,其对应的灰度值范围为[0,2L-1]之间的整数,则每一个灰度值作为信号S,通过公式计算每个灰度值的信噪比,最终得到整个灰度值范围内的信噪比曲线。
例如,对于量化级数为8级的遥感图像,可得到灰度范围为[0,255]区间上的信噪比曲线;对于量化级数为10级的遥感图像,可得到灰度范围为[0,1023]区间上的信噪比曲线。
实施例
本实施例利用大小为2048×2048个像素的真实高分辨率遥感图像(图4a所示)。为得到不含噪声的遥感图像,需要进行如下处理:首先,对图像通过MTF卷积核模糊(频域的MTF卷积核如图4b所示),卷积后的图像如图4c所示,大小为2048×2048个像素;然后,对图4c采用4倍的降采样处理,得到512×512个像素的降采样图像,此时的图像可看作是不含噪声的图像。
为得到实际噪声的大小,本实施例在降采样图像中人为地加入泊松分布的信号噪声,再加入高斯分布的加性白噪声,得到的图像即为带有噪声的图像(大小为512×512个像素),并给出理论的信噪比曲线。
1、对上述人为加入噪声的图像进行傅里叶变换。
2、设定R1=0.52,R2=0.60,Rmax=0.68,则频谱掩膜核函数一的非零区间为[0.52,0.68],频域中对应的图像如图5a所示,掩膜后的噪声图像如图5b所示。频谱掩膜核函数一的非零区间为[0.60,0.68],频域中对应的图像如图5c所示,掩膜后的噪声图像如图5d所示。
3、对人为加入噪声的图像图像进行分块,每块大小为8×8个像素,共得到64×64块图像,计算所有块图像的灰度平均值。
对上述得到的噪声图像分别进行分块,每块噪声图像大小也为8×8个像素,计算所有块噪声图像的均值和方差。所有块图像的信号均值与噪声方差的关系如图6所示。
4、根据所有图像块的灰度平均值与标准方差剔除不满足条件的图像块,通过逐步迭代拟合得到灰度平均值与噪声方差的直线方程,即得到信号与噪声方差的线性关系。本实施例取阈值δ=0.5,图7是本发明基于最小二乘法剔除不满足条件的图像块后的信号的灰度均值与噪声方差,以及信号的灰度均值与噪声方差拟合得到的直线方程:σ2=0.00168×S+0.36。
5、本实施例量化级数为8级的遥感图像,根据拟合得到的信号与噪声方差的关系式σ2=0.00168×S+0.36,进一步得到[0,255]灰度范围内的信噪比曲线,图8a是本发明实施例中加入随机噪声的100次试验,每次试验的测量值与理论值的2阶范数误差图,图8b与图8c给出了第一次、第二次试验的测量值与理论值曲线图。
用以上的方案进行了测量值与理论值的对比分析,测量的曲线与理论上的曲线非常逼近,获取的噪声大小的精度优于90%。
本发明未详细说明部分属本领域技术人员公知常识。
Claims (1)
1.一种高分辨率遥感图像信噪比曲线自适应获取方法,其特征在于步骤如下:
1)在高分辨率遥感图像g(x,y)上建立二维空间直角坐标系;以g(x,y)的左上角为原点、X轴正方向向下、Y轴正方向向右建立二维空间直角坐标系,x=0,1,…,M-1,y=0,1,…,N-1,g(x,y)的大小为M×N个像素;M、N为正整数;
2)对高分辨率遥感图像g(x,y)进行傅里叶变换,在得到的傅里叶变换上建立二维频率直角坐标系;以傅里叶变换的左上角为原点、U轴正方向向下、V轴正方向向右建立二维频率直角坐标系,将g(x,y)的傅里叶变换的零频中心变换到二维频率直角坐标系上的(M/2,N/2),得到零频中心平移后的傅里叶变换,表示为G(u,v),u,v为频率变量,u=0,1,…,M-1,v=0,1,…,N-1;
3)在二维频率直角坐标系中构造频谱掩膜核函数一K1(u,v)和频谱掩膜核函数二K2(u,v);其中,K1(u,v)表示为:
该函数在区间[R1,Rmax]范围内的值为1,称为非零区间一,在其它区间该函数的值为0;
K2(u,v)表示为:
该函数在区间[R2,Rmax]范围内的值为1,称为非零区间二,在其它区间该函数的值为0;
非零区间一[R1,Rmax]和非零区间二[R2,Rmax]中的R1、R2、Rmax的关系为:
4)对上述零频中心平移后的傅里叶变换G(u,v)分别与K1(u,v)和K2(u,v)进行点乘运算,得到掩膜频谱一和掩膜频谱二,分别将掩膜频谱一和掩膜频谱二的高频中心平移到二维频率直角坐标系的(M/2,N/2),得到高频中心平移后的掩膜频谱一N1(u,v)和掩膜频谱二N2(u,v),分别对N1(u,v)和N2(u,v)进行傅里叶反变换,将得到的二维空间直角坐标系上的反变换结果分别作为掩膜后的噪声图像一n1(x,y)和掩膜后的噪声图像二n2(x,y);
5)对高分辨率遥感图像g(x,y)进行分块,块大小为m×m个像素,共得到个图像块,计算所有图像块的灰度平均值mg(x,y)(i)和灰度标准方差σg(x,y)(i), 表示上取整;
分别对掩膜后的噪声图像一n1(x,y)和掩膜后的噪声图像二n2(x,y)进行分块,块大小为m×m个像素,都得到个噪声图像块,计算所有噪声图像块的噪声均值和噪声方差
将噪声方差与分别乘以系数λ1、λ2,得到转化的噪声方差和表示为:
其中,λ1=(M×N)/Γ1,λ2=(M×N)/Γ2,Γ1、Γ2分别表示掩膜频谱一和掩膜频谱二中的非零个数;其中m取值为8或16;
根据转化得到的噪声方差和则图像块的噪声方差计算公式如下:
所述R1、R2、Rmax的取值范围为0.52≤R1≤0.56,0.60≤R2≤0.64,0.66≤Rmax≤0.7;
6)根据上述所得到的所有图像块的灰度平均值mg(x,y)(i)与灰度标准方差σg(x,y)(i)剔除不满足条件的图像块,通过逐步迭代拟合得到信号与噪声方差的直线方程,进而得到整个灰度范围内的信噪比曲线,具体步骤如下:
(6.1)遍历所有图像块,判断每个图像块的灰度平均值mg(x,y)(i)与标准方差σg(x,y)(i)是否满足关系式如果不满足则剔除该图像块,否则保留该图像块;
(6.2)将保留下来的每个图像块的灰度平均值与其对应的噪声图像块的噪声方差作为一个点,所有图像块的灰度平均值与其对应的噪声图像块的噪声方差构成点集其中,点集中的点个数小于划分的图像块数即
(6.3)利用最小二乘法对点集进行直线拟合,并计算点集中的各个点到直线的距离,如果距离小于设定的阈值,则保留点集中的该点,否则剔除点集中的该点,该步骤中剔除的点集中的点个数记为η;
(6.4)重复步骤(6.2)、步骤(6.3),直到剔除的点集中的点个数η=0,利用点集中的剩余的点拟合直线方程,p、q为拟合系数,则得到信号S与噪声方差σ2的关系式可表示为σ2=p×S+q,信噪比即信号S与噪声标准方差σ的比值,其计算公式如下:
7)对于量化级数为L的遥感图像,其对应的灰度值范围为[0,2L-1]之间的整数,将每一个灰度值作为信号S,通过公式计算每个灰度值的信噪比,最终得到整个灰度值范围内的信噪比曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510106719.8A CN104732531B (zh) | 2015-03-11 | 2015-03-11 | 一种高分辨率遥感图像信噪比曲线自适应获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510106719.8A CN104732531B (zh) | 2015-03-11 | 2015-03-11 | 一种高分辨率遥感图像信噪比曲线自适应获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104732531A CN104732531A (zh) | 2015-06-24 |
CN104732531B true CN104732531B (zh) | 2017-07-07 |
Family
ID=53456399
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510106719.8A Expired - Fee Related CN104732531B (zh) | 2015-03-11 | 2015-03-11 | 一种高分辨率遥感图像信噪比曲线自适应获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104732531B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105096319B (zh) * | 2015-09-10 | 2017-11-07 | 北京空间机电研究所 | 一种基于凝视成像的卫星在轨信噪比测试方法 |
CN106528899B (zh) * | 2015-09-10 | 2019-09-03 | 中芯国际集成电路制造(上海)有限公司 | 用于光源-掩膜优化的图形选择方法 |
CN106023175B (zh) * | 2016-05-13 | 2018-11-23 | 哈尔滨工业大学(威海) | 基于立体模式判别水中航行体热尾流的方法 |
CN107133958B (zh) * | 2017-04-26 | 2021-04-23 | 北方工业大学 | 基于块粒度预判平衡直方图的光学遥感船舶切片分割方法 |
CN108257113B (zh) * | 2017-12-28 | 2021-06-11 | 北京空间机电研究所 | 一种基于全链路的噪声分析方法 |
CN109829899B (zh) * | 2019-01-18 | 2020-08-07 | 创新奇智(广州)科技有限公司 | 一种针对钢卷端面缺陷检测的背景抑制算法 |
CN111223095B (zh) * | 2020-03-13 | 2024-03-12 | 中冶长天国际工程有限责任公司 | 一种烧结机的台车箅条间距检测方法及*** |
CN117011193B (zh) * | 2023-09-28 | 2023-12-05 | 生态环境部长江流域生态环境监督管理局生态环境监测与科学研究中心 | 一种轻量化凝视卫星视频去噪方法及去噪*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101441764A (zh) * | 2008-12-31 | 2009-05-27 | 中国资源卫星应用中心 | 一种mtfc遥感图像复原方法 |
EP2429214A2 (en) * | 2010-09-08 | 2012-03-14 | AKG Acoustics GmbH | Method for acoustic signal tracking |
CN103679652A (zh) * | 2013-11-29 | 2014-03-26 | 北京空间机电研究所 | 一种大幅提高成像质量的图像复原*** |
CN104083177A (zh) * | 2014-07-08 | 2014-10-08 | 天津大学 | 数字x射线成像***的分辨率性能评估方法 |
-
2015
- 2015-03-11 CN CN201510106719.8A patent/CN104732531B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101441764A (zh) * | 2008-12-31 | 2009-05-27 | 中国资源卫星应用中心 | 一种mtfc遥感图像复原方法 |
EP2429214A2 (en) * | 2010-09-08 | 2012-03-14 | AKG Acoustics GmbH | Method for acoustic signal tracking |
CN103679652A (zh) * | 2013-11-29 | 2014-03-26 | 北京空间机电研究所 | 一种大幅提高成像质量的图像复原*** |
CN104083177A (zh) * | 2014-07-08 | 2014-10-08 | 天津大学 | 数字x射线成像***的分辨率性能评估方法 |
Non-Patent Citations (5)
Title |
---|
A Highly Accurate and Low Bias SNR Estimator:Algorithm and Implementation;Chao GONG 等;《RADIOENGINEERING》;20111231;第976-981页 * |
A Signal-to-Noise Ratio Estimator for Generalized Linear Model Systems;Gabriela Czanner 等;《World Congress on Engineering》;20080704;正文第1-7页 * |
MTFC 在光学遥感成像***优化设计中的应用研究;满益云 等;《航天返回与遥感》;20071231;第39-47页 * |
The modulation transfer function and signal-to-noise ratio of different digital filters: a technical approach;DD Brullmann,B d’Hoedt;《Dentomaxillofacial Radiology》;20111231;第222-229页 * |
光学遥感器像质预估与评价技术研究;刘兆军 等;《红外与激光工程》;20061031;第222-226页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104732531A (zh) | 2015-06-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104732531B (zh) | 一种高分辨率遥感图像信噪比曲线自适应获取方法 | |
CN104123464B (zh) | 一种高分辨率InSAR时序分析反演地物高程与地面沉降量的方法 | |
CN101980293B (zh) | 一种基于刃边图像的高光谱遥感***mtf检测方法 | |
CN104535981B (zh) | 海杂波Pareto幅度分布参数的双分位点估计方法 | |
CN105651263B (zh) | 浅海水深多源遥感融合反演方法 | |
CN113534159B (zh) | 一种基于sar卫星遥感数据的海面风速反演方法及*** | |
CN106772342A (zh) | 一种适用于大梯度地表沉降监测的时序差分雷达干涉方法 | |
CN110110654B (zh) | 一种针对下降型海洋内孤立波的振幅反演方法及装置 | |
CN103645476B (zh) | 一种合成孔径雷达差分干涉图序列的时空同质滤波方法 | |
CN103822598B (zh) | 地基sar在时间去相关严重区域的形变监测方法 | |
CN103646395B (zh) | 一种基于格网法的高精度影像配准方法 | |
CN105678757A (zh) | 一种物***移测量方法 | |
CN102779333A (zh) | 一种基于卡尔曼滤波的光学影像复原方法 | |
CN103839234A (zh) | 一种基于可控核的双几何非局部均值图像去噪方法 | |
CN103236067B (zh) | 一种像素级sar影像时间序列构建的局部自适应配准方法 | |
CN106530334A (zh) | 一种机载干涉合成孔径雷达复图像配准方法及复图像配准*** | |
CN104574424B (zh) | 基于多分辨率dct边缘梯度统计的无参照图像模糊度评价方法 | |
CN105678716B (zh) | 一种地基sar大气干扰相位校正方法及装置 | |
CN103325104B (zh) | 基于迭代稀疏表达的人脸图像超分辨率重建方法 | |
CN106683084A (zh) | 一种基于行间像偏移估计的tdi图像形变程度客观评价方法 | |
CN102928872B (zh) | 一种提高滩坝砂储层地震识别性能与描述精度的方法 | |
CN104157009A (zh) | 一种多源遥感影像质量定量比选方法 | |
CN105389466A (zh) | 一种校正尺度效应的中低分辨率遥感产品真值获取方法 | |
CN115830476A (zh) | 地形因子空间降尺度方法 | |
CN115980697A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170707 Termination date: 20180311 |