/// <summary> /// /// </summary> /// <param name="coordinateSystem"></param> /// <param name="tokenizer"></param> /// <returns></returns> private static GeoAPI.CoordinateSystems.ICoordinateSystem ReadCoordinateSystem(string coordinateSystem, GisSharpBlog.NetTopologySuite.IO.WktStreamTokenizer tokenizer) { GeoAPI.CoordinateSystems.ICoordinateSystem returnCS = null; switch (coordinateSystem) { /*case "VERT_CS": * IVerticalCoordinateSystem verticalCS = ReadVerticalCoordinateSystem(tokenizer); * returnCS = verticalCS; * break;*/ case "GEOGCS": GeoAPI.CoordinateSystems.IGeographicCoordinateSystem geographicCS = ReadGeographicCoordinateSystem(tokenizer); returnCS = geographicCS; break; case "PROJCS": GeoAPI.CoordinateSystems.IProjectedCoordinateSystem projectedCS = ReadProjectedCoordinateSystem(tokenizer); returnCS = projectedCS; break; case "COMPD_CS": /* ICompoundCoordinateSystem compoundCS = ReadCompoundCoordinateSystem(tokenizer); * returnCS = compoundCS; * break;*/ case "GEOCCS": case "FITTED_CS": case "LOCAL_CS": throw new InvalidOperationException(String.Format("{0} coordinate system is not recongized.", coordinateSystem)); } return(returnCS); }
private void Initialize() { GisSharpBlog.NetTopologySuite.Geometries.PrecisionModel precisionModel = new GisSharpBlog.NetTopologySuite.Geometries.PrecisionModel(GeoAPI.Geometries.PrecisionModels.Floating); SharpMap.CoordinateSystems.CoordinateSystem wgs84 = SharpMap.CoordinateSystems.GeographicCoordinateSystem.WGS84; GeoAPI.CoordinateSystems.ICoordinateSystemFactory cFac = new SharpMap.CoordinateSystems.CoordinateSystemFactory(); //Create geographic coordinate system based on the WGS84 datum GeoAPI.CoordinateSystems.IEllipsoid ellipsoid = cFac.CreateFlattenedSphere("WGS 84", 6378137, 298.257223563, SharpMap.CoordinateSystems.LinearUnit.Metre); GeoAPI.CoordinateSystems.IHorizontalDatum datum = cFac.CreateHorizontalDatum("WGS_1984", GeoAPI.CoordinateSystems.DatumType.HD_Geocentric, ellipsoid, null); GeoAPI.CoordinateSystems.IGeographicCoordinateSystem gcsWGS84 = cFac.CreateGeographicCoordinateSystem("WGS 84", SharpMap.CoordinateSystems.AngularUnit.Degrees, datum, SharpMap.CoordinateSystems.PrimeMeridian.Greenwich, new GeoAPI.CoordinateSystems.AxisInfo("Lon", GeoAPI.CoordinateSystems.AxisOrientationEnum.East), new GeoAPI.CoordinateSystems.AxisInfo("Lat", GeoAPI.CoordinateSystems.AxisOrientationEnum.North)); List <GeoAPI.CoordinateSystems.ProjectionParameter> parameters = new List <GeoAPI.CoordinateSystems.ProjectionParameter>(); parameters.Add(new GeoAPI.CoordinateSystems.ProjectionParameter("latitude_of_origin", 0)); parameters.Add(new GeoAPI.CoordinateSystems.ProjectionParameter("central_meridian", 117)); parameters.Add(new GeoAPI.CoordinateSystems.ProjectionParameter("scale_factor", 0.9996)); parameters.Add(new GeoAPI.CoordinateSystems.ProjectionParameter("false_easting", 500000)); parameters.Add(new GeoAPI.CoordinateSystems.ProjectionParameter("false_northing", 10000000)); GeoAPI.CoordinateSystems.IProjection projection = cFac.CreateProjection("Transverse Mercator", "Transverse Mercator", parameters); GeoAPI.CoordinateSystems.IProjectedCoordinateSystem utmWGS84 = cFac.CreateProjectedCoordinateSystem("WGS84 UTM Zone 50S", gcsWGS84, projection, SharpMap.CoordinateSystems.LinearUnit.Metre, new GeoAPI.CoordinateSystems.AxisInfo("East", GeoAPI.CoordinateSystems.AxisOrientationEnum.East), new GeoAPI.CoordinateSystems.AxisInfo("North", GeoAPI.CoordinateSystems.AxisOrientationEnum.North)); int SRID_utm50 = Convert.ToInt32(utmWGS84.AuthorityCode); //UTM50 SRID SharpMap.CoordinateSystems.Transformations.CoordinateTransformationFactory ctFact = new SharpMap.CoordinateSystems.Transformations.CoordinateTransformationFactory(); _Transform = ctFact.CreateFromCoordinateSystems(utmWGS84, wgs84); _UTM50Factory = new GisSharpBlog.NetTopologySuite.Geometries.GeometryFactory(precisionModel, SRID_utm50); }