法務省登記所備付地図データから古道・水路を抽出

こんにちは。
「法務省登記所備付地図データのG空間情報センター配信；地図データ変換済」というデータを見つけました。

それから、（古）道（濃黄色）および水路（水色）を抽出し、地図上に拡大表示した例です：

• ただし今回の方法ではかえって最近の新しい道が抽出できない場合があるようです。

データ

上記例は、香川県小豆郡小豆島町・土庄町 2024年のGeojsonデータを用いたものです。

• https://www.geospatial.jp/ckan/dataset/aigid-moj-37324/resource/37191a00-2a5c-47d2-b0ff-38f7ef5ddbcb
• https://www.geospatial.jp/ckan/dataset/aigid-moj-37322/resource/81f191f8-09ae-4774-9e7d-0072b637c680

シェルスクリプト

抽出には下記を用いました。

$ ./extract_path.sh *.geojson > path_extracted.jsonl

extract_path.sh

#!/bin/sh

# function
jq_f() {
  case "$2" in
    道) color="#909000";;
    水) color="#00ffff";;
     *) return;;
 esac
 jq -c '.features[]' "$1" | grep -e "$2"- | sed -r 's/(},\"geometry\":)/,\ \"color\":\ \"'"$color"'\"\1/'
}

# main
for file_geojson in "$@"; do
  jq_f "$file_geojson" 道
  jq_f "$file_geojson" 水
done