用ChatGPT结合GIS计算植被归一化指数NDVI出图教程
用ENVI计算比较繁琐,如今AI的盛行,我们可以轻松解决计算问题,只需1一分钟变可以出图。
详细教学请看上方视频步骤。
更多ChatGPT教学内容请见:ChatGPT结合GIS:一分钟完成植被归一化指数NDVI计算并出图
----
以下是我们通常用ENVI计算NDVI的流程和方法,共参考。
需要准备2个文件:1、植被影像;2、研究区域的边界图层。(如果你影像和边界都有,就直接看第三步)
一、准备工作(2个文件的准备)
1、首先植被影像,我们在"地理空间数据云"下载
看着箭头进行操作,选择Landsat8OLI_TIRS卫星数字产品——行政区——青海省(研究区域所在的省份)——海西(研究区域所在的市)(如果你研究的区域是一个市,就选择在这里就结束了,后面就不用选择了)——都兰县(如果是一个县或者区你就往下再选择)(如果你研究的是某个保护区或者更小的地方,看你研究的区域被包含在哪个区或者县下载就好了,后面再掩膜就行)——时间范围(最好选择你研究需要的时间范围)——月份(选择你研究区域植被茂盛的时候就行,一般78月份都行;但如果你研究的是冬季下雪植被的覆盖那肯定选下雪的日子嘛那就是1月丫对吧)——云量选择5——检索
1.1检索出来之后的界面:红色框里面的箭头,可以看到一个是灰色,一个是深色的,我们选择深色的,灰色是不能下载的;选中深色箭头,你据可以看到图里面绿色框起来的部分,你要在检索的结果里面选中能全部覆盖住你研究区域的图层,他可能是一个图层就包含了你的研究区域,也可能是好几个,我这里下载的天津市区的他就需要4张才能完全覆盖演剧区域,那么我们下载这4张就好。(我用天津给你们演示的)
OK我们的影像就准备好了。
2、研究区域边界下载
2.1打开全国乡镇边界矢量数据——在GIS中提取你的研究区域步骤如下:
首先打开gis——加载进去“全国乡镇边界矢量数据”/县级边界矢量图(看你研究的区域是县还是镇了),我这里提取县的,我就选“全国县级边界矢量图”
选中右键单击打开属性表——选中你的目标区域就可以了
选中右击,数据——导出数据——确定——是——关掉其他的图层显示我们要的图层就好啦,我命名为“都兰”
OK到这里我们研究区域边界就下载好啦
影像和边界我们都有啦,我们开始今天的重头戏,开始计算NDVI
三、计算详细步骤
Tips:如果你的研究区域区域如果涉及多景,后面就需要进行镶嵌哈。但是你不涉及的话,只有1景那么恭喜你太爽啦,我还是给大家讲一下涉及多景应该怎么弄。
多景(就是你下载了好几块):先辐射定标——大气校正(可做可不做,最好做,他们都这么说,但我这次选择不做哈哈哈哈)——镶嵌——掩膜——提取NDVI——(计算FVC:植被覆盖度,大家不需要这个就不用看计算FVC这一步哈)。
1景(你下载的只有1块):辐射定标——大气校正——掩膜——提取NDVI——计算FVC
或者掩膜——辐射定标——大气校正——提取NDVI——(计算FVC)
3.1 打开ENVI这个软件,界面就长下面那样。
3.2 加载进去我之前1.1让大家下载的影像;步骤:File——open as——optical sensors——Landsat——GEO-Tif 加载进去我们下载的影像,注意加载进去的是:MTL.txt格式的文件,你有几个就加载几个。
加载进去之后是这样的,我就加载了两景(你们有多少景加载多少景)
3.3 辐射定标
(OK接下来先进行辐射定标,如果你是多景一定要先进行辐射定标,要不然你先镶嵌完再辐射定标是不可以的,报错会告诉你图层的什么波段啥的不对什么什么的,我已经踩过坑了,你们避雷一下)
输入:rad——选中radiom,,,双击它:出现下面这个界面,选中红框里面框起来的那种格式的,选中——OK
好的,我们将刚才辐射定标完成的图命名为fsdb1.dat(你有多景的,你要重复操作刚才的过程,将你的每一景都要辐射定标)
辐射定标结束,我得到了fsdb1和fsdb2两个图层
3.4 镶嵌
在工具框里面输入:sea——选中seamless mosaic
出现下面的界面:(发现旁边东西没了是不是,不要着急正常的啦)
接下来如图操作:
加载完毕之后,如图操作:
不要点finish!!!我刚在图上箭头多指了一步,我已经截掉了;
在接着继续设置,如图操作:
接着等他搞完,我们就导出:
然后就是漫长的等待啦,镶嵌这一步我们就完成啦
3.5 掩膜
也就是裁剪啦,将镶嵌好的图层保留,其他的就可以移除啦,OK我们再加载进去研究区域,也就是2.1里的下载的边界啦
如图:
我们就加载进来啦,开始掩膜,如图操作:
就掩膜成功啦
3.6 提取NDVI
如图操作:(注意:里面输入公式的时候一定要在英文的状态下哈)
良心发现我还是把公式放这里吧:(float(b1)-b2)/(b1+b2)
float表示浮点,b1表示近红外,b2表示远红外
NDVI=(b1)-b2)/(b1+b2)
OK,NDVI我们就成功啦
4、植被覆盖度计算(大家需要我再出,不需要我就不出啦)(PS:本来是打算继续讲完这个的,但是我要干饭去啦,干饭人干饭去咯,大家需要就留言,后面可以补)