private static Ellipsoid ReadEllipsoid(WktStreamTokenizer tokenizer, bool includeAuthority)
{
//SPHEROID["Airy 1830",6377563.396,299.3249646,AUTHORITY["EPSG","7001"]]
tokenizer.ReadToken("[");
string name = tokenizer.ReadDoubleQuotedWord();
tokenizer.ReadToken(",");
tokenizer.NextToken();
double majorAxis = tokenizer.GetNumericValue();
tokenizer.ReadToken(",");
tokenizer.NextToken();
double e = tokenizer.GetNumericValue();
string authority = String.Empty;
string authorityCode = String.Empty;
if (includeAuthority)
{
tokenizer.ReadToken(",");
tokenizer.ReadAuthority(ref authority, ref authorityCode);
}
tokenizer.ReadToken("]");
Ellipsoid ellipsoid = new Ellipsoid(majorAxis, 0.0, e, true, LinearUnit.Meters, String.Empty, authority, authorityCode, name, String.Empty, String.Empty);
return(ellipsoid);
}
private static Projection ReadProjection(WktStreamTokenizer tokenizer, bool includeAuthority)
{
//tokenizer.NextToken();// PROJECTION
tokenizer.ReadToken("PROJECTION");
tokenizer.ReadToken("[");//[
string projectionName = tokenizer.ReadDoubleQuotedWord();
tokenizer.ReadToken("]"); //]
tokenizer.ReadToken(","); //,
tokenizer.ReadToken("PARAMETER");
ParameterList paramList = new ParameterList();
while (tokenizer.GetStringValue() == "PARAMETER")
{
tokenizer.ReadToken("[");
string paramName = tokenizer.ReadDoubleQuotedWord();
//added by monoGIS team: parameters names may be capitalized, but the code works only with lower-case!
paramName = paramName.ToLower();
tokenizer.ReadToken(",");
tokenizer.NextToken();
double paramValue = tokenizer.GetNumericValue();
tokenizer.ReadToken("]");
paramList.Add(paramName, paramValue);
if (!tokenizer.TryReadToken(","))
{
break;
}
if (!tokenizer.TryReadToken("PARAMETER"))
{
break;
}
}
ProjectionParameter[] paramArray = new ProjectionParameter[paramList.Count];
int i = 0;
foreach (string key in paramList.Keys)
{
ProjectionParameter param = new ProjectionParameter();
param.Name = key;
param.Value = (double)paramList[key];
paramArray[i] = param;
i++;
}
string authority = String.Empty;
string authorityCode = String.Empty;
Projection projection = new Projection(projectionName, paramArray, String.Empty, String.Empty, authority, authorityCode);
return(projection);
}
private static PrimeMeridian ReadPrimeMeridian(WktStreamTokenizer tokenizer, bool includeAuthority)
{
//PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]]
tokenizer.ReadToken("[");
string name = tokenizer.ReadDoubleQuotedWord();
tokenizer.ReadToken(",");
tokenizer.NextToken();
double longitude = tokenizer.GetNumericValue();
string authority = String.Empty;
string authorityCode = String.Empty;
if (includeAuthority)
{
tokenizer.ReadToken(",");
tokenizer.ReadAuthority(ref authority, ref authorityCode);
}
// make an assumption about the Angular units - degrees.
PrimeMeridian primeMeridian = new PrimeMeridian(name, new AngularUnit(180 / Math.PI), longitude, String.Empty, authority, authorityCode, String.Empty, String.Empty);
tokenizer.ReadToken("]");
return(primeMeridian);
}
private static WGS84ConversionInfo ReadWGS84ConversionInfo(WktStreamTokenizer tokenizer, bool includeAuthority)
{
//TOWGS84[0,0,0,0,0,0,0]
tokenizer.ReadToken("[");
WGS84ConversionInfo info = new WGS84ConversionInfo();
tokenizer.NextToken();
info.Dx = tokenizer.GetNumericValue();
tokenizer.ReadToken(",");
tokenizer.NextToken();
info.Dy = tokenizer.GetNumericValue();
tokenizer.ReadToken(",");
tokenizer.NextToken();
info.Dz = tokenizer.GetNumericValue();
tokenizer.ReadToken(",");
tokenizer.NextToken();
info.Ex = tokenizer.GetNumericValue();
tokenizer.ReadToken(",");
tokenizer.NextToken();
info.Ey = tokenizer.GetNumericValue();
tokenizer.ReadToken(",");
tokenizer.NextToken();
info.Ez = tokenizer.GetNumericValue();
tokenizer.ReadToken(",");
tokenizer.NextToken();
info.Ppm = tokenizer.GetNumericValue();
tokenizer.ReadToken("]");
info.IsInUse = true;
return(info);
}
private static Unit ReadUnit(WktStreamTokenizer tokenizer, bool includeAuthority)
{
//UNIT["degree",0.01745329251994433,AUTHORITY["EPSG","9102"]]
Unit unit = null;
tokenizer.ReadToken("[");
string unitName = tokenizer.ReadDoubleQuotedWord();
tokenizer.ReadToken(",");
tokenizer.NextToken();
double unitsPerUnit = tokenizer.GetNumericValue();
string authority = String.Empty;
string authorityCode = String.Empty;
if (includeAuthority)
{
tokenizer.ReadToken(",");
tokenizer.ReadAuthority(ref authority, ref authorityCode);
}
tokenizer.ReadToken("]");
switch (unitName.ToUpper())
{
// take into account the different spellings of the word meter/metre.
case "METER":
case "METRE":
case "KILOMETRE":
case "KILOMETER":
unit = new LinearUnit(unitsPerUnit, String.Empty, authority, authorityCode, unitName, String.Empty, String.Empty);
break;
case "DEGREE":
case "RADIAN":
case "GRAD":
unit = new AngularUnit(unitsPerUnit, String.Empty, authority, authorityCode, unitName, String.Empty, String.Empty);
break;
//case "CLARKE'S LINK":
//case "GOLD COAST FOOT":
//case "US SURVEY FOOT":
//case "CLARKE'S FOOT":
//case "FOOT":
//case "LINK":
//case "INDIAN YARD":
//case "GERMAN LEGAL METRE":
//case "BRITISH CHAIN (SEARS 1922)":
//case "BRITISH FOOT (SEARS 1922)":
//case "BRITISH CHAIN (SEARS 1922 TRUNCATED)":
//case "BRITISH CHAIN (BENOIT 1922 TRUNCATED)":
// unit = new LinearUnit(unitsPerUnit, String.Empty, authority, authorityCode, unitName, String.Empty, String.Empty);
// break;
default:
string u = unitName.ToUpper();
if (u.Contains("YARD") || u.Contains("CHAIN") ||
u.Contains("FOOT") || u.Contains("LINK") ||
u.Contains("METRE") || u.Contains("METER") ||
u.Contains("KILOMETER") || u.Contains("KILOMETRE"))
{
unit = new LinearUnit(unitsPerUnit, String.Empty, authority, authorityCode, unitName, String.Empty, String.Empty);
}
else
{
throw new NotImplementedException(String.Format("{0} is not recognized a unit of measure.", unitName));
}
break;
}
return(unit);
}
private static Geotransformation ReadGeotransformation(WktStreamTokenizer tokenizer, bool includeAuthority)
{
// GEOGTRAN[
// "MGI_To_WGS_1984_2",
// GEOGCS[
// "GCS_MGI",
// DATUM["D_MGI",SPHEROID["Bessel_1841",6377397.155,299.1528128]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],
// METHOD["Position_Vector"],
// PARAMETER["X_Axis_Translation",577.326],
// PARAMETER["Y_Axis_Translation",90.129],
// PARAMETER["Z_Axis_Translation",463.919],
// PARAMETER["X_Axis_Rotation",5.1365988],
// PARAMETER["Y_Axis_Rotation",1.4742],
// PARAMETER["Z_Axis_Rotation",5.2970436],
// PARAMETER["Scale_Difference",2.4232]
//]
tokenizer.ReadToken("[");
string name = tokenizer.ReadDoubleQuotedWord();
tokenizer.ReadToken(",");
tokenizer.ReadToken("GEOGCS");
GeographicCoordinateSystem geographicCS1 = WktCoordinateSystemReader.ReadGeographicCoordinateSystem(tokenizer, includeAuthority);
tokenizer.ReadToken(",");
tokenizer.ReadToken("GEOGCS");
GeographicCoordinateSystem geographicCS2 = WktCoordinateSystemReader.ReadGeographicCoordinateSystem(tokenizer, includeAuthority);
tokenizer.ReadToken(",");
string method = String.Empty;
ParameterList paramList = new ParameterList();
if (tokenizer.TryReadToken("METHOD"))
{
tokenizer.ReadToken("[");
method = tokenizer.ReadDoubleQuotedWord();
tokenizer.ReadToken("]");
}
if (tokenizer.TryReadToken(","))
{
tokenizer.ReadToken("PARAMETER");
while (tokenizer.GetStringValue() == "PARAMETER")
{
tokenizer.ReadToken("[");
string paramName = tokenizer.ReadDoubleQuotedWord();
paramName = paramName.ToLower();
tokenizer.ReadToken(",");
tokenizer.NextToken();
double paramValue = tokenizer.GetNumericValue();
tokenizer.ReadToken("]");
paramList.Add(paramName, paramValue);
if (!tokenizer.TryReadToken(","))
{
break;
}
if (!tokenizer.TryReadToken("PARAMETER"))
{
break;
}
}
}
return(new Geotransformation(name, method, geographicCS1, geographicCS2, paramList, "", "", "", ""));
}
