Lab上でデータを取り扱う方法を紹介してきました。Labで扱う際に役立つテクニックをまとめて紹介します。OSでのデータの閲覧方法など使い方が一部変更になっております。OSの基本操作は以下のリンクをご参照ください。Lab上でデータを取り扱う方法を紹介してきました。Labで扱う際に役立つテクニックをまとめて紹介します。画像を取得したい地点の座標を簡単に求めたり、1枚の大きな画像を取得するにはどんなコードを書いたらいいか、サンプルを動かしながら覚えていきましょう。Notebook)を使う方法は、「Tellusの開発環境でAPI引っ張ってみた」をご覧ください。OSからも取得できます。「データ選択」末尾の「グリッド」から「タイル座標」を選択すると、下記のようにタイル座標が表示されます。取得したいタイル座標をクリックしたのち、「タイル座標」右に記載されているコピーアイコンを押すことで値をコピーすることも可能です。1.地図タイルとは地図タイルとはWebメルカトル地図を256×256ピクセルの正方形に切り分けたものです。y)が決められており、世界地図(ズーム率0)を元に、ズーム率が1つ上がるごとに4分割して作られます。Webメルカトル地図は、社会や地理の教科書などでお馴染みのメルカトル図法で投影された世界地図を、東経180°から西経180°、北緯85.0511°から南緯85.0511°の範囲で切り取ったもので、ネットの地図サービスでよく使われています。y=1613)2.緯度経度からタイル座標を求めるTellusは主に地図タイル形式で衛星データを提供しています。ある地点の画像を取得したいのにその地点の緯度経度はわかるけどタイル座標がわからない、なんてことがよくあります。地図タイル確認のツールなどを使えば調べることはできますが、毎回調べるのは少々手間です。そこでPythonを使って経度緯度からタイル座標を計算するサンプルコードを紹介します。ytile)※コードの一部にopenstreetmapのサンプルをお借りしています。lat)とズーム率(zoom)を引数に与えると、その点を含んだ地図タイルの座標を返します。1617)を求めている例です。3.タイル座標から緯度経度求めるまた逆に、タイル画像を使ってデータ処理するとき画像に含まれる範囲の緯度経度が知りたくなることがあります。これもPythonを使ってタイル座標から緯度経度を計算するサンプルコードを紹介します。right_top[1])※コードの一部にopenstreetmapのサンプルをお借りしています。y)を引数に与えると、その画像の左下と右上の緯度経度を返します。バウンディングボックスがわかると、左上/右下の緯度経度も知ることができます。35.3890°)を求めている例です。4.タイル画像の各ピクセルの緯度経度を求めるタイル画像の四隅の緯度経度はget_tile_bboxメソッドを使って求めることができました。画像の各ピクセルの緯度経度は求めるにはどうすればいいでしょうか?簡易的にであれば、四隅の緯度経度の差を256で等分すれば求めることができます。lonlats※コード内でget_tile_bboxメソッドを利用しています。y)を引数に与えると、その画像の各ピクセルの左上の緯度経度を計算し256*256*2の配列として返します。y=1617)の各ピクセルの座標を求めると以下のようになります。緯度]の順で返ることに注意してください。これはあくまでも近似値です。地図タイルはメルカトル図法という投影法を使っているため、同じタイル内でも場所によって緯度の変化率がわずかに異なります(北半球なら緯度が高くなるほど1ピクセルあたりの緯度の変化が小さくなる)。そのため256で等分すると実際の座標からはずれてしまいます。変化率の差は微々たるものなので、ちょっとした計算であればこれで問題ありませんが、解析によっては厳密に座標を求めたいケースもあるでしょう。その場合は次のサンプルコードを使ってみてください。lonlats※コード内でget_tile_bboxメソッドを利用しています。y)を引数に与えると、メルカトル図法においてより高精度な緯度経度の値を256*256*2の配列として返します。5.タイルを連結して大きな画像を作る地図タイルの画像は1枚あたり256×256ピクセルとサイズが決められています。解析をするとき、画像に収まる範囲が狭くて困るケースがよくあります。ズームアウトをすれば1枚の中に広い範囲を収めることはできますが、解像度が下がってしまい見たいものがよくわからなくなっては本末転倒です。ここでは複数のタイル画像をつなぎ合わせて、解像度を落とさず広い範囲を収めた画像を作るサンプルを紹介します。params)を与えると、結合されたタイル画像を取得することができます。APIで画像を取得するメソッドは、APIごとに個別で作る必要があります。APIでASNARO-1の画像を取得するメソッドは次のように作ります。このメソッドを使って256×256ピクセルよりも大きい画像を取得してみましょう。512ピクセルの画像を取得しています。ASNARO-1以外の画像も同様に取得できます。以下は、APIでPalsar-2の画像を取得するメソッドです。このメソッドを使って今度は縦長の衛星画像を取得してみましょう。768ピクセルの画像を取得しています。Labで地図タイルを扱う際に役立つテクニックをまとめて紹介しました。今回紹介したのはライブラリへの依存を最小限にした初歩的なテクニックばかりです。タイル画像を扱っていく上でこんなことができたらいいのにといったアイデアや、こうするともっと便利だぜといった実例があれば、ぜひ宙畑編集部までお寄せください。2020/2/28【ゼロからのTellusの使い方】提供元オリジナルデータを取得する本記事では、提供元オリジナルデータ(標準処理データや高次処理プロダクトなど)とは何か、実際にデータを取得する方法をご紹介します。通説について、AMeDASのデータを用いて考察してみました。