rasterioを使用してTIFの概要から読み取る


9

rasterioを使用してtifのウィンドウを読んでいます。利用可能な概要を利用して、実際の解像度で読み取る代わりに、データのダウンサンプリングされたバージョンを取得するにはどうすればよいですか?

今私はこのようなウィンドウを読んでいます:

rmin, cmin = source.index(xmin, ymin)
rmax, cmax = source.index(xmax, ymax)
data = source.read(band, window=((rmax, rmin), (cmin, cmax)))

私がやりたいこと:

rmin, cmin = source.index(xmin, ymin, use_overview=True)
rmax, cmax = source.index(xmax, ymax, use_overview=True)
data = source.read(band, window=((rmax, rmin), (cmin, cmax)), use_overview=True)

何か案は?私の場合、概要は内部的なものですが、概要を読み取るソリューションはどれも優れています。

回答:


6

https://github.com/mapbox/rasterio/issues/710で回答。

30x30のウィンドウを3x3の配列に読み込んで、概要(可能な場合)が表示される例。

arr = np.empty(shape=(3, 3)).astype(src.profile['dtype'])
arr = src.read(1, out=arr, window=((0, 30), (0, 30)))
arr array([[9195, 9116, 9134],
           [9158, 9144, 9085],
           [9010, 8935, 9002]], dtype=uint16)

それがあなたのために働くかどうか私に知らせてください。


1
可能な場合、rasterioは常に概要を使用しますか?完全に関連する概要が利用できない場合、rasterioはリサンプリングにどのアルゴリズムを使用しますか?それは選択可能ですか?(著しく異なる出力例えば最近傍とエリア平均化弾力に注意してください。)..Perhapsは..推論rasterioのデフォルトの動作にGDALのドキュメントを参照してください
benjimin
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.