Stress2Grid 的 Matlab 脚本挺实用的,专门搞应力场那块的朋友可以留意下。是想S_Hmax方向和它的空间分布模式的时候,这个工具用起来会比较顺手。它不是插值,而是基于前人工作做了升级优化,算速度快、功能还多。
固定搜索半径和可变搜索范围这两个方案都能跑,按你的需求灵活切换。有时候你数据点不多,它也能给你一个大致靠谱的结果,不至于全靠猜。
还有一点挺方便的,就是它对规则网格支持比较好。做地图展示或者跟板块运动方向对比的时候,坐标统一了多事就清晰了。脚本结构清晰,变量命名也比较直白,改起来没啥负担。
如果你之前用过 Heidbach 那套插值脚本,转过来会蛮快上手的,逻辑差不多,但这版扩展性强一些。应力波长或者搞图像输出都还不错,配合其他工具比如 ArcGIS 来用,效果也挺好。
想进一步玩主应力计算的,也可以看看这个相关资源:3D 主应力计算脚本,两个脚本一起搭配起来用,深度一下子就上来了。
如果你做地应力场建模、搞插值可视化或是需要跟板块运动方向对比的项目,这脚本还挺值得试试的。用 Matlab 跑起来也流畅,依赖少,上手快。