I am looking for a way to detect shearlines or confluence regions such as that shown in the figure below (red dashed line).
I was able to derive the stretching, shear, and total deformations.
The shearline coincides with the region of maximum stretching deformation.
Currently, I am drawing the shearline manually. But I am wondering if there is a way to automatically detect this line or the deformation zone.
I am thinking of two ways
Looking at the dynamic variables (i.e., stretching deformation)
Using Python's pattern recognition function to detect this kind of pattern.
Problems
Any suggestions on how I can detect this automatically. I don't know how to implement the above methods. I am not sure what criteria to set in order to draw the shearline.
I am also new to Python so I don't know how to implement the second one.
I have wind data and calculated deformations in netcdf format. I uploaded the files to the following links:
Link to calculated deformation parameters
The deformation parameters are calculated this way:
where A, B, and V are the shearing, stretching, and total deformation parameters.
In the deform.nc file, shear_sph(shearing deformation; A in the equation above); stretch_sph(stretching deformation; B in the equation B); deform_sph(total deformation; V in the equation above).
I'll appreciate any help about this matter.
SUMMARY
I am looking for the NCL Shear, stretch, deform functionality in Python. Are there any existing software packages or do I have to write my own ?