遥感里通过夜间灯光数据提取城市建成区步骤(保姆级教学)
目录
一、对夜间灯光数据进行预处理
<1>图像裁剪(在ENVI处理)
<2>量纲转换(在ENVI中处理)
<3>去噪处理(这里我们在ENVI中进行中值滤波处理)
<4>重投影(在Arcmap中进行投影栅格)
二、提取城市建成区的方法(在Arcmap中操作)
<1>对上述的数据进行按掩膜提取
<2>转为整型
<3>栅格转面
<4>添加“面积”字段
<5>融合
<6>表转EXCEL
<7>通过筛选提取城市建成区
一、对夜间灯光数据进行预处理
<1>图像裁剪(在ENVI处理)
由于获取到的是全球影像数据,为了节约时间、增加效率,我们首先进行了研究区域的裁剪。
<2>量纲转换(在ENVI中处理)
由于NPP_VIIRS数据单位量级非常小,为了避免数值显示不完全的问题,方便后续的操作,我们利用Band Math对数据进行量纲转换,这里我们将DN值放大10的6次方(根据大家需要)。
<3>去噪处理(这里我们在ENVI中进行中值滤波处理)
利用中值滤波进行平滑处理。中值滤波会选取数字图像或数字序列中像素点及其周围临近像素点的像素值,将这些像素值排序,然后将位于中间位置的像素值作为当前像素点的像素值,让周围的像素值接近真实值,从而消除孤立的噪声点。通过相关研究, 在进行过滤时选用5×5的中值滤波栅格模块处理结果最优且对城市灯光影响较小。
<4>重投影(在Arcmap中进行投影栅格)
为使数据更精确,避免网格变形带来的影响以及便于计算影像中灯光亮值像元的面积,通常将夜光遥感影像数据转换成兰伯特方位角等积投影(Asia_Lambert_Conformal_Conic),该投影比较适合中国地形。同时我们将将栅格大小设置成10×10。
二、提取城市建成区的方法(在Arcmap中操作)
<1>对上述的数据进行按掩膜提取
得到上海市的夜间灯光数据,必须做这步。
<2>转为整型
通过截断将栅格的每个像元值转换为整型,防止后期报错
<3>栅格转面
将栅格数据集转换为面要素。字段一般默认为Value字段。
<4>添加“面积”字段
添加“面积”字段(默认即可)后通过计算几何自动计算出即可,注意单位选择。
<5>融合
融合字段选择亮度值的字段,统计字段选择面积,统计类型选择SUM。
<6>表转EXCEL
在EXCEL表中通过对GRIDCODE倒序(从大到小)排列,累加面积大于或等于建成区面积的亮度值即为阈值。表中的阈值即为29824596。
<7>通过筛选提取城市建成区
通过SQL表达式选择大于或等于GRIDCODE的像元即可。
最终即可提取出城市建成区。
效果图如下:
遥感里通过夜间灯光数据提取城市建成区步骤(保姆级教学)相关推荐
- 一、DMSP/OLS、NPP/VIIRS夜间灯光数据之城市建成区提取之理论介绍
图1 环渤海城市群城市建成区扩张图 一.前言 对于上述图1,很多小伙伴可能见过类似的,那么这种图怎么来通过GIS或者通过其他软件实现.早期由于遥感影像种类有限,大部分城市建成区范围通过Landsat等 ...
- 二、DMSP/OLS、NPP/VIIRS夜间灯光数据之建成区提取——阈值确定
一.DMSP/OLS夜间灯光数据之建成区提取--理论https://mp.csdn.net/mp_blog/creation/editor/112169085 一.前言 前文提到基于DMSP/OLS夜 ...
- 十三、基于DMSP/OLS、NPP/VIIRS夜间灯光数据和兴趣点(POI)城市建成区提取——理论
一.前言 夜间灯光数据是卫星传感器对夜间建筑物,道路和交通工具等发出灯光的记录[1,被广泛应用于城市建成区提取中.利用夜间灯光数据提取城市的建成区,能在很大程度上避免传统多光谱遥感的光谱混淆 .但由于 ...
- 二、DMSP/OLS、NPP/VIIRS等夜间灯光数据能源碳排放空间化——灯光指数计算(一)
一.前言 图1 某地区能源碳排放空间化分布图 其实能源碳排放数据直接找[中国能源统计年鉴]就可以,这里与GDP空间化有一点不同的就是,就是在计算夜间灯光指数之前,需要将研究区域所属年份建成区范围提取出 ...
- 十四、基于DMSP/OLS、NPP/VIIRS夜间灯光数据和兴趣点(POI)城市建成区提取——POI数据建成区提取
一.前言 POI核密度估计法建成区提取核心思路就是利用POI核密度分析,然后绘制出整个区域核密度等值线,定义Sd为以核密度值为d的等值线包围的环状曲线内部的面积,定义Sd开根号为Sd的理论半径,得到这 ...
- 三、DMSP/OLS、NPP/VIIRS夜间灯光数据之建成区提取
一.前言 在正式开始提取之前,还是要继续啰嗦一下,其实基于夜间灯光数据建成区提取,可以分为好几个层次,基于县一级.市一级.省一级,很多论文的做法都是以市一级为准,也就是统计好市一级城市的每一年的统计年 ...
- 三、DMSP/OLS、NPP/VIIRS夜间灯光数据之建成区提取——阈值确定(2)
https://blog.csdn.net/weixin_44725365/article/details/112169151https://blog.csdn.net/weixin_44725365 ...
- Python 分析中国城市夜间灯光数据
在城市化(城镇扩展).环境和能源等社会科学研究中常使用夜间灯光影像作为人类活动的表征.研究中使用最多的夜间灯光数据来自美国国防气象卫星计划(Defense Meteorological Satelli ...
- Python读取夜间灯光数据看中国城市发展
本文转载自气象学渣,详情可以扫描下方二维码: 在城市化(城镇扩展).环境和能源等社会科学研究中常使用夜间灯光影像作为人类活动的表征.研究中使用最多的夜间灯光数据来自美国国防气象卫星计划(Defense ...
最新文章
- linux硬盘拷贝系统,教你在第二块硬盘上拷贝安装Linux系统
- 高效搭建Storm全然分布式集群
- 查看SAP CRM和C4C的UI technical信息
- IP地址的网络位与主机位
- 10_python基础—数据序列(推导式)
- 数据结构思维 第一章 接口
- 医药/医疗/互联网医疗平台/问诊/挂号/开药/处方/复诊/患者管理/开药问诊/视频问诊/电话问诊/图文问诊/医生端/医师认证/医院/药品/续方/常用处方/电子处方/抢单/接诊/退诊/预约/科室/开方
- 通过的镜像源安装python包
- MyMathLib系列(行列式计算2)
- ASP.NET 2.0 中的SqlCacheDependency特性
- 动态链接库的设计(DLL)
- ActiveX控件进行注册
- 罗永浩演讲有趣的PPT
- 1.机器学习入门-初学者的机器学习
- win7——win server 2012 iis中使用asp程序出现Microsoft OLE DB Provider for ODBC Drivers 错误 '80004005'
- 微信或QQ屏蔽域名,爆红域名如何在微信打开,如何进行微信域名防封?
- 认识kubernetes(k8s),k8s单节点etcd部署
- 画图软件visio安装
- 怎么把移动硬盘挂载到linux,linux 怎么挂载移动硬盘
- 【Codeforces】764A Compote
热门文章
- Spring中的RowMapper接口的作用
- android view 虚线分割线,RecyclerView设置分割线---DividerItemDecoration
- 使用chkrootkit检测rootkit
- 九层妖塔SCALA单层循环实现
- 阿里云负载均衡SLB网络型NLB负载均衡架构性能详解
- 韩沉船事件幸存者忆惊魂时刻:船就像在喝水
- Mysql中长连接和短连接的区别
- 刷题 | top k问题
- 我的世界宅霸服务器文件在哪里,宅霸我的世界怎么联机?我的世界联机攻略
- oracle数据库报错12154,PL/SQL登录Oracle数据库报错ORA-12154:TNS:无法解析指定的连接标识符解决方法...