Idrisi

Idrisi32 Release Two. © 1987 – 2002 J. Ronald Eastman & Clark University.

Idrisi je akademický, modulární a silně rastrový projekt Clarkovy univerzity. K výpočtům všech níže uvedených charakteristik či operací používá modifikovaný přístup Jensona a Dominguea (1988) – v čem však tkví ona modifikace prameny Idrisi neříkají.

Příkaz Flow, který je určen k výpočtu flow direction, najdeme v sekci GIS Analysis – Surface Analysis – Feature Extraction. Do procedury vstupuje výškový grid. Dále lze zvolit možnost odstranění sinků – tato operace je modifikací přístupu Jensona a Dominguea (1988) a během ní jsou vytvořeny tři dočasné soubory, které jsou po ukončení procesu smazány – dva jsou typu integer a třetí může být real, ale o jejich funkci více help Idrisi neříká. Také lze definovat pozadí (místa s Null hodnotou), které bude při výpočtu pominuto. Výstupem příkazu Flow je grid celých čísel – kódování je zřejmé z obrázku. Jedná se tedy o algoritmus D8. Pixely, kterým nelze přiřadit flow direction, nesou hodnotu -1. Hodnota 0 je pak spojena s pixely na okrajích a rovněž s pixely tvořící pozadí.

Příkaz Runoff zde provádí výpočet akumulace srážkových jednotek na pixel. Vstupní vrstvou je zde opět výškový grid a lze uplatnit odstranění sinků. Z toho je zřejmé, že opět bude vypočítán grid flow direction, který je nutnou podmínkou pro výpočet FA, a následně bude jako meziprodukt smazán. Ve výstupním gridu nemají pixely na okrajích gridu platné hodnoty FA, ale pouze hodnotu 1, která je počáteční hodnotou pro všechny pixely. Je to způsobeno algoritmem výpočtu flow direction, jelikož tyto okrajové pixely nemají všech osm sousedů.

Stejně jako v ArcView GIS i zde lze uplatnit alternativní výpočet FD pomocí sklonu a orientace pixelu. Příkaz Surface, resp. Aspect má na vstupu výškový grid a na výstupu grid reálných čísel se spojitým rozmezím hodnot od 0º do 360º. Aby mohly být řádně ohodnoceny i okrajové pixely, příkaz Surface dočasně rozšíří vstupní data o jednu řadu pixelů na všech stranách a tyto pixely pak mají hodnotu stejnou jako jejich nejbližší soused.

U příkazu Watershed, který vykresluje povodí, je opět vstupní vrstvou výškový grid a při průběhu výpočtu je vytvořeno až sedm nových dočasných gridů – je jasné, že se jedná mimo jiné o grid FD, grid FA a pomocné gridy při odstranění sinků. Takový postup je nanejvýš zbytečný zvláště pak, když při hydrologickém modelování je často zapotřebí všech zmíněných vrstev. Řešením by bylo buď změnit požadavky na vstupní vrstvu, nebo pak vytvořit jediný příkaz (Watershed) s možností uchování dílčích meziproduktů.