Version 5.0.1: Dropped MapProjection.UnitsPerDegree

This commit is contained in:
ClemensF 2020-04-02 18:06:39 +02:00
parent 2b41d298f4
commit 43c33c2564
12 changed files with 29 additions and 50 deletions

View file

@ -18,19 +18,19 @@ namespace MapControl
public override Point LocationToMap(Location location) public override Point LocationToMap(Location location)
{ {
var xScale = UnitsPerDegree * Math.Cos(Center.Latitude * Math.PI / 180d); var xScale = Wgs84MetersPerDegree * Math.Cos(Center.Latitude * Math.PI / 180d);
return new Point( return new Point(
xScale * (location.Longitude - Center.Longitude), xScale * (location.Longitude - Center.Longitude),
UnitsPerDegree * location.Latitude); Wgs84MetersPerDegree * location.Latitude);
} }
public override Location MapToLocation(Point point) public override Location MapToLocation(Point point)
{ {
var xScale = UnitsPerDegree * Math.Cos(Center.Latitude * Math.PI / 180d); var xScale = Wgs84MetersPerDegree * Math.Cos(Center.Latitude * Math.PI / 180d);
return new Location( return new Location(
point.Y / UnitsPerDegree, point.Y / Wgs84MetersPerDegree,
point.X / xScale + Center.Longitude); point.X / xScale + Center.Longitude);
} }
} }

View file

@ -27,7 +27,7 @@ namespace MapControl
GetAzimuthDistance(Center, location, out azimuth, out distance); GetAzimuthDistance(Center, location, out azimuth, out distance);
var mapDistance = distance * UnitsPerDegree * 180d / Math.PI; var mapDistance = distance * Wgs84EquatorialRadius;
return new Point(mapDistance * Math.Sin(azimuth), mapDistance * Math.Cos(azimuth)); return new Point(mapDistance * Math.Sin(azimuth), mapDistance * Math.Cos(azimuth));
} }
@ -42,7 +42,7 @@ namespace MapControl
var azimuth = Math.Atan2(point.X, point.Y); var azimuth = Math.Atan2(point.X, point.Y);
var mapDistance = Math.Sqrt(point.X * point.X + point.Y * point.Y); var mapDistance = Math.Sqrt(point.X * point.X + point.Y * point.Y);
var distance = mapDistance / (UnitsPerDegree * 180d / Math.PI); var distance = mapDistance / Wgs84EquatorialRadius;
return GetLocation(Center, azimuth, distance); return GetLocation(Center, azimuth, distance);
} }

View file

@ -33,23 +33,23 @@ namespace MapControl
public override Point LocationToMap(Location location) public override Point LocationToMap(Location location)
{ {
return new Point( return new Point(
location.Longitude * UnitsPerDegree, Wgs84MetersPerDegree * location.Longitude,
location.Latitude * UnitsPerDegree); Wgs84MetersPerDegree * location.Latitude);
} }
public override Location MapToLocation(Point point) public override Location MapToLocation(Point point)
{ {
return new Location( return new Location(
point.Y / UnitsPerDegree, point.Y / Wgs84MetersPerDegree,
point.X / UnitsPerDegree); point.X / Wgs84MetersPerDegree);
} }
public override string GetBboxValue(Rect rect) public override string GetBboxValue(Rect rect)
{ {
return string.Format(CultureInfo.InvariantCulture, return string.Format(CultureInfo.InvariantCulture,
CrsId == "CRS:84" ? "{0},{1},{2},{3}" : "{1},{0},{3},{2}", CrsId == "CRS:84" ? "{0},{1},{2},{3}" : "{1},{0},{3},{2}",
rect.X / UnitsPerDegree, rect.Y / UnitsPerDegree, rect.X / Wgs84MetersPerDegree, rect.Y / Wgs84MetersPerDegree,
(rect.X + rect.Width) / UnitsPerDegree, (rect.Y + rect.Height) / UnitsPerDegree); (rect.X + rect.Width) / Wgs84MetersPerDegree, (rect.Y + rect.Height) / Wgs84MetersPerDegree);
} }
} }
} }

View file

@ -31,7 +31,7 @@ namespace MapControl
GetAzimuthDistance(Center, location, out azimuth, out distance); GetAzimuthDistance(Center, location, out azimuth, out distance);
var mapDistance = distance < Math.PI / 2d var mapDistance = distance < Math.PI / 2d
? Math.Tan(distance) * UnitsPerDegree * 180d / Math.PI ? Math.Tan(distance) * Wgs84EquatorialRadius
: double.PositiveInfinity; : double.PositiveInfinity;
return new Point(mapDistance * Math.Sin(azimuth), mapDistance * Math.Cos(azimuth)); return new Point(mapDistance * Math.Sin(azimuth), mapDistance * Math.Cos(azimuth));
@ -47,7 +47,7 @@ namespace MapControl
var azimuth = Math.Atan2(point.X, point.Y); var azimuth = Math.Atan2(point.X, point.Y);
var mapDistance = Math.Sqrt(point.X * point.X + point.Y * point.Y); var mapDistance = Math.Sqrt(point.X * point.X + point.Y * point.Y);
var distance = Math.Atan(mapDistance / (UnitsPerDegree * 180d / Math.PI)); var distance = Math.Atan(mapDistance / Wgs84EquatorialRadius);
return GetLocation(Center, azimuth, distance); return GetLocation(Center, azimuth, distance);
} }

View file

@ -30,7 +30,7 @@ namespace MapControl
private double GetLineDistance() private double GetLineDistance()
{ {
var pixelPerDegree = ParentMap.ViewTransform.Scale * ParentMap.MapProjection.UnitsPerDegree; var pixelPerDegree = ParentMap.ViewTransform.Scale * MapProjection.Wgs84MetersPerDegree;
var minDistance = MinLineDistance / pixelPerDegree; var minDistance = MinLineDistance / pixelPerDegree;
var scale = 1d; var scale = 1d;

View file

@ -56,15 +56,6 @@ namespace MapControl
get { return 90d; } get { return 90d; }
} }
/// <summary>
/// Gets the scale factor from geographic to cartesian coordinates, on the line of true scale of a
/// cylindrical projection (usually the equator), or at the projection center of an azimuthal projection.
/// </summary>
public virtual double UnitsPerDegree
{
get { return Wgs84MetersPerDegree; }
}
/// <summary> /// <summary>
/// Gets the relative map scale at the specified Location. /// Gets the relative map scale at the specified Location.
/// </summary> /// </summary>

View file

@ -29,11 +29,10 @@ namespace MapControl
var lat0 = Center.Latitude * Math.PI / 180d; var lat0 = Center.Latitude * Math.PI / 180d;
var lat = location.Latitude * Math.PI / 180d; var lat = location.Latitude * Math.PI / 180d;
var dLon = (location.Longitude - Center.Longitude) * Math.PI / 180d; var dLon = (location.Longitude - Center.Longitude) * Math.PI / 180d;
var s = UnitsPerDegree * 180d / Math.PI;
return new Point( return new Point(
s * Math.Cos(lat) * Math.Sin(dLon), Wgs84EquatorialRadius * Math.Cos(lat) * Math.Sin(dLon),
s * (Math.Cos(lat0) * Math.Sin(lat) - Math.Sin(lat0) * Math.Cos(lat) * Math.Cos(dLon))); Wgs84EquatorialRadius * (Math.Cos(lat0) * Math.Sin(lat) - Math.Sin(lat0) * Math.Cos(lat) * Math.Cos(dLon)));
} }
public override Location MapToLocation(Point point) public override Location MapToLocation(Point point)
@ -43,9 +42,8 @@ namespace MapControl
return new Location(Center.Latitude, Center.Longitude); return new Location(Center.Latitude, Center.Longitude);
} }
var s = UnitsPerDegree * 180d / Math.PI; var x = point.X / Wgs84EquatorialRadius;
var x = point.X / s; var y = point.Y / Wgs84EquatorialRadius;
var y = point.Y / s;
var r2 = x * x + y * y; var r2 = x * x + y * y;
if (r2 > 1d) if (r2 > 1d)

View file

@ -30,7 +30,7 @@ namespace MapControl
GetAzimuthDistance(Center, location, out azimuth, out distance); GetAzimuthDistance(Center, location, out azimuth, out distance);
var mapDistance = Math.Tan(distance / 2d) * 2d * UnitsPerDegree * 180d / Math.PI; var mapDistance = Math.Tan(distance / 2d) * 2d * Wgs84EquatorialRadius;
return new Point(mapDistance * Math.Sin(azimuth), mapDistance * Math.Cos(azimuth)); return new Point(mapDistance * Math.Sin(azimuth), mapDistance * Math.Cos(azimuth));
} }
@ -45,7 +45,7 @@ namespace MapControl
var azimuth = Math.Atan2(point.X, point.Y); var azimuth = Math.Atan2(point.X, point.Y);
var mapDistance = Math.Sqrt(point.X * point.X + point.Y * point.Y); var mapDistance = Math.Sqrt(point.X * point.X + point.Y * point.Y);
var distance = 2d * Math.Atan(mapDistance / (2d * UnitsPerDegree * 180d / Math.PI)); var distance = 2d * Math.Atan(mapDistance / (2d * Wgs84EquatorialRadius));
return GetLocation(Center, azimuth, distance); return GetLocation(Center, azimuth, distance);
} }

View file

@ -42,15 +42,15 @@ namespace MapControl
public override Point LocationToMap(Location location) public override Point LocationToMap(Location location)
{ {
return new Point( return new Point(
UnitsPerDegree * location.Longitude, Wgs84MetersPerDegree * location.Longitude,
UnitsPerDegree * LatitudeToY(location.Latitude)); Wgs84MetersPerDegree * LatitudeToY(location.Latitude));
} }
public override Location MapToLocation(Point point) public override Location MapToLocation(Point point)
{ {
return new Location( return new Location(
YToLatitude(point.Y / UnitsPerDegree), YToLatitude(point.Y / Wgs84MetersPerDegree),
point.X / UnitsPerDegree); point.X / Wgs84MetersPerDegree);
} }
public static double LatitudeToY(double latitude) public static double LatitudeToY(double latitude)

View file

@ -42,15 +42,15 @@ namespace MapControl
public override Point LocationToMap(Location location) public override Point LocationToMap(Location location)
{ {
return new Point( return new Point(
UnitsPerDegree * location.Longitude, Wgs84MetersPerDegree * location.Longitude,
UnitsPerDegree * LatitudeToY(location.Latitude)); Wgs84MetersPerDegree * LatitudeToY(location.Latitude));
} }
public override Location MapToLocation(Point point) public override Location MapToLocation(Point point)
{ {
return new Location( return new Location(
YToLatitude(point.Y / UnitsPerDegree), YToLatitude(point.Y / Wgs84MetersPerDegree),
point.X / UnitsPerDegree); point.X / Wgs84MetersPerDegree);
} }
public static double LatitudeToY(double latitude) public static double LatitudeToY(double latitude)

View file

@ -110,11 +110,6 @@ namespace MapControl.Projections
get { return isWebMercator; } get { return isWebMercator; }
} }
public override double UnitsPerDegree
{
get { return unitsPerDegree; }
}
public override Point LocationToMap(Location location) public override Point LocationToMap(Location location)
{ {
if (LocationToMapTransform == null) if (LocationToMapTransform == null)

View file

@ -33,11 +33,6 @@ namespace MapControl.Projections
this.falseNorthing = falseNorthing; this.falseNorthing = falseNorthing;
} }
public override double UnitsPerDegree
{
get { return scaleFactor * Wgs84MetersPerDegree; }
}
public override Vector GetRelativeScale(Location location) public override Vector GetRelativeScale(Location location)
{ {
var lat = (north ? location.Latitude : -location.Latitude) * Math.PI / 180d; var lat = (north ? location.Latitude : -location.Latitude) * Math.PI / 180d;