CN106651799B - 一种基于地物光谱属性的坏线修复方法 - Google Patents
一种基于地物光谱属性的坏线修复方法 Download PDFInfo
- Publication number
- CN106651799B CN106651799B CN201611189755.6A CN201611189755A CN106651799B CN 106651799 B CN106651799 B CN 106651799B CN 201611189755 A CN201611189755 A CN 201611189755A CN 106651799 B CN106651799 B CN 106651799B
- Authority
- CN
- China
- Prior art keywords
- bad
- pixel
- image
- bad line
- spectrum
- 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
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000001228 spectrum Methods 0.000 title claims abstract description 25
- 230000003595 spectral effect Effects 0.000 claims abstract description 11
- 238000007796 conventional method Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 abstract description 6
- 230000008439 repair process Effects 0.000 abstract description 6
- 201000004569 Blindness Diseases 0.000 abstract description 2
- 239000000284 extract Substances 0.000 abstract description 2
- 238000003384 imaging method Methods 0.000 abstract 1
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000013179 statistical model Methods 0.000 description 2
- 238000000701 chemical imaging Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/77—Retouching; Inpainting; Scratch removal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于地物光谱属性的坏线修复方法，该方法首先提取坏线像元非坏线部分对应光谱曲线，以该光谱曲线片段为依据，通过全影像搜索，寻找与之相似度最高的地物光谱曲线，以此光谱曲线对应地物类型作为坏线像元地物类型，并用对应光谱值修复坏线像元的像元值；该方法不再盲目的利用统计关系修复坏线，而是利用先验知识判断坏线像元的地物类型，再依靠同时成像的其他相同像元的光谱特征，填补坏线区域。这种坏线修复方法，不仅修复了破损的遥感数据，提高数据的利用率。相比较于传统的方法，该方法不盲目的依赖于统计相关性，而是依靠地物本身的光谱特征来实现坏线区域的修复，科学可靠，精度更好。
Description
技术领域
本发明涉及一种坏线修复方法，特别涉及一种基于地物光谱属性的坏线修复方法，属于图像处理研究领域。
背景技术
坏线是指影像上一整列或一行明显的暗纹，它的像元值一般为很小的非典型的像元值，如0或者很小的值。坏线存在于很多类型的影像上，特别是高光谱影像。其形成的主要原因是传感器的定标问题。因为面阵列CCD器件上万个探测单元的定标十分困难。坏线的存在严重的影响了遥感影像的使用，必须加以纠正。而现有的坏线校正方法多是依靠坏线像元与左右相邻像元的相关性，建立统计模型，利用相邻像元对其进行恢复。这些方法对相邻像元的依赖性很强，缺乏理论依据，修复像元值为统计意义的估计值，而并非真实的地物的像元值。
发明内容
本发明要解决的技术问题是克服现有的坏线校正方法多是依靠坏线像元与左右相邻像元的相关性，建立统计模型，利用相邻像元对其进行恢复，这些方法对相邻像元的依赖性很强，缺乏理论依据，修复像元值为统计意义的估计值，而并非真实的地物的像元值的缺陷，提供一种基于地物光谱属性的坏线修复方法。
为了解决上述技术问题，本发明提供了如下的技术方案：
本发明提供了一种基于地物光谱属性的坏线修复方法，包括以下步骤：
1)定义多波段图像有p行，q列和n个波段，其图像大小为p*q*n，并假设该图像第b个波段的第i行第j列为坏线像元；
2)提取该影像除坏线波段外n-1个波段第i行第j列像元值，也就是坏线像元的地物光谱特征，用fbad表示，其数学表达式为：
fbad＝[DN1(i，j)，DN2(i，j)，…，DNb-1(i，j)，DNb+1(i，j)，…，DNn(i，j)]
影像上任意像元(x，y)除波段b之外地物光谱特征fref可以表示为：
fref＝[DN1(x，y)，DN2(x，y)，…，DNb-1(x，y)，DNb+1(x，y)，…，DNn(x，y)]
利用代价函数计算fbad与影像上其他像元光谱曲线fref的相似度，
ε＝fref-fbad
对于整幅影像就可以得到ε的集合E，那么E中ε的最小值就被认为该像元与fbad的相似度最高；
3)利用该像元第b个波段的像元值代替第b个波段的第i行第j列为坏线像元，以此实现其坏线的修复。
作为本发明的一种优选技术方案，选取TM图像的一块，并人为制造出坏线，通过比较传统方法和本方法的校正结果，分析所述象元值的修复精度和准确度。
本发明所达到的有益效果是：这种坏线修复方法，不仅修复了破损的遥感数据，提高数据的利用率。相比较于传统的方法，该方法不盲目的依赖于统计相关性，而是依靠地物本身的光谱特征来实现坏线区域的修复，科学可靠，精度更好。
附图说明
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。在附图中：
图1是本发明的Hyperion第15波段第114列的坏线示意图；
图2是本发明的Hyperion坏线修复后结果图；
图3是本发明的原始图像；
图4是本发明的增加坏线后的图像；
图5是本发明的基于地物光谱属性法修复图像。
具体实施方式
以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
实施例1
如图1-5所示，本发明提供了一种基于地物光谱属性的坏线修复方法，首先定义多波段图像有p行，q列和n个波段，其图像大小为p*q*n，并假设该图像第b个波段的第i行第j列为坏线像元；再提取该影像除坏线波段外n-1个波段第i行第j列像元值，也就是坏线像元的地物光谱特征，用fbad表示，其数学表达式为：
fbad＝[DN1(i，j)，DN2(i，j)，…，DNb-1(i，j)，DNb+1(i，j)…，DNn(i，j)]
影像上任意像元(x，y)除波段b之外地物光谱特征fref可以表示为：
fref＝[DN1(x，y)，DN2(x，y)，…，DNb-1(x，y)，DNb+1(x，y)，…，DNn(x，y)]
利用代价函数计算fbad与影像上其他像元光谱曲线fref的相似度，
ε＝fref-fbad
然后，对于整幅影像就可以得到ε的集合E，那么E中ε的最小值就被认为该像元与fbad的相似度最高；并利用该像元第b个波段的像元值代替第b个波段的第i行第j列为坏线像元，以此实现其坏线的修复。
在坏线修复好之后，选取TM图像的一块(如图3)，并人为制造出坏线(如图4)，通过比较传统方法和本方法的校正结果，分析其修复精度和准确度。例如使用了邻域均值法和本发明所述的基于地物光谱属性的坏线修复方法对TM数据的人造坏线进行修复，并通过比较两种方法的计算结果和真值之间的差异来评价不同方法的优劣。分别均方根误差(RMSE)，和平均精度作为评价不同方法的指标。通过计算，传统方法修复结果的平均精度为96.2％，均方根误差1.3。而基于地物光谱属性坏线修复法修复结果平均精度为99.6％，均方根误差1.1。精度提高了3.4个百分点。
方法原理：该方法首先提取坏线像元非坏线部分对应光谱曲线，以该光谱曲线片段为依据，通过全影像搜索，寻找与之相似度最高的地物光谱曲线，以此光谱曲线对应地物类型作为坏线像元地物类型，并用对应光谱值修复坏线像元的像元值；该方法不再盲目的利用统计关系修复坏线，而是利用先验知识判断坏线像元的地物类型，再依靠同时成像其他相同像元的光谱特征，填补坏线区域。
本发明所达到的有益效果是：这种坏线修复方法，不仅修复了破损的遥感数据，提高数据的利用率。相比较于传统的方法，该方法不盲目的依赖于统计相关性，而是依靠地物本身的光谱特征来实现坏线区域的修复，科学可靠，精度更好。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (2)
1.一种基于地物光谱属性的坏线修复方法，其特征在于，包括以下步骤：
1)定义多波段图像有p行，q列和n个波段，其图像大小为p*q*n，并假设该图像第b个波段的第i行第j列为坏线像元；
2)提取影像除坏线波段外n-1个波段第i行第j列像元值，也就是坏线像元的地物光谱特征，用fbad表示，其数学表达式为：
fbad＝[DN1(i，j)，DN2(i，j)，…，DNb-1(i，j)，DNb+1(i，j)，…，DNn(i，j)]
影像上任意像元(x，y)除波段b之外地物光谱特征fref可以表示为：
fref＝[DN1(x，y)，DN2(x，y)，…，DNb-1(x，y)，DNb+1(x，y)，…，DNn(x，y)]
利用代价函数计算fbad与影像上其他像元光谱曲线fref的相似度，
ε＝fref-fbad
对于整幅影像就可以得到ε的集合E，那么E中ε的最小值就被认为该像元与fbad的相似度最高；
3)利用该像元第b个波段的像元值代替第b个波段的第i行第j列为坏线像元，以此实现其坏线的修复。
2.根据权利要求1所述的一种基于地物光谱属性的坏线修复方法，其特征在于，选取TM图像的一块，并人为制造出坏线，通过比较传统方法和本方法的校正结果，分析所述像 元值的修复精度和准确度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611189755.6A CN106651799B (zh) | 2016-12-16 | 2016-12-16 | 一种基于地物光谱属性的坏线修复方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611189755.6A CN106651799B (zh) | 2016-12-16 | 2016-12-16 | 一种基于地物光谱属性的坏线修复方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106651799A CN106651799A (zh) | 2017-05-10 |
| CN106651799B true CN106651799B (zh) | 2019-05-28 |
Family
ID=58833398
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611189755.6A Expired - Fee Related CN106651799B (zh) | 2016-12-16 | 2016-12-16 | 一种基于地物光谱属性的坏线修复方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106651799B (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104380713A (zh) * | 2014-03-12 | 2015-02-25 | 深圳市大疆创新科技有限公司 | 数码影像的坏点矫正方法和系统 |
| CN105787925A (zh) * | 2016-02-03 | 2016-07-20 | 中国科学院空间应用工程与技术中心 | 推扫型光学遥感载荷原始图像坏线的自动检测方法和系统 |
-
2016
- 2016-12-16 CN CN201611189755.6A patent/CN106651799B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104380713A (zh) * | 2014-03-12 | 2015-02-25 | 深圳市大疆创新科技有限公司 | 数码影像的坏点矫正方法和系统 |
| CN105787925A (zh) * | 2016-02-03 | 2016-07-20 | 中国科学院空间应用工程与技术中心 | 推扫型光学遥感载荷原始图像坏线的自动检测方法和系统 |
Non-Patent Citations (2)
| Title |
|---|
| Detection and correction of abnormal picels in Hyperion images;T.Han et al.;《IEEE》;20021107;第1327-1330页 |
| EO．1 Hyperion高光谱数据坏线修复方法的对比分析;于颖 等;《东北林业大学学报》;20091130;第37卷(第11期);第123-124页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106651799A (zh) | 2017-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106023185B (zh) | 一种输电设备故障诊断方法 | |
| WO2017031716A1 (zh) | 自然场景图像中手写体数学公式结构分析与识别方法 | |
| Horé et al. | Is there a relationship between peak‐signal‐to‐noise ratio and structural similarity index measure? | |
| CN108960404B (zh) | 一种基于图像的人群计数方法及设备 | |
| CN112254656B (zh) | 一种基于结构表面点特征的立体视觉三维位移测量方法 | |
| CN107622239B (zh) | 一种层次化局部结构约束的遥感图像指定建筑区检测方法 | |
| CN107230197B (zh) | 基于卫星云图和rvm的热带气旋客观定强方法 | |
| CN110490082A (zh) | 一种有效融合神经网络特征的道路场景语义分割方法 | |
| CN105095896A (zh) | 一种基于查找表的图像畸变校正方法 | |
| CN111178392B (zh) | 基于深度神经网络的航空发动机孔探图像损伤分割方法 | |
| CN105092607A (zh) | 球面光学元件表面缺陷评价方法 | |
| CN104809724A (zh) | 多波段遥感影像的自动精配准方法 | |
| CN114494821B (zh) | 基于特征多尺度感知与自适应聚合的遥感影像云检测方法 | |
| CN114723750A (zh) | 基于改进yolox算法的输电线路耐张线夹缺陷检测方法 | |
| CN113326846A (zh) | 一种基于机器视觉的桥梁表观病害快速检测方法 | |
| Sun et al. | Analysis of non-local Euclidean medians and its improvement | |
| Mustafa | A modified Hamming distance measure for quick detection of dissimilar binary images | |
| CN111027542A (zh) | 一种基于Faster RCNN算法改进的目标检测方法 | |
| CN106651799B (zh) | 一种基于地物光谱属性的坏线修复方法 | |
| CN107392948B (zh) | 一种分振幅实时偏振成像系统的图像配准方法 | |
| CN112950562A (zh) | 一种基于线结构光的扣件检测算法 | |
| Wu et al. | Deep learning‐based super‐resolution with feature coordinators preservation for vision‐based measurement | |
| CN105354833A (zh) | 一种阴影检测的方法和装置 | |
| Zhang et al. | The modified gradient edge detection method for the color filter array image of the CMOS image sensor | |
| Kim et al. | Building footprints extraction of dense residential areas from Lidar data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for 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: 20190528 Termination date: 20191216 |