envi单波段阈值法提取水体_等吸收双波长法

envi单波段阈值法提取水体_等吸收双波长法水体提取方法之单波段阈值法升级版,利用MNDWI,NDWI,INDVI结合分析,单波段提取水体

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE稳定放心使用

本操作执行环境:Windows系统下Erdas 2015。

操作对象:巴尔喀什湖部分行列号影像,此教程使用行列号为:2017年3月28日行列号为152028影像(巴尔喀什湖:地处哈萨克斯坦共和国的东南部,在一个大型盆地内,主要有伊犁河流入。在世界众多的湖泊中,它因湖水一半为咸水湖一半是淡水湖、东段与西段颜色不同而独具特色。1864年中国清朝和俄罗斯帝国签订不平等的《勘分西北界约记》以后,巴尔喀什湖正式脱离中国。苏联时期在哈萨克斯坦加盟共和国境内,苏联解体以后就在哈萨克斯坦共和国。)。

1:影像如何下载?

大家应该都知道的地理空间数据云,然而地理空间数据云影像绝大多数为国内遥感影像,那么如何下载巴尔喀什湖影像尼?具体操作过程见我的个人网站http://www.shaofee.com/archives/375.html,或者我的博客:http://blog.csdn.net/qq_36305327/article/details/71154996

2:影像预处理

一般从地理空间空间数据云或者其他网站一般都是处理过的影像,所以我们不需要做太多的处理。

正常的预处理过程包括:(一)几何精校正与图像配准(二) 图像融合(三)图像镶嵌与裁剪(四)大气校正,在这里我就不多说了,对于区域或者精度要求比较高的朋友,这些过程必不可少。

在这里我需要对下载影像进行融合一下,以便后续操作。

envi单波段阈值法提取水体_等吸收双波长法

landsat8 OIL有很多波段,但是我们做水体提取,并不需要那么多。这里我们对前7个波段进行合成,选择影像时候注意按照顺序输入。合成后影像为1520282017mix.img.

envi单波段阈值法提取水体_等吸收双波长法

3:利用水体指数对水体信息进行增强。

NDWI=(绿光-近红外)/(绿光+近红外),

MNDWI=(绿光-中红外)/(绿光+中红外),

INDVI=(红光–近红外)/(红光+近红外),

注意:这里的绿光、红光、近红外和中红外波段,分别一一对应于OLI影像中的第3、4、5、6波段;对于TM(或ETM+)影像,则分别一一对应于影像中的第2、3、4、5波段;

这里以NDWI为例。

下面我们以MNDWI为举例:

打开Erdas2015的Model maker,在模板里寻找NDVI模板,对它进行改写。

envi单波段阈值法提取水体_等吸收双波长法

如果你使用的不是OLI影像则不要要改写3处,只需要改写1,2处就可以。

输入影像

1处由:$n1_1520282017mix(4) – $n1_1520282017mix(3) 改写为:$n1_1520282017mix(3) – $n1_1520282017mix(6)

2处由:$n1_1520282017mix(4) +$n1_1520282017mix(3) 改写为:$n1_1520282017mix(3) +$n1_1520282017mix(6)

3处由:($n11_memory – GLOBAL MIN ( $n11_memory )) / ( GLOBAL MAX ( $n11_memory ) – GLOBAL MIN ( $n11_memory )) * 255 改写为 $n11_memory * 10000+100001。

为什么要拉伸?因为正常MNDWI,NDWI,INDVI的值位于(-1,1),这里3处主要对其进行拉伸为后面我们进行阈值提取做准备。

上述NDWI,INDVI,计算过程略。过程同MNDWI。

4:如何采集阈值?

在操作之前我们利用将源波段3,4,5与MNDWI,NDWI,INDWI波段输入顺序进行融合合成后为1520282017.img。融合过程同2:影像预处理中合成过程。

下面我们来确定水体阈值。下图为影像在1520282017.img ,3,2,1波段下显示效果,即是相当于1520282017mix.img的5,4,3波段下效果图。

envi单波段阈值法提取水体_等吸收双波长法

在1520282017.img的视图下,自由调节波段组合顺序,直到你自己可以分清水体与其他地物。

例如:envi单波段阈值法提取水体_等吸收双波长法

下面打开Spectral Profile工具进行集点。

envi单波段阈值法提取水体_等吸收双波长法

envi单波段阈值法提取水体_等吸收双波长法

当水体采集点完成后采集,其他地物点。这里修改下线条颜色。我们将水体修改为红色。

envi单波段阈值法提取水体_等吸收双波长法

envi单波段阈值法提取水体_等吸收双波长法

选中所有点然后点击Color列,选择Color修改为红色。点击应用。如下图所示,再次点击应用。

envi单波段阈值法提取水体_等吸收双波长法

envi单波段阈值法提取水体_等吸收双波长法

envi单波段阈值法提取水体_等吸收双波长法

采集尽可能多的点,分不清陆地水体可以结合Google Erath。上图为采集后点,红色线为水体,绿色线为其他地物。

从表中我们可以看出在4波段,水体与其他地物有很明显的分界。那么我们就可以利用4号波段来进行阈值提取。

下面我们导出数据。点击右上角文件——- 导出数据。导出后的文件为.Sif格式。

利用文本编辑工具直接打开,全选复制粘贴到EXCEL中,下面对数据进行处理。

envi单波段阈值法提取水体_等吸收双波长法

选择4号波段,对数据进行整理。得到下图:

envi单波段阈值法提取水体_等吸收双波长法

下面得出阈值我们就可以通过Erdas2015的Model maker,进行处理,提取水体。

envi单波段阈值法提取水体_等吸收双波长法

通过此模板提取水体,采集的点数据比较多,提取效果十分好。下面附一张效果图。

envi单波段阈值法提取水体_等吸收双波长法

 

envi单波段阈值法提取水体_等吸收双波长法

 

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/179423.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • 从Java培训班出来后是怎么找到工作的?

    从Java培训班出来后是怎么找到工作的?经常听到有人抱怨,我参加过IT培训,公司一看我的简历就直接pass我了,把培训经历去掉后,面试通知是多了,但是往往面试完,就不了了之了,难道我参加过培训,公司就不会选择我了吗?一些企业有点排斥培训班出来的,他们排斥的原因相信大家都知道,主要是简历过度包装,面试都会一上手就废而且现在相当多的公司,会选择在培训机构培训过的程序员,很多公司会去到相关培训机构做双选会,直接挑选适合本公司职位的学员,即能快速了解此人的技术水平,还能很直观的了解到他的学习水平。没有一个企业喜欢学习能力差的人,你可以在进公司的时

    2022年7月7日
    21
  • android短信验证码代码,Android短信验证码自动填写实现代码

    android短信验证码代码,Android短信验证码自动填写实现代码今天给大家分享一个Android短信验证码自动填写的功能。先看下效果图,我发送了一条短信到手机,自动填写验证码。如图:这个小功能运用到了观察者模式,什么是观察者模式?观察者模式:定义对象间的一种一个(Subject)对多(Observer)的依赖关系,当一个对象的状态发送改变时,所以依赖于它的对象都得到通知并被自动更新。在本例中,我们在短信中注册一个观察者,当短信功能(被观察者)收到信息时,就会…

    2022年7月25日
    7
  • springboot+websocket实现服务端、客户端

    springboot+websocket实现服务端、客户端一、引言小编最近一直在使用springboot框架开发项目,毕竟现在很多公司都在采用此框架,之后小编也会陆续写关于springboot开发常用功能的文章。什么场景下会要使用到websocket的呢?websocket主要功能就是实现网络通讯,比如说最经典的客服聊天窗口、您有新的消息通知,或者是项目与项目之间的通讯,都可以采用websocket来实现。二、websocket介…

    2022年7月11日
    16
  • JVM – 内存模型

    JVM – 内存模型JVM-内存模型

    2022年5月24日
    30
  • QMap与QHash

    QMap与QHash转载:http://newfaction.net/2010/11/17/qt-qhash-and-qmap-difference.html关联容器可以保存任意多个具有相同类型的项,且它们由一个键索引。Qt提供两个主要的关联容器类:QMap<K,T>和QHash<K,T>。     QMap<K,T>是一个以升序键顺序存储键值对的数据结构…

    2022年5月29日
    31
  • linux时间戳转换成日期格式_数据库时间戳转换成时间

    linux时间戳转换成日期格式_数据库时间戳转换成时间//ConvertUnixtimestamptonormaldatestylepublicStringTimeStamp2Date(StringtimestampString){Longtimestamp=Long.parseLong(timestampString)*1000;Stringdate=newjava.text.S…

    2022年10月2日
    0

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号