Kar sınıflamasını otomatikleştirmek için işleme alet kutusunu kullanma

Snowy Hydro Limited (SHL) is a peak electricity generator based in the Snowy Mountains of New South Wales, Australia. Despite owning a variety of generation assets across the country, Snowy Hydro Limited is best known for the custodianship of an iconic piece of Australian infrastructure, the Snowy Mountains Hydro Electric Scheme.

The Snowy Mountains Hydro Electric Scheme, is made possible by seasonal snow falls in the New South Wales alpine region. Water from the spring snow melt is used to generate electricity addressing peak demand in the national electricity market, while also diverting water to the drier inland areas of South Eastern Australia for irrigation and agriculture.

Classifying snow cover over the winter months allows SHL to monitor the life cycle of the snow covered area, archive coverage areas for historical records, and provide coverage statistics to the hydrographic services team for water yield and storage forecasts.

Harita 1.1 Karlı Su Tutma Harita

Harita 1.1 Karlı Su Tutma Harita

Arkaplan

Each year during the winter months (1st June to the 31st October) SHL obtains as many cloud free images as possible from the MODIS Terra (EOS AM) satellite. These images are then ingested into an algorithm known as the Normalised Difference Snow Index (NDSI), which is used to indicate areas where snow is present. While this algorithm is not unique to any one particular sensor, the MODIS instrument was chosen due to its adequate spatial resolution and daily coverage of the project area (Map 1.1 Snowy Water Catchment).

Hedef

Bu projenin amacı NDSI akışını model olarak QGIS işleme aracında otomatikleştirmektir. İstenen çıktı karlı su tutma bünyesinde artırılmış, kilometre kare de çevirili toplam alanın tanımlanması ile beraber karla kaplı vektör poligondur.

Yöntem

NDSI algoritması oldukça basittir(Formül 1.1 de gösterilen Normalleştirilmiş Fark Kar İndeksi(NFKİ)) ve izleyen adımlar vektör formatında kaplayan karı almayı gerektirir hatta model araç işlemeye izin verir.

Formül 1.1 - Normalleştirilmiş Fark Kar İndeksi(NFKİ)

Formül 1.1 - Normalleştirilmiş Fark Kar İndeksi(NFKİ)

This NDSI algorithm will return results between -1 and 1, with the threshold for snow typically being values greater than 0.4. An example of the NDSI output is shown below (Map 1.3 - NDSI output), alongside a natural colour image (Map 1.2 - MODIS true colour) for the 13th July 2014. As seen from these images the NDSI output is very good at discriminating snow covered areas from other land cover types and clouds. It can however misclassify portions of large water bodies as snow, caused by water having absorption characteristics in the Near Infrared part of the spectrum. This characteristic of the NDSI indicates the need to incorporate a water mask at some stage in the model.

Harita 1.2 - MODIS gerçek renk

Harita 1.2 - MODIS gerçek renk

Harita 1.3 -NFKİ çıkış

Harita 1.3 -NFKİ çıkış

Harita 1.4 - Kaydedilen "Karlı" vs "Karsız"

Harita 1.4 - Kaydedilen “Karlı” vs “Karsız”

The NDSI output is created with some very simple band map by using the GRASS module “r.mapcalculator”. Once the NDSI grid was available the NDSI threshold of 0.4 needed to be applied and this would result in a boolean grid depicting “Snow” and “Non snow” areas. To get that classification another GRASS module “r.recode” was used and produced the result shown above in Map 1.4 Recoded ‘Snow’ vs ‘Non Snow’.

The recoded grid was then converted to a vector polygon dataset and the snow covered areas were extracted (shown below, Map 1.5 - Extracted snow). The snow covered areas were then clipped to the Snowy Water Catchment (shown below, Map - 1.6 Snow inside the SWC), then the large water bodies are masked out (shown below, Map 1.7 - Final snow covered area), and finally the area statistics for the snow covered area is generated.

Harita 1.5 - Çıkarılan Kar

Harita 1.5 - Çıkarılan Kar

Harita 1.6 - SWC içindeki kar

Harita 1.6 - SWC içindeki kar

Harita 1.7 - Son kar kaplı alan

Harita 1.7 - Son kar kaplı alan

The output from the model is a vector polygon showing the snow covered area within the Snowy Water Catchment. The attributes of the polygon include the area covered by the snow and the percentage of the catchment that is covered. The model in its entirety can be seen below in “Image 1.1 - Normalised Difference Snow Index (NDSI) QGIS processing toolbox model”.

Resim 1.1 - Normalleştirilmiş Farklı Kar İndeksi(NDKİ) QGIS model araç işleme

Resim 1.1 - Normalleştirilmiş Farklı Kar İndeksi(NDKİ) QGIS model araç işleme

Netice/Sonuç

The QGIS processing toolbox is a powerful way to automate a workflow such as this snow classification. The model shown above successfully automates the snow classification process, and delivers a GIS dataset that can be used for visualisation and input into further statistical applications. It successfully removes subjective user error associated with previous methods used by SHL, while providing a repeatable and consistent classification tool. Some seasonal results can be viewed here.

Yazar

Andrew Jeffrey Cooma NSW Avusturalya temelli Snowy Hydro Limited de GIS Analistidir. Wagga Wagga daki Charles Strut Üniversitesinden 2004 yılında Uzaysal Bilgi Sisteminde uzmanlaşarak mezun olmuştur. Snowy Hydro Limited deyken Cloud Seeding deneyi gibi mevcut projelere uzaysal destek sağladı ve günden güne kurumun her bölümünde çalıştı.