// XAML Map Control - https://github.com/ClemensFischer/XAML-Map-Control
// © 2017 Clemens Fischer
// Licensed under the Microsoft Public License (Ms-PL)
using System;
using System.Globalization;
#if NETFX_CORE
using Windows.Foundation;
using Windows.UI.Xaml.Media;
#else
using System.Windows;
using System.Windows.Media;
#endif
namespace MapControl
{
///
/// Defines a map projection between geographic coordinates and cartesian map coordinates
/// and viewport coordinates, i.e. pixels.
///
public abstract partial class MapProjection
{
public const double Wgs84EquatorialRadius = 6378137d;
public const double MetersPerDegree = Wgs84EquatorialRadius * Math.PI / 180d;
///
/// Gets or sets the WMS 1.3.0 CRS Identifier.
///
public string CrsId { get; set; }
///
/// Indicates if this is a web mercator projection, i.e. compatible with map tile layers.
///
public bool IsWebMercator { get; protected set; } = false;
///
/// Indicates if this is an azimuthal projection.
///
public bool IsAzimuthal { get; protected set; } = false;
///
/// Gets the scale factor from longitude to x values of a normal cylindrical projection.
/// Returns NaN if this is not a normal cylindrical projection.
///
public double LongitudeScale { get; protected set; } = 1d;
///
/// Gets the absolute value of the minimum and maximum latitude that can be transformed.
///
public double MaxLatitude { get; protected set; } = 90d;
///
/// Gets the scaling factor from cartesian map coordinates to viewport coordinates.
///
public double ViewportScale { get; protected set; }
///
/// Gets the transformation from cartesian map coordinates to viewport coordinates (pixels).
///
public MatrixTransform ViewportTransform { get; } = new MatrixTransform();
///
/// Gets the scaling factor from cartesian map coordinates to viewport coordinates for the specified zoom level.
///
public virtual double GetViewportScale(double zoomLevel)
{
return Math.Pow(2d, zoomLevel) * TileSource.TileSize / 360d;
}
///
/// Gets the map scale at the specified Location as viewport coordinate units per meter (px/m).
///
public abstract Point GetMapScale(Location location);
///
/// Transforms a Location in geographic coordinates to a Point in cartesian map coordinates.
///
public abstract Point LocationToPoint(Location location);
///
/// Transforms a Point in cartesian map coordinates to a Location in geographic coordinates.
///
public abstract Location PointToLocation(Point point);
///
/// Translates a Location in geographic coordinates by the specified small amount in viewport coordinates.
///
public abstract Location TranslateLocation(Location location, Point translation);
///
/// Transforms a BoundingBox in geographic coordinates to a Rect in cartesian map coordinates.
///
public virtual Rect BoundingBoxToRect(BoundingBox boundingBox)
{
return new Rect(
LocationToPoint(new Location(boundingBox.South, boundingBox.West)),
LocationToPoint(new Location(boundingBox.North, boundingBox.East)));
}
///
/// Transforms a Rect in cartesian map coordinates to a BoundingBox in geographic coordinates.
///
public virtual BoundingBox RectToBoundingBox(Rect rect)
{
var sw = PointToLocation(new Point(rect.X, rect.Y));
var ne = PointToLocation(new Point(rect.X + rect.Width, rect.Y + rect.Height));
return new BoundingBox(sw.Latitude, sw.Longitude, ne.Latitude, ne.Longitude);
}
///
/// Transforms a Location in geographic coordinates to a Point in viewport coordinates.
///
public Point LocationToViewportPoint(Location location)
{
return ViewportTransform.Transform(LocationToPoint(location));
}
///
/// Transforms a Point in viewport coordinates to a Location in geographic coordinates.
///
public Location ViewportPointToLocation(Point point)
{
return PointToLocation(ViewportTransform.Inverse.Transform(point));
}
///
/// Transforms a Rect in viewport coordinates to a BoundingBox in geographic coordinates.
///
public BoundingBox ViewportRectToBoundingBox(Rect rect)
{
return RectToBoundingBox(ViewportTransform.Inverse.TransformBounds(rect));
}
///
/// Sets ViewportScale and ViewportTransform values.
///
public virtual void SetViewportTransform(Location projectionCenter, Location mapCenter, Point viewportCenter, double zoomLevel, double heading)
{
ViewportScale = GetViewportScale(zoomLevel);
ViewportTransform.Matrix = MatrixEx.TranslateScaleRotateTranslate(
LocationToPoint(mapCenter), ViewportScale, -ViewportScale, heading, viewportCenter);
}
///
/// Gets a WMS 1.3.0 query parameter string from the specified bounding box,
/// e.g. "CRS=...&BBOX=...&WIDTH=...&HEIGHT=..."
///
public virtual string WmsQueryParameters(BoundingBox boundingBox, string version = "1.3.0")
{
if (string.IsNullOrEmpty(CrsId))
{
return null;
}
var format = "CRS={0}&BBOX={1},{2},{3},{4}&WIDTH={5}&HEIGHT={6}";
if (version.StartsWith("1.1."))
{
format = "SRS={0}&BBOX={1},{2},{3},{4}&WIDTH={5}&HEIGHT={6}";
}
else if (CrsId == "EPSG:4326")
{
format = "CRS={0}&BBOX={2},{1},{4},{3}&WIDTH={5}&HEIGHT={6}";
}
var rect = BoundingBoxToRect(boundingBox);
var width = (int)Math.Round(ViewportScale * rect.Width);
var height = (int)Math.Round(ViewportScale * rect.Height);
return string.Format(CultureInfo.InvariantCulture, format, CrsId,
rect.X, rect.Y, (rect.X + rect.Width), (rect.Y + rect.Height), width, height);
}
}
}