CN102176741A - 星空图像中帧转移型ccd传感器拖尾亮线消除方法 - Google Patents
星空图像中帧转移型ccd传感器拖尾亮线消除方法 Download PDFInfo
- Publication number
- CN102176741A CN102176741A CN 201110005161 CN201110005161A CN102176741A CN 102176741 A CN102176741 A CN 102176741A CN 201110005161 CN201110005161 CN 201110005161 CN 201110005161 A CN201110005161 A CN 201110005161A CN 102176741 A CN102176741 A CN 102176741A
- Authority
- CN
- China
- Prior art keywords
- sigma
- bright line
- row
- subgraph
- formula
- 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.)
- Pending
Links
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种星空图像中帧转移型CCD传感器拖尾亮线消除方法,用于解决现有的拖尾亮线检测方法由于增加的噪声点而降低了弱小目标的信噪比的技术问题。技术方案是采用二维小波变换的方法把原图分解为低频子图与垂直、水平、对角高频子图,而CCD拖尾为垂直或水平亮线,因此只有低频子图、垂直或水平子图包含亮线,只对这些子图进行亮线的消除,有效的保留其他细节信息;在检测时使用基于列灰度和奇异值检测方法,复杂度低;消除时充分考虑到过亮恒星点的影响,对亮线灰度强度的估计更加准确。由于不使用插值方法,不会削减弱小目标的信噪比,相反会使得拖尾亮线上目标的信噪比平均增加5dB以上。
Description
技术领域
本发明涉及一种CCD传感器拖尾亮线消除方法,特别是一种星空图像中帧转移型CCD传感器拖尾亮线消除方法。
背景技术
帧转移型CCD传感器在拍摄星空背景图像时,当视场有过亮恒星时会产生拖尾亮线,这种拖尾亮线对目标检测有严重影响,现有的拖尾亮线消除方法主要有:改变CCD传感器硬件设计和图像处理方法。
文献“星图中亮线的自动检测算法,光电技术应用,2009,Vol.24(4),p41-44”公开了一种星图拖尾亮线的自动去除算法。该方法利用恒星成像的频域各向同性性质和Fourier变换的自配准性质,以及Fourier谱对噪声和恒星不敏感性对亮线进行检测,并分别在空域使用中值滤波方法或在频域使用线性、三次样条插值、多项式拟合插值等方法对亮线进行滤除。文献所述方法在检测拖尾亮线时对图像进行Fourier变换,实际上由于星空图像大部分为黑色背景,恒星及目标的白点占图像少部分,只需要对原图统计列的灰度和极值即可检测出拖尾亮线的位置,Fourier变换增加了检测的复杂度;在消除拖尾时,文献中使用了插值的方法,这会导致增加新的噪声点;另外,插值方法会使得拖尾亮线上的弱小目标变得更弱,目标信噪比平均下降3dB。
发明内容
为了克服现有的拖尾亮线检测方法由于增加的噪声点而降低了弱小目标的信噪比的不足,本发明提供一种星空图像中帧转移型CCD传感器拖尾亮线消除方法。该方法使用二维小波分解,将拖尾亮线分离到低频分量及垂直分量中分别进行检测消除处理,可以避免现有的使用插值方法而导致增加新的噪声点问题,可以增加弱小目标的信噪比。
本发明解决其技术问题所采用的技术方案:一种星空图像中帧转移型CCD传感器拖尾亮线消除方法,其特点包括以下步骤:
(a)用S(i,j)表示像素(i,j)的测量信号强度,I(i,j)表示没有拖尾亮线影响的信号强度,建立帧转移型CCD传感器拖尾亮线模型:
式中,a=ηΔτ/τ,其中i,j是图像像素坐标,1≤i≤M,1≤j≤N,令根据公式(4),把消除拖尾亮线转换为从S(i,j)中减除附加的部分,还原输入图像I(i,j),即估计Tj的值;
(b)二维小波将图像S(i,j)分解为:
式中,j,k表示图像上两个方向的位移,是对应尺度空间Vj的展开系数; 分别对应小波空间的小波展开系数,其中Vj强调低频特征,分别强调各自系数是用以描述图像的水平、垂直和对角特征;公式(5)中的系数计算:
各列灰度和求得灰度和的均值μB与方差σ2如下:
式中,N为总列数。根据均值与方差确定阈值T:
T=μB+λ·σB (15)
列灰度和大于阈值的列即存在亮线的列,RL(j)为1表示第j列存在拖尾亮线;
各包含亮线的列减去该列的亮线灰度强度估计值Tj,如公式(17)所示,即在子图上消除了拖尾亮线;
I(i,j)=S(i,j)-Tj (17)
式中,S(i,j)是像素(i,j)的测量信号强度,I(i,j)是没有拖尾亮线影响的信号强度;
(d)将小波分解的第一层到第三层的水平分量子图和对角分量子图、经过修改的第一层到第三层的垂直分量子图、第三层低频分量子图进行小波的逆变换,使用二维mallat对小波分解的低频、各高频子图进行重构,即可得到消除拖尾后的图像;公式如下:
本发明的有益结果是:由于采用二维小波变换的方法把原图分解为低频子图与垂直、水平、对角高频子图,而CCD拖尾为垂直或水平亮线,因此只有低频子图、垂直或水平子图包含亮线,只对这些子图进行亮线的消除,有效的保留其他细节信息;在检测时使用基于列灰度和奇异值检测方法,复杂度低;消除时充分考虑到过亮恒星点的影响,对亮线灰度强度的估计更加准确。由于不使用插值方法,不会削减弱小目标的信噪比,相反会使得拖尾亮线上目标的信噪比平均增加5dB以上。
下面结合具体实施方式对本发明作详细说明。
具体实施方式
1.建立帧转移型CCD传感器拖尾亮线的模型。
星空图像视场内恒星亮度太高造成CCD传感器整行或整列电子溢出,拖尾区域通常为过某颗高亮恒星的中心,且贯穿整幅图像的行或列方向的直线。以下以垂直亮线为例阐述消除方法;原图像分为两部分,一部分为没有拖尾的图像,另一部分为拖尾亮线。i为第j列像素的行坐标,恒星光源在曝光周期以及转移时间内均恒定。在曝光周期τ内,像素i积累的光子数为:
N′i,j=kτAi,j (1)
式中,Ai,j是i像素上的辐射强度,k是比例常数。光子帧转移时间为Δτ,则因第i行以前的像元光子曝光附加在第i行的光子数为:
ΔNi,j=kηΔτAi,j (2)
式中,η是电荷转移单位时间内曝光光子生成数目与曝光单位时间光子生成数目之比。
当CCD所有像素都被照明的情形,像素i的总电荷为:
S(i,j)表示像素(i,j)的测量信号强度,I(i,j)表示没有拖尾亮线影响的信号强度。S(i,j)=k1Ni,j,k1为比例常数,是反映CCD光电转移的因子,则I(i,j)=k1kτAi,j,得图像S(i,j):
根据公式(4),把消除拖尾亮线问题转换为从S(i,j)中减除附加的部分,还原输入图像I(i,j)问题,即估计出Tj的值,并用原图减去Tj。
2.对原始图像进行二维小波分解。
对原始星空图像进行三层二维Haar小波Mallat快速小波分解,得到低频子图像、高频子图像;高频子图像由水平高频子图像、垂直高频子图像和对角高频子图像组成。
原图像为S(i,j),二维小波分解为:
式中,j,k表示图像上两个方向的位移,是对应尺度空间Vj的展开系数; 分别对应小波空间的小波展开系数,其中Vj强调低频特征,分别强调各自系数是用以描述图像的水平、垂直和对角特征。公式(5)中的系数通过二维的Mallat快速级联算法计算:
进行二维小波分解时选择Haar小波基函数,Haar小波函数是一个差分函数,它的两尺度方程φ(x)为:
3.对低频子图与垂直高频子图上的亮线进行检测与消除。
各列灰度和求得灰度和的均值μB与方差σ2如公式(13)和(14)所示:
式中,N为总列数。根据均值与方差确定阈值T,列灰度和大于阈值的列即存在亮线的列,RL(j)为1表示第j列存在拖尾亮线情况。
T=μB+λ·σB (15)
找到包含亮线的列之后,估计出该列的亮线灰度强度,即Tj的值。Tj可由该列背景灰度均值近似逼近,即但过亮恒星点的像素灰度将会影响逼近的效果,过亮点的去除方法:统计该列上每行像素的灰度值,过亮点会出现奇异峰值,根据该列中各行灰度值的均值与方差设置适当的阈值,把大于阈值的点去除后再求灰度均值即Tj的近似逼近值。
最后,各包含亮线的列减去该列的亮线灰度强度估计值Tj,如公式(17)所示,即在子图上消除了拖尾亮线。
I(i,j)=S(i,j)-Tj (17)
式中,S(i,j)为像素(i,j)的测量信号强度,I(i,j)为没有拖尾亮线影响的信号强度。
4.利用各子图进行小波逆变换重构原图。
将小波分解的第一层到第三层的水平分量子图和对角分量子图、经过修改的第一层到第三层的垂直分量子图、第三层低频分量子图进行小波的逆变换,使用二维mallat算法重构出消除拖尾亮线后的初始输入图像。公式如下:
Claims (1)
1.一种星空图像中帧转移型CCD传感器拖尾亮线消除方法,其特征在于包括下述步骤:
(a)用S(i,j)表示像素(i,j)的测量信号强度,I(i,j)表示没有拖尾亮线影响的信号强度,建立帧转移型CCD传感器拖尾亮线模型:
(b)二维小波将图像S(i,j)分解为:
式中,j,k表示图像上两个方向的位移,是对应尺度空间Vj的展开系数; 分别对应小波空间的小波展开系数,其中Vj强调低频特征,分别强调各自系数是用以描述图像的水平、垂直和对角特征;公式(5)中的系数计算:
各列灰度和求得灰度和的均值μB与方差σ2如下:
式中,N为总列数。根据均值与方差确定阈值T:
T=μB+λ·σB (15)
列灰度和大于阈值的列即存在亮线的列,RL(j)为1表示第j列存在拖尾亮线;
各包含亮线的列减去该列的亮线灰度强度估计值Tj,如公式(17)所示,即在子图上消除了拖尾亮线;
I(i,j)=S(i,j)-Tj (17)
式中,S(i,j)是像素(i,j)的测量信号强度,I(i,j)是没有拖尾亮线影响的信号强度;
(d)将小波分解的第一层到第三层的水平分量子图和对角分量子图、经过修改的第一层到第三层的垂直分量子图、第三层低频分量子图进行小波的逆变换,使用二维mallat对小波分解的低频、各高频子图进行重构,即可得到消除拖尾后的图像;公式如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110005161 CN102176741A (zh) | 2011-01-06 | 2011-01-06 | 星空图像中帧转移型ccd传感器拖尾亮线消除方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110005161 CN102176741A (zh) | 2011-01-06 | 2011-01-06 | 星空图像中帧转移型ccd传感器拖尾亮线消除方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102176741A true CN102176741A (zh) | 2011-09-07 |
Family
ID=44519835
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110005161 Pending CN102176741A (zh) | 2011-01-06 | 2011-01-06 | 星空图像中帧转移型ccd传感器拖尾亮线消除方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102176741A (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104052904A (zh) * | 2014-05-28 | 2014-09-17 | 刘远明 | 多分辨率反锐化掩膜usm的图像处理方法及装置 |
CN104159035A (zh) * | 2014-08-22 | 2014-11-19 | 深圳市中兴移动通信有限公司 | 一种星轨拍摄效果自动优化的方法及装置 |
CN104394334A (zh) * | 2014-11-18 | 2015-03-04 | 中国科学院西安光学精密机械研究所 | 图像传感器在轨噪声智能抑制方法 |
CN105894477A (zh) * | 2016-06-03 | 2016-08-24 | 深圳市樊溪电子有限公司 | 天文图像噪声去除方法 |
CN107360379A (zh) * | 2012-12-27 | 2017-11-17 | 松下电器(美国)知识产权公司 | 信息通信方法 |
CN107449416A (zh) * | 2017-06-20 | 2017-12-08 | 中国人民解放军国防科学技术大学 | 基于矢量累积的恒星拖尾星点提取方法 |
CN108780571A (zh) * | 2015-12-31 | 2018-11-09 | 上海联影医疗科技有限公司 | 一种图像处理方法和*** |
CN111405207A (zh) * | 2018-12-27 | 2020-07-10 | 天津大学青岛海洋技术研究院 | 一种应用于高速脉冲传感器的重构图像拖尾消除算法 |
CN113160080A (zh) * | 2021-04-16 | 2021-07-23 | 桂林市啄木鸟医疗器械有限公司 | 一种cr影像降噪方法、装置、设备及介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101853501A (zh) * | 2010-05-13 | 2010-10-06 | 西北工业大学 | CCD传感器Smear效应处理方法 |
-
2011
- 2011-01-06 CN CN 201110005161 patent/CN102176741A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101853501A (zh) * | 2010-05-13 | 2010-10-06 | 西北工业大学 | CCD传感器Smear效应处理方法 |
Non-Patent Citations (4)
Title |
---|
《中国体视学与图像分析》 20101130 高建伟等 一种帧转移型CCD传感器Smear效应消除方法 全文 1 第15卷, 第4期 * |
《中国科技信息》 20070630 韩采芹等 帧转移型CCD相机图像拖影研究 全文 1 , 第11期 * |
《信号处理》 20100430 黄宗福等 天文全帧CCD图像拖尾的快速去除方法 全文 1 第26卷, 第4期 * |
《科学技术与工程》 20090831 韩采芹等 科学级帧转移型CCD相机拖影问题的处理 全文 1 第9卷, 第15期 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107395977B (zh) * | 2012-12-27 | 2019-12-17 | 松下电器(美国)知识产权公司 | 信息通信方法 |
CN107360379A (zh) * | 2012-12-27 | 2017-11-17 | 松下电器(美国)知识产权公司 | 信息通信方法 |
CN107395977A (zh) * | 2012-12-27 | 2017-11-24 | 松下电器(美国)知识产权公司 | 信息通信方法 |
CN107360379B (zh) * | 2012-12-27 | 2020-06-30 | 松下电器(美国)知识产权公司 | 信息通信方法 |
CN104052904A (zh) * | 2014-05-28 | 2014-09-17 | 刘远明 | 多分辨率反锐化掩膜usm的图像处理方法及装置 |
CN104159035A (zh) * | 2014-08-22 | 2014-11-19 | 深圳市中兴移动通信有限公司 | 一种星轨拍摄效果自动优化的方法及装置 |
CN104159035B (zh) * | 2014-08-22 | 2019-08-16 | 努比亚技术有限公司 | 一种星轨拍摄效果自动优化的方法及装置 |
CN104394334A (zh) * | 2014-11-18 | 2015-03-04 | 中国科学院西安光学精密机械研究所 | 图像传感器在轨噪声智能抑制方法 |
CN104394334B (zh) * | 2014-11-18 | 2017-06-27 | 西北工业大学 | 图像传感器在轨噪声智能抑制方法 |
US11880978B2 (en) | 2015-12-31 | 2024-01-23 | Shanghai United Imaging Healthcare Co., Ltd. | Methods and systems for image processing |
CN108780571B (zh) * | 2015-12-31 | 2022-05-31 | 上海联影医疗科技股份有限公司 | 一种图像处理方法和*** |
CN108780571A (zh) * | 2015-12-31 | 2018-11-09 | 上海联影医疗科技有限公司 | 一种图像处理方法和*** |
CN105894477B (zh) * | 2016-06-03 | 2019-04-02 | 深圳市樊溪电子有限公司 | 天文图像噪声去除方法 |
WO2017206322A1 (zh) * | 2016-06-03 | 2017-12-07 | 深圳市樊溪电子有限公司 | 天文图像噪声去除方法 |
CN105894477A (zh) * | 2016-06-03 | 2016-08-24 | 深圳市樊溪电子有限公司 | 天文图像噪声去除方法 |
CN107449416B (zh) * | 2017-06-20 | 2019-09-03 | 中国人民解放军国防科学技术大学 | 基于矢量累积的恒星拖尾星点提取方法 |
CN107449416A (zh) * | 2017-06-20 | 2017-12-08 | 中国人民解放军国防科学技术大学 | 基于矢量累积的恒星拖尾星点提取方法 |
CN111405207A (zh) * | 2018-12-27 | 2020-07-10 | 天津大学青岛海洋技术研究院 | 一种应用于高速脉冲传感器的重构图像拖尾消除算法 |
CN111405207B (zh) * | 2018-12-27 | 2022-07-26 | 天津大学青岛海洋技术研究院 | 一种应用于高速脉冲传感器的重构图像拖尾消除方法 |
CN113160080A (zh) * | 2021-04-16 | 2021-07-23 | 桂林市啄木鸟医疗器械有限公司 | 一种cr影像降噪方法、装置、设备及介质 |
CN113160080B (zh) * | 2021-04-16 | 2023-09-22 | 桂林市啄木鸟医疗器械有限公司 | 一种cr影像降噪方法、装置、设备及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102176741A (zh) | 星空图像中帧转移型ccd传感器拖尾亮线消除方法 | |
US11215700B2 (en) | Method and system for real-time motion artifact handling and noise removal for ToF sensor images | |
Yang et al. | Depth recovery using an adaptive color-guided auto-regressive model | |
CN101727662B (zh) | Sar图像非局部均值去斑方法 | |
CN102768760B (zh) | 一种基于图像纹理的图像快速去雾方法 | |
CN103369209A (zh) | 视频降噪装置及方法 | |
CN102682442B (zh) | 基于光流场的运动目标超分辨率图像重建方法 | |
CN103198463A (zh) | 基于整体结构和空间细节信息融合的光谱图像全色锐化方法 | |
CN101661611A (zh) | 基于贝叶斯非局部均值滤波器的实现方法 | |
CN102800086A (zh) | 一种海上场景显著性检测方法 | |
CN105913392A (zh) | 复杂环境下退化图像综合质量提升方法 | |
CN103136734A (zh) | POCS超分辨率图像重建时边缘Halo效应的抑制方法 | |
CN100417191C (zh) | 一种图像组合降噪方法 | |
CN101491080A (zh) | 用于图像去噪和可控去模糊的保曲率滤波器 | |
CN104272347A (zh) | 去除包含在静止图像中的雾的图像处理装置及其方法 | |
CN105631898A (zh) | 基于空时显著性融合的红外运动目标检测方法 | |
CN104680485A (zh) | 一种基于多分辨率的图像去噪方法及装置 | |
Sharma et al. | Satellite image contrast and resolution enhancement using discrete wavelet transform and singular value decomposition | |
CN112712467A (zh) | 基于计算机视觉与色彩滤波阵列的图像处理方法 | |
CN105184761A (zh) | 一种基于小波分析的图像去雨方法及*** | |
CN102789634B (zh) | 一种获取光照均一化图像的方法 | |
JP6819794B2 (ja) | レーダ画像処理装置、レーダ画像処理方法およびレーダ画像処理プログラム | |
CN102547363A (zh) | 基于轮廓波变换域图像特征的无参考图像质量评测方法 | |
CN102254306B (zh) | 一种基于图像简化分层模型的实时图像去雾方法 | |
CN101908206B (zh) | 一种基于mca的sar图像抑噪方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20110907 |