Pythonで衛星画像を任意の範囲・解像度に処理してGeoTiffで保存する - Qiita

はじめに 衛星画像を処理する際には、解像度の違いが問題となります。例えば異なる解像度（ピクセル数）の衛星プロダクトをNumpy arrayとして読み込んだ場合、それらは要素数異なる行列として読み込まれるため、行列についての多くの演算を行う事ができません。そのため、複数の衛星画像プロダクトを用いるのに先立ち解像度や投影を統一しておく必要があります。この課題を解決する為の処理コードを書きましたので、簡単にご紹介したいと思います。 コードができる事 ①基準となる画像に合わせてリサンプリング及びクリップ処理を行う ②任意の領域に入力画像をクリップする ③任意の解像度へリサンプリングする 前提条件 以下の状態にフォルダが構成されていることが前提となっています。例としましては、Target folderにプロダクトの名前、Sub folderにある年、もしくは月の画像を入れてある様な状態です。投影法は