Mosaic to New Raster ile minimum ve maksimum değerleri anlama?

Birkaç DEM rasterini bu kaynaktan bir araya getirmem gerekiyor: http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp , bu yüzden Mosaic to New Raster aracını kullanıyorum.

Aracı kurdum:

Ben rasterleri girdim (hepsi aynı kaynaktan, aynı boyutta ve projeksiyon yok);
Piksel Türünü ayarlayın (orijinal rasterlerdeki gibi 32 bit kayan nokta);
Hücre boyutunu orijinal rasterler olarak ayarlayın;
Bant Sayısı = orijinal raster olarak 1; Mozaik Operatörü, BLEND ve MEAN ile yaptım (aynı sonucu alıyorum).

Sahip olduğum sorun, elde edilen mozaiğin, raster 1 (-5123.8, 23.25), raster 2 (-5974.6, 40.09), raster 3'ün maksimum ve minimum değerlerinden farklı bir maksimum ve minimum değer aralığı göstermesidir. (-57770.2, 38), raster 4 (-2534.3, 23.55) ve son mozaik raster (-5975.8, 81.1).

Sanırım bu çözüm doğru değil, en azından ben bunu beklemiyordum. Herkesin iyi olup olmadığı hakkında bir fikri vardır ve eğer değilse, nasıl çözülür ve doğru maksimum ve minimum değerlerle uygun bir mozaik raster elde edilir?

Masaüstü için ArcGIS 10.2.2 kullanıyorum.

— tom
kaynak

Girdi ve çıktı rasterlerinin istatistiklerini nasıl bulduğunuzu söyleyebilir misiniz? Bazen, verimlilik açısından ArcGIS, istatistikleri tahmin etmek için verilerin yalnızca bir alt kümesini aldatır ve kullanır. Özellikle bu, minima ve maxima raporlamasında yoldan çıkabileceği anlamına gelir.

— whuber

İlginiz için teşekkürler. Katmanı sağ tıklıyorum (hem ilk rasterler hem de mozaik) ve Özellikler'e bakıyorum. Sembolojideki her bir rasterin değerlerini Streched'den Unique değerlerine geçirmek için kontrol ettiğimi belirtmek gerekir ve tüm katmanlar için Streched değerlerine geçersem, bir katmandaki min Streched değeriyle eşleşir. max Başka bir değerden mozaik oluşturulduğunda gördüğüm min ve max değerlerine Streched değeri. Prosedürün hala doğru olup olmadığını ve Streched ile Uniques değerleri arasındaki bu farkları kimse açıklayabilir mi? Teşekkürler!

— tom

"Özellikler" sekmesinin hile yaptığından şüpheleniyorum. Bununla ilgili daha fazla bilgi için gis.stackexchange.com/questions/8780 adresine bakın .

— whuber

Teşekkürler. 'Özellikler' sekmesinin hile yapma olasılığını hiç düşünmedim. Belki de, ArcGIS'in arkasındaki nedenin, tek bir raster katmanı için Sembolojide Streched ve Unique değerleri için farklı değerler verdiğini biliyoruz. Nedenini bilen var mı?

— tom

Whuber'ın belirttiği gibi, raster özelliklerinde bulunan istatistikler bazen yaklaşıktır veya güncel değildir. Bunlar, gerçek raster değerlerine yanıltıcı olabilecek önceden belirlenmiş özelliklerdir.

NumPy dizilerini kullanarak gerçek verilerin% 100'ünden kendi min / maks değerlerinizi hesapladı. Bkz ArcGIS numpy ile Çalışma ve RasterToNumPyArray (arcpy) . Örneğin:

import arcpy

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

Eksik değerleriniz (NODATA) varsa, doğru istatistikleri almak için maskelenmiş bir dizi gerekir:

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()))

Ayrıca, sayıları NumPy dizileri olarak okumak için ArcGIS'e ihtiyacınız yoktur; örneğin GDAL veya rasterio benzer olabilir .

— Mike T
kaynak

Diğer yanıtların söylediği gibi, istatistiklerin güncelliğini yitirmesi muhtemeldir. ArcGIS kullanmayı tercih ediyorsanız, Veri Yönetimi araç kutusundaki İstatistikleri Hesapla aracını deneyin. Bu sizin için istatistikleri güncellemelidir.

— kırlangıç
kaynak