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 <Point>(geom);
Assert.AreEqual(point, geom);
Assert.IsTrue(point.EqualsExact(geom));
Assert.AreEqual(point.Coordinate, geom.Coordinate);
Assert.IsTrue(point.Coordinate.Equals(geom.Coordinate));
Assert.IsInstanceOf <CoordinateZ>(geom.Coordinate);
var coordZ = (CoordinateZ)geom.Coordinate;
Assert.IsTrue(((CoordinateZ)point.Coordinate).Equals3D(coordZ));
Assert.AreEqual(point.Z, ((Point)geom).Z);
Assert.AreEqual(point.M, ((Point)geom).M);
}
}
});
}
public virtual void New_pointZ_should_be_written()
{
var coord = new Coordinate(11.11, 22.22, 33.33);
var point = GeometryFactory.Default.CreatePoint(coord);
DoTest(conn =>
{
using (var cmd = conn.CreateCommand())
{
const string sql = @"
DELETE FROM [sample_feature_table_z];
INSERT INTO [sample_feature_table_z] ([id], [geometry])
VALUES (1, @pbytes);";
cmd.CommandText = sql;
var writer = new GeoPackageGeoWriter
{
HandleOrdinates = Ordinates.XYZ
};
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_z]
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
{
HandleOrdinates = Ordinates.XYZ
};
var geom = reader.Read(blob);
Assert.AreEqual(point, geom);
Assert.AreEqual(point.AsBinary(), geom.AsBinary());
Assert.IsTrue(rdr.GetBoolean(1));
Assert.IsFalse(rdr.GetBoolean(2));
string wkt = rdr.GetString(3);
Assert.IsFalse(string.IsNullOrEmpty(wkt));
Assert.AreEqual("POINT Z(11.11 22.22 33.33)", wkt);
}
}
});
}
public virtual void Srid_should_be_read_only_if_configured()
{
var coord = new Coordinate(11.11, 22.22);
var point = GeometryFactory.Default.CreatePoint(coord);
point.SRID = 3004;
DoTest(conn =>
{
using (var cmd = conn.CreateCommand())
{
const string sql = @"
DELETE FROM [sample_feature_table];
INSERT INTO [sample_feature_table] ([id], [geometry])
VALUES (1, @pbytes);";
cmd.CommandText = sql;
var writer = new GeoPackageGeoWriter
{
HandleOrdinates = Ordinates.XY
};
Assert.IsTrue(writer.HandleSRID);
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]
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
{
HandleOrdinates = Ordinates.XY,
HandleSRID = true
};
var geom = reader.Read(blob);
Assert.AreEqual(point, geom);
Assert.AreEqual(point.SRID, geom.SRID);
}
}
});
}
public virtual void Existing_pointZ_should_be_read()
{
var coord = new CoordinateZ(11.11, 22.22, 33.33);
var point = GeometryFactory.Default.CreatePoint(coord);
DoTest(conn =>
{
using (var cmd = conn.CreateCommand())
{
const string sql = @"
DELETE FROM [sample_feature_table_z];
INSERT INTO [sample_feature_table_z] ([id], [geometry])
VALUES (1, gpkgMakePointZ(@px, @py, @pz, 4326));";
cmd.CommandText = sql;
cmd.Parameters.AddWithValue("px", coord.X);
cmd.Parameters.AddWithValue("py", coord.Y);
cmd.Parameters.AddWithValue("pz", coord.Z);
int ret = cmd.ExecuteNonQuery();
Assert.AreEqual(1, ret);
}
using (var cmd = conn.CreateCommand())
{
const string sql = @"
SELECT geometry
FROM [sample_feature_table_z]
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();
var geom = reader.Read(blob);
Assert.IsNotNull(geom);
Assert.IsInstanceOf <Point>(geom);
Assert.AreEqual(point, geom);
Assert.IsTrue(point.EqualsExact(geom));
Assert.AreEqual(coord, geom.Coordinate);
Assert.IsTrue(coord.Equals(geom.Coordinate));
Assert.IsInstanceOf <CoordinateZ>(geom.Coordinate);
var coordZ = (CoordinateZ)geom.Coordinate;
Assert.IsTrue(coord.Equals3D(coordZ));
Assert.AreEqual(coord.Z, geom.Coordinate.Z);
}
}
});
}
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);
}
}
});
}
protected override Geometry Read(byte[] bytes)
{
var reader = new GeoPackageGeoReader();
return(reader.Read(bytes));
}