(475224.0, 4186282.0) のような座標値を持つ地理参照された画像があります。私の画像の寸法は (647, 2180) です。つまり、647 列と 2180 行があります。各ピクセルの座標を配列として取得できるように、座標値をサイズ (647, 2180) の numpy 配列に取りたいと思います。以下のようにコーディングします。
rr = rasterio.open(fname) #fname is the georefered image
col = rr.width
row = rr.height
coord = np.empty(shape=(col,row),dtype=rr.dtypes[0])
for i in range(0,col):
for j in range(0,row):
coord[i,j] = rr.transform*(i,j)
問題は、rr.transform*(i,j) が (475224.0, 4186282.0) のような値を与えることです。セルに保存する方法。上記のプログラムの場合、以下のようなエラーが発生します
トレースバック (最新の呼び出しが最後): ファイル "/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/IPython/core/interactiveshell.py"、2881 行目、run_code exec(code_obj 、self.user_global_ns、self.user_ns) ファイル ""、3 行目、coord[i,j] = rr.transform*(i,j) 内 ValueError: 配列要素をシーケンスで設定しています。