Pythonを使用してラスターの各ピクセルのXY座標とセル値を取得する方法は?
私は本当にPythonが初めてで、ArcGIS 10のPythonを使用して、ピクセルごとのラスターのセル値と座標(各ピクセルの中心のXY座標をマップ)を取得する簡単な方法があるかどうかを知りたいですか? これをさらに説明するには、最初のピクセルのマップX、マップY、およびセル値を取得し、それらの3つの値を3つの変数に割り当て、他のピクセルの残りに対してこの手順を繰り返す必要があります(ラスター全体をループします)。 私の質問をもっと説明する必要があると思います。問題は、最初のラスタのピクセルのXY位置を取得し、そのXY位置に対応するいくつかの他のラスタのセル値を取得する必要があることです。このプロセスは、中間ポイントシェープファイルを作成せずに、最初のラスターのすべてのピクセルをループ処理する必要があります。80億ピクセル近くのラスターを処理する必要があるため、本当に時間がかかります。また、ArcGIS 10でPythonを使用してこれを行う必要があります。 @JamesS:ご提案ありがとうございます。はい、これは1つのラスターで機能しますが、他のいくつかのラスターのセル値も収集する必要があります。問題は、最初のラスタの最初のピクセルのX座標とY座標を取得した後、最初のラスタのX、Y位置、3番目のラスタなどに対応する2番目のラスタのセル値を取得する必要があることです。だから、最初のラスターをループするとき、ピクセルのXとYの位置を取得し、その位置に対応する他のラスターのセル値を取得することは同時に行われるべきだと思いますが、わかりません。これは、最初のラスターをポイントシェープファイルに変換し、ArcGIS 10でポイント値への複数値の抽出を実行することで実行できますが、 @hmfly:ありがとう、はい、このメソッド(RastertoNumpyarray)は、配列の既知の行と列の値の座標を取得できれば機能します。 @whuber:計算を実行したくないので、テキストファイルにXY座標とセル値を書き込むだけで十分です。