この地図が答えていること/いないこと
Landsat 8/9 の熱赤外で 2013〜2026年の夏(7〜8月)を合成し、福岡市の建物に載せたもの。使ったデータはすべて無料で、アカウント登録も課金もありません。
これは気温ではありません
測っているのは地表面の温度です。アスファルトや屋根の表面が何度になっているかで、天気予報の気温とは別物です。真夏の昼は舗装が50度を超える一方、同じ時刻の海面は29度前後になります。体感や熱中症に効くのはこちらの面の分布ですが、気温として読まないでください。
晴れた日の昼にしか測れていない
光学・熱赤外の衛星は雲を見通せないので、曇りと雨の日は欠測です。撮影は太陽同期軌道で午前10時半ごろに固定されており、夜間や夕方の暑さは一切入っていません。「夜まで熱が残る場所」を知りたい場合、この地図では答えられません。
屋根1枚より画素のほうが大きい
熱赤外センサーの分解能は約100mで、配布時に30mグリッドへ再配分されています。建物1棟の屋根は画素より小さいので、この地図の値は「その建物の屋根の温度」ではなく「その建物の周囲90mの地面がどれくらい暑いか」です。隣り合う建物が同じ値になるのはそのためで、細かい差は意味を持ちません。
検証できたことと、できないこと
型としてこの地図は「面を測る」もので、正解が存在しません。誰も福岡市全域の地表面温度を実測していないので、答え合わせができない。当てられる3つを当てた結果が下です。3件のうち1件は、検定そのものが成立しませんでした。
作っている途中で間違えた測り方
どれもエラーを出さず、それらしい数字を返していました。同じ罠は他の人も踏むはずなので残しておきます。
- 雲を外していなかった。(以下の数字は、この地図を作る前に1枚のシーンで測っていたときのもので、上に出ている建物どうしの振れ幅とは別の量です。)面の振れ幅が 34.0 degC と出て、最低が 16.6 degC になった。真夏にありえない。雲は上空にあって冷たいので、熱赤外では「とても冷たい地面」として写ります。QA_PIXEL で外したら 22.8 degC、最低は 28.9 degC で海面水温と一致した
- 絶対温度で切った。「45度以上が高温」で切ったら対象域の 53.8% が該当した。真夏の都市面はふつうに45度を超えます。閾値は絶対でなく、その画像の中での相対(分位)で置くのが正しい。この地図に「市内の順位」モードがあるのはそのためです
- 年ごとの枚数の差を無視しかけた。使えるシーン数は年によって1〜8枚とばらつきます。全部まとめて中央値を取ると枚数の多い年だけが効くので、年内の中央値を取ってから年をまたいで中央値を取る2段にしました
- 対象範囲が市域を覆っていなかった。最初の枠だと市の南端の185棟がはみ出していた。「福岡市の」と名乗るなら市域を覆っていないといけないので枠を直しました。現在も小呂島の201棟だけは40km沖合のため対象外です(欠測ではなく対象外)
データと再現
- 熱赤外: Landsat Collection 2 Level-2(USGS/Microsoft Planetary Computer 経由)。アカウント不要で署名URLが返る
- 換算: DN × 0.00341802 + 149.0 − 273.15、雲マスクは QA_PIXEL の bit1〜4
- Landsat 7 は除外。2003年以降 SLC-off で縞状の欠損があり、8/9 とは較正も違うため、混ぜると「平年」が機器の差を測ることになります
- 建物: Project PLATEAU 福岡市2024 v4 LOD1(国土交通省・CC BY 4.0)
- 背景地図: 国土地理院タイル