/// <summary>
/// Находится ли точка в зоне ошибки?
/// </summary>
/// <param name="x0">Центр первого скана</param>
/// <param name="y0">Центр первого скана</param>
/// <param name="x1">Центр зоны ошибки (центр второго скана)</param>
/// <param name="y1">Центр зоны ошибки (центр второго скана)</param>
/// <param name="x2">Проверяемая точка</param>
/// <param name="y2">Проверяемая точка</param>
/// <returns></returns>
private bool isPointInPieZone(int x0, int y0, int x1, int y1, int x2, int y2)
{
//пока простая проверка: ровный разброс по углу и по длине
double l_sp2 = Parameters.getSquaredDistance(x0, y0, x2, y2);
double l_sp = Math.Pow(l_sp2, 0.5);
double psi_rl_rad = Parameters.getAngleRadian(x0, y0, x1, y1);
double psi_sp_rad = Parameters.getAngleRadian(x0, y0, x2, y2);
bool angleFlag = false; //входит ли по угловой зоне
double l_rl2 = Parameters.getSquaredDistance(x0, y0, x1, y1);
double l_rl = Math.Pow(l_rl2, 0.5);
double sgm_lrl = Parameters.getSgm_l(l_rl);
double l_rlPlus3sgm = l_rl + 3 * sgm_lrl;
double l_rlMinus3sgm = l_rl - 3 * sgm_lrl;
if ((psi_rl_rad > Math.PI / 2) && (psi_sp_rad < -Math.PI / 2))
{
if ((psi_sp_rad + 2 * Math.PI >= psi_rl_rad - 3 * Parameters.getSgm_psi_rad()) && (psi_sp_rad + 2 * Math.PI <= psi_rl_rad + 3 * Parameters.getSgm_psi_rad()))
{
angleFlag = true;
}
else
{
angleFlag = false;
}
}
else
{
if ((psi_rl_rad < -Math.PI / 2) && (psi_sp_rad > Math.PI / 2))
{
if ((psi_sp_rad - 2 * Math.PI >= psi_rl_rad - 3 * Parameters.getSgm_psi_rad()) && (psi_sp_rad - 2 * Math.PI <= psi_rl_rad + 3 * Parameters.getSgm_psi_rad()))
{
angleFlag = true;
}
else
{
angleFlag = false;
}
}
else
{
if ((psi_sp_rad >= psi_rl_rad - 3 * Parameters.getSgm_psi_rad()) && (psi_sp_rad <= psi_rl_rad + 3 * Parameters.getSgm_psi_rad()))
{
angleFlag = true;
}
}
}
if (((l_sp >= l_rlMinus3sgm) && (l_sp <= l_rlPlus3sgm)) && angleFlag)
{
return(true);
}
else
{
return(false);
}
}
/// <summary>
///
/// </summary>
/// <param name="x1"></param>
/// <param name="y1"></param>
/// <param name="x2"></param>
/// <param name="y2"></param>
/// <returns></returns>
private List <int[]> pieErrorZoneSearch(int x1, int y1, int x2, int y2)
{
List <int[]> pieErrorZone = new List <int[]>();
double l2 = Parameters.getSquaredDistance(x1, y1, x2, y2);
double l = Math.Pow(l2, 0.5);
double lminus3sgm = l - 3 * Parameters.getSgm_l(l);
double lplus3sgm = l + 3 * Parameters.getSgm_l(l);
double psi_rad = Parameters.getAngleRadian(x1, y1, x2, y2);
//половина стороны квадрата, в который точно вписан искомый сектор
int searchingSquareHalfSide = (int)Math.Ceiling(Math.Pow(Math.Pow(lplus3sgm * 3 * Parameters.getSgm_psi_rad(), 2) + Math.Pow(3 * Parameters.getSgm_l(l), 2), 0.5));
double l_rl2 = Parameters.getSquaredDistance(X0, Y0, X1, Y1);
double l_rl = Math.Pow(l_rl2, 0.5);
double sgm_lrl = Parameters.getSgm_l(l_rl);
double l_rlPlus3sgm = l_rl + 3 * sgm_lrl;
double l_rlMinus3sgm = l_rl - 3 * sgm_lrl;
for (int x = -searchingSquareHalfSide; x <= searchingSquareHalfSide; x++)
{
for (int y = -searchingSquareHalfSide; y <= searchingSquareHalfSide; y++)
{
double l_sp2 = Parameters.getSquaredDistance(x1, y1, x2 + x, y2 + y);
double l_sp = Math.Pow(l_sp2, 0.5);
double psi_sp_rad = Parameters.getAngleRadian(x1, y1, x2 + x, y2 + y);
if (isPointInPieZoneExpress(ref l_sp, ref psi_sp_rad, ref psi_rad, ref l_rlMinus3sgm, ref l_rlPlus3sgm))
{
pieErrorZone.Add(new int[2] {
x2 + x, y2 + y
});
}
}
}
return(pieErrorZone);
}