cost733超简单安装运行
本帖最后由 SunJiaming 于 2024-3-31 14:01 编辑最近无意间发现有人将cost733class 1.4上传到了anaconda，可以仿照conda安装NCL的方式，安装并运行cost733class原网址：Cost733Class :: Anaconda.org
我的系统是Ubuntu 22.04 LTS
(1) 创建环境，命名为cost733，这步不用管python版本
conda create -n cost733
(2) 创建好环境后，启动环境cost733并运行如下代码安装cost733
conda activate cost733
conda install -c conda-forge cost733class
如果这步出现问题，最好换一下conda的源，不要用镜像源。
(3) 等待安装完成后，我们拿hgt.1999.nc做个测试，空间分辨率2.5度，时间分辨率6小时，算法使用PCT（T-mode PCA）
cost733class -dat pth:/home/sunjm/data/ncep_r2/daily/hgt.1999.nc var:hgt sle:500 slo:80:130:2.5 sla:10:60:2.5 -met PCT -ncl 4 -claoutput-PCT.txtpth换成自己的文件地址，sle代表高度层，slo代表经度范围，sla代表纬度范围，2.5表示文件分辨率，-met表示算法名称（PCT），-ncl表示分型数量，-cla表示输出文件名称，自己起个名字就行，如果不指定输出位置，就会在根目录输出（我的是/home/sunjm/)
(4) output-PCT.txt中长这样，前4列为年月日小时，最后一列代表每个时次对应的天气类型
(5) 再计算一下ECV （explain cluster variance）
cost733class -dat pth:/home/sunjm/data/ncep_r2/daily/hgt.1999.nc-clain pth:/home/sunjm/output-PCT.txt dtc:4 ddt:6h -met ECV -idx clain-ECV
或者：cost733class -dat pth:/home/sunjm/data/ncep_r2/daily/hgt.1999.nc var:hgt sle:500 slo:80:130:2.5 sla:10:60:2.5 -clain pth:/home/sunjm/output-PCT.txt dtc:4 ddt:6h -met ECV -idx clain-ECV
-dat是原始数据的位置，-clain是输出结果的位置，dtc代表文件的日期数据有4列，这个必须要加；ddt代表时间分辨率6小时，-met代表计算ECV
(6) ECV最后会在屏幕输出，我的结果是0.0556
static/image/hrline/line3.png
2024-03-31 更新
在计算PCT时，最好用var声明下变量名，不然cost733有可能识别不出来参与运算的变量而报错，并且对netcdf文件还有个要求，那就是nc文件的纬度（lat）只能从高纬度排序至低纬度（即从90°~-90°），但是在计算时sla的范围要求从小到大填写（如例子中10:60，而不能是60:10，否则也会报错）
点赞{:eb502:} 太棒了！感谢楼主分享！ 你好，博主大大，这个分型数量是根据什么弄出来的哇 南极的七色鹿 发表于 2024-3-7 12:31
你好，博主大大，这个分型数量是根据什么弄出来的哇
Ning, G., Yim, S.H.L., Wang, S. et al. Synergistic effects of synoptic weather patterns and topography on air quality: a case of the Sichuan Basin of China. Clim Dyn 53, 6729–6744 (2019). https://doi.org/10.1007/s00382-019-04954-3
可以参考下这篇论文，简单来说，每一种分型数都对应一个ECV，取前后两个ECV的差值最大时，对应的分形数就是最佳的分型 SunJiaming 发表于 2024-3-7 14:53
Ning, G., Yim, S.H.L., Wang, S. et al. Synergistic effects of synoptic weather patterns and topogr ...
好嘞博主大大！超级感谢诶 南极的七色鹿 发表于 2024-3-7 16:35
好嘞博主大大！超级感谢诶
客气客气~一起学习！ SunJiaming 发表于 2024-3-7 22:01
客气客气~一起学习！
博主大大你好~可以麻烦再问您一个问题吗？就是如果我对某个省份进行天气分型，那么用PCT方法进行命令操作时，经纬度的范围是不是要比该省的实际经纬度更大点哇？ 南极的七色鹿 发表于 2024-3-7 23:43
博主大大你好~可以麻烦再问您一个问题吗？就是如果我对某个省份进行天气分型，那么用PCT方法进行命令操作 ...
最好能大一些，因为范围太小的话数据点太少，可能会报错，并且分型出来的结果可能太相似，没法分析 SunJiaming 发表于 2024-3-8 10:16
最好能大一些，因为范围太小的话数据点太少，可能会报错，并且分型出来的结果可能太相似，没法分析
好嘞~太感谢博主大大啦
页:
[1]
2