最終目標は
各セルごとに気温データのあるラスタファイルから、任意の緯度経度に対応したデータを取り出したいのです。
しかし、ラスタファイルを触るのははじめてのためまずは任意のセルに埋め込まれた数値データを取得することを目的とします。
データは
world Climhttp://worldclim.org/version2
のaverage temperature (°C)(tavg 10m)を用いました。
このファイルセットには12枚のtifファイルが有り、気温帯ごとに?ファイルが分かれているようです。
これをQGISに読み込み、QGISのpythonコンソールから
https://gis-oer.github.io/gitbook/book/materials/python/07/07.html
に習って、バンドを行列データに変換してみました。
結果
array([[-3.4000000e+38, -3.4000000e+38, -3.4000000e+38, ...,
-3.4000000e+38, -3.4000000e+38, -3.4000000e+38],
[-3.4000000e+38, -3.4000000e+38, -3.4000000e+38, ...,
-3.4000000e+38, -3.4000000e+38, -3.4000000e+38],
[-3.4000000e+38, -3.4000000e+38, -3.4000000e+38, ...,
-3.4000000e+38, -3.4000000e+38, -3.4000000e+38],
...,
[-2.6385750e+01, -2.6817751e+01, -2.6819750e+01, ...,
-2.5207001e+01, -2.5195000e+01, -2.4301750e+01],
[-2.4037750e+01, -2.4944750e+01, -2.4939999e+01, ...,
-2.3297501e+01, -2.3284500e+01, -2.2035500e+01],
[-2.4511250e+01, -2.4910000e+01, -2.4910000e+01, ...,
-2.3329250e+01, -2.3312000e+01, -2.2505249e+01]], dtype=float32)
というarrayを取得できました。
しかし、埋め込まれいるのは気温データのはずなのに
-3.4000000e+38 というのは明らかに違う情報であるようです。
これではセルごとのデータを取得することはできないのでしょうか。
完全に行き詰まっており、何をすれば良いかもわかりません。助けていただけませんでしょうか。
回答1件
あなたの回答
tips
プレビュー
バッドをするには、ログインかつ
こちらの条件を満たす必要があります。