编辑-地球科学江南体育网页版堆栈交换江南电子竞技平台

你不是登录。你的编辑将被放置在一个队列,直到同行评议。

我们欢迎编辑,使读者更容易理解和更有价值。因为社区成员审查编辑,请尽量使文章更好比你怎么发现,例如,通过修正语法或添加额外的资源和超链接。

牧师

定位shearlines或变形区域的地图

我正在寻找一种方式来检测shearlines或融合等地区,如下图所示(红色虚线)。

我能够推导出拉伸、剪切和总变形。

shearline正值地区最大拉伸变形。

目前,我画shearline手动。但是我想知道如果有一种方法可以自动检测这条线或变形区。

我想两种方式

看(即动态变量。、拉伸变形)
使用Python的模式识别功能来检测这种模式。

问题

任何建议我可以自动检测到。我不知道如何实现上面的方法。我不确定什么标准设置为了画出shearline。

我也是新到Python所以我不知道如何实现第二个。

我有风netcdf格式的数据和计算变形。我上传文件到以下链接:

UWND

VWND

链接到变形计算参数

变形参数计算:

在A、B和V是剪切、拉伸,和总变形参数。

在变形。数控文件,shear_sph(剪切变形;在上面的方程中);stretch_sph(拉伸变形;B在方程中B);deform_sph(总变形;在上面的方程V)。

我非常感谢关于这件事的任何帮助。

总结

我在寻找NCL剪切、拉伸变形在Python的功能。有没有现有的软件包或我必须编写自己的吗?

回答

我不能够充分实际代码在Python中但我希望我可以点你在正确的方向上如果我计算分流和汇合。可以先做一个[流线阴谋使用matplotlib] [1]。一旦你有了,你需要实际的简化数据中描述这个链接- [https://matplotlib.org/api/_as_gen/matplotlib.pyplot.streamplot.html] [2], LineCollection给出数据结构。每简化LineCollection需要计算(切向和法向量)[3]基于局部笛卡尔坐标系统的起源的想法集中在某处沿着流线。一旦你有了切向和法向量您将需要使用这个方程即散度在自然坐标系$ $ \微分算符。v_h = (e_n \压裂{\部分e_s}{\部分n} * v_h) + \压裂{\部分v_h}{\部分年代}$ $使用的组合(Frenet艾史蒂夫公式][4]这可以显示等于$ $ \压裂{| v_h |} {R_n} + \压裂{\部分v_h}{\部分年代}$ $从你的角度你只要第一项方程计算分流和汇合。所以R_n美元的曲率半径法向量(之前你已经计算)。为了计算曲率半径可以使用这个(python代码的最小二乘拟合)[5](使用最小二乘的方法——获得椭圆然后计算半径)。如果半径是积极的我. .e R_n > 0美元然后你有分流和消极的融合。 **UPDATE** If you are looking for the exact [NCL shear stretch deform][6] in Python what you can do is to use this package in python [pyspharm][7]. The data has to be global and then you need to orient the grid in a south to north direction. Once you have that you can use that package to calculate the shear, stretch and deformation parameters by using the [pyspharm's **getgrad()**][8] method to calculate the vector gradient of the wind vectors. If you pass in the components of the wind vectors(both u and v) to getgrad() you will get the zonal and meridional components of the gradient of the wind vectors. Adding and subtracting them will get you the shear and stretch. Sample code from spharm import Spharmt uspectral = Spharmt.grdtospec(u,ntruncation) vspectral = Spharmt.grdtospec(v,ntruncation) dudx,dudy = Spharmt.getgrad(uspectral) dvdx,dvdy = Spharmt.getgrad(vspectral) #shear and stretch shear = dvdx + dudy stretch = dudx - dvdy [1]: https://matplotlib.org/gallery/images_contours_and_fields/plot_streamplot.html [2]: https://matplotlib.org/api/_as_gen/matplotlib.pyplot.streamplot.html [3]: https://scicomp.stackexchange.com/questions/16444/determine-unit-outward-normal-vector-for-a-curve/16445#16445 [4]: https://en.wikipedia.org/wiki/Frenet%E2%80%93Serret_formulas [5]: https://github.com/bdhammel/least-squares-ellipse-fitting [6]: https://www.ncl.ucar.edu/Document/Functions/Contributed/shear_stretch_deform.shtml [7]: https://pypi.org/project/pyspharm/ [8]: https://github.com/jswhit/pyspharm/blob/master/Lib/spharm.py#L845 [9]: https://ajdawson.github.io/windspharm/latest/

编辑总结

取消

\ begingroup美元嗨@gunsub,非常感谢。我将看更多。我猜没有直截了当的解决方案。我计算deformatiion参数在NCL使用球面谐波。这是更精确的比为中心有限差分。 \ endgroup美元
- - - - - -Lyndz
2019年5月20日在十二24吗
1

\ begingroup美元适当@Lyndz我已经更新了我的答案。 \ endgroup美元
- - - - - -gansub
2019年5月20日12:26

添加一个评论 |

正确的小错误或错误
澄清的意思,而不更改它
添加相关资源或链接
总是尊重作者的意图
不使用编辑与作者回复

创建代码用引号”或式子~篱笆
' ' '
像这样
' ' '
添加语言标识符来突出显示的代码
”“python
def函数(foo):
打印(foo)
' ' '
把段落之间的回报
在结束linebreak添加2个空间
_italic_或* *的* *
报价将>开始的
做链接(尽量使用https)
< https://example.com >
(例子)(https://example.com)
< a href = " https://example.com " > < / >示例
MathJax方程罪\ \ ^ 2θ美元

格式帮助»
回答帮助»

MathJax帮助»