/// <summary>
///
/// </summary>
/// <param name="tokenizer"></param>
/// <returns></returns>
private static IProjectedCoordinateSystem ReadProjectedCoordinateSystem(WktStreamTokenizer tokenizer)
{
/*PROJCS[
* "OSGB 1936 / British National Grid",
* GEOGCS[
* "OSGB 1936",
* DATUM[...]
* PRIMEM[...]
* AXIS["Geodetic latitude","NORTH"]
* AXIS["Geodetic longitude","EAST"]
* AUTHORITY["EPSG","4277"]
* ],
* PROJECTION["Transverse Mercator"],
* PARAMETER["latitude_of_natural_origin",49],
* PARAMETER["longitude_of_natural_origin",-2],
* PARAMETER["scale_factor_at_natural_origin",0.999601272],
* PARAMETER["false_easting",400000],
* PARAMETER["false_northing",-100000],
* AXIS["Easting","EAST"],
* AXIS["Northing","NORTH"],
* AUTHORITY["EPSG","27700"]
* ]
*/
tokenizer.ReadToken("[");
string name = tokenizer.ReadDoubleQuotedWord();
tokenizer.ReadToken(",");
tokenizer.ReadToken("GEOGCS");
IGeographicCoordinateSystem geographicCS = ReadGeographicCoordinateSystem(tokenizer);
tokenizer.ReadToken(",");
IProjection projection = null;
IUnit unit = null;
List <AxisInfo> axes = new List <AxisInfo>(2);
string authority = String.Empty;
long authorityCode = -1;
TokenType ct = tokenizer.NextToken();
while (ct != TokenType.Eol && ct != TokenType.Eof)
{
switch (tokenizer.GetStringValue())
{
case ",":
case "]":
break;
case "PROJECTION":
projection = ReadProjection(tokenizer);
ct = tokenizer.GetTokenType();
continue;
//break;
case "UNIT":
unit = ReadLinearUnit(tokenizer);
break;
case "AXIS":
axes.Add(ReadAxis(tokenizer));
tokenizer.NextToken();
break;
case "AUTHORITY":
tokenizer.ReadAuthority(ref authority, ref authorityCode);
//tokenizer.ReadToken("]");
break;
}
ct = tokenizer.NextToken();
}
//This is default axis values if not specified.
if (axes.Count == 0)
{
axes.Add(new AxisInfo("X", AxisOrientationEnum.East));
axes.Add(new AxisInfo("Y", AxisOrientationEnum.North));
}
IProjectedCoordinateSystem projectedCS = new ProjectedCoordinateSystem(geographicCS.HorizontalDatum, geographicCS, unit as LinearUnit, projection, axes, name, authority, authorityCode, String.Empty, String.Empty, String.Empty);
return(projectedCS);
}
/// <summary>
///
/// </summary>
/// <param name="tokenizer"></param>
/// <returns></returns>
private static IProjectedCoordinateSystem ReadProjectedCoordinateSystem(WktStreamTokenizer tokenizer)
{
/*PROJCS[
"OSGB 1936 / British National Grid",
GEOGCS[
"OSGB 1936",
DATUM[...]
PRIMEM[...]
AXIS["Geodetic latitude","NORTH"]
AXIS["Geodetic longitude","EAST"]
AUTHORITY["EPSG","4277"]
],
PROJECTION["Transverse Mercator"],
PARAMETER["latitude_of_natural_origin",49],
PARAMETER["longitude_of_natural_origin",-2],
PARAMETER["scale_factor_at_natural_origin",0.999601272],
PARAMETER["false_easting",400000],
PARAMETER["false_northing",-100000],
AXIS["Easting","EAST"],
AXIS["Northing","NORTH"],
AUTHORITY["EPSG","27700"]
]
*/
tokenizer.ReadToken("[");
string name = tokenizer.ReadDoubleQuotedWord();
tokenizer.ReadToken(",");
tokenizer.ReadToken("GEOGCS");
IGeographicCoordinateSystem geographicCS = ReadGeographicCoordinateSystem(tokenizer);
tokenizer.ReadToken(",");
IProjection projection = null;
IUnit unit = null;
List<AxisInfo> axes = new List<AxisInfo>(2);
string authority = String.Empty;
long authorityCode = -1;
TokenType ct = tokenizer.NextToken();
while (ct != TokenType.Eol && ct != TokenType.Eof)
{
switch (tokenizer.GetStringValue())
{
case ",":
case "]":
break;
case "PROJECTION":
projection = ReadProjection(tokenizer);
ct = tokenizer.GetTokenType();
continue;
//break;
case "UNIT":
unit = ReadLinearUnit(tokenizer);
break;
case "AXIS":
axes.Add(ReadAxis(tokenizer));
tokenizer.NextToken();
break;
case "AUTHORITY":
tokenizer.ReadAuthority(ref authority, ref authorityCode);
//tokenizer.ReadToken("]");
break;
}
ct = tokenizer.NextToken();
}
//This is default axis values if not specified.
if (axes.Count == 0)
{
axes.Add(new AxisInfo("X", AxisOrientationEnum.East));
axes.Add(new AxisInfo("Y", AxisOrientationEnum.North));
}
IProjectedCoordinateSystem projectedCS = new ProjectedCoordinateSystem(geographicCS.HorizontalDatum, geographicCS, unit as LinearUnit, projection, axes, name, authority, authorityCode, String.Empty, String.Empty, String.Empty);
return projectedCS;
}
