如何在GRASS中全景显示LANDSAT 8?


9

使用QGIS已有相当长的时间进行地图绘制后,我被迫深入研究GRASS,因为全景投影不能直接在QGIS中执行。在大致了解结构之后,我对i.fusion.brovey进行了一些试验。我能够生成15m分辨率的栅格,但是当加载到QGIS中时,我无法实现/恢复原始30m栅格中存在的良好的颜色对比度。

较早的帖子中,有人证实了这一点,并建议对landsat使用另一个功能:“ i.his.rgb”(在使用i.landsat.toar之后,我没有使用)。在使用i.his.rgb之前,他开始将最小值/最大值转换为1-255。我的原始栅格的范围是0-65535。它对我没有用...

经过几天努力奋斗只是为了进行全景拍摄,我退出了,我感到被困住了...为什么仅仅应用i.pansharp不能像其他人建议的那样对Landsat起作用?有人可以建议一个更简单的工作流程来执行Landsat 8图像的全清晰度吗?期待一些建议,在此先感谢...

Answers:


7

看一下 i.landsat.rgb-可能在运行全景图像之前对LANDSAT图像执行颜色自动平衡。您也可以考虑使用i.landsat.toar将各个频道的数字转换为大气顶辐射或反射率。另请参见http://grasswiki.osgeo.org/wiki/LANDSAT

顺便说一句:原始LANDSAT数据的范围是0-65535,这很奇怪,因为它们是作为8位通道(0-255)交付的。因此,请首先检查您的输入数据。

更新:10/2013

有关使用当前GRASS GIS 7的简化过程,请参见GRASS GIS 7中的处理Landsat 8数据在此处输入图片说明:RGB复合和泛锐化


我最初从neteler.org上看到您要链接的那篇文章,但这些说明根本不起作用(由于使用生成的栅格的浮点格式i.landsat.toar,因此我最终遵循以下建议:gis.stackexchange.com / a / 39239/9857
胡安

3
请注意,该建议是多年以前提出的,并不反映最新的GRASS GIS代码。顺便说一句:新的附加组件是一个非常出色的全锐化模块:i.fusion.hpf-基于高通滤波器附加技术融合高分辨率全色和低分辨率多光谱数据,请参阅grass.osgeo.org/grass70/manuals/ addons / i.fusion.hpf.html
markusN 2013年

4

感谢Michal和Markus。最终,我可以使用所示的功能brov进行全景拍摄。在GRASS中执行完之后,我将RGB栅格导出为数据类型为Uint16的gdal GTiff。马库斯(Markus),我一直在下载和使用数张Landsat 8图像(11条带中的每条带的色差),以自然/假色获得不错的视觉效果。据我所知,原始的Landat Tiff是Uint16,其值介于0到65535之间(或者我错过了什么?)...感谢GRASS LANDSAT链接...


3
是的,Landsat 8图像存储为Uint16类型。
webrian

3

如果遵循以下文章中的公式,则可以使FIHS或Braovery转换变得非常简单:http : //www.sensorsportal.com/HTML/DIGEST/august_2014/Vol_177/P_RP_0183.pdf

假设您正在使用GRASS GIS并使Landsat 8锐化。

使用平移8 tif设置您的区域-重要的是要设置15m的分辨率,如果执行此操作将被设置。

对于没有权重的Brovery,分别将其输入到“地图计算器”中分别用于蓝色,绿色和红色:

band8 /(band2 + band3 + band4)/ 3 * band2

band8 /(band2 + band3 + band4)/ 3 * band3

band8 /(band2 + band3 + band4)/ 3 * band4

对上述每个输出运行r.colors并将直方图均衡和颜色表设置为灰色。

运行r.composite,选择适当的波段(红色= 4,绿色= 3,蓝色= 2)

检出图像。

显然,如上文所述,可以在考虑大气校正,颜色拉伸和色带权重的情况下完成此过程。



2

在解决了这个问题之后,我找到了一个简单的解决方案,可以让您使用SCP在QGIS(2.18.21是我的版本)中进行快速处理。https://fromgistors.blogspot.com/2015/07/major-update-semi-automatic-44.html上介绍了采用此方法的SCP更新。

这很容易做到,而且就像魅力一样!我检查了最终输出,像素的分辨率为15 m,与全色8波段匹配。如果有人好奇-它会按照手册中的描述使用Brovey变换方法(请参阅https://media.readthedocs.org/pdf/semiautomaticclassificationmanual/latest/semiautomaticclassificationmanual.pdf)。

SCP Pansharpen工具的步骤1 SCP Pansharpen的步骤3-4



By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.