高斯投影正反算 c#代码

* Math.Cos(Bf)); Nf=c/ Vf; double ynf = pcrdCar.y / Nf; pcrdGeo.dLatitude = Bf- 1.0 / 2.0 * Vf* Vf* tf* (ynf * ynf - 1.0 / 12.0 * Math.Pow(ynf, 4)* (5+3* tf* tf+ itaf * itaf -9* Math.Pow(itaf * tf, 2)) + 1.0 / 360.0 * (61 + 90* tf* tf+ 45* Math.Pow(tf, 4)) * Math.Pow(ynf, 6)); pcrdGeo.dLongitude = (ynf / Math.Cos(Bf) - (1+2* tf* tf+ itaf * itaf) * Math.Pow(ynf, 3)/ 6.0 / Math.Cos(Bf) + (5+ 28* tf* tf+ 24* Math.Pow(tf, 4)+6* itaf * itaf +8* Math.Pow(itaf * tf, 2)) * Math.Pow(ynf, 5)/ 120.0 / Math.Cos(Bf)); pcrdGeo.dLongitude = pcrdGeo.dLongitude + midlong; //pcrdGeo.dLatitude = pcrdGeo.dLatitude; richTextBox2.Text = "Results:\nLatitude: "+ Convert.ToString(pcrdGeo.dLatitude) + "\nLongtitude: "+ Convert.ToString(pcrdGeo.dLongitude); } private void label13_Click(object sender, EventArgs e) {}