Подсчет гектаров по спутниковому снимку

[Alexey13]

Всем добрый вечер!

Помогите пожалуйста, как можно подсчитать количество гектаров в растровом снимке GeoTiff?

Снимки проективного покрытия примерно по 400 га,(около 300 шт.) с 10 классами 0-10%, 10-20%, 20-30%, 30-40%, и.т.д. до 100%. В каждой карте по 4-5 классов. Необходимо узнать площадь каждого класса, можно даже в процентах от общай площади.


Из инструментов QGIS, Global Mapper, Erdas в процессе установки...

Заранее спасибо!

[novia]

В эрдасе делаете модель в model maker, состоящую из одного снимка на входе, условного оператора и одного снимка на выходе. В условном оператора применяете границы классов к пикселам снимка и записываете номер класса в новых файл.

Вот тут написано как запустить модель для 100500 файлов.

А тут писали про то, как посчитать количество пикселов каждого класса.
Соответственно "гектары" - это количество пикс. * площадь пикс. в кв. м /10000

[Сообщение с мобильного устройства]

[Alexey13]

Спасибо большое. Но способ описанный в Эрдасе слишком сложный для моего текущего уровня .

Существует ли более простой способ выделения контура по цвету? Данная операция легко проделывается в Гимпе и фотошопе, но площадь определить нельзя.

Спасибо!

[novia]

Все таки 300 снимков "не руками" обработать не совсем тривиальная задача.
Попробуйте разобраться с model maker и запуском пакетной обработки, они только с виду страшные и непонятные. Задавайте тут вопросы, у вас все получится.

[Александр Мурый]

В QGIS можно исп-ть "Инструменты GRASS": меню Растр --> Отчёты и статистика --> Статистика --> r.report
Вот ссылка на мануал. Работа с "Инструментами GRASS" описана в русскоязычном руководстве (лежит здесь на сайте).

[Максим Дубинин]

так вы пробовали вышеприведенный http://gis-lab.info/qa/im-countpixels.html ?

[Alexey13]

Второй день мучаю модуль Grass в Qgis. Выходит ошибка - неверный слой "Слой не является действительным и не может быть добавлен на карту". Ошибка выходит при импорте растра Gdal в Grass, хотя все делаю по руководству. И кстати почему то слетает кодировке в модуле Grass в путях к файлам. Возможно причина в этом?

[Alexey13]

Спасибо большое. Попробовал быстро не получилось и пока отложил, увидев возможности Grass. Но если с не получится вернусь к упомянутому вами способу.

[Александр Мурый]

Alexey13, опишите последовательность действий с GRASS через QGIS. Главное правило здесь: по возможности никаких пробелов и других "нехороших" символов в путях в файлам и в именах файлов и минимум кириллицы.
С кодировками -- это давняя (так и не решённая проблема в грасс на Windows), в системах с нормальным UTF-8 всё намного лучше.
P.S. не зазорно иногда и на грассовском форуме спросить, если что, кстати

[Alexey13]

Amuriy, Для начала создаю каталог grassdata в корне HDD, дабы избавится от кириллицы. Модуль Grass включен.
* Создаю новый набор путь GRASS c базjq данных C:/grassdata/.
* Создаю район 123
* Выбираю проекцию WGS84
* Регион Kazakhstan, жму кнопку установить, называю сам набор и все, набор готов.
* Пробую r.in.gdal, выбираю 30.TIF с HDD, вписываю имя выходного файла 30.
* Жму на кнопку выполнить в некоторых случаях появляется ошибка "G_set_window(): ������������ ������ ��� ������
Завершено с ошибкой".
В большинстве случаев конвертация проходит успешно но после нажатия кнопки открыть вывод появляется ошибка "Слой не является действительным и не может быть добавлен на карту".

Я тоже заметил, что Grass очень требовательна к мелочам, поэтому сразу извиняюсь если вдруг ошибка окажется очевидной:)

[Александр Мурый]

* Выбираю проекцию WGS84

Растры точно в географич. координатах?

* Пробую r.in.gdal, выбираю 30.TIF с HDD, вписываю имя выходного файла 30. * Жму на кнопку выполнить в некоторых случаях появляется ошибка "G_set_window(): ������������ ������ ��� ������

Если работать в "родной" GRASS, то в таких случаях есть флаг "-l" для импорта растров в географ. координатах. Пинок в сторону "Инструментов GRASS": там нет почти никаких полезных опций из исходного модуля, в том числе и "-l".

Я тоже заметил, что Grass очень требовательна к мелочам, поэтому сразу извиняюсь если вдруг ошибка окажется очевидной:)

Ну не знаю, по-моему, кодировка и (тем более) системы координат -- это довольно сурьёзные штуки.