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);
}
}
});
}
