Null neznamená nula?

Anonymous User

<\/HEAD>

Jak sčítat rastry různých velikostí s hodnotami NoData<\/SPAN><\/EM> <\/H2>
Jedním z nejčastějších témat na fórech je, jak sečíst dva rastry<\/A> s různými prostorovými rozsahy a získat výsledek nejen tam, kde se oba rastry překrývají. Důvod, proč je výsledek pouze v těchto oblastech, souvisí s hodnotami null.

Existují proměnné prostředí, které lze upravit, aby se zajistilo zarovnání rastrů (snap raster) a zachování rozsahu během procesu (Spatial Extent). Pouhé opravení nastavení prostředí však stále nezajistí správný výstup procesu. Pokud upravíte prostorový rozsah tak, aby odpovídal "sjednocení vstupů", výstupní raster bude mít správný prostorový rozsah, ale hodnoty budou pouze tam, kde se oba rastry překrývaly. Tak co je špatně?

Očividně jsou to hodnoty pixelů ... nebo počkat, asi to není tak očividné. Vyzkoušejte tuto moudrost<\/A>, kterou poskytl uživatel fóra a Product Engineer Eric R.

Hodnota pixelu + NoData pixel = NoData pixel
Hodnota pixelu + 0 = Hodnota pixelu
Hodnota pixelu + Nic (protože vaše rozsahy jsou různé) = Nic.

Dva různé rastry použité ve výpočtu (NoData je červená).<\/P><\/DIV>

Součet rastrů<\/P><\/DIV>

Možná se ptáte, "Co se vlastně děje? Jsou hodnoty v analýze použity nebo ne?" Odpověď je, že budou použity, pokud část rovnice neobsahuje hodnotu null. Pokud je přítomna hodnota null, která způsobí, že výstup bude null, pak výstup nebude správný. 

Další otázka je, jak to opravit. Mějte na paměti, že obě řešení vyžadují, aby rastry používané jako vstupy měly shodu v prostorovém rozsahu, registraci pixelů a velikosti buněk. K tomu můžete použít nástroj Copy Raster a nastavit příslušná nastavení prostředí. Jakmile máte rastry opravené, existují dvě možnosti v závislosti na požadovaném výstupu,1.<\/STRONG>

Nastavit hodnoty NoData<\/A> na hodnotu nula, takže když sečtete rastry dohromady, výstup bude mít správný prostorový rozsah a hodnoty z obou rastrů.

Výstup po nastavení hodnoty NoData (modré pixely = 0)<\/P><\/DIV>2.<\/STRONG> Využít workflow

Con<\/A>( IsNull<\/A> ) k zahrnutí všech pixelů a jejich přičtení tam, kde je to vhodné.

Con(IsNull("RasterSMALL"),"RasterBIG", ("RasterSMALL" + "RasterBIG"))

Výstup z nástroje Con<\/P><\/DIV>

Shrnuto, existují možnosti jak pracovat s buňkami NoData. V závislosti na tom, jaký výstup chcete dosáhnout, existují různé techniky. Je také důležité pochopit vliv NoData při použití ve výpočtu rastrů. 

Existuje související příspěvek na fóru, který diskutuje použití Pythonu pro součet více rastrů dohromady:

Vypočítat součet 2334 rastrových vrstev<\/A>.Jeff S. - Geodata Raster Support Analyst<\/EM>



<\/>BODY><\/>HTML>