一、实现方法
底层(100m栅格)
ConvLSTM + 膨胀卷积每个100m×100m栅格作为图节点,输入闪电密度、反射率及地形高程数据。采用 ConvLSTM (3×3核) 捕捉闪电发展时序特征,叠加膨胀率=3的膨胀卷积扩大感受野至3km范围 [^6]。
特征提取:123# Python代码示例 - 底层特征编码器low_level_feat = ConvLSTM(input_data, filters=128, kernel_size=(3,3))low_level_feat = DilatedConv2D(low_level_feat, dilation_rate=3)
中层(1km超节点)
图池化聚类按1km²区域将底层栅格聚类为超节点,特征由反射率加权聚合:$h_{mid} = \sum_{v \in Cluster} w_v \cdot h_v ,权重 w_v = \text{sigmoid}(\text{Reflectivity}_v)$ ^2。
门控GCN跨邻域(r=5km)传递信息,门控公式:$g_{gate} ...
Code logic analysis
model: Sonnet 4.5
time: 2025/11/6 13:14:35
conversationUUID: 4145c238-4a65-48d7-ad86-487e424e7182
human:分析这段代码的逻辑
file: 未命名附件
assistant:代码逻辑分析这是一个气象数据处理脚本,主要处理预测数据的CSV整合和概率图的文件管理。
📋 核心功能代码实现了两种时间跨度的数据处理:
6小时数据(h0-h5):短期预测
48小时数据(h0-h47):长期预测
🔄 执行流程12345678910111213141516171819202122输入: start_time="202511042000", end_time="202511061900" (2025年11月4日20:00 ~ 11月6日19:00,共48小时)┌─────────────────────────────────────┐│ 第一阶段:6小时处理 ...
2025 11 03使用长短期记忆LSTM ,预测得到的0到12小时效果好 POD 93.4 FAR 17.4 CSI为78.1
提出的知识和想法
雷暴早期需要满足静力不稳定、水汽和抬升触发 3 个基本条件
对流有效位能、沙氏指数、抬升指数和 K 指数与雷暴强度有较好的相关性大风指数、修正指数、静力能对流稳定度等23 个探空物理量参数与雷电活动关联度大于 0. 8,可作为雷电潜势预报指标因子。
分析超参数的选取对 LSTM 雷电活动预测模型的影响,确定学习率为 0. 000 1、批次样本量为32、序列长度为 10 用于预测,有利于模型的快速拟合,提高模型泛化能力。
LSTM 模型在(0,12] h 的雷电潜势预测效果明显优于(12,24〛h,其正确率为 95.8%,被试工作特征 AUC 为 0.989,POD 为 93.4%,FAR 为 17.4%,CSI为 78.1%。
随着预测时限的增加,(12,24]h 的雷电潜势预测存在 POD 较低、FAR 较高,下一步将提高模型深度、优化模型参数,结合雷达、卫星等观测资料增加空间维度信息,开展基于深度学习网络的格点化雷电预警技术研究及应用 ...
Linux常用命令速查手册nohup python train_new.py > train.log 2>&1 & 离线也可运行指令
📊 磁盘空间管理查看磁盘使用情况1234df -h # 查看所有磁盘分区使用情况(人类可读格式)df -h /home # 查看指定目录所在分区du -sh /path/to/dir # 查看指定目录总大小du -h --max-depth=1 # 查看当前目录下各子目录大小
查找大文件12du -ah /path | sort -rh | head -n 20 # 找出最大的20个文件/目录find / -type f -size +100M # 查找大于100M的文件
📁 文件和目录操作查看文件列表123456ls -l # 详细列表格式ls -lh # 以可读格式显示文件大小ls -lht # 按修改时 ...
平台错误一览当虚拟环境没有初始化时,conda init 然后开一个新的终端即可解决如果说某个包使用了默认环境的,使用指令unset PYTHONPATH 解决雷达数据处理完后的数据位置 /home/nmic/work/project/results/processedFile_radar工程代码路径 /home/nmic/project
main
清理旧文件(cleanup_intermediate_files)
删除上次遗留的csv lock等文件,以免程序中断导致错误
采集新气象数据(acquire_data)
获取到数据的起止时间:st前推时间 , et 当前时间 st_new 当前时间 et_new 往后5小时时间
数据处理与整合(process_all_data)
数据格式转换(convert_data_to_final_format)
调用AI模型进行预测(run_prediction_with_preloaded / create_fallback_pred ...
terrific 极好的 了不起的magnificent 壮观的 宏伟的 值得表扬的grand 壮丽的 堂皇的 极好的 精彩的
magnitude 巨大 重大 重要性 = importantmagnify 扩大 放大
verge 边缘 小草地emerge 出现 浮现 暴露emergency 紧急情况 突发事件 应急的 紧急的submerge 掩膜 湮没 掩盖submarine 潜艇 海底的 海平面下的marine 海洋的 海产的 海事的swamp 淹没 淹overwhelm 充溢 压倒 击败 征服compound 复合物 混合物 化合物 复合的
10-27hybrid 杂种动物,杂交植物 混合物compound 混合物chronic 慢性的,长久的benign 善良的 和善的 慈祥的 良性的acute 严重的 急性的convenience 方便 便利 方便的东西diagnose 诊断 判断(问题的原因)remedy 治 ...
public每24小时记录一次,一次预测240H,基准时间是当天的
2025/11/26问题一天内我只发现了两个文化,12点和0点的文件。但是拿到的nc数据中是每一小时一个。怎么转换。
气温分析文件: ECMFC1D_TEM_100_2025082112_ANEA_0.grib1简单读取结果:每12小时一个文件,包含所有level的温度可视化的图片
露点温度需要使用TEM和RHU但是给出的grib中只有表层和2m高的露点温度,所以使用TEM_100和RHU_100模仿 dewtemp_trh运算
至此 已成功计算出TT,K总结一下grib文件分析的思路和过程
使用ecCodes 的grib_ls 查看grib的基础信息grib_ls ECMFC1D_RHU_100_2025082112_ANEA_1.grib1 > RHU_100_1.log 2>&1 –重定向到当前文件文件
都是用_1的文件,包含各个层次
问题:刚获取的ECMFC1D数据中缺少我们需要的数据(Cape,Dcape等)。缺少前置知识,无法使用这些数据得到我们需要的数据(Dc ...


