public void TestReadValueZFromLineStringZShapeFile()
{
var file = GetTestDataFilePath("LineString_With_Z.shp");
var sh = new SharpMap.Data.Providers.ShapeFile(file, true);
sh.Open();
var fc = sh.GetFeatureCount();
NUnit.Framework.Assert.AreEqual(1221, fc);
NUnit.Framework.Assert.AreEqual(0, sh.GetObjectIDsInView(sh.GetExtents())[0]);
var featsInView = sh.GetGeometriesInView(new GeoAPI.Geometries.Envelope(sh.GetExtents()));
NUnit.Framework.Assert.AreEqual(1221, featsInView.Count);
foreach (var item in featsInView)
{
NUnit.Framework.Assert.IsNotNull(item.Coordinate.Z);
}
NUnit.Framework.Assert.AreEqual(featsInView[0].Coordinates[0].Z, 35.865);
NUnit.Framework.Assert.AreEqual(featsInView[0].Coordinates[1].Z, 35.743);
NUnit.Framework.Assert.AreEqual(featsInView[1].Coordinates[0].Z, 35.518);
NUnit.Framework.Assert.AreEqual(featsInView[1].Coordinates[1].Z, 35.518);
NUnit.Framework.Assert.AreEqual(featsInView[2].Coordinates[0].Z, 37.438);
NUnit.Framework.Assert.AreEqual(featsInView[2].Coordinates[1].Z, 37.441);
NUnit.Framework.Assert.AreEqual(featsInView[3].Coordinates[0].Z, 37.441);
NUnit.Framework.Assert.AreEqual(featsInView[3].Coordinates[1].Z, 37.441);
sh.Close();
}
public void TestReadValueZFromPointZShapeFile()
{
var file = GetTestDataFilePath("Point_With_Z.shp");
var sh = new SharpMap.Data.Providers.ShapeFile(file, true);
sh.Open();
var fc = sh.GetFeatureCount();
NUnit.Framework.Assert.AreEqual(1149, fc);
NUnit.Framework.Assert.AreEqual(0, sh.GetObjectIDsInView(sh.GetExtents())[0]);
var featsInView = sh.GetGeometriesInView(new GeoAPI.Geometries.Envelope(sh.GetExtents()));
NUnit.Framework.Assert.AreEqual(1149, featsInView.Count);
foreach (var item in featsInView)
{
NUnit.Framework.Assert.IsNotNull(item.Coordinate.Z);
}
NUnit.Framework.Assert.AreEqual(featsInView[0].Coordinate.Z, 146.473);
NUnit.Framework.Assert.AreEqual(featsInView[1].Coordinate.Z, 181.374);
NUnit.Framework.Assert.AreEqual(featsInView[2].Coordinate.Z, 146.676);
NUnit.Framework.Assert.AreEqual(featsInView[3].Coordinate.Z, 181.087);
NUnit.Framework.Assert.AreEqual(featsInView[4].Coordinate.Z, 169.948);
NUnit.Framework.Assert.AreEqual(featsInView[5].Coordinate.Z, 169.916);
sh.Close();
}
public void SetupFixture()
{
SqlConnectionStringBuilder connStrBuilder = new SqlConnectionStringBuilder(UnitTests.Properties.Settings.Default.SqlServer2008);
if (string.IsNullOrEmpty(connStrBuilder.DataSource) || string.IsNullOrEmpty(connStrBuilder.InitialCatalog))
{
Assert.Ignore("Requires SQL Server connectionstring");
}
GeoAPI.GeometryServiceProvider.Instance = new NetTopologySuite.NtsGeometryServices();
// Set up sample table
using (SqlConnection conn = new SqlConnection(UnitTests.Properties.Settings.Default.SqlServer2008))
{
conn.Open();
using (SqlCommand cmd = conn.CreateCommand())
{
// The ID column cannot simply be int, because that would cause GetObjectIDsInView to fail. The provider internally works with uint
cmd.CommandText = "CREATE TABLE roads_ugl(ID decimal(10,0) identity(1,1) PRIMARY KEY, NAME nvarchar(100), GEOM geometry)";
cmd.ExecuteNonQuery();
}
// Load data
using (SharpMap.Data.Providers.ShapeFile shapeFile = new SharpMap.Data.Providers.ShapeFile(GetTestFile()))
{
shapeFile.Open();
IEnumerable <uint> indexes = shapeFile.GetObjectIDsInView(shapeFile.GetExtents());
indexes = indexes.Take(100);
foreach (uint idx in indexes)
{
SharpMap.Data.FeatureDataRow feature = shapeFile.GetFeature(idx);
using (SqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText = "INSERT INTO roads_ugl(NAME, GEOM) VALUES (@Name, geometry::STGeomFromText(@Geom, @Srid))";
cmd.Parameters.AddWithValue("@Geom", feature.Geometry.AsText());
cmd.Parameters.AddWithValue("@Name", feature["NAME"]);
cmd.Parameters.AddWithValue("@Srid", shapeFile.SRID);
cmd.ExecuteNonQuery();
}
}
}
// Create spatial index
using (SqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText = "CREATE SPATIAL INDEX [IX_roads_ugl_GEOM] ON [dbo].[roads_ugl](GEOM)USING GEOMETRY_GRID WITH (BOUNDING_BOX =(-98, 40, -82, 50), GRIDS =(LEVEL_1 = MEDIUM,LEVEL_2 = MEDIUM,LEVEL_3 = MEDIUM,LEVEL_4 = MEDIUM))";
cmd.ExecuteNonQuery();
}
}
}
public void SetupFixture()
{
SqlConnectionStringBuilder connStrBuilder = new SqlConnectionStringBuilder(UnitTests.Properties.Settings.Default.SqlServer2008);
if (string.IsNullOrEmpty(connStrBuilder.DataSource) || string.IsNullOrEmpty(connStrBuilder.InitialCatalog))
{
Assert.Ignore("Requires SQL Server connectionstring");
}
GeoAPI.GeometryServiceProvider.Instance = new NetTopologySuite.NtsGeometryServices();
// Set up sample table
using (SqlConnection conn = new SqlConnection(UnitTests.Properties.Settings.Default.SqlServer2008))
{
conn.Open();
using(SqlCommand cmd = conn.CreateCommand())
{
// The ID column cannot simply be int, because that would cause GetObjectIDsInView to fail. The provider internally works with uint
cmd.CommandText = "CREATE TABLE roads_ugl(ID decimal(10,0) identity(1,1) PRIMARY KEY, NAME nvarchar(100), GEOM geometry)";
cmd.ExecuteNonQuery();
}
// Load data
using (SharpMap.Data.Providers.ShapeFile shapeFile = new SharpMap.Data.Providers.ShapeFile(GetTestFile()))
{
shapeFile.Open();
IEnumerable<uint> indexes = shapeFile.GetObjectIDsInView(shapeFile.GetExtents());
indexes = indexes.Take(100);
foreach (uint idx in indexes)
{
SharpMap.Data.FeatureDataRow feature = shapeFile.GetFeature(idx);
using (SqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText = "INSERT INTO roads_ugl(NAME, GEOM) VALUES (@Name, geometry::STGeomFromText(@Geom, @Srid))";
cmd.Parameters.AddWithValue("@Geom", feature.Geometry.AsText());
cmd.Parameters.AddWithValue("@Name", feature["NAME"]);
cmd.Parameters.AddWithValue("@Srid", shapeFile.SRID);
cmd.ExecuteNonQuery();
}
}
}
// Create spatial index
using(SqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText = "CREATE SPATIAL INDEX [IX_roads_ugl_GEOM] ON [dbo].[roads_ugl](GEOM)USING GEOMETRY_GRID WITH (BOUNDING_BOX =(-98, 40, -82, 50), GRIDS =(LEVEL_1 = MEDIUM,LEVEL_2 = MEDIUM,LEVEL_3 = MEDIUM,LEVEL_4 = MEDIUM))";
cmd.ExecuteNonQuery();
}
}
}
public void TestReadPointZShapeFile()
{
var file = GetTestDataFilePath("SPATIAL_F_SKARVMUFF.shp");
var sh = new SharpMap.Data.Providers.ShapeFile(file, true);
sh.Open();
var fc = sh.GetFeatureCount();
NUnit.Framework.Assert.AreEqual(4342, fc);
NUnit.Framework.Assert.AreEqual(0, sh.GetObjectIDsInView(sh.GetExtents())[0]);
var featsInView = sh.GetGeometriesInView(new GeoAPI.Geometries.Envelope(sh.GetExtents()));
NUnit.Framework.Assert.AreEqual(4342, featsInView.Count);
sh.Close();
}
public void OneTimeSetUp()
{
try
{
GeoAPI.GeometryServiceProvider.Instance = new NetTopologySuite.NtsGeometryServices();
var connStrBuilder = new NpgsqlConnectionStringBuilder(Properties.Settings.Default.PostGis);
if (string.IsNullOrEmpty(connStrBuilder.Host) || string.IsNullOrEmpty(connStrBuilder.Database))
{
Assert.Ignore("Requires PostgreSQL connectionstring");
}
// Set up sample table
using (var conn = new NpgsqlConnection(Properties.Settings.Default.PostGis))
{
conn.Open();
// Load data
using (var shapeFile = new SharpMap.Data.Providers.ShapeFile(TestUtility.GetPathToTestFile("roads_ugl.shp"), false, false, 4326))
{
shapeFile.Open();
using (var cmd = conn.CreateCommand())
{
cmd.CommandText =
"SELECT COUNT(*) FROM \"geometry_columns\" WHERE \"f_table_name\" = 'roads_ugl_g' AND \"f_geometry_column\"='geom';";
var count = cmd.ExecuteScalar();
if (Convert.ToInt32(count) > 0)
{
cmd.CommandText = "SELECT DropGeometryColumn('roads_ugl_g', 'geom');";
cmd.ExecuteNonQuery();
cmd.CommandText = "DROP TABLE roads_ugl_g";
cmd.ExecuteNonQuery();
}
// The ID column cannot simply be int, because that would cause GetObjectIDsInView to fail. The provider internally works with uint
cmd.CommandText =
"CREATE TABLE roads_ugl_g(id integer primary key, name character varying(100));";
cmd.ExecuteNonQuery();
cmd.CommandText = "SELECT AddGeometryColumn('roads_ugl_g', 'geom', " + shapeFile.SRID +
", 'GEOMETRY', 2);";
cmd.ExecuteNonQuery();
}
IEnumerable <uint> indexes = shapeFile.GetObjectIDsInView(shapeFile.GetExtents());
_insertedIds = new List <uint>(indexes.Take(100));
using (NpgsqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText =
"INSERT INTO roads_ugl_g(id, name, geom) VALUES (@PId, @PName, @PGeom);";
var @params = cmd.Parameters;
@params.AddRange(
new[]
{
new NpgsqlParameter("PId", NpgsqlDbType.Integer),
new NpgsqlParameter("PName", NpgsqlDbType.Varchar, 100),
new NpgsqlParameter("PGeom", NpgsqlDbType.Bytea)
});
var writer = new PostGisWriter();
foreach (var idx in _insertedIds)
{
var feature = shapeFile.GetFeature(idx);
@params["PId"].NpgsqlValue = (int)idx;
@params["PName"].NpgsqlValue = feature["NAME"];
@params["PGeom"].NpgsqlValue = writer.Write(feature.Geometry);
cmd.ExecuteNonQuery();
}
}
}
}
}
catch
{
Assert.Ignore("Failed to connect to PostgreSQL/PostGIS Server");
}
}
public void FixtureSetup()
{
var connStrBuilder = new NpgsqlConnectionStringBuilder(Properties.Settings.Default.PostGis);
if (string.IsNullOrEmpty(connStrBuilder.Host) || string.IsNullOrEmpty(connStrBuilder.Database))
{
Assert.Ignore("Requires PostgreSQL connectionstring");
}
GeoAPI.GeometryServiceProvider.Instance = new NetTopologySuite.NtsGeometryServices();
try
{
// Set up sample table
using (var conn = new NpgsqlConnection(Properties.Settings.Default.PostGis))
{
conn.Open();
// Load data
using (var shapeFile = new SharpMap.Data.Providers.ShapeFile(GetTestFile()))
{
shapeFile.SRID = 4326;
shapeFile.Open();
using (var cmd = conn.CreateCommand())
{
cmd.CommandText = "DROP TABLE IF EXISTS roads_ugl";
cmd.ExecuteNonQuery();
// The ID column cannot simply be int, because that would cause GetObjectIDsInView to fail. The provider internally works with uint
cmd.CommandText =
"CREATE TABLE roads_ugl(id integer primary key, name character varying(100), geog geography);";
cmd.ExecuteNonQuery();
}
IEnumerable<uint> indexes = shapeFile.GetObjectIDsInView(shapeFile.GetExtents());
_insertedIds = new List<uint>(indexes.Take(100));
using (NpgsqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText =
"INSERT INTO roads_ugl(id, name, geog) VALUES (@PId, @PName, ST_GeogFromWKB(@PGeom));";
var @params = cmd.Parameters;
@params.AddRange(
new[]
{
new NpgsqlParameter("PId", NpgsqlDbType.Integer),
new NpgsqlParameter("PName", NpgsqlDbType.Varchar, 100),
new NpgsqlParameter("PGeom", NpgsqlDbType.Bytea)
});
var writer = new PostGisWriter();
foreach (var idx in _insertedIds)
{
var feature = shapeFile.GetFeature(idx);
@params["PId"].NpgsqlValue = (int)idx;
@params["PName"].NpgsqlValue = feature["NAME"];
@params["PGeom"].NpgsqlValue = writer.Write(feature.Geometry);
cmd.ExecuteNonQuery();
}
}
// Verify
foreach (var pgp in GetTestProvider())
{
foreach (var idx in _insertedIds)
{
var g1 = pgp.GetGeometryByID(idx);
var g2 = shapeFile.GetGeometryByID(idx);
Assert.AreEqual(g1, g2);
}
}
}
}
}
catch (Exception ee)
{
Assert.Ignore("Failed to connect to PostgreSQL/PostGIS Server");
}
}
public void FixtureSetup()
{
var connStrBuilder = new NpgsqlConnectionStringBuilder(Properties.Settings.Default.PostGis);
if (string.IsNullOrEmpty(connStrBuilder.Host) || string.IsNullOrEmpty(connStrBuilder.Database))
{
Assert.Ignore("Requires PostgreSQL connectionstring");
}
GeoAPI.GeometryServiceProvider.Instance = new NetTopologySuite.NtsGeometryServices();
try
{
// Set up sample table
using (var conn = new NpgsqlConnection(Properties.Settings.Default.PostGis))
{
conn.Open();
// Load data
using (var shapeFile = new SharpMap.Data.Providers.ShapeFile(GetTestFile(), false, false, 4326))
{
shapeFile.Open();
using (var cmd = conn.CreateCommand())
{
cmd.CommandText = "DROP TABLE IF EXISTS roads_ugl";
cmd.ExecuteNonQuery();
// The ID column cannot simply be int, because that would cause GetObjectIDsInView to fail. The provider internally works with uint
cmd.CommandText =
"CREATE TABLE roads_ugl(id integer primary key, name character varying(100), geog geography);";
cmd.ExecuteNonQuery();
}
IEnumerable <uint> indexes = shapeFile.GetObjectIDsInView(shapeFile.GetExtents());
_insertedIds = new List <uint>(indexes.Take(100));
using (NpgsqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText =
"INSERT INTO roads_ugl(id, name, geog) VALUES (@PId, @PName, ST_GeogFromWKB(@PGeom));";
var @params = cmd.Parameters;
@params.AddRange(
new[]
{
new NpgsqlParameter("PId", NpgsqlDbType.Integer),
new NpgsqlParameter("PName", NpgsqlDbType.Varchar, 100),
new NpgsqlParameter("PGeom", NpgsqlDbType.Bytea)
});
var writer = new PostGisWriter();
foreach (var idx in _insertedIds)
{
var feature = shapeFile.GetFeature(idx);
@params["PId"].NpgsqlValue = (int)idx;
@params["PName"].NpgsqlValue = feature["NAME"];
@params["PGeom"].NpgsqlValue = writer.Write(feature.Geometry);
cmd.ExecuteNonQuery();
}
}
// Verify
foreach (var pgp in GetTestProvider())
{
foreach (var idx in _insertedIds)
{
var g1 = pgp.GetGeometryByID(idx);
var g2 = shapeFile.GetGeometryByID(idx);
Assert.AreEqual(g1, g2);
}
}
}
}
}
catch (Exception ee)
{
Assert.Ignore("Failed to connect to PostgreSQL/PostGIS Server");
}
}
public void OneTimeSetUp()
{
SqlConnectionStringBuilder connStrBuilder = new SqlConnectionStringBuilder(Properties.Settings.Default.SqlServer2008);
if (string.IsNullOrEmpty(connStrBuilder.DataSource) || string.IsNullOrEmpty(connStrBuilder.InitialCatalog))
{
NUnit.Framework.Assert.Ignore("Requires SQL Server connectionstring");
}
GeoAPI.GeometryServiceProvider.Instance = new NetTopologySuite.NtsGeometryServices();
//SqlServerTypes.Utilities.LoadNativeAssemblies(AppDomain.CurrentDomain.BaseDirectory);
// Set up sample tables (Geometry + Geography)
using (SqlConnection conn = new SqlConnection(UnitTests.Properties.Settings.Default.SqlServer2008))
{
conn.Open();
using (SqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText = "CREATE TABLE roads_ugl_geom(ID int identity(1,1) PRIMARY KEY, NAME nvarchar(100), GEOM geometry)";
cmd.ExecuteNonQuery();
}
using (SqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText = "CREATE TABLE roads_ugl_geog(ID int identity(1,1) PRIMARY KEY, NAME nvarchar(100), GEOG geography)";
cmd.ExecuteNonQuery();
}
// Load data
using (SharpMap.Data.Providers.ShapeFile shapeFile = new SharpMap.Data.Providers.ShapeFile(TestUtility.GetPathToTestFile("roads.ugl")))
{
shapeFile.Open();
_geometrySrid = shapeFile.SRID;
IEnumerable <uint> indexes = shapeFile.GetObjectIDsInView(shapeFile.GetExtents());
// Note that spatial indexes may only kick in at certain number of records
// so for thorough testing comment out next line and load all features (approx 3500)
indexes = indexes.Take(100);
var cmdGeom = new SqlCommand("INSERT INTO roads_ugl_geom(NAME, GEOM) VALUES (@Name, geometry::STGeomFromText(@Geom, @Srid))", conn);
var cmdGeog = new SqlCommand("INSERT INTO roads_ugl_geog(NAME, GEOG) VALUES (@Name, geography::STGeomFromText(@Geog, @Srid))", conn);
foreach (uint idx in indexes)
{
SharpMap.Data.FeatureDataRow feature = shapeFile.GetFeature(idx);
string wkt;
if (feature.Geometry == null || feature.Geometry.IsEmpty)
{
wkt = "LINESTRING EMPTY";
}
else
{
wkt = feature.Geometry.AsText();
}
if (cmdGeom.Parameters.Count == 0)
{
cmdGeom.Parameters.AddWithValue("@Geom", wkt);
cmdGeom.Parameters.AddWithValue("@Name", feature["NAME"]);
cmdGeom.Parameters.AddWithValue("@Srid", _geometrySrid);
}
else
{
cmdGeom.Parameters[0].Value = wkt;
cmdGeom.Parameters[1].Value = feature["NAME"];
}
cmdGeom.ExecuteNonQuery();
if (cmdGeog.Parameters.Count == 0)
{
cmdGeog.Parameters.AddWithValue("@Geog", wkt);
cmdGeog.Parameters.AddWithValue("@Name", feature["NAME"]);
cmdGeog.Parameters.AddWithValue("@Srid", GeographySrid);
}
else
{
cmdGeog.Parameters[0].Value = wkt;
cmdGeog.Parameters[1].Value = feature["NAME"];
}
cmdGeog.ExecuteNonQuery();
}
cmdGeom.Dispose();
cmdGeog.Dispose();
}
// ensure we have some features with NULL and EMPTY geometries
using (var cmd = conn.CreateCommand())
{
// To find invalid geometries:
// SELECT {OidColumn}, {GeometryColumn}.STIsValid() AS STIsValid, {GeometryColumn}.IsValidDetailed() AS IsValidDetailed FROM {QualifiedTableName}
// NULL
cmd.CommandText = "INSERT INTO roads_ugl_geom(NAME, GEOM) VALUES ('Test null wkt', NULL)";
cmd.ExecuteNonQuery();
cmd.CommandText = "INSERT INTO roads_ugl_geog(NAME, GEOG) VALUES ('Test null wkt', NULL)";
cmd.ExecuteNonQuery();
// EMPTY
cmd.CommandText = "INSERT INTO roads_ugl_geom(NAME, GEOM) VALUES ('Test empty wkt', 'LINESTRING EMPTY')";
cmd.ExecuteNonQuery();
cmd.CommandText = "INSERT INTO roads_ugl_geog(NAME, GEOG) VALUES ('Test empty wkt', 'LINESTRING EMPTY')";
cmd.ExecuteNonQuery();
// INVALID (ID 144 from shape file; see also ID 2055)
cmd.CommandText = $"INSERT INTO roads_ugl_geom(NAME, GEOM) VALUES ('Test invalid wkt', geometry::STGeomFromText('LINESTRING (-84.652756071629071 42.676743004284312, -84.652924071615374 42.676624004283632, -84.652756071629071 42.676743004284312, -84.652512071649028 42.676922004285323, -84.641022072594438 42.685478004332808, -84.638779072781034 42.687271004342172, -84.636932072941363 42.689831004350026, -84.634491073153043 42.693100004360424, -84.62387107404335 42.701092004405112, -84.603256075794022 42.715752004493233, -84.603142075803831 42.715832004493734, -84.599823076091937 42.718651004508146, -84.588676077031693 42.722431004556235, -84.586021077270672 42.725533004568049)', {_geometrySrid}))";
cmd.ExecuteNonQuery();
cmd.CommandText = $"INSERT INTO roads_ugl_geog(NAME, GEOG) VALUES ('Test invalid wkt', geography::STGeomFromText('LINESTRING (-84.652756071629071 42.676743004284312, -84.652924071615374 42.676624004283632, -84.652756071629071 42.676743004284312, -84.652512071649028 42.676922004285323, -84.641022072594438 42.685478004332808, -84.638779072781034 42.687271004342172, -84.636932072941363 42.689831004350026, -84.634491073153043 42.693100004360424, -84.62387107404335 42.701092004405112, -84.603256075794022 42.715752004493233, -84.603142075803831 42.715832004493734, -84.599823076091937 42.718651004508146, -84.588676077031693 42.722431004556235, -84.586021077270672 42.725533004568049)', {GeographySrid}))";
cmd.ExecuteNonQuery();
}
// Create GEOM spatial index
using (SqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText = $"CREATE SPATIAL INDEX [{GeometrySpatialIndex}] ON [dbo].[roads_ugl_geom](GEOM) USING GEOMETRY_GRID WITH (BOUNDING_BOX =(-98, 40, -82, 50), GRIDS =(LEVEL_1 = MEDIUM,LEVEL_2 = MEDIUM,LEVEL_3 = MEDIUM,LEVEL_4 = MEDIUM))";
cmd.CommandTimeout = 300;
cmd.ExecuteNonQuery();
}
// Create GEOG spatial index
using (SqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText = $"CREATE SPATIAL INDEX [{GeographySpatialIndex}] ON [dbo].[roads_ugl_geog](GEOG)";
cmd.CommandTimeout = 300;
cmd.ExecuteNonQuery();
}
// initialise counts and test IDs
using (SqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText = $"SELECT COUNT(ID) FROM roads_ugl_geom WHERE Geom.STIsEmpty() = 0 AND Geom.STIsValid() = 1";
_numValidGeoms = (int)cmd.ExecuteScalar();
cmd.CommandText = $"SELECT COUNT(ID) FROM roads_ugl_geom WHERE GEOM IS NOT NULL AND Geom.STIsEmpty() = 0 AND Geom.STIsValid() = 0";
_numInvalidGeoms = (int)cmd.ExecuteScalar();
_numValidatedGeoms = _numValidGeoms + _numInvalidGeoms;
cmd.CommandText = $"SELECT COUNT(ID) FROM roads_ugl_geom";
_numFeatures = (int)cmd.ExecuteScalar();
_idNullGeom = (uint)(_numFeatures - 2);
_idEmptyGeom = (uint)(_numFeatures - 1);
_idInvalidGeom = (uint)(_numFeatures);
}
}
}
public void FixtureSetup()
{
var connStrBuilder = new NpgsqlConnectionStringBuilder(Properties.Settings.Default.PostGis);
if (string.IsNullOrEmpty(connStrBuilder.Host) || string.IsNullOrEmpty(connStrBuilder.Database))
{
Assert.Ignore("Requires PostgreSQL connectionstring");
}
GeoAPI.GeometryServiceProvider.Instance = new NetTopologySuite.NtsGeometryServices();
try
{
// Set up sample table
using (var conn = new NpgsqlConnection(Properties.Settings.Default.PostGis))
{
conn.Open();
// Load data
using (var shapeFile = new SharpMap.Data.Providers.ShapeFile(GetTestFile(), false, false, 4326))
{
shapeFile.Open();
using (var cmd = conn.CreateCommand())
{
cmd.CommandText =
"SELECT COUNT(*) FROM \"geometry_columns\" WHERE \"f_table_name\" = 'roads_ugl_g' AND \"f_geometry_column\"='geom';";
var count = cmd.ExecuteScalar();
if (Convert.ToInt32(count) > 0)
{
cmd.CommandText = "SELECT DropGeometryColumn('roads_ugl_g', 'geom');";
cmd.ExecuteNonQuery();
cmd.CommandText = "DROP TABLE roads_ugl_g";
cmd.ExecuteNonQuery();
}
// The ID column cannot simply be int, because that would cause GetObjectIDsInView to fail. The provider internally works with uint
cmd.CommandText =
"CREATE TABLE roads_ugl_g(id integer primary key, name character varying(100));";
cmd.ExecuteNonQuery();
cmd.CommandText = "SELECT AddGeometryColumn('roads_ugl_g', 'geom', " + shapeFile.SRID +
", 'GEOMETRY', 2);";
cmd.ExecuteNonQuery();
}
IEnumerable<uint> indexes = shapeFile.GetObjectIDsInView(shapeFile.GetExtents());
_insertedIds = new List<uint>(indexes.Take(100));
using (NpgsqlCommand cmd = conn.CreateCommand())
{
cmd.CommandText =
"INSERT INTO roads_ugl_g(id, name, geom) VALUES (@PId, @PName, @PGeom);";
var @params = cmd.Parameters;
@params.AddRange(
new[]
{
new NpgsqlParameter("PId", NpgsqlDbType.Integer),
new NpgsqlParameter("PName", NpgsqlDbType.Varchar, 100),
new NpgsqlParameter("PGeom", NpgsqlDbType.Bytea)
});
var writer = new PostGisWriter();
foreach (var idx in _insertedIds)
{
var feature = shapeFile.GetFeature(idx);
@params["PId"].NpgsqlValue = (int) idx;
@params["PName"].NpgsqlValue = feature["NAME"];
@params["PGeom"].NpgsqlValue = writer.Write(feature.Geometry);
cmd.ExecuteNonQuery();
}
}
}
}
}
catch
{
Assert.Ignore("Failed to connect to PostgreSQL/PostGIS Server");
}
}
public void TestReadPointZShapeFile()
{
var file = GetTestDataFilePath("SPATIAL_F_SKARVMUFF.shp");
var sh = new SharpMap.Data.Providers.ShapeFile(file, true);
sh.Open();
var fc = sh.GetFeatureCount();
NUnit.Framework.Assert.AreEqual(4342, fc);
NUnit.Framework.Assert.AreEqual(0, sh.GetObjectIDsInView(sh.GetExtents())[0]);
var featsInView = sh.GetGeometriesInView(new GeoAPI.Geometries.Envelope(sh.GetExtents()));
NUnit.Framework.Assert.AreEqual(4342, featsInView.Count);
sh.Close();
}