+ double r = \r
+ Math.Sqrt(1 - Const.FirstEccentricitySquared * sin * sin);\r
+ return Const.SemiMajorAxisA / r;\r
+ }\r
+\r
+ public static System GetSystemWithNearestOriginTo(\r
+ GeodeticCoordinate coord)\r
+ {\r
+ System nearestSystem = System.Unknown;\r
+ double min = Double.MaxValue;\r
+ for (int i = 1; i <= (int)System.XIX; ++i)\r
+ {\r
+ System s = (System)i;\r
+ Jpc c = new Jpc(coord, s);\r
+ double d = c.X * c.X + c.Y * c.Y;\r
+ if (d < min)\r
+ {\r
+ min = d;\r
+ nearestSystem = s;\r
+ }\r
+ }\r
+ return nearestSystem;\r