/// <summary>
/// Returns the average curvature of the substroke
/// </summary>
public double getAvgCurvature(Substroke sub)
{
this.arclength = new ArcLength(sub.Points);
this.slope = new Slope(sub.Points);
this.curve = new Curvature(sub.Points, arclength.Profile, slope.TanProfile);
return(curve.AverageCurvature);
}
private void btnSCal_Click(object sender, EventArgs e)
{
//获取所选椭球参数值
double a = Convert.ToDouble(txtA.Text);
//由于double精度有限，所以选择通过扁率的倒数β来计算e的平方
double e2 = (2 * Convert.ToDouble(txtB.Text) - 1) /
Convert.ToDouble(txtB.Text) / Convert.ToDouble(txtB.Text);
double n = 0; //迭代次数或者迭代精度
if (txtSE.Enabled)
{
n = Convert.ToDouble(txtSE.Text);
}
//存储从ListView中获取到弧长
List <string> list = new List <string>();
Data.FromListView(this.listLS, list);
ArcLength arc = new ArcLength(a, e2);
List <string> listRes = new List <string>();
//经典迭代算法结果
List <string> list_c = new List <string>();
foreach (var item in list)
{
double rad = arc.B_standard(Convert.ToDouble(item), 1e-10);
list_c.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
//根据数值积分计算弧长
switch (this.cboxSMethod.SelectedIndex)
{
case 2: //单点迭代法
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_single(Convert.ToDouble(item), (int)n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 1: //二分法
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_half(Convert.ToDouble(item), n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 3: //牛顿迭代法
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_newton(Convert.ToDouble(item), (int)n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 4: //割线法
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_sec(Convert.ToDouble(item), (int)n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 6: //欧拉迭代公式
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_Euler_ex(Convert.ToDouble(item), (int)n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 7: //欧拉预估-校正算法
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_Euler_prex(Convert.ToDouble(item), (int)n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 8: //2阶龙哥格库塔算法
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_RK2x(Convert.ToDouble(item), (int)n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 9: //4阶龙格-库塔算法
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_RK4x(Convert.ToDouble(item), (int)n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 11: //欧拉与牛顿迭代
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_Euler_EX(Convert.ToDouble(item), (int)n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 12: //欧拉预估-校正与牛顿迭代
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_Euler_preX(Convert.ToDouble(item), (int)n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 13: //2阶龙格-库塔算法与牛顿迭代
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_RK2X(Convert.ToDouble(item), (int)n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 14: //4阶龙格-库塔算法与牛顿迭代
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_RK4X(Convert.ToDouble(item), (int)n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
default: break;
}
//在ListView中显示与结果
if (cboxSMethod.SelectedIndex == 0)
{
Data.AddDataList(this.listLS, list, list_c);
}
else
{
List <string> listDif = XZ_dif(listRes, list_c);
Data.AddDataList(this.listLS, list, listRes, listDif);
}
}
//计算子午线弧长
private void btnLZCal_Click(object sender, EventArgs e)
{
//获取所选椭球参数值
double a = Convert.ToDouble(txtA.Text);
//由于double精度有限，所以选择通过扁率的倒数β来计算e的平方
double e2 = (2 * Convert.ToDouble(txtB.Text) - 1) /
Convert.ToDouble(txtB.Text) / Convert.ToDouble(txtB.Text);
//获取步长或者精度
double n = 0;
if (txtZE.Enabled)
{
n = Convert.ToDouble(txtZE.Text);
}
//存储从ListView中获取到的“纬度B”
List <string> list = new List <string>();
Data.FromListView(this.listLZ, list);
List <string> listStr = new List <string>();
foreach (var item in list)
{
listStr.Add((Convert.ToDouble(item)).ToString("0.0000"));
}
//将list中的“纬度”转换为弧度制
//根据所选算法计算弧长，将结果存放在listRes集合中。
List <string> listRes = new List <string>();
ArcLength arc = new ArcLength(a, e2);
//经典算法结果
List <string> list_c = new List <string>();
foreach (var item in list)
{
list_c.Add(arc.X_standard(DmsRad.dms2rad(item)).ToString("0.000"));
}
//根据数值积分计算弧长
switch (this.cboxZMethod.SelectedIndex)
{
case 1: //复化梯形算法
listRes.Clear();
foreach (var item in list)
{
double B_rad = DmsRad.dms2rad(item.ToString());
listRes.Add(arc.X_ft((int)(B_rad / (n / DmsRad.p)), B_rad).ToString("0.000"));
}
break;
case 2: //复化辛普森算法
listRes.Clear();
foreach (var item in list)
{
double B_rad = DmsRad.dms2rad(item.ToString());
listRes.Add(arc.X_fs((int)(B_rad / (n / DmsRad.p)), B_rad).ToString("0.000"));
}
break;
case 3: //复化科特斯算法
listRes.Clear();
foreach (var item in list)
{
double B_rad = DmsRad.dms2rad(item.ToString());
listRes.Add(arc.X_fc((int)(B_rad / (n / DmsRad.p)), B_rad).ToString("0.000"));
}
break;
case 4: //龙贝格算法
listRes.Clear();
foreach (var item in list)
{
double B_rad = DmsRad.dms2rad(item.ToString());
listRes.Add(arc.X_r(n, B_rad).ToString("0.000"));
}
break;
//case 5：
case 6: //欧拉迭代公式
listRes.Clear();
foreach (var item in list)
{
double B_rad = DmsRad.dms2rad(item.ToString());
listRes.Add(arc.X_Ee(Convert.ToInt32(n), B_rad).ToString("0.000"));
}
break;
case 7: //欧拉矫正公式
listRes.Clear();
foreach (var item in list)
{
double B_rad = DmsRad.dms2rad(item.ToString());
listRes.Add(arc.X_Epre(Convert.ToInt32(n), B_rad).ToString("0.000"));
}
break;
case 8: //二阶龙格-库塔算法
listRes.Clear();
foreach (var item in list)
{
double B_rad = DmsRad.dms2rad(item.ToString());
listRes.Add(arc.X_RK2(Convert.ToInt32(n), B_rad).ToString("0.000"));
}
break;
case 9: //龙贝格算法
listRes.Clear();
foreach (var item in list)
{
double B_rad = DmsRad.dms2rad(item.ToString());
listRes.Add(arc.X_RK4(Convert.ToInt32(n), B_rad).ToString("0.000"));
}
break;
default: break;
}
//在ListView中显示与结果
if (cboxZMethod.SelectedIndex == 0)
{
Data.AddDataList(this.listLZ, listStr, list_c);
}
else
{
List <string> listDif = XZ_dif(listRes, list_c);
Data.AddDataList(this.listLZ, listStr, listRes, listDif);
}
}
private void btnSCal_Click(object sender, EventArgs e)
{
//获取所选椭球参数值
double a = Convert.ToDouble(txtA.Text);
//由于double精度有限，所以选择通过扁率的倒数β来计算e的平方
double e2 = (2 * Convert.ToDouble(txtB.Text) - 1) /
Convert.ToDouble(txtB.Text) / Convert.ToDouble(txtB.Text);
int n = 0; //迭代次数
if (txtSE.Enabled)
{
n = Convert.ToInt32(txtSE.Text);
}
//存储从ListView中获取到弧长
List <string> list = new List <string>();
Data.FromListView(this.listLS, list);
ArcLength arc = new ArcLength(a, e2);
List <string> listRes = new List <string>();
//经典迭代算法结果
List <string> list_c = new List <string>();
foreach (var item in list)
{
double rad = arc.B_standard(Convert.ToDouble(item), 1e-10);
list_c.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
//根据数值积分计算弧长
switch (this.cboxSMethod.SelectedIndex)
{
case 1: //单点迭代法
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_single(Convert.ToDouble(item), n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 2: //牛顿迭代法
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_newton(Convert.ToDouble(item), n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
case 3: //割线法
listRes.Clear();
foreach (var item in list)
{
double rad = arc.B_sec(Convert.ToDouble(item), n);
listRes.Add(Convert.ToDouble(DmsRad.rad2dms(rad)).ToString("0.0000"));
}
break;
default: break;
}
//在ListView中显示与结果
if (cboxSMethod.SelectedIndex == 0)
{
Data.AddDataList(this.listLS, list, list_c);
}
else
{
List <string> listDif = XZ_dif(listRes, list_c);
Data.AddDataList(this.listLS, list, listRes, listDif);
}
}
