.Netでシェープファイルの属性の型を変更したい。

前提・実現したいこと

.Net Frameworkでシェープファイルの属性の型を変更したい。

以下の質問でpythonで実施していたことを、.Net Frameworkでやるにはどうしたらいいかという質問です。
シェープファイルのdbfのみ修正する方法

元のシェープファイルを変更するか、新たに作成するのでも構いません。

想定している修正内容は
・実数型の小数点以下桁数の変更
・文字列型の桁数の変更
・列順序の変更
・列名の変更

やり方や使えるライブラリなど教えて欲しいです。

試したこと

型を変更しないで内容をコピーは
DotSpatialをNugetでインストールして以下で出来ました。

C#
1            var fromlayer = DotSpatial.Data.DataManager.DefaultDataManager.OpenFile(@"D:\study\sample.shp")
2                as DotSpatial.Data.FeatureSet;
3
4            var tolayer = DotSpatial.Data.DataManager.DefaultDataManager.CreateVector(@"D:\study\sampleto.shp",fromlayer.FeatureType)
5                as DotSpatial.Data.FeatureSet;
6
7            foreach (System.Data.DataColumn column in fromlayer.DataTable.Columns) {
8                tolayer.DataTable.Columns.Add(column.Caption,column.DataType);
9            }
10
11            foreach (System.Data.DataRow row in fromlayer.DataTable.Rows)
12            {
13                System.Data.DataRow newrow = tolayer.DataTable.NewRow();
14                newrow.ItemArray = row.ItemArray;
15                tolayer.DataTable.Rows.Add(newrow);
16            }
17
18            tolayer.Save();

ただ、型の変更が出来ません。
列名の変更くらいなら、出来そうな気もしますが
小数点以下桁数の変更や文字列型の桁数の変更は、どうやればいいかわかりません。

追記

VB.NETでgdal/ogrを使ってみる。 [Chapter 1]

を参考にして
・実数型の小数点以下桁数の変更
・文字列型の桁数の変更
・列名の変更
は以下のようにして出来たのではないかと思います。

VBNet
1Imports System.IO
2Imports OSGeo.OGR
3
4Public Class Form1
5    ' この部分はリンク先をほぼコピペです。 START
6    Private Function SetupGDALOGR() As Boolean
7
8        SetupGDALOGR = False
9
10        Dim strNewPath As String = Environment.GetEnvironmentVariable("PATH")
11        Dim gdalPath As String = "C:\Program Files"
12        Dim tmpPath As String = "c:\tmp"
13        Dim GDAL_CACHEMAX As String = "100000"
14
15        Dim info As New FileInfo(gdalPath)
16        If (info.Attributes And FileAttributes.Directory) > 0 Then
17            Dim oldPath As [String] = Environment.GetEnvironmentVariable("PATH")
18            Dim PATH2 As [String] = gdalPath & "\GDAL;" &
19                                      gdalPath & "\GDAL\gdalplugins;" &
20                                      gdalPath & "\GDAL\csharp;"
21
22            If Not strNewPath.Contains(PATH2) Then
23                strNewPath = PATH2 & ";" & strNewPath
24                Environment.SetEnvironmentVariable("PATH", strNewPath)
25            End If
26
27            Dim GDAL_DATA As [String] = gdalPath & "\GDAL\gdal-data"
28            Environment.SetEnvironmentVariable("\GDAL_DATA", GDAL_DATA)
29            OSGeo.GDAL.Gdal.SetConfigOption("\GDAL_DATA", GDAL_DATA)
30
31            Dim GDAL_DRIVER_PATH As [String] = gdalPath & "\GDAL\gdalplugins"
32            Environment.SetEnvironmentVariable("GDAL_DRIVER_PATH", GDAL_DRIVER_PATH)
33            OSGeo.GDAL.Gdal.SetConfigOption("GDAL_DRIVER_PATH", GDAL_DRIVER_PATH)
34
35
36            OSGeo.OGR.Ogr.RegisterAll()
37            OSGeo.GDAL.Gdal.AllRegister()
38
39            SetupGDALOGR = True
40
41        End If
42
43    End Function
44    ' この部分はリンク先をほぼコピペです。 END
45
46    Private Sub Form1_Load(sender As Object, e As EventArgs) Handles MyBase.Load
47        SetupGDALOGR()
48    End Sub
49
50    Private Sub Button1_Click(sender As Object, e As EventArgs) Handles Button1.Click
51
52        Dim p_Ds_From As DataSource = Ogr.Open("C:\Work\Shape\sample.shp", 0)
53        Dim p_Layer_From As Layer = p_Ds_From.GetLayerByIndex(0)
54
55        Dim p_Drv As Driver = Ogr.GetDriverByName("ESRI Shapefile")
56        Dim p_Ds_To As DataSource = p_Drv.CreateDataSource("C:\Work\Shape\sample2.shp", Nothing)
57
58        Dim p_Layer_To As Layer = p_Ds_To.CreateLayer("layer", p_Layer_From.GetSpatialRef, p_Layer_From.GetGeomType, Nothing)
59
60
61        Dim p_Layer_Defn_From As FeatureDefn = p_Layer_From.GetLayerDefn()
62
63
64        For p_Ind As Long = 0 To p_Layer_Defn_From.GetFieldCount - 1
65            Dim p_Field_Defn As FieldDefn
66            p_Field_Defn = p_Layer_Defn_From.GetFieldDefn(p_Ind)
67            ' ここで列の定義を変更。名前の変更や文字桁数の変更、小数点以下桁数の変更が実現出来ました。
68            If p_Field_Defn.GetName = "XXXX" Then
69                p_Field_Defn.SetWidth(10)
70            ElseIf p_Field_Defn.GetName = "XXXX2" Then
71                p_Field_Defn.SetName("YYYY2")
72                p_Field_Defn.SetWidth(10)
73                p_Field_Defn.SetPrecision(1)
74
75            End If
76
77            p_Layer_To.CreateField(p_Field_Defn, 1)
78
79
80        Next
81
82        Dim p_Feature_From As Feature
83        p_Feature_From = p_Layer_From.GetNextFeature()
84
85
86        Do While p_Feature_From IsNot Nothing
87            p_Layer_To.CreateFeature(p_Feature_From)
88            p_Feature_From = p_Layer_From.GetNextFeature()
89        Loop
90
91        p_Ds_From.Dispose()
92        p_Ds_To.Dispose()
93
94    End Sub
95End Class
96

尚、環境構築としては、上記のリンクに書いてある通り
gisinternals Support Site
のHPにあるgdal-303-1928-x64-core.msiをインストールして
参照設定にgdal_csharp.dll,ogr_csharp.dll,osr_csharp.dllを設定しました。

後は列順序の変更なのですが、どのようにするかまだ不明です。
gdalでやる方法や、他の方法をご存知であれば教えて下さい。

行動規範の内容に同意します

回答1件

自己解決

順序の変更も出来たので自己解決とします。

追記の通り
FieldDefnを変更することで、フィールド名、桁数、小数点以下桁数は変更できました。

あと新たに作成した、Layerに対してCreateFieldの順番を変えるとフィールドの順番の変更になります。

データの移行の所はSetFromWithMapを使用することで、属性の入れ替えが出来ました。

以下に追記と似たような感じですが、サンプルコードを記載しておきます。

VBNET
1    ' SetupGDALOGR部分は省略
2
3    Private Sub Button1_Click(sender As Object, e As EventArgs) Handles Button1.Click
4
5        Dim p_Ds_From As DataSource = Ogr.Open("C:\Work\Shape\sample.shp", 0)
6        Dim p_Layer_From As Layer = p_Ds_From.GetLayerByIndex(0)
7
8        Dim p_Drv As Driver = Ogr.GetDriverByName("ESRI Shapefile")
9        Dim p_Ds_To As DataSource = p_Drv.CreateDataSource("C:\Work\Shape\sample2.shp", Nothing)
10
11        Dim p_Layer_To As Layer = p_Ds_To.CreateLayer("layer", p_Layer_From.GetSpatialRef, p_Layer_From.GetGeomType, Nothing)
12
13
14        Dim p_Layer_Defn_From As FeatureDefn = p_Layer_From.GetLayerDefn()
15
16        Dim p_Dict_Field_Defn As New Dictionary(Of String, FieldDefn)
17
18        For p_Ind As Long = p_Layer_Defn_From.GetFieldCount - 1 To 0 Step -1
19            Dim p_Field_Defn As FieldDefn
20            p_Field_Defn = p_Layer_Defn_From.GetFieldDefn(p_Ind)
21
22            If p_Field_Defn.GetName = "XXXX" Then
23                p_Field_Defn.SetType(FieldType.OFTInteger)
24                p_Field_Defn.SetWidth(10)
25            ElseIf p_Field_Defn.GetName = "YYYY" Then
26                p_Field_Defn.SetName("YYYY2")
27                p_Field_Defn.SetWidth(10)
28                p_Field_Defn.SetPrecision(1)
29            Else
30                p_Field_Defn.SetWidth(10)
31            End If
32
33            p_Dict_Field_Defn.Add(p_Field_Defn.GetName, p_Field_Defn)
34
35        Next
36
37        p_Layer_To.CreateField(p_Dict_Field_Defn.Item("XXXXX"), 0)
38        p_Layer_To.CreateField(p_Dict_Field_Defn.Item("YYYY2"), 0)
39        p_Layer_To.CreateField(p_Dict_Field_Defn.Item("test1"), 0)
40        p_Layer_To.CreateField(p_Dict_Field_Defn.Item("test2"), 0)
41
42        Dim p_Feature_From As Feature
43        p_Feature_From = p_Layer_From.GetNextFeature()
44
45        Dim p_Feature_To As Feature
46        p_Feature_To = New Feature(p_Layer_To.GetLayerDefn())
47
48
49        Do While p_Feature_From IsNot Nothing
50            p_Feature_To.SetFromWithMap(p_Feature_From, 0, 4, {0, 2, 3, 1})
51            p_Layer_To.CreateFeature(p_Feature_To)
52            p_Feature_From = p_Layer_From.GetNextFeature()
53        Loop
54
55        p_Ds_From.Dispose()
56        p_Ds_To.Dispose()
57
58    End Sub