日本列島 自然環境・温泉GIS
地形・地質・火山・温泉を一枚の地図に重ね、公開データの範囲でその関係を実測する
教育・研究・技術実験用の Web GIS。
本番: https://japan-geo-onsen.vercel.app
- 地図: MapLibre GL JS
- 外側: Next.js（静的書き出し）— サーバ側の処理はひとつも無い
- データ加工: Python（ETL は手動実行）
できること
| 画面 | 中身 |
|---|
/ | 地形・陰影・色別標高・地質・等高線・火山・温泉のレイヤー地図。クリックで属性 |
/stats/ | 温泉と対照群を、地質・標高・傾斜・火山距離で比べた集計 |
/onsen-stats/ | 環境省の温泉利用状況（源泉数・温度別源泉数・湧出量）を都道府県別に |
/about/ | 出典・利用条件・加工方法と、探して見つからなかったデータの記録 |
等高線は事前に作った線ではなく、表示範囲の標高タイルをその場で読んで
marching squares で引いている。
この地図が出している数
| 対象 | 件数 | 出典 |
|---|
| 活火山 | 111 | 気象庁 |
| 温泉（観光資源） | 1,320 | 国土交通省 国土数値情報 P12(2014) |
| 対照点 | 1,320 | 同上（温泉以外の観光地点を県ごとに同数抽出） |
| 地質凡例 | 2,416 | 産総研 シームレス地質図V2 |
分かったこと（実測 2026-09-07）
- 気象庁の火山一覧 JSON は 120 件あるが、活火山は 111 である。
120 = 111 + 3（座標を持たない一覧表示用の項目）+ 6（親火山の火口細分）。
この分解が崩れたらデータ生成は止まる。
- 温泉の泉質・泉温・湧出量を地点ごとに載せた全国の公開データは存在しない。
国土数値情報の全 157 データセットに温泉は無く、環境省は都道府県別の集計 PDF のみ。
よってこれらの欄は null のまま出荷し、画面には「公開データに無い」と書く。
- 環境省が数える温泉地 2,839 の一覧は公開されていない。 集計 PDF に載るのは数だけで、
名前も座標も無い。だから 2,839 そのものを点にすることはできない。
被覆を広げるため Wikidata(CC0)の温泉 1,477 件を第二の層として足した
（重複を除いた地点 2,455＝温泉地の 86%）。ただし悉皆調査ではないので統計には使わない。
- 重なった点は一覧から選ばせる。 画素 1 点の問い合わせで先頭だけを開くと、上に描かれた層の点しか選べない。
層を優先せず、カーソルの下に重なる点だけを一覧にする(近いだけの点は入れない・押し損ねは最寄りの 1 点)。
判定の半径は描画と同じ段から作る(
lib/pick.ts)。
- AI の説明(§35)は、その点を学習していないモデルで出す。 学習に使ったモデルだと RF は点を覚えて 1 に近い値を出す。
出力は確率と呼ばない。SHAP の加法性と、分割外 ROC-AUC が合否のモデルと一致することをテストで押さえる。
- AI(Phase 4〜5)は合否を学習の前に書いた。 当てる問いは「P12 に登録された地点が温泉の分類か」(温泉の存在ではない)。
都道府県分割で検証(無作為分割は空間の漏れで高く出る)。全国ポテンシャル地図は §36 の読まれ方を招くので塗らない。
- 地域比較(§58)は境界を描かない。 基準点(Wikidata の項目・気象庁の活火山)+ 半径 10/25 km の円で定める。
同名の別物は分類で落とし、1 つに決まらなければ止める。P12 の点が 5 つ未満の地域(別府 3 点)は要約に注記。
- 温泉比較(§57)は数を作らない。 設計書の表には泉温・湧出量・pH が並ぶが地点別の公開データは無いので、
両列とも「公開データに無い」と出す。所在都道府県の集計は表の外に「この温泉の値ではない参考」として置く。
- 周辺環境分析(5/10/25/50km)は数え切れるものだけ数える。 火山・温泉・湖沼(縁まで)・植生自然度。
河川と標高の起伏は計算量の理由で数えない。候補の絞り込みは総当たりと完全一致を確かめてある。
数は点ファイルに入れず、点を開いたときに読む
surroundings/<層>.json に分けた(起動時の取得を増やさない)。
- 辿れない理由は三つある。 ①分類が違う(ほったらかし温泉: 日帰り入浴施設)、
②座標が無い(大田区の黒湯温泉: 分類は温泉だが P625 が無い)、
③記事と項目が紐付いていない(山中湖温泉: 項目は在るがサイトリンクが無い)。
②③は地図に出せないので、名前の一覧として
/missing/ に出す(518+432 件)。
- 0 件は「無い」の証拠にならない。
Category:日本の温泉 は存在せず(実体は都道府県別)、
存在しない入口を叩いて 0 件を得ていた。入口が在ることを確かめてから数える。
- 分類の食い違いで、どちらの出典にも入らない温泉がある。 ほったらかし温泉は
観光資源データに 1 件も無く、Wikidata では「日帰り入浴施設」に分類されていて
温泉クラスに当たらない。同じ理由で外れる 108 件を第三の層にした
（ただし銭湯・公衆浴場も入るので温泉とは限らない。統計には使わない）。
- 地形断面は縦を引き伸ばして描くので、倍率を必ず出す。 上り・下りの合計は刻みの
細かさで変わる量なので、そのことも書く（起伏は刻みが 1 画素に近ければ落ち着く）。
- 植生自然度で差が出るのは中ほどの段階だけ。 温泉は市街地・農耕地に少なく（−14.2 pt）、
植林地・二次林に多い（+7.1 / +8.5 pt）一方、自然林・自然草原では +0.09 pt でほぼ差が無い。
- 温泉点の被覆には穴がある。 10 都府県に 1 点も無い。「点が無い＝温泉が無い」ではない。
環境省の公表値と並べると、地図の点 1,320 は温泉地 2,839 の 46%、源泉 27,899 の 5% にすぎない。
- 泉温は都道府県別になら実データがある。 環境省の集計で、温度の記録がある源泉 24,668 のうち
51.5% が 42 度以上。ただし源泉全体の 11.6% は温度が測られていない。
- 環境省 PDF は抽出テキストでは数字が割れており、書式では列を決められない割れ方がある
（「4」+「7」= 47）。文字の座標で切ると全 47 行が例外なく 20 列になった。
- OpenStreetMap の
amenity=public_bath 等は全国 5,803 件で被覆は広いが、
源泉ではなく入浴施設であり、ODbL の継承条件が政府標準利用規約のデータと混ざるため採らなかった。
作り直す
python -m venv .venv
./.venv/Scripts/python.exe -m pip install pyshp Pillow pytest
npm install
# 生データの取得（raw/ にキャッシュされ、二度目からは外部を叩かない）
./.venv/Scripts/python.exe etl/build_volcanoes.py
./.venv/Scripts/python.exe etl/build_onsen_base.py
./.venv/Scripts/python.exe etl/enrich_onsen.py
./.venv/Scripts/python.exe etl/build_onsen_wikidata.py # 第二の点レイヤー(CC0)
./.venv/Scripts/python.exe etl/enrich_onsen.py onsen_wikidata.geojson
./.venv/Scripts/python.exe etl/build_control.py
./.venv/Scripts/python.exe etl/enrich_onsen.py control.geojson
./.venv/Scripts/python.exe etl/build_water.py # 河川・湖沼(下記の zip が要る)
./.venv/Scripts/python.exe etl/parse_onsen_stats.py # 環境省 PDF(都道府県別)
./.venv/Scripts/python.exe etl/build_onsen_facility.py # 入浴施設(第三の層・Wikidata)
./.venv/Scripts/python.exe etl/build_onsen_wikipedia.py # 温泉記事(第四の層)と未掲載一覧
./.venv/Scripts/python.exe etl/build_surroundings.py # まわり 5/10/25/50km の数(全 4 層)
./.venv/Scripts/python.exe etl/build_regions.py # 地域比較(5 地域 + 活火山 111 を基準点に)
./.venv/Scripts/python.exe etl/build_ai_dataset.py # AI のデータセット(P12 温泉 × 対照)
./.venv/Scripts/python.exe etl/train_baselines.py # LR・RF・XGBoost の比較(ローカルで学習・合否は docs/ai_preregistration.md)
./.venv/Scripts/python.exe etl/explain_ai.py # 温泉ごとの見立て(分割外の出力 + SHAP の寄与)
./.venv/Scripts/python.exe etl/build_cnn_patches.py # CNN に見せる地形の画像(z14 64×64・raw/ai/ に置き積まない)
./.venv/Scripts/python.exe etl/train_cnn.py # CNN の比較(約 4 時間・途中から再開可・合否は docs/ai_preregistration_cnn.md)
./.venv/Scripts/python.exe etl/build_vegetation.py # 植生(下記の lzh が要る)
./.venv/Scripts/python.exe etl/make_profile_fixture.py # 断面の二実装照合フィクスチャ
./.venv/Scripts/python.exe etl/check_elevation.py # 標高を地理院 標高 API と突き合わせる
./.venv/Scripts/python.exe etl/build_stats.py
./.venv/Scripts/python.exe etl/build_manifest.py
npm run verify # 型検査 → vitest → next build
node harness/smoke.mjs --shot # 実ブラウザ検品(手元)
node harness/smoke.mjs --url https://japan-geo-onsen.vercel.app # 本番に対する検品
植生の生データ
curl -L -o raw/biodic/veg_c02.lzh https://www.biodic.go.jp/dload/veg_c02.lzh
curl -L -o raw/biodic/veg05m01.lzh https://www.biodic.go.jp/dload/veg05m01.lzh
"/c/Program Files/7-Zip/7z.exe" x -y -oraw/biodic raw/biodic/veg_c02.lzh
"/c/Program Files/7-Zip/7z.exe" x -y -oraw/biodic raw/biodic/veg05m01.lzh
LZH 形式なので 7-Zip が要る（Python の標準ライブラリでは開けない）。
展開物はリポジトリに入れていない（配布形の lzh だけを残す）。
河川・湖沼の生データ
リポジトリには入れていない（河川は 47 県で 296 MB）。次で raw/ に置く。
# 河川 W05（都道府県別・年版が県ごとに違うので一覧ページから拾う）
# https://nlftp.mlit.go.jp/ksj/gml/datalist/KsjTmplt-W05.html
# 湖沼 W09（全国 1 ファイル・7.7 MB）
curl -L -o raw/ksj/W09/W09-05_GML.zip https://nlftp.mlit.go.jp/ksj/gml/data/W09/W09-05/W09-05_GML.zip
河川データ(W05)の利用条件は非商用限定で、国土数値情報の他のデータとは異なる。
配る
raw/ は 570 MB あり、そのまま送るとアップロードが中断したうえ
Vercel free のファイル枠(5,000)を焼く。.vercelignore で除外する。
vercel deploy --prod --yes
--archive=tgz を使ってはならない。 この経路は .vercelignore を
ローカルで適用せず、raw/ を丸ごと詰めて 530 MB を送ろうとする
(実測 2026-09-08)。archive はファイル数の上限を避ける道具であって、
容量を減らす道具ではない。
デプロイ中に npm run build を走らせない。 送信元の out/ が
書き換わると ENOENT で落ちる。
配ったあと、出た版が意図した版かを実測する。 デプロイの終了コードを読み、
node harness/smoke.mjs --url <本番> を通し、さらに
配られている JS が意図した実装かを確かめること —— デプロイが失敗していても
古い版が動いているので、画面を見ただけでは気づけない。
生成はビルドより前に置くこと。 成果物は public/data/ に直接書いており、
data/ から写す段は作っていない（写す段があると、写した後に元を変えたときに
成果物だけが古く残る — HC-199）。
データ出典
| 提供元 | 使っているもの | 利用条件 |
|---|
| 国土地理院 | 地図・陰影起伏・色別標高・標高タイル | 国土地理院コンテンツ利用規約 |
| 産業技術総合研究所 地質調査総合センター | 20万分の1日本シームレス地質図V2 | 政府標準利用規約(第2.0版) |
| 気象庁 | 活火山の一覧と位置 | 気象庁ホームページの利用について |
| 国土交通省 | 国土数値情報 観光資源データ(P12) | 国土数値情報 利用約款 |
免責
この地図と集計は公開データをそのまま重ねたものです。
温泉の存在・泉質・湧出量を保証するものではありません。
「ここを掘れば温泉が出る」という主張はしません。
ライセンス
MIT License © 2026 坂田哲朗（コードについて。データは各提供元の条件に従う）
