Оценка показателей значений в ячейках регулярной сети

[Ремесленник]

Приветствую
Вас
Господа!
Вопрос в следующем. В одной из статей по оценке состояния лесов нашел способ проведения оценки с использованием регулярной сети с ячейками произвольного размера (в частности в статье использовалась ячейка 10 на 10 км). Подскажите, как это делается? Подробно насколько это возможно. Как такое можно сделать в ERDAS?
к примеру, разбить территорию на ячейки 5 на 5 км и посчитать сколько пикселов интересующего класса попадает в каждую такую ячейку?

Заранее благодарен!

ссылка на статью :http://zhurnal.ape.relarn.ru/articles/2006/100.pdf

[Saimon]

Добрый день,

Это можно сделать в ERDAS с помощью зональной статистики
Какая версия ERDAS используется?

С уважением

[Ремесленник]

Используется ERDAS 9.2 30-дневный.

[Ремесленник]

и вообще как это делается? Может есть какой - либо литературный источник, ссылка, или статья? А если делать это не в ERDAS?

[Ремесленник]

как затем сделать из этого карту подобную той что на рисунке в первом моем посте?

[SS_Rebelious]

В зависимости от того, какае у вас исходные данные, это можно сделать и в QGIS, к примеру, и в R.

Сам алгоритм прост - территория разбивается на полигоны(в данном случае квадраты), затем считается, сколько точек с определённым значением попадает в каждый полигон, затем полигоны раскрашиваются в различные цвета в зависимости от количества попавших в них точек.

[Saimon]

Как сказали уже выше это можно сделать. ну практически, в любой ГИС.
В свое время когда мне надо было посчитать плотность километража сейсмопрофилей на единицу площади я использовал ArcGIS.

[Ремесленник]

В зависимости от того, какае у вас исходные данные, это можно сделать и в QGIS, к примеру, и в R.

Сам алгоритм прост - территория разбивается на полигоны(в данном случае квадраты), затем считается, сколько точек с определённым значением попадает в каждый полигон, затем полигоны раскрашиваются в различные цвета в зависимости от количества попавших в них точек.

Это то что мне нужно!!!

Что касается данных то это - спутниковый снимок прошедший неконтролируемую классификацию. Сохранены только 3 класса. Обнажения горных пород (доломитов), лес и все остальное. Площадь обнажений не велика в пределах 10-15 пиксел на снимке landsat 5 TM.

Добрый день,

Это можно сделать в ERDAS с помощью зональной статистики
Какая версия ERDAS используется?

С уважением

Господа! Я искренне прошу у Вас прощенья, но я совершенно не знаю что мне делать. Скачав Qgis и R я понимаю что не успею их освоить. Я заблуждался по поводу другой гис.


И я вынужден спросить:
А как использовать зональную статистику в ERDAS (грубо говоря, где это и как это? как применить эту процедуру к моим данным?) и как получить тот результат который описал "SS_Rebelious"

... считается, сколько точек с определённым значением попадает в каждый полигон, затем полигоны раскрашиваются в различные цвета в зависимости от количества попавших в них точек.

Это то самое , что я ищу!!!

[Григорий]

Про зональную статистику смотрите здесь:
viewtopic.php?t=1141

[SS_Rebelious]

как получить тот результат который описал "SS_Rebelious"

SS_Rebelious писал(а):... считается, сколько точек с определённым значением попадает в каждый полигон, затем полигоны раскрашиваются в различные цвета в зависимости от количества попавших в них точек.

Всё довольно просто.
1. Вычисляете необходимые индексы с помощью калькулятора растров (Raster->RasterCalculator) - получаете результирующий растр.
2. Устанавливаете модуль Point sampling tool (мануал QGIS вам в помощь)
3. С помощью Vector->ResearchTools->RegularPoints создаёте сетку в которой каждая точка соответствует отдельному пикселю растра и при помощи модуля Point sampling tool задаёте значения пикселей значениям точек. Подробнее здесь
4. Разделяете точечный слой на несколько точечных слоёв в соответствии с требуемыми значениями (выделяете точки с нужными значениями расширенным поиском в таблице атрибутов, затем правой кнопкой мыши кликаете по слою и выбираете "сохранить выделенное как").
5. Теперь идёте Vector->ResearchTools->VectorGrid и создаёте полигональную сетку нужного размера.
6. Несколько раз выполняете Vector->AnalysisTools->Points in polygon - с использованием полигональной сетки и слоёв, полученных в пункте 4.
7. Открываете получившиеся слои полигональных сеток (с записанным в них количеством попавших точек) заходите в свойства слоя и настраиваете легенду (дважды щлкаете по слою - выбираете вкладку Stile, нажимаете new symbology, выбираете Graduated (градуированный знак) задаёте количество классов, классифицируете, редактируете классификацию если надо). Опять же здесь поможет мануал.
8. В процессе работы кратко записываете каждый шаг, делаете скриншоты, верстаете необльшую статейку-руководство и публикуете здесь, на сайте ГИС-Лаб, благо рыба перед вами

P.S. про статью я серьёзно. сделайте хотя бы скриншоты.

[Ремесленник]

Уважаемый Господин SS_Rebelious! Благодарю Вас за оказанное внимание и Ваш дельный совет. Инструкция которую Вы подготовили для меня разрешила все мои трудности.
С Уважением , Ремесленник.

[Ремесленник]

Также я не забываю прочих откликнувшихся. Господа, я искренне Вам благодарен.

С Уважением , Ремесленник.