请教高手剖面曲线拐点(曲率最大点)的判定方法
一剖面曲线(如附件图1所示),横轴为测量值X,纵轴为深度,为了提取其跃层,需先判定剖面曲线的拐点,请教该情况拐点的判定方法?相关定义如下:某要素垂直分布曲线(如附件图2所示)上曲率最大的点A、B被称为拐点,当ΔX/ΔZ>0.2时,AB段即为跃层,A、B分别为跃层的顶界和底界。垂直分布曲线上可以有单跃层(如附件图3所示),也可以有双跃层(如附件图4所示)。
根据定义知“拐点”是整体趋势上曲率最大的点,针对图1所示的实测曲线(数据源如附件data1.txt说示,第一列为深度,第二列为测量值X),请教如何判定拐点?我先对曲线进行简单的平滑,然后进行圆的拟合,求曲率最大点,但结果不理想,很难找到整体趋势上的曲率最大点。
恳请高手赐教,帮忙找出解决的方法,若能给出C/C++程序更是感激不尽!!!
图1234和数据源.rar
(91.37 KB)