应用开源软件去掉影像中的云
参考《Open Source Geospatial Tools Applications in Earth Observation》—Chapter 17 Case Study on Multispectral Land Cover Classification
Depending on the time of the year and the location on Earth, the satellite image
scenes can contain some clouds. The images used within this case study contain
some cloud cover. Clouds directly affect image quality and can severely hamper the
image analysis and classification. To obtain a scene with minimal clouds, we combine
the downloaded scenes and create a “cloud free” composite. If you are using your own
satellite imagery acquired for another location and it is cloud free than you can skip
the next processing step. The BQA dataset is in a bit-packed format2 and contains
information on water and clouds. This means we need to decode the information
in order to obtain, for instance, a cloud mask. More information on the bit-packing
encoding scheme is explained in Sect. 14.3, where we also implement a tool to create
a cloud mask from the BQA dataset. Alternatively, you can use the official tool, the
L-LDOPE Toolbelt, that is provided by the USGS.3 As an example, we create a cloud
mask with the tool bqa2cloud from Sect. 14.3:
bqa2 cloud -i LC 82070232013160 LGN 00_ BQA .TIF -o LC 82070232013160 LGN 00_ CLD . TIF
bqa 2 cloud -i LC 82070232013192 LGN 00_ BQA . TIF -o LC 82070232013192 LGN 00_ CLD . TIF
We then stack all bands into a single multi-band raster dataset in GeoTIFF format
using the utility gdal_merge.py. The panchromatic band contains little extra
spectral information and the thermal bands have been acquired with a coarser spatial
resolution. Therefore, we do not include bands 8, 10 and 11.
gdal_ merge.py -o LC82070232013160LGN00_ merge.tif -of GTiff -ps → 30 30 - separate -co INTERLEAVE = BAND -co COMPRESS = LZW → LC82070232013160LGN00_B1.TIF LC82070232013160 LGN00_ B2.TIF → LC82070232013160LGN00_B3.TIF LC82070232013160 LGN00_ B4.TIF → LC82070232013160LGN00_B5.TIF LC82070232013160 LGN00_ B6.TIF → LC82070232013160LGN00_B7.TIF LC82070232013160 LGN00_ B9.TIF → LC82070232013160LGN00_CLD.TIF
//
-o LC82070232013160LGN00_merge.tif Name of the output multi-band raster dataset. -of GTiff We explicitly set the GeoTIFF image format for the output (not required here as this is the default). -ps 30 30 Set a common spatial resolution for all inputs. Most bands are already at this resolution, but the panchromatic band needs to be resampled from 15 to 30 m. -separate Place each input file into a separate stacked band. -co INTERLEAVE=BAND Creation option for band sequential encoding (see Sect. 3.5). -co COMPRESS=LZW Creation option for LZW compressed output (see Sect. 5.5.4).
To create the cloud free composite in Fig. 17.2c based on the two multi-band
raster datasets in Fig. 17.2a, b, we can use the utility pkcomposite from pktools
(see Sect. 12.2).
pkcomposite -i LC82070232013160LGN00. tif -i → LC82070232013192LGN00.tif -o LC82070232013_ composite. tif → -bndnodata 8 -srcnodata 1 -dstnodata 0 -cr maxndvi -cb 3 -cb 4
//
-i LC82070232013160LGN00.tif -i ... Sequence of input raster datasets. -o LC82070232013_composite.tif Name of the output composite based on the two input datasets. myinline-bndnodata 8 Index of the cloud band defining nodata values in the input raster datasets (9th band as the index starts from 0). -srcnodata 1 Value indicating nodata, i.e. cloudy values in this case. -dstnodata 0 Value written in output raster dataset in case of nodata (i.e. if both scenes are cloudy). -cr maxndvi In case more than one input is cloud free, we select the pixel with the maximum NDVI value. -cb 3 -cb 4 The composite bands: the respective indices for the red and near infrared spectral bands in the multi-band input raster to calculate the NDVI on the fly (see Table 17.1)
相关文章
- 使用Spring MVC构建REST风格WEB应用
- .NET Core实战项目之CMS 第二章 入门篇-快速入门ASP.NET Core看这篇就够了 [WPF]为旧版本的应用添加触控支持 js css等静态文件版本控制,一处配置多处更新.net版【原创】 C#图片处理,缩略图制作 SqlDataReader 结果集 转成 DataTable DataTable转成List集合 如何使用linq读取DataTable集合?
- 第十四届蓝桥杯软件类省赛备赛分享帖(Web应用开发组)
- [译]深度学习:模仿人脑还是完善应用
- Android应用性能优化最佳实践.
- 《SAFe 4.0参考指南:精益软件与系统工程的规模化敏捷框架》一2.7 原则#7——应用节奏,通过跨领域计划进行同步
- web渗透测试----31、检测到绕过web应用防火墙
- 【快应用】原生广告下载状态监听案例
- 当今商业中使用的三种十分重要的IT应用系统
- RSA加密和数字签名在Java中常见应用【原创】
- iPhone应用程序开发基础之一: IBOutlet与IBAction
- 黑客教父龚蔚:扫码应用要警惕 公共WiFi攻不破支付软件
- 《Core Data应用开发实践指南》一2.3 实体
- 《OpenGL ES应用开发实践指南:Android卷》—— 1.4 创建Renderer类
- 《应用时间序列分析:R软件陪同》——1.2 时间序列例子
- 《应用时间序列分析:R软件陪同》——2.3 随机游走
- 《应用时间序列分析:R软件陪同》——2.4 趋势平稳过程
- 《应用时间序列分析:R软件陪同》——2.5 一般线性模型
- 《应用时间序列分析:R软件陪同》——2.11 习题
- 如何卸载和安装windows自带软件,内置应用,右键无法直接卸载
- LwIP应用开发笔记之一:LwIP无操作系统基本移植
- 应用Python处理空间关系数据
- 应用Canary文件类型阻击勒索软件
- 【SoC FPGA学习】六、基于虚拟地址映射的 UART 编程应用
- 用于IoT应用程序开发的10大开源软件
- UVa 12657 Boxes in a Line(应用双链表)
- Gartner:2017年应用和基础设施中间件软件市场收入将突破270亿美元
- Gartner:2017年应用和基础设施中间件软件市场收入将突破270亿美元
- 关于Web应用和容器的指纹收集以及自动化软件的制作