My Blog

欢迎访问我的博客

使用 MATLAB 绘制中国各省的煤炭产量和山西的煤矿分布

Plotting Coal Production in Chinese Provinces and Coal Mining Distribution in Shanxi Using MATLAB

这里记录并分享下我在写上一篇论文(Unveiling unprecedented methane hotspots in China’s leading coal production hub: A satellite mapping revelation)时作图的代码。效果如图 左边这个图上有两个图层,分别是各省的煤炭产量(来自能源局)及煤矿甲烷排放(来自EDGAR清单),并用了不同...

使用MATLAB保存NC文件

Save NC file Using MATLAB

很多卫星产品和气象模型喜欢使用nc文件类型,其优点是存储高效。NC文件的结构大致为变量和属性,属性又分为全局属性和变量属性。 一、常见函数的基本用法 创建NC文件需要三步 第一步:创建变量,定义这个变量的维度,以及类型(如double or float) 在名为 myncclassic.nc 的文件中创建一个名为 peaks 的二维变量。使用 Dimensions 名称-值参量指...

使用 MATLAB 绘制拟合线周围的 95% 置信区间

Plotting 95% Confidence Intervals Around Fitted Line Using MATLAB

一、绘图 读论文的时候发现一个好看的图(如下),之前知道这种图用R语言来画非常方便,这次想用MATLAB来实现一下。 clc;clear; x=[1 2 3 4 5 6 7 8 9 10]; % in_situ_emission_rate y=[11 13 15 14 17 14 18 16 19 20] - 10; % estimated_emission_rate e = 0...

个人学术主页搭建详细教程-基于HugoBlox搭建

一、概念和搭建思路 1.1 GitHub Pages、 Hugo 是什么 GitHub Pages 是一组静态网页集合(Static Web Page),这些静态网页由 GitHub 托管(host)和发布,所以是 GitHub + Pages。 Hugo 是用 Go 语言写的静态网站生成器(Static Site Generator)。可以把 Markdown 文件转化成 H...

总结Ueff和U10的关系

根据WRF-LES模拟的nc文件,拟合出总结Ueff和U10的关系式

根据WRF-LES生成的逐分钟级的Plume,可以用来总结Ueff和U10的关系式,也可以用来模拟给定排放速率下的羽流 ​ 未完待续…. 引用 如果你使用了这些代码,请对以下论文施加引用: Pei Z, Han G, Mao H, et al. Improving quantification of methane point source emissions from ...

Plume Gif 制作

根据WRF-LES模拟的nc文件生成不同排放流速率下、不同时间戳的柱浓度图像

一、初步测试 先写个简单的代码测试下单帧的柱浓度图像。 clc;clear % 可能需要修改的量 pixel_res = 30; % 空间分辨率 emission_per_hour = 2000; % kg/h M_gas = 16.04; % CO2:44 CH4:16.04 %% 第一步,得到缩放因子 n_total_air_1m2 = 3.6055e+05; % 底面积是 1m2...

干空气柱积分质量和有效风速的自动获取

Automatic acquisition of integral mass and effective wind speed for dry air columns

IME算法中,需要用到的气象参数有: 干空气柱积分质量 vertical_mass_of_atmosphere 有效风速 Effective wind speed 一、干空气柱积分质量 ERA5官网中:single level 提供了空气柱积分质量,尽管不是干空气的,但差异不大;pressure_level 提供了不同高度处的温湿压等参数,我们可以手动来...

理解误差传播定律

在统计学中,不确定性的传播(或误差的传播)是变量的不确定性(或误差,更具体地说是随机误差)对基于它们的函数的不确定性的影响。当变量是实验测量值时,由于测量限制(如仪器精度),它们具有不确定性,这种不确定性由于函数中变量的组合而传播。 一个比较简单的例子:使用温度计测量温度,单位为摄氏度,假设测得20度,仪器误差为10%,即测量结果为20 +- 0.2 度。如何用华氏度表示这次测量的不确...

Matlab 自定义Legend

模仿SA论文-在卫星底图上绘制点

起初是看到Sci. Adv.上一篇论文里的一张图挺好看的,想模仿画一下。 当然,如果作图量不大的话,也可以用Arcgis来出图,这篇博客主要介绍用MATLAB如何来实现绘制这种图。 首先要用高分辨率卫星作为底图,所以我们采用了geobasemap satelite来实现。 用Marker size来映射排放量的大小,这个可以通过geoscatter或者geoplot来实现(类似于...

利用爬虫获取山西能源局提供的煤矿级数据

本来就不怎么会用Python,再加上好久好久没用了,花了两天时间解决了科研过程遇到的实际问题。 总体思路 使用爬虫获取山西能源局提供的煤矿级的数据(包括煤矿名称、所属城市、开采深度、瓦斯类型等) 根据煤矿名称,数量少的可以直接通过坐标拾取系统,数量多的可以利用百度地图api,获得煤矿对应的经纬度(百度坐标系下) ...