全国版：国土基本図1/2500図郭をPythonで作成する
JGD2011対応：国土基本図1/2500図郭をPythonで作成する JGD2011対応：国土基本図1/2500図郭をPythonで作成する 平面直角座標系ごとに正式な図郭コード（HD***形式）付きの 1/2500 図郭データを自動生成する手順メモ。 1. 参考にした資料 本記事の図郭コード計算ロジックは、公共測量標準図式および下記記事の内容を参考にしています。 @shiba54 (Y S) さん： 「国土基本図の図郭」を作成する 実装は上記アイデアをもとに、”全国対応版”JGD2011の各平面直角座標系に対応するよう再構成したものです。 2. 前提知識のざっくり整理 2.1 平面直角座標系（JGD2011） 日本の平面直角座標系（Gauss–Krüger 投影）は I〜XIX の 19 系。 各系には中央経線と EPSG コードが割り当てられています。 VII系：EPSG 6675 （石川・富山・岐阜・愛知 ほか） X系：EPSG 6678 （青森・秋田・岩手・山形・宮城 ほか） …といった形で全国をカバー 2.2 国土基本図の図郭体系（レベル2500） 公共測量標準図式に基づき、図郭は階層的に構成されています。 1/50,000 ：40km × 30km（2文字コード：AA〜TT） 1/5,000 ：4km × 3km（上記 50,000 を 10×10 分割 → 2桁追加） 1/2,500 ：2km × 1.5km（5,000 を 2×2 分割 → 1桁（1〜4）追加） 例： HD433 のようなコードが 1/2500 の図郭記号となります。 3. 必要な環境 Python 3.x と、次のライブラリを使用します。 pip install geopandas sha...