Python读取夜间灯光数据看中国城市发展
本文转载自气象学渣,详情可以扫描下方二维码:
在城市化(城镇扩展)、环境和能源等社会科学研究中常使用夜间灯光影像作为人类活动的表征。研究中使用最多的夜间灯光数据来自美国国防气象卫星计划(Defense Meteorological Satellite Program,DMSP)一系列气象卫星的观测。处理夜间灯光数据的方法和软件有很多,常用的有Arcgis等。在这里小渣将尝试用Python来读取夜间灯光数据。
一、夜间灯光数据下载
DMSP/OLS夜间灯光数据下载地址:
https://www.ngdc.noaa.gov/eog/dmsp/downloadV4composites.html
小渣以2013年为例,下载后得到压缩文件F182013.v4.tar,解压之后:
目前发布的灯光数据产品中使用较多的是稳定灯光数据(stable lights),它剔除了由火灾、爆炸等造成的短暂亮光,只包含了相对稳定的灯光,且均为无云数据。让我们再次解压稳定灯光数据:F182013.v4c_web.stable_lights.avg_vis.tif.gz,获得的是1张标签图像文件格式(tiff)的图片,我们需要的灯光数据就在这张图里啦。
1、安装libtiff库
小渣的实践经验表明,直接用pip/conda install在线安装libtiff库可能会报错(具体原因不明)。建议先下载好本机适配的libtiff-0.4.2-xxxx-xxxx-xxxxx.whl文件,采用pip install 路径+libtiff库的离线方式安装最靠谱。
下载地址:https://www.lfd.uci.edu/~gohlke/pythonlibs/
安装命令:pip install 路径+libtiff-0.4.2-xxxx-xxxx-xxxxx.whl
2、Python提取指定区域的夜间灯光
安装好Python的libtiff库后,接下来用Python提取并显示中国东部的夜间灯光,附上源代码(在这里小渣不禁想感叹一下Python的简洁与高效):
from libtiff import TIFF
import matplotlib.pyplot as plt# 读取夜间灯光数据
tif = TIFF.open('F182013.v4c_web.stable_lights.avg_vis.tif', mode='r')
img = tif.read_image()#提取包含中国东部范围的夜间灯光数据画图
#夜灯数据覆盖范围为-65~75oN,-180-180oN,分辨率是30’
lons=100;lone=137;lats=15;late=52
lons_grid=int((lons+180.0)/(30.0/3600))
lone_grid=int((lone+180.0)/(30.0/3600))
lats_grid=int((75.0-lats)/(30.0/3600))
late_grid=int((75.0-late)/(30.0/3600))
img2=img[late_grid:lats_grid,lons_grid:lone_grid]
plt.imshow(img2)
plt.show()
采用我们上期介绍的绘制动图的方法(戳这里复习:
◆ ◆ ◆ ◆ ◆
长按二维码关注我们
数据森麟公众号的交流群已经建立,许多小伙伴已经加入其中,感谢大家的支持。大家可以在群里交流关于数据分析&数据挖掘的相关内容,还没有加入的小伙伴可以扫描下方管理员二维码,进群前一定要关注公众号奥,关注后让管理员帮忙拉进群,期待大家的加入。
管理员二维码:
Python读取夜间灯光数据看中国城市发展相关推荐
- Python 分析中国城市夜间灯光数据
在城市化(城镇扩展).环境和能源等社会科学研究中常使用夜间灯光影像作为人类活动的表征.研究中使用最多的夜间灯光数据来自美国国防气象卫星计划(Defense Meteorological Satelli ...
- 夜间灯光数据dn值_中国区域DMSP/OLS夜间灯光影像的校正
DMSP/OLS夜间灯光影像由于传感器未经过星上定标,导致长时间序列的各期影像之间缺乏连续性和可比性,且由于光谱分辨率的限制,使得大量的城市中心像元出现DN值为63的集聚饱和现象等问题,因此学者们对以 ...
- 一、DMSP/OLS、NPP/VIIRS等夜间灯光数据之GDP空间化的理论介绍
一.前言 自上世纪90年代以来,DMSP/OLS夜间灯光数据问世,由于和其他遥感数据不同,该数据可以探测到地球夜间微弱的灯光.随着数据的积累,可以获取长时间序列的数据之后,随后常用于人类社会经济指标的 ...
- DMSP夜间灯光数据
数据和详细信息参见https://ngdc.noaa.gov/eog/dmsp/dmsp.html) 1.美国国防气象卫星计划(Defense Meteorological Satellite Pro ...
- 遥感里通过夜间灯光数据提取城市建成区步骤(保姆级教学)
目录 一.对夜间灯光数据进行预处理 <1>图像裁剪(在ENVI处理) <2>量纲转换(在ENVI中处理) <3>去噪处理(这里我们在ENVI中进行中值滤波处理) & ...
- 三、DMSP/OLS夜间灯光数据校正一相互校正、饱和校正
DMSP由于多代传感器,每一代传感器之间存在重叠年份,且未完成星上辐射校正,DN值存在异常情况.对其的预处理主要解决年份重叠.DN值异常.灯光饱和等问题,要完成相互校正.连续性校正.饱和校正. 当从N ...
- 二、DMSP/OLS夜间灯光数据校正一理论介绍
夜间灯光数据介绍: 1.美国国防气象卫星计划(Defense Meteorological Satellite Program,DMSP)由美国空军航天与导弹系统中心运作,卫星运行的线性扫描系统(Op ...
- 二、DMSP/OLS、NPP/VIIRS等夜间灯光数据能源碳排放空间化——灯光指数计算(一)
一.前言 图1 某地区能源碳排放空间化分布图 其实能源碳排放数据直接找[中国能源统计年鉴]就可以,这里与GDP空间化有一点不同的就是,就是在计算夜间灯光指数之前,需要将研究区域所属年份建成区范围提取出 ...
- 四、NPP/VIIRS夜间灯光数据校正之连续性校正
三.NPP/VIRRS夜间灯光数据校正之去除负值和异常值https://blog.csdn.net/weixin_44725365/article/details/123782077一.前言 上文对N ...
最新文章
- 【转】登录后自动显示未读快件
- JDK源码分析:Byte.java
- Linux(centos7)安装nginx并配置代理
- 20175305张天钰 《java程序设计》第四周课下测试总结
- Java之static的内容
- 选择多级分类_商用车齿轮油的选择与运用!!
- 判断两个时间段是否有交集_判断两个人是否处于暧昧关系,就看四点,特明显...
- 小米平板2wifi驱动下载_小米平板2 win10驱动
- 挖金子修订版开发进行中
- Elasticsearch(三)使用 Kibana 操作 ES
- ubuntu创建桌面快捷方式
- mysql非常规的中文乱码问题
- 在小百合注册了一个帐号
- 重装系统要花钱?教你免费一键重装win10系统
- ESP通过乐为物联控制灯,微信发送数值,ESP上传传感器数据
- 过来人的经验:自学编程的6个建议
- Matlab 斜率和曲率,曲率_与闪电共舞_新浪博客
- AI ProCon圆满落幕,五大技术专场精彩瞬间不容错过
- python学习⑥|基于tesseract的视频字幕提取
- python P51-60