IDL代码实现湖泊水体范围遥感提取
1 基本原理与方法
最 常 用 的 水 体 提 取 方 法 是 归 一 化 植 被 指 数 法 NDVI，计 算 公 式 为（Rnir-Rred）/（Rnir+Rred）。有许多研究使用了 NDWI 或它的改进形式，用以在遥感影像上区分陆地和水体。
2 水体提取方法
在辐射传输模拟的基础上，Hu（2009）提出了一种浮藻指数 FAI（Floating Algae Index）。利用 MODIS 遥感影像获取 FAI的计算公式如下：
FAI=Rrc.859 - R'rc.859. (1)
R'rc.859=Rrc.645+（Rrc.1 240-Rrc.645）（859-645）/（1 240-645）(2)
在公式(1),(2)中，数字为 MODIS 的各个波段的中心波长,500 m 分辨率的 1 240 nm 波段先需要用锐化的方法重采样到 500 m（跟 1 真彩色合成影像类似）。因为气溶胶反射率在 645～1 240 nm之间，随着波长的增加呈近似线性衰减趋势，基线减法可以视作是一种有效的大气校正方法。
3 主要代码 ：
；645，859 和 1240 三个波段数据读取 read_hdf , hdf_file=MODIS_file,tag_name='CorrRefl_01' ,data=Ref_645read_hdf,hdf_file=MODIS_file ,tag_name='CorrRefl_02' ，data=Ref_859 read_hdf,hdf_file=MODIS_file,tag_name='CorrRefl_03',data=Ref_1240 width=（size（Ref_645))[1] & length=（size（Ref_645))[2]；梯度 gradient 计算 radius=1 strucElem=SHIFT（DIST（2*radius+1）,radius，radius） LE radius morphImg=MORPH_GRADIENT（FAI_img,strucElem）compute_hist，morphimg，1，hist_x，hist_y ；水体范围提取 hist_x=hist_x［idx］ hist_y=hist_y［idx］ Water_idx=where（FAI_imgltThreshold，count） water_area=bytarr（width，length） water_area［idx］=255 ；结果以 png 格式输出 write_png,’water_area.png’,water_area，/order 4 精度验证
HJ-1A/1B CCD 数据具有高空间分辨率（30 m），可以将其提取的水体范围视作近似“真实值”来验证 MODIS 数据的提取结果。HJCCD 的水体提取方法与 MODIS 类似，但由于其没有短波红外的设置，无法估算 FAI，因此此处选用 NDVI 和梯度的方法来提取水边界线。对相同日期 MODIS 和 HJCCD 影像提取的结果进行了叠加分析。HJ-1A/1B CCD 数据较高的空间分辨率（30 m）使其获取的结果能表达更多的细节特征。
更多遥感和GIS知识，关注我的个人微信公众号：遥感加油站
IDL代码实现湖泊水体范围遥感提取相关推荐
- 免疫性属于计算机病毒特征码,基于免疫和代码重定位的计算机病毒特征码提取与检测方法.doc...
基于免疫和代码重定位的计算机病毒特征码提取与检测方法 基于免疫和代码重定位的计算机病毒特征码提取与检测方法 张瑜 LIU Qing-zhong 宋丽萍 罗自强 曹均阔 海南师范大学信息学院 Depar ...
- lego-loam代码分析（1）-地面提取和点云类聚
lego-loam代码分析(1)-地面提取和点云类聚 概述 imageProjecion.cpp 获取点云角度范围 (findStartEndAngle()) 无序的点云变为有序(projectPoi ...
- 【IDL代码模块】提取文件名中的部分字符串
先知道这些: 1.file_basename():不要文件路径(d:/-/),只要文件名本身(如mod.20220305.hdf). 2.strmid((Expression, First_Chara ...
- 代码实战 | 用LeGO-LOAM实现地面提取
编辑丨计算机视觉life 作者介绍:Zach,移动机器人从业者,热爱移动机器人行业,立志于科技助力美好生活. LeGO-LOAM框架设计思路的第一步就是提取并分离地面.本篇文章就来详细说明LeGO-L ...
- 基于MODIS数据的NDVI与LST相关性分析（IDL代码实现）
1 数据预处理 (1)数据提取 我们可以选取2018年5月初华东地区MODIS中的MOD11A2和MOD13A2的16天合成LST和NDVI产品数据,下载地址:MODIS数据下载 网站下载数据需要注册 ...
- matlab人体轮廓识别,人体动作识别代码，求问这是提取的什么轮廓信息
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 目地是进行视频中人体动作的识别,这个代码好像是计算人体边界轮廓特征参数,与模版进行比对. 我知道里面有面积 周长 宽 高 ,可是那个旋转计算交点个数那个参 ...
- 基于高分辨率影像的潍坊大棚遥感提取
1 背景 朋友,当你坐在经过潍坊的火车上的时候,是否被窗外的风景震惊过?那一望无际的大棚,像一片海洋,又像一片草原,连接天际,一眼望不到边.在卫星地图上,白茫茫一片,蔚为壮观. 本文基于0.6米的高分 ...
- 使用代码从抓包文件中提取H264码流
1 从抓包文件中过滤出包含单条流的RTP包 使用ssrc 或者 payload type过滤皆可,具体如下图: 2 使用wireshark可以提前分析下本条码流到底有没有丢包 这样就可以预先知道这条流 ...
- PHP代码去快手视频水印实现提取快手水印视频播放地址
原理分析 从快手App拿到视频分享链接在浏览器中访问,比如地址 https://f.kuaishou.com/1mqmxe 然后会重定向,打开浏览器控制台点击左上角的手机图片,刷新模拟手机请求,观察控 ...
- 遥感提取土地利用（从下载到ENVI到ACRGIS 【简单粗暴版】）
一.影像下载 下载地址:地理空间数据云http://www.gscloud.cn/search,可以用邮箱注册一个账户. 选择影像:数据集根据自己所需要的影像精度选择 影像选择注意点: 可选择的数据集 ...
最新文章
- 终端(terminal)、tty、shell、控制台（console）、bash之间的区别与联系
- python不能处理excel文件-python处理Excel文件
- shell中if 变量里包含字符串的判断
- 工程图样中粗实线的用途_图纸天天画，粗实线和细实线的线宽比例是多少？2:1还是3:1？...
- linux配置sftp-server,Ubuntu Server如何配置SFTP（建立用户监狱）
- What is the difference between LINQ to SQL and LINQ to Entities?
- ubintu yum装mysql_Ubuntu 下安装Mysql centos8 下安装mysql8.0
- 事件分发机制 事件拦截 滑动冲突 MD
- TopOn的两种测试方法
- java单例的生命周期_Spring Bean生命周期详解
- office未响应（2010、2013……）崩溃
- 8K HDR!｜为 Chromium 实现 HEVC 硬解 - 原理/实测指南
- 卸载wps后安装office图标异常的问题
- 感谢以色列，站在对抗邪恶的最前线
- 广东省工科赛 智能终端配送机器人：创作心路历程，踩过的坑和解决方案
- vue样式初始化_前端Vue项目——初始化及导航栏
- 最新计算机台式机小机箱,几款热门台式机主机机箱的参数配置【图文】
- 鸿蒙系统网络连接设置ip,网络ip地址错误解决方法
- 买学生台灯应该怎么选择？挑选国AA级照度的护眼台灯
- 你离拍出好的雪景作品，还有多远？ @教摄影