自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+

GIS

GIS小白进阶之路

  • 博客(17)
  • 收藏
  • 关注

原创 WRF替换GVF数据

GVF(Green Vegetation Fraction) 是决定WRF模拟局地/区域气候研究的关键参数。研究表明GVF对模式模拟温度、湿度和潜热误差校正贡献率分别为62%,87%和92%。因此本文提供对WRF中的GVF进行替换的具体方法。GVF的替换和LAI的替换几乎一模一样,区别就是将wrflowinp_d01中的LAI参数替换成VEGFRA即可。LAI的替换见我的另外一个帖子。同时如果有更好的方法,欢迎评论区留言讨论!执行完代码后就可以将生成的tiff格式的GVF数据下载到本地进行处理。

2023-07-31 12:19:44 367 1

原创 WRF中替换LAI数据

读取替换后的geo_em.d02_LAI.nc,可以看第二张图,只有一层数据(目前本人还不知道如何实现如模型默认的生成12层LAI数据,留作讨论和后续摸索)也符合index文件中的设定,但是这里并没有将像元值乘以index中设置的scale_factor=0.01,结果如第三张图所示。最后一步就是,将geo_em.d02_LAI.nc中的2018年第257天的LAI数据替换掉geo_em.d02.nc中的变量LAI12M中的9月份数据,即可。这样就替换完成了一个月的LAI数据,其他月份也是如此。

2023-04-14 15:37:06 1146 5

原创 中国1km分辨率的DEM数据以及合并后的中国行政区划数据

   本人在进行将中科院土地覆被类型分类转换为IGBP类型分类时,需要用到中国1km公里分辨率的dem数据。数据的来源和制作来源于如何下载拼接全球陆地1km分辨率DEM,非常感谢他的分享,我这里只是将他全球的数据裁剪出中国的部分,以便需要的小伙伴可以直接使用,不用再次浪费时间进行裁剪!   行政区划数据为shp格式,我将全国合并为一个面。如果需要省级、市级以及县级的行政区划数据,可以私聊我。百度网盘链接:https://pan.baidu.com/s/1mhI3HaA-kbyd6hycbBbBhQ提取

2022-05-28 16:20:57 1035 2

原创 WRF Output Fields描述

Lu_index = land use category 土地利用类型Znu= eta values on half(mass) levels 一半(质量)水平上的eta值Znw=eta values on full(w) levels 全(w)级eta值ZS=depths of centers of soil layers 土层中心深度,单位mDzs=thicknesses of soil layers 土层厚度,单位mVar_SSO =variance of subgrid-sc

2022-05-21 00:38:46 1407

原创 WRF中namelist.wps与namelist.input参数详解

   作为一位的WRF初学者,在设置namelist.wps和namelist.input时,我总是感觉在设置namelist的时候较为模糊。于是,我通过参考官网方文档、网络上相关的资料和自己模拟时设置的经验,对namelist中的一些参数进行解析,以便我们更好的设置和理解这些参数。1、namelist.wps&shareWrf_core=‘ARW’, #告诉WPS输入数据准备用哪个动态核心,可选‘ARW’或‘NMM’Max_dom=2, #最大嵌套层数,包括粗域start_date =

2022-05-15 23:09:33 19051 11

原创 WRF中使用ERA5再分析数据

   本人在做一些实例时,需要用到最新的ERA5再分析数据集,其中我下载的是ERA5 hourly data on pressure levels form 1979 to present和ERA5 hourly data on single levels from 1979 to present中的3天数据。1、注册climate copernicus账号下载ERA5数据前先要注册其账号2、选择变量   WRF中对于气象数据的字段有具体的要求,......

2022-05-13 11:51:50 4733 18

原创 WRF替换静态地理数据中的土地利用数据(WRF替换下垫面数据)

   在进行自己的WRF项目模拟的时候,因为需要用中国土地利用数据替换WRF内嵌老旧的土地利用数据,故需要将土地利用进行替换。我使用的中国的土地利用数据是来自中国土地利用数据(1980-2015)。其中我以下载的lucc1980为例进行演示。1、坐标转换   (1)下载的lucc1980数据如下图所示:   (2)在ArcMap中打开如下图所示:打开属性表可以看到,lucc1980的坐标系为投影坐标系Krasovsky_1940_Albers,数据格式为GRID。这里需要将投影坐标系转换成地理坐

2022-05-13 01:38:46 7231 68

原创 WRF模型模拟时所遇到的问题及解决方法

在运行real.exe时出现的问题-------------- FATAL CALLED ---------------FATAL CALLED FROM FILE: LINE: 401error opening met_em.d01.2016-10-06_00:00:00.nc for input; bad date in namelist or file not in directory-------------- FATAL CALLED ---------------FATA

2022-05-10 20:34:16 3997 21

原创 WRF案例运行-Nested Model Runs

1、说明   本次案例是利用WRF进行Nested Model Runs的嵌套模拟,下面的例子来自官网Nested Model Runs。如果还没安装和配置好WRF,可以先看我这篇博客WRF模型配置。如果之前没有做过单个域WRF模拟的同学,建议可以先试试单个域的WRF模拟WRF案例运行-Single Domain Case。2、Nested Model Runs的基本流程...

2022-04-27 17:36:07 1860 2

原创 WRF案例运行-Single Domain Case

1、相关说明   本次案例运行是利用WRF进行Single Domain Case的模拟,模拟的区域是中国西北和青藏地区。如果还没安装和配置好WRF,可以先看我这篇博客WRF模型配置。   参考资料有:1、Single Domain Case - Hurricane Matthew 2、WRF官网案例运行2、知识准备   2.1 WRF案例运行基本流程,如下图所示:   2.2 WPS中三个模块的作用:   (1)geogrid:确定模式区域并把静态地理数据插值到格点。   (2)ungri

2022-04-12 10:58:35 3080 7

原创 WRF模型配置

WRF模型配置笔记1.相关说明2.System Environment Tests3.Building Libraries4.Library Compatibility Tests5.Building WRF6.Building WPS7.Static geography data8.Post processing1.相关说明VMware:VMware Workstation 16 ProUbuntu:ubuntu-14.04.6-desktop-amd64.iso官方文档:How to Compi

2022-02-23 01:27:19 3144 8

原创 通量数据处理(Python)——NC格式转为CSV格式

通量数据处理(Python)——NC格式转为CSV格式1.数据结构  NetCDF(network Common Data Form)网络通用数据格式包括变量、维和属性。通量数据RDMF_2011_L3.nc可利用软件Panoply进行可视化,如下图所示:2.数据读取  Pyton读取nc数据,现在大部分的方法利用netCDF4包的Dataset方法读取文件,但nc格式的通量数据则无法利用Dataset读取出变量的值。import pandas as pdimport osfrom netCD

2021-06-27 20:40:32 5110 3

原创 Matlab——数据拟合与插值

  前言:本次使用的Matlab版本为R2016b,使用的拟合方式是cftool工具拟合和fittype与fit自定义函数拟合。文章目录  一、cftool工具  1.二维数据  2.三维数据  二、fittype与fit自定义函数拟合  三、参考文献  一、cftool工具  打开Matlab,在命令行窗口中输入cftool或者在菜单栏中的APP中点击Curve Fitting按钮都可以打开cftool工具界面,如图1所示:图1  1.二维数据  这里的实验二维数据是:distance=[0

2020-05-19 16:17:07 7379

原创 室内定位——卡尔曼滤波原理、Matlab与javascript实现

前言:在基于Beacon的室内定位中,由于获取的RSSI值中含有大量噪声,故对RSSI值进行去噪是非常有必要的。

2020-05-18 10:59:20 4841 3

原创 室内定位——如何在微信小程序中获取Beacon的RSSI值

  前言:本文使用的Beacon是智石科技的Smart Beacon,手机的操作系统是Android,其中实现室内定位最重要的参数是Beacon发射的蓝牙信号中的RSSI值。废话不多说,请看下面分解。1.打开微信小程序开发文档:单击这个链接微信官方文档,找到蓝牙部分的API。2.使用的API依次是:  1) wx.openBluetoothAdapter():打开手机蓝牙适配器  2) wx.startBluetoothDevicesDiscovery():开始搜索附近的蓝牙设备  3) wx.

2020-05-16 12:23:27 3776 3

原创 Matlab------在Matlab中如何画圆

1、利用rectangle函数画圆&nbsp&nbsprectangle函数功能

2020-03-10 18:34:17 49313 4

原创 微信小程序------如何在其他页面获取和改变全局变量

最近自己在做毕业论文是利用iBeacon来进行室内定位,因为需求刚好使用到了微信小程序。由于自己是微信小程序这方面的小白,所以这个问题困扰了自己很久,其实这个问题很简单。结合自己入的坑和上网查的资料,谈谈自己经验,希望能够帮助到刚自学微信小程序的小白!

2020-01-13 10:07:40 5769 5

空空如也

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

提示
确定要删除当前文章?
取消 删除