JPG / PNG / BMP / GIFをGeoTIFF（TFW付き）に一括変換する方法

JPG / PNG / BMP / GIFをGeoTIFF（TFW付き）に一括変換する方法

空中写真やスクリーンショット等の JPG / PNG / BMP / GIF を、座標付きの GeoTIFF（.tif） に一括変換する手順です。
同名のワールドファイル（.jgw, .pgw, .bpw, .gfw, .wld）が隣にあれば、自動的に取り込まれます。変換時に TFW（.tfw） も同時出力します。

対象と前提

  • 対象画像: .jpg, .jpeg, .png, .bmp, .gif
  • オプションで圧縮 (COMPRESS=LZW)、タイル化 (TILED=YES)、TFW出力 (TFW=YES) を有効化
  • ワールドファイルがあれば座標を自動解釈（例: photo.jpg と photo.jgw）
  • 座標系（EPSG）はワールドファイルに含まれないため、必要に応じて後述の方法で付与

方法①：コマンド一発（単発変換）

GDAL が入っていれば、1枚だけなら以下でOK。

gdal_translate -of GTiff ^
  -co COMPRESS=LZW -co TILED=YES -co TFW=YES ^
  input.jpg output.tif
  • TFW=YES により output.tfw も生成されます。
  • 同名ワールドファイル（input.jgw など）があれば自動で反映。

方法②：フォルダ内を一括変換（Pythonスクリプト）

フォルダを選ぶだけで、JPG/PNG/BMP/GIFをまとめてGeoTIFF化します。
出力先は「同じフォルダ」または「サブフォルダ（_tif）」に変更可能です。

import os
import tkinter as tk
from tkinter import filedialog
from osgeo import gdal

# ===== 設定（必要なら変更）=====================================
VALID_EXTS = {".jpg", ".jpeg", ".png", ".bmp", ".gif"}
OUTPUT_TO_SUBFOLDER = True   # True: 入力フォルダ下に "_tif" を作ってそこへ出力 / False: 同じ場所に出力
ASSIGN_EPSG = None           # 例: 6668 など。None のままなら付与しない
CREATE_OPTIONS = [
    "COMPRESS=LZW",          # 可逆圧縮
    "TILED=YES",             # タイル化
    "TFW=YES"                # .tfw も出力
]
# ============================================================

def ensure_outdir(input_dir):
    if not OUTPUT_TO_SUBFOLDER:
        return input_dir
    outdir = os.path.join(input_dir, "_tif")
    os.makedirs(outdir, exist_ok=True)
    return outdir

def main():
    root = tk.Tk(); root.withdraw()
    folder = filedialog.askdirectory(title="画像フォルダを選択（JPG/PNG/BMP/GIF）")
    if not folder:
        print("キャンセルしました"); return

    outdir = ensure_outdir(folder)
    print(f"[INFO] 入力: {folder}")
    print(f"[INFO] 出力: {outdir}")

    converted, failed = [], []

    for name in os.listdir(folder):
        in_path = os.path.join(folder, name)
        if not os.path.isfile(in_path):
            continue
        ext = os.path.splitext(name)[1].lower()
        if ext not in VALID_EXTS:
            continue

        stem = os.path.splitext(name)[0]
        out_tif = os.path.join(outdir, stem + ".tif")

        print(f"[RUN ] {name} → {os.path.basename(out_tif)}")
        ds = gdal.Translate(
            out_tif,
            in_path,
            creationOptions=CREATE_OPTIONS,
            outputSRS=(f"EPSG:{ASSIGN_EPSG}" if ASSIGN_EPSG else None)
        )
        if ds is None:
            print(f"[ERR ] 変換失敗: {name}")
            failed.append(in_path)
        else:
            ds = None
            converted.append(out_tif)

    print("\\n===== 結果 =====")
    print(f"変換成功: {len(converted)} 件")
    for p in converted: print("  -", p)
    if failed:
        print(f"失敗: {len(failed)} 件")
        for p in failed: print("  -", p)

if __name__ == "__main__":
    main()

ポイント

  • GDALは同名ワールドファイルがあれば自動で参照します（.jgw/.pgw/.bpw/.gfw/.wld）。
  • ASSIGN_EPSG を設定すれば、変換と同時に座標系（EPSG）を付与できます。
  • COMPRESS=LZW + TILED=YES はQGISでも快適に動く定番設定です。

座標系（EPSG）の付与・確認

ワールドファイルには投影情報が無いため、必要に応じてEPSGを付与します。

# 既存の GeoTIFF に後付けする場合
gdal_edit.py -a_srs EPSG:6668 photo.tif

# 中身確認
gdalinfo photo.tif
  • 例：JGD2011 平面直角系（帯域によりEPSGコードは異なる）
  • 付与後はQGISの「レイヤプロパティ → 情報」で投影が入っているか確認

よくある質問

  • Q: .tfw が出てこない
    → 上のスクリプト・コマンドは TFW=YES で出力します。既存のGeoTIFFからTFWだけ作りたい場合は、
    gdal_translate -of GTiff -co TFW=YES input.tif output.tif として再書き出しします（中身は劣化しません）。
  • Q: 変換後に色がおかしい/サイズが大きい
    → 8bit RGBであればそのままOK。容量が大きい場合は COMPRESS=LZW の他に、16bit画像なら PREDICTOR=2 を追加すると圧縮効率が上がることがあります。
    例：-co COMPRESS=LZW -co PREDICTOR=2 -co TILED=YES
  • Q: 表示が重い
    → 変換後にピラミッド（オーバービュー）を付けると高速化します。
    gdaladdo -r average photo.tif 2 4 8 16

まとめ

  • JPG/PNG/BMP/GIFをワンクリックでGeoTIFF化（TFW同時出力）
  • 同名ワールドファイルがあれば自動で座標を取り込み
  • 必要に応じてEPSGを付与（gdal_edit.py -a_srs）
  • 圧縮・タイル化・オーバービューで実運用も軽快

コメント

コメントを投稿

このブログの人気の投稿

石川県：土砂災害（特別）警戒区域＋ＣＳ立体図

石川県の土砂災害（特別）警戒区域まとめ 2025 年 5 月 16 日 時点 石川県砂防課が公表する 土砂災害（特別）警戒区域指定状況 （PDF）を整理し、市町別に 土石流・地すべり・急傾斜地崩壊 の警戒区域（イエロー）と 特別警戒区域（レッド）を集計しました。 市町 警戒区域 (イエロー) 特別警戒区域 (レッド) 合計 土石流 地すべり 急傾斜 土石流 地すべり 急傾斜 輪島市 374 217 375 313 0 372 966 金沢市 230 82 247 170 0 235 559 能登町 240 68 218 210 0 209 526 七尾市 192 45 201 165 0 201 438 珠洲市 91 76 164 83 0 163 331 小松市 181 13 66 162 0 64 260 白山市 139 20 91 73 0 84 250 志賀町 130 15 100 106 0 99 245 穴水町 93 12 115 79 0 113 220 加賀市 78 21 101 65 0 96 200 津幡町 78 48 57 56 0 56 183 中能登町 130 4 36 99 0 35 170 能美市 64 1 41 59 0 40 106 宝達志水町 33 30 53 26 0 52 116 かほく市 12 2 30 5 0 30 44 川北町 0 0 0 0 0 0 0 野々市市 0 0 0 0 0 0 0 内灘町 0 0 5 0 0 5 5 合計 2065 654 1900 1671 0 1854 4619 ＣＳ立体図＋土砂災害警戒区域 MAP を開く 関連資料・データ ▸ 土砂災害（特別）警戒区域指定状況 PDF（R7.5.16） ▸ 石川県土砂災害（特別）警戒区域図 オープンデータ（SHP・GeoJSON 等） ▸ CS 立体図を使った 地形判読マニュアル（林野庁） ハザードマップの入手方法 お住まいの 市役所・町役場の防災担当窓口 や自...
続きを読む

日本でよく使う EPSG コード 一覧

1. 経緯度（地理）座標系 区分 名称 EPSG 主な用途・備考 世界標準（地理） WGS 84（緯度・経度） 4326 GNSS/GPS の測位やデータ交換で使用。単位は「度」。 Web地図（投影） WGS 84 / Pseudo-Mercator（Web Mercator） 3857 Google/OSM/国土地理院タイルなどの 表示 に使用。単位は「メートル」。別名：Popular Visualization / Google Mercator。 日本最新 JGD2011（緯度・経度） 6668 測地成果2011。公共測量の基準。平面直角座標系は下表参照。 過去の全国基準 JGD2000（緯度・経度） 4612 2002–2011 年に使用。JGD2011 と混在時は変換推奨。 旧日本測地系 Tokyo Datum 4301 1940年代〜2002年まで使用。WGS84とは数百mのズレあり。 2. 平面直角座標系（Gauss–Krüger 投影） 系 (ローマ数字) 系 (数字) 中央経線 JGD2011 EPSG JGD2000 EPSG 適用都道府県（概要） I 1 129°30′E 6669 2443 長崎、鹿児島西部 II 2 131°00′E 6670 2444 福岡・佐賀・熊本・大分・宮崎、鹿児島...
続きを読む

設計定数を求めるための代表N値について

代表 N 値（設計用採用値）の決め方 ― 要点まとめ 1．実務でよく尋ねられる「出典」の悩み（改訂） 発注者から「この整理手法の出典はどこか」と問われることが多い一方で、道路土工・道路橋示方書・地盤調査法など主要文献には 「平均値±1/2σで採用する」と明記した箇所はほとんど見当たりません。 しかし、 『報告書作成マニュアル［土質編］』（社団法人 全国地質調査業協会連合会） には、 ばらつきデータから設計用値を定める 三手法 （ a： 平均、 b： 最小／最大、 c： 平均と標準偏差による補正）が明記されています。 特に c) として次式の例が示され、 設計値 ＝ 平均値 − (1/2) × 標準偏差 （平均値を標準偏差で補正して設計値とする） と記載されています（第2章 2.7.3、p.74）。 したがって、本稿で用いる 「平均 − 1/2σ」 という整理は単なる経験則ではなく、 公的マニュアルに裏づけのある実務手法 です。 なお同マニュアルは、標準偏差の重み係数（1/2 など）の選び方には議論があること、すなわち k の選定はデータ量・ばらつき・重要度等に応じて調整する余地があることも併記しています。 出典： 『報告書作成マニュアル［土質編］』第2章 2.7.3「設計用土質定数の設定方法」p.74（a〜c の三手法と「平均 − (1/2) 標準偏差」の例示） 2．建築基礎構造設計規準・同解説（1974）の記述 地盤の場合、土質試験などによって得られる諸数値にかなりのばらつきがあるのがむしろ普通であり、このようにばらつきのある測定値の取り扱いについては、ばらつきの原因とその程度、基礎底面と各試料採取位置との相対的な関係など、複雑な要素が絡むため一概には言えない。しかし通常、1地点で同一土質層について得られた測定値は、その平均値を用いて差し支えないとする。そのうえで、とくにばらつきが著しい場合などは、平均値よりさらに安全側を見込んで 設計用採用値 ＝ （平均値） ± （標準偏差）／2 とする。ここで複号 ± は、算定結果が安全側となる符号をとる...
続きを読む