Страница 1 из 1
Подсчет гектаров по спутниковому снимку
Добавлено: 13 авг 2011, 19:08
Alexey13
Всем добрый вечер!
Помогите пожалуйста, как можно подсчитать количество гектаров в растровом снимке GeoTiff?
Снимки проективного покрытия примерно по 400 га,(около 300 шт.) с 10 классами 0-10%, 10-20%, 20-30%, 30-40%, и.т.д. до 100%. В каждой карте по 4-5 классов. Необходимо узнать площадь каждого класса, можно даже в процентах от общай площади.
Из инструментов QGIS, Global Mapper, Erdas в процессе установки...
Заранее спасибо!
Re: Подсчет гектаров по спутниковому снимку
Добавлено: 13 авг 2011, 20:59
novia
В эрдасе делаете модель в model maker, состоящую из одного снимка на входе, условного оператора и одного снимка на выходе. В условном оператора применяете границы классов к пикселам снимка и записываете номер класса в новых файл.
Вот
тут написано как запустить модель для 100500 файлов.
А
тут писали про то, как посчитать количество пикселов каждого класса.
Соответственно "гектары" - это количество пикс. * площадь пикс. в кв. м /10000
[Сообщение с мобильного устройства] 
Re: Подсчет гектаров по спутниковому снимку
Добавлено: 14 авг 2011, 16:48
Alexey13
Спасибо большое. Но способ описанный в Эрдасе слишком сложный для моего текущего уровня

.
Существует ли более простой способ выделения контура по цвету? Данная операция легко проделывается в Гимпе и фотошопе, но площадь определить нельзя.
Спасибо!
Re: Подсчет гектаров по спутниковому снимку
Добавлено: 14 авг 2011, 21:21
novia
Все таки 300 снимков "не руками" обработать не совсем тривиальная задача.
Попробуйте разобраться с model maker и запуском пакетной обработки, они только с виду страшные и непонятные. Задавайте тут вопросы, у вас все получится.
Re: Подсчет гектаров по спутниковому снимку
Добавлено: 14 авг 2011, 22:11
Александр Мурый
В QGIS можно исп-ть "Инструменты GRASS": меню Растр --> Отчёты и статистика --> Статистика --> r.report
Вот ссылка на мануал. Работа с "Инструментами GRASS" описана в русскоязычном руководстве (лежит здесь на сайте).
Re: Подсчет гектаров по спутниковому снимку
Добавлено: 14 авг 2011, 23:02
Максим Дубинин
Re: Подсчет гектаров по спутниковому снимку
Добавлено: 16 авг 2011, 21:25
Alexey13
Второй день мучаю модуль Grass в Qgis. Выходит ошибка - неверный слой "Слой не является действительным и не может быть добавлен на карту". Ошибка выходит при импорте растра Gdal в Grass, хотя все делаю по руководству. И кстати почему то слетает кодировке в модуле Grass в путях к файлам. Возможно причина в этом?
Re: Подсчет гектаров по спутниковому снимку
Добавлено: 16 авг 2011, 21:28
Alexey13
Спасибо большое. Попробовал быстро не получилось и пока отложил, увидев возможности Grass. Но если с не получится вернусь к упомянутому вами способу.
Re: Подсчет гектаров по спутниковому снимку
Добавлено: 16 авг 2011, 21:43
Александр Мурый
Alexey13, опишите последовательность действий с GRASS через QGIS. Главное правило здесь: по возможности никаких пробелов и других "нехороших" символов в путях в файлам и в именах файлов и минимум кириллицы.
С кодировками -- это давняя (так и не решённая проблема в грасс на Windows), в системах с нормальным UTF-8 всё намного лучше.
P.S. не зазорно иногда и на грассовском форуме спросить, если что, кстати

Re: Подсчет гектаров по спутниковому снимку
Добавлено: 16 авг 2011, 22:19
Alexey13
Amuriy, Для начала создаю каталог grassdata в корне HDD, дабы избавится от кириллицы. Модуль Grass включен.
* Создаю новый набор путь GRASS c базjq данных C:/grassdata/.
* Создаю район 123
* Выбираю проекцию WGS84
* Регион Kazakhstan, жму кнопку установить, называю сам набор и все, набор готов.
* Пробую r.in.gdal, выбираю 30.TIF с HDD, вписываю имя выходного файла 30.
* Жму на кнопку выполнить в некоторых случаях появляется ошибка "G_set_window(): ������������ ������ ��� ������
Завершено с ошибкой".
В большинстве случаев конвертация проходит успешно но после нажатия кнопки открыть вывод появляется ошибка "Слой не является действительным и не может быть добавлен на карту".
Я тоже заметил, что Grass очень требовательна к мелочам, поэтому сразу извиняюсь если вдруг ошибка окажется очевидной:)
Re: Подсчет гектаров по спутниковому снимку
Добавлено: 16 авг 2011, 23:12
Александр Мурый
Alexey13 писал(а):
* Выбираю проекцию WGS84
Растры точно в географич. координатах?
Alexey13 писал(а):
* Пробую r.in.gdal, выбираю 30.TIF с HDD, вписываю имя выходного файла 30. * Жму на кнопку выполнить в некоторых случаях появляется ошибка "G_set_window(): ������������ ������ ��� ������
Если работать в "родной" GRASS, то в таких случаях есть флаг "-l" для импорта растров в географ. координатах. Пинок в сторону "Инструментов GRASS": там нет почти никаких полезных опций из исходного модуля, в том числе и "-l".
Alexey13 писал(а):
Я тоже заметил, что Grass очень требовательна к мелочам, поэтому сразу извиняюсь если вдруг ошибка окажется очевидной:)
Ну не знаю, по-моему, кодировка и (тем более) системы координат -- это довольно сурьёзные штуки.