緯度経度を一括で平面直角座標に変換するPythonプログラム

本ページでは、緯度・経度を平面直角座標系に一括で変換するPythonプログラムを掲載しています。以下の操作が可能です：

• 緯度経度（EPSG:4326 など）から平面直角座標系（EPSG:6669〜6686）への座標変換を行うPythonコードを生成
• 生成されたPythonコードを .py ファイルとしてダウンロード
• EPSGコードの変更に応じて、リアルタイムでコード内容が更新

使い方解説

１．変換前の座標系を選ぶ

「変換前の座標系」ドロップダウンでは、以下のような選択肢があります：

• EPSG:4326（WGS84、緯度経度）
• EPSG:3857（Webメルカトル）

通常、緯度経度が与えられている場合は「EPSG:4326」を選択します。

２．変換後の座標系を選ぶ

「変換後の座標系」では、平面直角座標系 第1〜18系 をEPSGコードで選択できます：

地域に応じた系を選ぶことで、正確な座標変換が可能になります。

３．変換コードを確認

座標系を選択すると、ページ内のコードブロックに座標コードがリアルタイムで反映されます：

from pyproj import Transformer
transformer = Transformer.from_crs("EPSG:4326", "EPSG:6677", always_xy=True)

このコードは、CSVファイルを読み込み、緯度経度から指定した座標系へ変換し、変換結果をCSVに出力するスクリプトになっています。

４．Pythonコードをダウンロード

「⬇ Pythonコードをダウンロード」ボタンをクリックすると、表示中のPythonコードが .py ファイルとして自動的にダウンロードされます。

• ファイル名：coordinate_transform.py
• 中身は、現在選択している座標系に合わせてカスタマイズされたコード
• そのままPython実行環境で動かすことができます（必要なライブラリ：pandas, pyproj）

３．読み込みファイル名の変更

ダウンロードしたPythonプログラムの読み取りファイル名（ファイルパス）を変更して実行してください。

fname = “input.csv”

読み込むCSVファイルのフォーマットはヘッダー付きの緯度,経度です。

💡 利用上の注意

Pythonスクリプトを実行するには、Python環境（3.x）と、pandas および pyproj ライブラリが必要です。

pip install pandas pyproj

• 入力ファイル（input.csv）は、少なくとも lat と lon の列を含んでいる必要があります。
• 出力ファイル（output.csv）には、元のデータに加えて X, Y 列が追加されます。