「NTS」を利用して、ポリゴンの中心点を抽出すると3次元に設定される問題

「NTS」を利用して、ポリゴンの中心点を抽出し、shpに保存したところポイントの図形がPOINTではなくPOINTZ(3次元)に設定されていました。

１．原因分析１

ポリゴンから中心点の抽出は以下のプロパティを使っています。PointOnSurfaceは、図形上の
中心点を取得します。それでも取得できない場合は、Centroidを使用します。Centroidはポリゴンの形状により図形上にない場合があります。

var center = feature.Geometry.PointOnSurface;
if(center == null)
{
  center = feature.Geometry.Centroid;
}

var center = feature.Geometry.PointOnSurface;

if(center == null)

{

center = feature.Geometry.Centroid;

}

２．原因分析２

これらのプロパティの処理をを追っていくと、図形生成クラス(GeometryFactory)でポイントを生成するメソッドCreatePoint()を読んでします。

public IPoint CreatePoint(Coordinate coordinate)
{
    return CreatePoint(coordinate != null ? CoordinateSequenceFactory.Create(new[] { coordinate }) : null);
}

public IPoint CreatePoint(Coordinate coordinate)

{

return CreatePoint(coordinate != null ? CoordinateSequenceFactory.Create(new[] { coordinate }) : null);

}

３．原因分析３

さらに追っていくと、CoordinateArraySequenceクラスのコンストラクタで、dimensionが4に設定されていました。
dimensionは次元を意味し、デフォルトの場合は何故か4が設定されています。これが原因でした。

/// dimensionなしのコンストラクタ
public CoordinateArraySequence(Coordinate[] coordinates): this(coordinates, 4) { }

/// dimension付きのコンストラクタ
public CoordinateArraySequence(int size, int dimension)
{
    Coordinates = new Coordinate[size];
    _dimension = dimension;
    for (var i = 0; i < size; i++)
        Coordinates[i] = new Coordinate();
}

/// dimensionなしのコンストラクタ

public CoordinateArraySequence(Coordinate[] coordinates): this(coordinates, 4) { }

/// dimension付きのコンストラクタ

public CoordinateArraySequence(int size, int dimension)

{

Coordinates = new Coordinate[size];

_dimension = dimension;

for (var i = 0; i < size; i++)

Coordinates[i] = new Coordinate();

}

４．対策１

図形生成クラス(GeometryFactory)では、CreatePointを使用する場合は、Coordinateではなく、CoordinateArraySequenceを引数に使用します。CoordinateArraySequenceでは、dimensionの指定ができるので2次元を設定します。

var point = feature.Geometry.Centroid;

// 2次元を指定
var coordinates =  new CoordinateArraySequence(
    new Coordinate[] { new Coordinate(point.X, point.Y) }, 2)

center = GeometryFactory.Default.CreatePoint(coordinates));

var point = feature.Geometry.Centroid;

// 2次元を指定

var coordinates = new CoordinateArraySequence(

new Coordinate[] { new Coordinate(point.X, point.Y) }, 2)

center = GeometryFactory.Default.CreatePoint(coordinates));

５．対策２

もう一つの対策は、WKTを利用します。WKTの場合は、デフォルトの場合でも2次元に設定されます。
以下例

var wktReader = new WKTReader();
var feature = new Feature();
var wkt = new StringBuilder();
wkt.AppendFormat("POINT(");
wkt.AppendFormat("{0} {1})",point.X,point.Y);
feature.Geometry = wktReader.Read(wkt.ToString());