2017/03/14 海洋内波和另外关于opencv中的椭圆检测方法

最近终于在内波的问题上有一点点的突破了。

之前在查阅文献的时候，总是了许多的公式，从KdV方程的，到Bragg散射和组合模型的，自己也明白内波到海洋表面需要一个调制过程，但是却没有注意到一点，内波到海表面的调制过程是可以归结为两种表达式：

第一种是使用得到的内波的幅度表达式Ux=−Cpahseηhupper，其中η为内波的幅度，h为上层水体的厚度，C为内波的相速度，结合内波传播的相速度以及上层液体的厚度，使用简单的分式就可以求出海表面的流场，有了流场之后，可以利用M4S软件进行后面的波高谱计算，求解后向散射系数和成像也都集成在M4S软件中；

另外一种则是通过求解流函数，然后对流函数进行三个维度的求偏导，得到三个维度的速度Ux=∂Φ∂x、Ux=∂Φ∂y，之后的步骤也是通过M4S软件实现。

以下内容均针对一维内波进行考虑，暂时不考虑二维内波

前面对于内波幅度的获得就是通过求解KdV方程。这一部分属于动力学的内容，而后面部分，尤其是在组合模型上面，是直接与雷达参数相关的。当获取了内波的振幅之后，求得海表面的流场是一个简单的表达式问题。后面部分首先应该能够利用公式来解释各项参数对雷达成像的影响，同时要提出利用公式完成的仿真来显示这种影响，同时利用M4S的结果来支持仿真的结果和解释的合理性。

思路这样一捋顺之后就知道后面的工作应该是什么了！

内波RCS的表达式大致可以表示为如下形式δσσ0=±4+γμA，其中A为内波的的流速偏导，是与内波相关的量，在进行雷达参数影响衡量的过程中，可以认定为恒定值，不予考虑。而γ=1+3k2τ/gρ1+k2τ/gρ，μ为张弛率，是张弛时间τ的倒数，剩下的字母k表示雷达波矢量的模值。这个波矢量在之前的σ0中也有涉及

今天完成了 Practical OpenCV书本中关于椭圆检测的代码誊写，该部分主要使用的方法是随机采样审查方法，其大致流程可以表述为下：

随机的从数据集中抽样

找出满足抽样出来的点的方程

找到“内嵌代码”—即满足这个方程的那些点。通过测量距离，来确定一个点是不是和这个模型匹配。通常需要一个距离门限来决定一个点是不是满足

迭代直到你确定一个模型，在该模型下面，“inlier”的数目超过你设定的门限