private static void ApplyParametersDefault(ProjectionInfo result, IEnumerable <INamedParameter> parameters)
{
Contract.Requires(result != null);
Contract.Requires(parameters != null);
var lon0Param = new MultiParameterSelector(
new CentralMeridianParameterSelector(),
new LongitudeOfNaturalOriginParameterSelector()
);
var loncParam = new LongitudeOfCenterParameterSelector();
var lat0Param = new MultiParameterSelector(
new LatitudeOfCenterParameterSelector(),
new LatitudeOfNaturalOriginParameterSelector());
var lat1Param = new StandardParallelParameterSelector(1);
var lat2Param = new StandardParallelParameterSelector(2);
var x0Param = new FalseEastingParameterSelector();
var y0Param = new FalseNorthingParameterSelector();
var k0Param = new ScaleFactorParameterSelector();
var alphaParam = new KeywordNamedParameterSelector("ALPHA", "AZIMUTH");
var angleSkewParam = new KeywordNamedParameterSelector("SKEW", "GAMMA");
var zoneParam = new FullMatchParameterSelector("ZONE");
var southParam = new FullMatchParameterSelector("SOUTH");
var standardParallel1Param = new MultiParameterSelector(
new StandardParallelParameterSelector(),
new StandardParallelParameterSelector(1));
var standardParallel2Param = new StandardParallelParameterSelector(2);
var paramLookup = new NamedParameterLookup(parameters);
paramLookup.Assign(lon0Param, loncParam, lat0Param, lat1Param, lat2Param, x0Param, y0Param, k0Param, alphaParam, angleSkewParam, zoneParam, southParam, standardParallel1Param, standardParallel2Param);
if (lon0Param.IsSelected)
{
result.CentralMeridian = lon0Param.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
}
if (loncParam.IsSelected)
{
result.LongitudeOfCenter = loncParam.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
}
if (lat0Param.IsSelected)
{
result.LatitudeOfOrigin = lat0Param.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
}
if (lat1Param.IsSelected)
{
result.StandardParallel1 = lat1Param.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
}
if (lat2Param.IsSelected)
{
result.StandardParallel2 = lat2Param.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
}
if (x0Param.IsSelected)
{
result.FalseEasting = x0Param.GetValueAsDouble(OgcLinearUnit.DefaultMeter);
}
if (y0Param.IsSelected)
{
result.FalseNorthing = y0Param.GetValueAsDouble(OgcLinearUnit.DefaultMeter);
}
if (k0Param.IsSelected)
{
result.ScaleFactor = k0Param.GetValueAsDouble(ScaleUnitUnity.Value) ?? 1.0;
}
if (alphaParam.IsSelected)
{
result.alpha = alphaParam.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
}
/*if (angleSkewParam.IsSelected)
* result.gamma = angleSkewParam.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);*/
if (southParam.IsSelected)
{
result.IsSouth = southParam.GetValueAsBoolean().GetValueOrDefault();
}
if (standardParallel1Param.IsSelected)
{
result.StandardParallel1 = standardParallel1Param.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
}
if (standardParallel2Param.IsSelected)
{
result.StandardParallel2 = standardParallel2Param.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
}
if (zoneParam.IsSelected)
{
result.Zone = zoneParam.GetValueAsInt32();
}
}
private SpheroidProjectionBase CreatePolarStereographic(ProjectionCompilationParams opData)
{
Contract.Requires(opData != null);
var latParam = new MultiParameterSelector(
new LatitudeOfNaturalOriginParameterSelector(),
new LatitudeOfTrueScaleParameterSelector());
var originLonParam = new MultiParameterSelector(
new LongitudeOfNaturalOriginParameterSelector(),
new FullMatchParameterSelector("LONDOWNFROMPOLE"));// new KeywordNamedParameterSelector("LON", "ORIGIN");
var offsetXParam = new FalseEastingParameterSelector();
var offsetYParam = new FalseNorthingParameterSelector();
var standardParallelParam = new StandardParallelParameterSelector();
var scaleFactorParam = new ScaleFactorParameterSelector();
opData.ParameterLookup.Assign(latParam, originLonParam, offsetXParam, offsetYParam, standardParallelParam, scaleFactorParam);
var linearUnit = opData.StepParams.RelatedOutputCrsUnit;
var fromSpheroid = opData.StepParams.ConvertRelatedOutputSpheroidUnit(linearUnit);
if (null == fromSpheroid)
return null;
double latitude, longitude;
if (!latParam.TryGetValueAsDouble(OgcAngularUnit.DefaultRadians, out latitude))
latitude = 0;
if (!originLonParam.TryGetValueAsDouble(OgcAngularUnit.DefaultRadians, out longitude))
longitude = 0;
double xOffset, yOffset;
var offset = (offsetXParam.TryGetValueAsDouble(linearUnit, out xOffset) && offsetYParam.TryGetValueAsDouble(linearUnit, out yOffset))
? new Vector2(xOffset, yOffset)
: Vector2.Zero;
double scaleFactor;
if (scaleFactorParam.TryGetValueAsDouble(ScaleUnitUnity.Value, out scaleFactor))
return new PolarStereographic(new GeographicCoordinate(latitude,longitude), scaleFactor, offset, fromSpheroid);
double standardParallel;
if (standardParallelParam.TryGetValueAsDouble(OgcAngularUnit.DefaultRadians, out standardParallel)) {
if (
(offsetXParam.IsSelected && opData.ParameterLookup.ParameterNameNormalizedComparer.Normalize(offsetXParam.Selection.Name).Contains("ORIGIN"))
|| (offsetYParam.IsSelected && opData.ParameterLookup.ParameterNameNormalizedComparer.Normalize(offsetYParam.Selection.Name).Contains("ORIGIN"))
) {
return PolarStereographic.CreateFromStandardParallelAndFalseOffsetAtOrigin(longitude, standardParallel, offset, fromSpheroid);
}
return PolarStereographic.CreateFromStandardParallel(longitude, standardParallel, offset, fromSpheroid);
}
throw new InvalidOperationException();
}
private static void ApplyParametersDefault(ProjectionInfo result, IEnumerable<INamedParameter> parameters)
{
Contract.Requires(result != null);
Contract.Requires(parameters != null);
var lon0Param = new MultiParameterSelector(
new CentralMeridianParameterSelector(),
new LongitudeOfNaturalOriginParameterSelector()
);
var loncParam = new LongitudeOfCenterParameterSelector();
var lat0Param = new MultiParameterSelector(
new LatitudeOfCenterParameterSelector(),
new LatitudeOfNaturalOriginParameterSelector());
var lat1Param = new StandardParallelParameterSelector(1);
var lat2Param = new StandardParallelParameterSelector(2);
var x0Param = new FalseEastingParameterSelector();
var y0Param = new FalseNorthingParameterSelector();
var k0Param = new ScaleFactorParameterSelector();
var alphaParam = new KeywordNamedParameterSelector("ALPHA", "AZIMUTH");
var angleSkewParam = new KeywordNamedParameterSelector("SKEW", "GAMMA");
var zoneParam = new FullMatchParameterSelector("ZONE");
var southParam = new FullMatchParameterSelector("SOUTH");
var standardParallel1Param = new MultiParameterSelector(
new StandardParallelParameterSelector(),
new StandardParallelParameterSelector(1));
var standardParallel2Param = new StandardParallelParameterSelector(2);
var paramLookup = new NamedParameterLookup(parameters);
paramLookup.Assign(lon0Param, loncParam, lat0Param, lat1Param, lat2Param, x0Param, y0Param, k0Param, alphaParam, angleSkewParam, zoneParam, southParam, standardParallel1Param, standardParallel2Param);
if (lon0Param.IsSelected)
result.CentralMeridian = lon0Param.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
if (loncParam.IsSelected)
result.LongitudeOfCenter = loncParam.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
if (lat0Param.IsSelected)
result.LatitudeOfOrigin = lat0Param.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
if (lat1Param.IsSelected)
result.StandardParallel1 = lat1Param.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
if (lat2Param.IsSelected)
result.StandardParallel2 = lat2Param.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
if (x0Param.IsSelected)
result.FalseEasting = x0Param.GetValueAsDouble(OgcLinearUnit.DefaultMeter);
if (y0Param.IsSelected)
result.FalseNorthing = y0Param.GetValueAsDouble(OgcLinearUnit.DefaultMeter);
if (k0Param.IsSelected)
result.ScaleFactor = k0Param.GetValueAsDouble(ScaleUnitUnity.Value) ?? 1.0;
if (alphaParam.IsSelected)
result.alpha = alphaParam.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
/*if (angleSkewParam.IsSelected)
result.gamma = angleSkewParam.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);*/
if (southParam.IsSelected)
result.IsSouth = southParam.GetValueAsBoolean().GetValueOrDefault();
if (standardParallel1Param.IsSelected)
result.StandardParallel1 = standardParallel1Param.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
if (standardParallel2Param.IsSelected)
result.StandardParallel2 = standardParallel2Param.GetValueAsDouble(OgcAngularUnit.DefaultDegrees);
if (zoneParam.IsSelected)
result.Zone = zoneParam.GetValueAsInt32();
}
private AlbersEqualArea CreateAlbersEqualArea(ProjectionCompilationParams opData)
{
Contract.Requires(opData != null);
var originLatParam = new MultiParameterSelector(
new LatitudeOfFalseOrigin(),
new LatitudeOfNaturalOriginParameterSelector());
var originLonParam = new MultiParameterSelector(
new LongitudeOfFalseOrigin(),
new CentralMeridianParameterSelector(),
new LongitudeOfNaturalOriginParameterSelector(),
new LongitudeOfCenterParameterSelector());
var offsetXParam = new FalseEastingParameterSelector();
var offsetYParam = new FalseNorthingParameterSelector();
var lat1SpParam = new StandardParallelParameterSelector(1);
var lat2SpParam = new StandardParallelParameterSelector(2);
opData.ParameterLookup.Assign(originLatParam, originLonParam, offsetXParam, offsetYParam, lat1SpParam, lat2SpParam);
var fromSpheroid = opData.StepParams.ConvertRelatedOutputSpheroidUnit(opData.StepParams.RelatedOutputCrsUnit);
if (null == fromSpheroid)
return null;
GeographicCoordinate origin;
if (!(originLatParam.IsSelected || originLonParam.IsSelected) || !TryCreateGeographicCoordinate(originLatParam.Selection, originLonParam.Selection, out origin))
origin = GeographicCoordinate.Zero;
Vector2 offset;
if (!(offsetXParam.IsSelected || offsetYParam.IsSelected) || !TryCreateVector2(offsetXParam.Selection, offsetYParam.Selection, opData.StepParams.RelatedOutputCrsUnit, out offset))
offset = Vector2.Zero;
double parallel1;
if (!lat1SpParam.IsSelected || !TryGetDouble(lat1SpParam.Selection, OgcAngularUnit.DefaultRadians, out parallel1))
parallel1 = 0;
double parallel2;
if (!lat2SpParam.IsSelected || !TryGetDouble(lat2SpParam.Selection, OgcAngularUnit.DefaultRadians, out parallel2))
parallel2 = 0;
return new AlbersEqualArea(origin, parallel1, parallel2, offset, fromSpheroid);
}
