Contents
はじめに
久しぶりにGISに触れることにした。恐れていたのであるが，まあ，諦めた。以下，時系列で書き殴ってゆく。
1. 位置ベクトルの準備 2025年11月10日
ひょんなことから，「丹の海」にのめり込むことに。桑田郡式内社などの分布と地形の関係を知る必要性が出て来た。ここでは，式内社の地球座標系上の位置情報をGoogle Earth Proで求めて，その平面直角座標系6の座標値を求めた。その過程はいわば次のMicrosoft Excel表示のスクリーンショットに示している。
自ら作成した十数年前の授業コンテンツ座標系と変換を参照して，国土地理院のリンクが切れていないことに驚き，感激した。ここに掲載しているMicrosoft Excelファイルを使ってみたが，当時の作成意図を忘れていた。敢えて問題を残して，授業で学生が修正するようにしていたものであった。
Google Earth Proで位置情報をどのような形で表示させる必要があったのかも忘れて，経緯度情報（緯度ならdd,mm, ss.sssss）をコピペした上でこのMicrosoft Excelファイルを使おうとしたが使えない。度表現をすべきであった。緯度であれば，dd.dddddというような形である。この形はそのまま算術計算ができるので最も優れた形であった。Google Earth Proで20点についてやり直した。その入力列がA,B列である。国土地理院の座標変換サイト 平面直角座標系への換算はかつて，経度そして緯度を入力する形であったので，そのようにMicrosoft Excelで入力欄を作成したのであるが，現在，緯度そして経度の順になっている。データ入力の途中から気づいたので，姑息に触っている。
一括処理をする。この国土地理院サイトの「入力ファイルの例」の中身を次に。
平面直角座標への変換入出力ファイル形式
「一括計算」処理で使用する入力ファイルや出力ファイルの形式を説明します。
入力ファイル
一度に複数の座標値を一括して計算する「一括計算」の際、読み込むファイルの形式は以下に従って下さい。 （入力ファイル例のダウンロード→bl2xy.in）
入力ファイル名形式
・入力ファイルは、テキストファイル（形式）です。
上記の「入力ファイル例のダウンロード」からダウンロードし、テキストエディタ(Windowsの場合メモ帳など)で数値を変更の上ご利用下さい。
・入力ファイル名に使えるのは半角数字、半角小文字アルファベット、半角大文字アルファベット、半角アンダースコア、半角ドットです。日本語は使用しないで下さい。
・文字コードは、Shift JIS又はUTF-8として下さい。UTF-8の場合は、BOMなしで入力ファイルを保存して下さい。
・入力ファイルの拡張子は「.in」に限定します。 (例) filename.in
入力ファイル
・1行に、緯度、経度の順番で入力して下さい。
・緯度と経度の間の区切りは、半角空白とします。
（2つ以上の連続した半角空白も、1つの半角空白として処理されます。)
・区切り文字として、全角空白やタブやカンマは使用しないで下さい。
・緯度、経度の形式は、度分秒を区切らずに記述して下さい。
・分および秒の整数部分は2桁とします。10の位がゼロのときは空白とせず0とします。
・ヘッダーの行頭は”#”（半角）とします。
・緯度、経度の数字は半角数字を使用して下さい。[例] # 緯度(dms) 経度(dms) 370856.12340 1395201.98760
[緯度、経度の良い例 37度8分56.1234秒の場合] 370856.12340 +370856.12340
［緯度、経度の悪い例 37度8分56.1234秒の場合］ 370856 ← 秒に小数点がない 37 08 56.12340 ← 度・分・秒の区切りが半角空白となっている 37:08:56.12340 ← 度・分・秒の区切りがコロンとなっている 37-08-56.12340 ← 度・分・秒の区切りがハイフンとなっている 37.085612340 ← 度と分の間に小数点が入っている 37 856.12340 ← 分が1桁になっている ３７８５６．１２３４０ ← 数字が全角となっている
出力ファイル
一度に複数の座標値を一括して計算する「一括計算」の際、作成される出力ファイルの形式は以下の例のとおりです。
出力ファイル形式
・出力ファイル名は、標準では入力ファイル名と同様です。
・ダウンロードダイアログで、任意のファイル名を付けることができます。
・出力ファイルのデフォルトの拡張子は「.out」になります。
・計算結果が範囲外の場合は、「-9999.」が出力されます。[例] filename.out
出力ファイル（例） 数字は実際の計算結果とは異なります
# このファイルは緯度経度から平面直角座標への換算結果（一括処理）です。 # 選択された測地系は「世界測地系」、平面直角座標系の系番号は「15」です。 # 出力項目は以下の通りです。 # 緯度(dms) 経度(dms) X座標(m) Y座標(m) 真北方向角(dms) 縮尺係数 コメント #---------------------------------------------------------------------------- 370856.12340 1395201.98761 127483.442 3010.007 -113.67 0.99990011 415200.00000 1431545.00000 -9999. -9999. -9999. 1.00817437
まあ，これで全部わかるのであるが，L列，M列での，国土地理院独特の位置情報形式に，経緯度を変換するするには，concatenate関数を使うのが適当だろう。
たとえば，135度32分5.90280秒の位置情報を、国土地理院の座標系変換ツールでは，1353205.90280という形で入力することがある。分と秒の整数部は二桁でなければならず，
Microsoft Excelのconcatenate関数を利用する場合，=concatenate(G9,H9,”0″,I9)と手作業で入力する必要がある。一括変換をしたいので，自動で「分と秒の整数部の二桁」を確保したい。その式が次のものである。
=CONCATENATE(G9, TEXT(H9, “00”), TEXT(I9, “00.#####”))
図1には散布図を二つ示している。左の散布図は，Google Earth Proでのぼくの読みが問題ないかを示している。GEProの表示との間に違和感を感じなかったので，データ取得に問題はないということだ。
図1の右手の散布図は平面直角座標系への変換後のもので，これも左の散布図とよく対応できている。変換過程も問題ないということだ。
GIS世界を構築するには，一つ一つ積み上げる必要があって，この程度ではあるが，確認をしていかないと，後に痛い目に会う。
ご参考になると思うので作業したエクセルファイル を次のように，アップロードした。
2. 国土地理院から実体視用カラー空中写真をダンロード 2025年11月12日
昨日，篠町山本の桑田神社をターゲットに，周辺地形などを観察した。トロッコ列車亀岡駅のそばの無料パーキングスペースをみつけて駐車した。何とも言いにくい観光地であったが，自然環境の観察には極めて有効な場所であった。もう寒かったが保津川下りの乗船客に手を振った。ひさしぶりの体験であった。
先延ばしにしていた空中写真判読をしないとどうしようもないことを実感した。桑田神社の後背はかなりの地滑りしているように見える。錯覚か？ 図3は桑田神社境内を遠望したところである。ケヤキだったかの紅葉後背に照葉樹林の幹群が見える。図4はG.E.Proでみたところだが，2025/4/12という撮影年月日がある。今年の春である。かなりの範囲で風化した泥岩の裸岩がみえる。本殿などが立地する平坦面高度はほぼ105mである。図に見えるレールは山陰本線で，トンネルに吸い込まれている。よく見るとこの左手の保津川左岸に接したレールはかつての山陰本線でいまは，トロッコ列車が嵐山から走る。
さて，亀岡盆地の水系が収束する保津川峡谷のはじまりは図4左手に見える。いくつかの課題を解決すべく，空中写真をダウンロードした。空中写真の無料ダウンロードは，ワンストップ，という表示があれば良い。400 dpiでできるので，今回の目的には十分である。一つ一つの写真について，全画像が提供されるものにこだわってきたが，マック画面上で実体視をするには縦二分割の方が便利ではないかと考えた。空中写真測量する訳でもないし。
CKK2020 2 – C6-1〜-22，C7-1〜-22の計44枚をダウンロードした。各写真は正方形をなすが，撮影コース沿い方向の2/3長が切り取られている。つまりコース沿いの両端1/3➗2が切り捨てられているのである。この両端は歪みがあるので，切り取ってしまうという発想は大胆である。思いもよらなかった。なお，提供されている空中写真は400dpiではあるが，眠い。とはいえ，画像アプリで彩度をあげると改善される。ぼくはGIMPを使っている。色 > 彩度，で，スケールの1.000の数字を選んで，2を入力してタブキーを押すと，彩度を2倍に上げることができる。そして，ファイル > ファイル名に上書きエキスポート，すると，同じパスのファイルに上書きされる。たとえば，1.6MBのファイルは1.8MBに，1.9MBのものは2.1または2.2MBに増大する。
以上，令和7年11月16日。
次の章では，Affinityアプリをつかって，実体視についてトライしてみる。
3. Affinityで実体視