Mosaic to New Rasterの最小値と最大値を理解していますか?


14

このソースから複数のDEMラスターを1つに結合する必要があります。 http //srtm.csi.cgiar.org/SELECTION/inputCoord.aspなので、Mosaic to New Rasterツールを使用します。

ツールをセットアップしました:

  1. ラスターを入力します(すべて同じソース、同じサイズ、投影なし)。
  2. ピクセルタイプを設定します(元のラスターのように32ビット浮動小数点)。
  3. セルサイズを元のラスターとして設定します。
  4. 元のラスターとしてのバンド数= 1。モザイクオペレーター、私はBLENDとMEANでそれを行いました(同じ結果が得られます)。

私が抱えている問題は、結果のモザイクが、個々のラスターの最大値と最小値とは異なる最大値と最小値の範囲を示すことです。たとえば、ラスター1(-5123.8、23.25)、ラスター2(-5974.6、40.09)、ラスター3 (-57770.2、38)、ラスター4(-2534.3、23.55)、および最終モザイクラスター(-5975.8、81.1)。

この解決策は正しくないと思います。少なくともそれを期待していませんでした。これが大丈夫かどうか、そしてそうでない場合、誰でもそれを解決し、正しい最大値と最小値を持つ適切なモザイクラスターを取得する方法を知っていますか?

ArcGIS 10.2.2 for Desktopを使用しています。


入力ラスターと出力ラスターの統計をどのように見つけたか教えてください。効率を高めるために、ArcGISはデータのサブセットのみをごまかして使用し、統計を推定する場合があります。特に、これは、最小値と最大値のレポートでかなり外れている可能性があることを意味します。
whuber

ご関心をお寄せいただきありがとうございます。レイヤー(初期ラスターとモザ​​イクの両方)を右クリックして、[プロパティ]を確認します。Symbologyの各ラスターの値をStrechedからUnique値に切り替えることを確認し、すべてのレイヤーでStreched値に切り替えると、1つのレイヤーの最小Streched値と一致することがわかりました。別のモザイクから作成されたときに表示される最小値と最大値への最大ストレッチ値。手順がまだ正しい場合、StretchedとUniquesの値のこれらの違いは誰でも説明できますか?ありがとう!
トム

2
「プロパティ」タブは不正行為であると思われます。詳細については、gis.stackexchange.com / questions / 8780を参照してください。
whuber

1
ありがとう。「プロパティ」タブが不正行為をしている可能性について考えたことはありません。ArcGISの背後にある理由が、1つのラスターレイヤーに対して、SymbologyのStreched値とUnique値に異なる値を提供していることがわかっている場合、これは明らかになる可能性があります。誰もがその理由を知っていますか?
トム

回答:


10

whuberが述べたように、多くの場合、ラスタープロパティで見つかった統計は近似値であるか、古くなっています。これらは、実際のラスター値を誤解させる可能性のある事前定義されたプロパティです。

NumPy配列を使用して、実際のデータの100%から独自の最小/最大値を計算しました。ArcGISでのNumPyの操作、およびRasterToNumPyArray(arcpy)を参照してください。例えば:

import arcpy

inrast = r'C:\data\inRaster.tif'
my_array = arcpy.RasterToNumPyArray(inrast)
print((my_array.min(), my_array.max()))

欠損値(NODATA)がある場合、正しい統計を取得するにはマスクされた配列が必要です。

import numpy as np
my_array = arcpy.RasterToNumPyArray(inrast)
my_masked_array = np.ma.masked_equal(my_array, arcpy.Raster(inrast).noDataValue)
print((my_masked_array.min(), my_masked_array.max()))

また、ラスターをNumPy配列として読み込むためにArcGISは必要ありません。例えば、GDALまたはrasterioは同様のことができます


6

他の回答では、統計は古くなっている可能性が高いと述べています。ArcGISを使用する場合は、[データ管理]ツールボックスの[統計の計算]ツールを試してください。これにより、統計が更新されます。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.