FAQ/ 左岸線・右岸線のshapefileを取り込みたい¶

河川の左岸線・右岸線を、shapefile に取り込む方法を教えてください。

DioVistaでは、shapefileファイルは表示のみに使われ、プロジェクトデータとして取り込む機能が提供されていません。

shapefile データをプロジェクトに取り込むための方法はありますでしょうか。なお、対象データは、.dbf, .prj, .shp, .shx の4つのファイルから構成されています。

回答¶

河心線・河岸線の要件

DioVISTAでは、左右岸線および河心線が次のようになっている必要があります。

下流から上流に向かっていること
1個のポリラインになっていること

以下に示す例では、shapefile 内で対象の河川名のポリラインが出現する順につないで、1個のポリラインにします。元データの shapefile で上記の要件が満たされていない場合、意図したとおりに河川を作ることができないこともあります。

次のステップで取り込んでください。

shapefile を geojson に変換します。QGIS を例に説明します。

対象の shapefile をQGISにドラッグ&ドロップします。
対象の shapefile を右クリックし、[エクスポート] > [地物の保存] を選択します。
[ベクタレイヤを名前を付けて保存] ダイアログで、geojson 形式でエクスポートします。[CRS] は、[EPSG: 4326 - WGS 84] にしてください。

河心線の geojson の例です。

{ "type": "Feature", "properties": { "W05_004": "OO川", }, "geometry": { "type": "MultiLineString", "coordinates": [ [ [ 133.000, 34.000 ], [ 133.001, 34.001 ], [ 133.002, 34.002 ] ] ] } }
...

左右岸線の geojson の例です。

{ "type": "Feature", "properties": { "河川名": "OO川", "左右岸": "左岸" }, "geometry": { "type": "MultiLineString", "coordinates": [ [ [ 133.000, 34.000 ], [ 133.001, 34.001 ], [ 133.002, 34.002  ] ] } }
...

PowerShell スクリプトを作ります。ImportFrom-GeojsonToFsxprojDRiver.ps1 として保存します。

スクリプトファイルのエンコーディング

PowerShell 5 をお使いの場合は UTF16 エンコーディングで保存してください。
PowerShell 7 をお使いの場合は UTF8 エンコーディングで保存してください。

# 設定 ここから
# 河川データを取り込む先のプロジェクトファイルを指定します。
$projPath = ".\sample.fsxproj"
# プロジェクトファイルの保存先を指定します。（フルパスで指定します）。
$projOutPath = "C:\Temp\sample1.fsxproj"
# 左右岸線の geojson を指定します。
$lines_lr = ".\lines_lr.json"
# 河心線の geojson を指定します。
$lines_center = ".\lines_center.json"
# 設定 ここまで

[xml]$xml = cat -Encoding utf8 $projPath
$lrLines = cat -Encoding utf8 $lines_lr | ConvertFrom-Json
$centerLines = cat -Encoding utf8 $lines_center | ConvertFrom-Json
$rivers = $lrLines.properties."河川名" | select -Unique
$riverId = $xml.floodSim.conditions.diffusionalRivers.diffusionalRiver.Length
foreach ($river in $rivers) {
    $xDiffusionalRiver = $xml.CreateElement("diffusionalRiver")
    [void]$xml.floodSim.conditions.diffusionalRivers.AppendChild($xDiffusionalRiver)
    # 河川の属性
    $attrs = @{
        "name"                          = $river
        "riverId"                       = $riverId
        "viewItems"                     = "8191"
        "inflow"                        = "1 0 0"
        "downstreamToeWaterLevel"       = "0"
        "closeDownstreamToe"            = "False"
        "enableConnectionToDrain"       = "True"
        "maxTransverseInterval"         = "100"
        "depthValidMin"                 = "0.01"
        "areaValidMin"                  = "0.003"
        "onewayRtoP"                    = "False"
        "prConnExpressionOfRiverLine"   = "0"
        "riverWallAllowReverseOverFlow" = "False"
        "valid"                         = "True"
        "centerLineStyle"               = "0 2 0 0 255 255"
        "leftLineStyle"                 = "0 2 51 102 255 255"
        "rightLineStyle"                = "0 2 255 153 204 255"
        "transverseLineStyle"           = "0 1 255 153 0 255"
        "riverWallLineStyle"            = "0 4 204 255 204 192"
        "flowGraphMin"                  = ""
        "flowGraphMax"                  = ""
        "upperRiverId"                  = "-1"
        "upperRiverDistance"            = "0"
        "confluenceIsLeft"              = "False"
        "upperRiverIsLeft"              = "False"
        "mainstreamRiverId"             = "-1"
        "mainstreamRiverDistance"       = "0"
        "distributaryIsLeft"            = "False"
        "mainstreamRiverIsLeft"         = "False"
        "distributaryRatio"             = "0"
        "distributaryLimit"             = ""
        "steadyFlowAlpha"               = "1"
        "steadyFlowBeta"                = "1"
        "initialPrecip"                 = ""
        "precipCorrectionFactor"        = ""
        "watershedDelineationOnly"      = "False"
        "conditionalInflow"             = "False"
    }
    $attrs.GetEnumerator() | % {
        $xAttr = $xml.CreateAttribute($_.Name)
        $xAttr.Value = $_.Value
        [void]$xDiffusionalRiver.Attributes.Append($xAttr)
    }
    # 左右岸線
    $sides = @{
        "elementName"       = "leftLine"
        "shapefileProperty" = "左岸"
    }, 
    @{
        "elementName"       = "rightLine"
        "shapefileProperty" = "右岸"
    }
    foreach ($side in $sides) {
        $xLine = $xml.CreateElement($side.elementName)
        [void]$xDiffusionalRiver.AppendChild($xLine)
        $lrLine = $lrLines | ? { $_.properties."河川名" -eq $river -and $_.properties."左右岸" -eq $side.shapefileProperty } 
        foreach ($coords  in $lrLine.geometry.coordinates) {
            foreach ($coord  in $coords) {
                $lon = ($coord[0] * [Math]::PI / 180.0)
                $lat = ($coord[1] * [Math]::PI / 180.0)
                $xCoordNode = $xml.CreateTextNode(("{0:0.000000000000000} {1:0.000000000000000}" -f $lat, $lon))
                $xCoord = $xml.CreateElement("coord")
                [void]$xCoord.AppendChild($xCoordNode)
                [void]$xLine.AppendChild($xCoord)
            }
        }
    }
    # 河心線
    $xCenterLine = $xml.CreateElement("centerLine")
    [void]$xDiffusionalRiver.AppendChild($xCenterLine)
    foreach ($centerLine in ($centerLines | ? { $_.properties."W05_004" -eq $river } )) {
        foreach ($coords  in $centerLine.geometry.coordinates) {
            foreach ($coord  in $coords) {
                $lon = ($coord[0] * [Math]::PI / 180.0)
                $lat = ($coord[1] * [Math]::PI / 180.0)
                $xCoordNode = $xml.CreateTextNode(("{0:0.000000000000000} {1:0.000000000000000}" -f $lat, $lon))
                $xCoord = $xml.CreateElement("coord")
                [void]$xCoord.AppendChild($xCoordNode)
                [void]$xCenterLine.AppendChild($xCoord)
            }
        }
    }
    $riverId++
}
$xml.Save($projOutPath)

PowerShell スクリプトを実行します。作成されたプロジェクトが "C:\Temp\sample1.fsxproj" として保存されました。

FAQ/ 左岸線・右岸線のshapefileを取り込みたい¶

回答¶

関連項目¶