public virtual void Existing_pointZM_should_be_read() { var sequenceFactory = new DotSpatialAffineCoordinateSequenceFactory(Ordinates.XYZM); var coordinateSequence = sequenceFactory.Create( new[] { 11.11, 22.22 }, new[] { 33.33 }, new[] { 44.44 }); var factory = new GeometryFactory(sequenceFactory); var point = factory.CreatePoint(coordinateSequence); Assert.AreEqual(33.33, point.Z); Assert.AreEqual(44.44, point.M); DoTest(conn => { using (var cmd = conn.CreateCommand()) { const string sql = @" DELETE FROM [sample_feature_table_zm]; INSERT INTO [sample_feature_table_zm] ([id], [geometry]) VALUES (1, gpkgMakePointZM(@px, @py, @pz, @pm, 4326));"; cmd.CommandText = sql; cmd.Parameters.AddWithValue("px", point.X); cmd.Parameters.AddWithValue("py", point.Y); cmd.Parameters.AddWithValue("pz", point.Z); cmd.Parameters.AddWithValue("pm", point.M); int ret = cmd.ExecuteNonQuery(); Assert.AreEqual(1, ret); } using (var cmd = conn.CreateCommand()) { const string sql = @" SELECT geometry FROM [sample_feature_table_zm] WHERE [id] = 1;"; cmd.CommandText = sql; using (var rdr = cmd.ExecuteReader()) { Assert.IsTrue(rdr.Read()); byte[] blob = rdr.GetFieldValue <byte[]>(0); var reader = new GeoPackageGeoReader(sequenceFactory, new PrecisionModel(PrecisionModels.Floating)); var geom = reader.Read(blob); Assert.IsNotNull(geom); Assert.IsInstanceOf <IPoint>(geom); Assert.AreEqual(point, geom); Assert.IsTrue(point.EqualsExact(geom)); Assert.AreEqual(point.Coordinate, geom.Coordinate); Assert.IsTrue(point.Coordinate.Equals(geom.Coordinate)); Assert.IsTrue(point.Coordinate.Equals3D(geom.Coordinate)); Assert.AreEqual(point.Z, ((Point)geom).Z); Assert.AreEqual(point.M, ((Point)geom).M); } } }); }
static void Main(string[] args) { //This sequence has to be initialized with ordinates XYZM. If you construct it //with just XYZ values, then you can't add M values to it later. var coordinateSequence = new DotSpatialAffineCoordinateSequence(new[] { 1.0, 2.0 }, new[] { 3.0 }, new[] { 4.0 }); var sequenceFactory = new DotSpatialAffineCoordinateSequenceFactory(GeoAPI.Geometries.Ordinates.XYZM); var geometryFactory = new GeometryFactory(sequenceFactory); var p3 = new Point(coordinateSequence, geometryFactory); p3.Y = 8; p3.M = 1; var m = p3.M; Debug.WriteLine(p3); Debug.WriteLine(p3.M); Debug.WriteLine(m); }
public virtual void New_pointZM_should_be_written() { var sequenceFactory = new DotSpatialAffineCoordinateSequenceFactory(Ordinates.XYZM); var coordinateSequence = sequenceFactory.Create( new[] { 11.11, 22.22 }, new[] { 33.33 }, new[] { 44.44 }); var factory = new GeometryFactory(sequenceFactory); var point = factory.CreatePoint(coordinateSequence); Assert.AreEqual(33.33, point.Z); Assert.AreEqual(44.44, point.M); DoTest(conn => { using (var cmd = conn.CreateCommand()) { const string sql = @" DELETE FROM [sample_feature_table_zm]; INSERT INTO [sample_feature_table_zm] ([id], [geometry]) VALUES (1, @pbytes);"; cmd.CommandText = sql; var writer = new GeoPackageGeoWriter { HandleOrdinates = Ordinates.XYZM }; byte[] bytes = writer.Write(point); cmd.Parameters.AddWithValue("pbytes", bytes); int ret = cmd.ExecuteNonQuery(); Assert.AreEqual(1, ret); } using (var cmd = conn.CreateCommand()) { const string sql = @" SELECT geometry ,ST_Is3D(geometry) ,ST_IsMeasured(geometry) ,ST_AsText(GeomFromGPB(geometry)) FROM [sample_feature_table_zm] WHERE [id] = 1;"; cmd.CommandText = sql; using (var rdr = cmd.ExecuteReader()) { Assert.IsTrue(rdr.Read()); byte[] blob = rdr.GetFieldValue <byte[]>(0); var reader = new GeoPackageGeoReader(sequenceFactory, new PrecisionModel(PrecisionModels.Floating)); var geom = reader.Read(blob); Assert.AreEqual(point, geom); Assert.AreEqual(point.AsBinary(), geom.AsBinary()); Assert.IsTrue(rdr.GetBoolean(1)); Assert.IsTrue(rdr.GetBoolean(2)); string wkt = rdr.GetString(3); Assert.IsFalse(string.IsNullOrEmpty(wkt)); Assert.AreEqual("POINT ZM(11.11 22.22 33.33 44.44)", wkt); } } }); }