ClemensFischer/XAML-Map-Control
1
0
Fork 0
mirror of https://github.com/ClemensFischer/XAML-Map-Control.git synced 2026-05-07 13:37:47 +00:00
Code Issues Projects Releases Packages Wiki Activity

Comments on map projections

Browse source
This commit is contained in:
Clemens 2022-02-23 23:19:32 +01:00
parent f6b4830f24
commit 6a14d740e8
9 changed files with 27 additions and 16 deletions
Download patch file Download diff file Expand all files Collapse all files

14
MapControl/Shared/AutoEquirectangularProjection.cs
View file

@ -9,6 +9,11 @@ using System.Windows;
namespace MapControl
{
/// <summary>
/// Auto-Equirectangular Projection - AUTO2:42004.
/// Equidistant cylindrical projection with standard parallel and central meridian set by the Center property.
/// See "Map Projections - A Working Manual" (https://pubs.usgs.gov/pp/1395/report.pdf), p.90-91.
/// </summary>
public class AutoEquirectangularProjection : MapProjection
{
public const string DefaultCrsId = "AUTO2:42004";
@ -16,24 +21,21 @@ namespace MapControl
public AutoEquirectangularProjection()
{
CrsId = DefaultCrsId;
IsNormalCylindrical = true;
}
public override Point LocationToMap(Location location)
{
var xScale = Wgs84MetersPerDegree * Math.Cos(Center.Latitude * Math.PI / 180d);
return new Point(
xScale * (location.Longitude - Center.Longitude),
Wgs84MetersPerDegree * (location.Longitude - Center.Longitude) * Math.Cos(Center.Latitude * Math.PI / 180d),
Wgs84MetersPerDegree * location.Latitude);
}
public override Location MapToLocation(Point point)
{
var xScale = Wgs84MetersPerDegree * Math.Cos(Center.Latitude * Math.PI / 180d);
return new Location(
point.Y / Wgs84MetersPerDegree,
point.X / xScale + Center.Longitude);
point.X / (Wgs84MetersPerDegree * Math.Cos(Center.Latitude * Math.PI / 180d)) + Center.Longitude);
}
}
}