Страница 1 из 1
Определить площадь растра одного цвета
Добавлено: 30 янв 2019, 20:15
_taras_
Доброго времени.
Я обрезал по нужной мне области, разбил снимок на одноканальное псевдоцветное изображение и записал отдельными слоями. Теперь нужно найти общую площадь всех пикселей слоя. Подскажите, пожалуйста, как это сделать?
Re: Определить площадь растра одного цвета
Добавлено: 30 янв 2019, 21:16
darsvid
Processing -> Raster analysis -> Raster layer unique values report
делает отчет по уникальным значениям растра - считает кол-во пикселей и площадь
Также в GRASS r.report
Re: Определить площадь растра одного цвета
Добавлено: 30 янв 2019, 21:42
_taras_
Спасибо за подсказку.
В результатах неясно
Код: Выделить всё
NODATA pixel count: 150319
Значение Pixel count Площадь (deg²)
0 77972 0.05417483069202151
1 9411 0.006538748930931802
Что значит NODATA pixel?
Площадь определяется умножением Pixel count на Площадь (deg²)? И в каких единицах измеряется?
Re: Определить площадь растра одного цвета
Добавлено: 30 янв 2019, 21:49
darsvid
Пиксели NoData используются для маскировки пикселей, которые не входят в вашу область интереса. На самом деле ваш растр прямоугольный, но именно благодаря тому, что все ненужное занято в нем пикселями NoData вы видите на картинке сложный контур Килийского устья.
Площадь считается в единицах проекции вашего растра (кол-во пикселей*разрешение растра), в данном случае - градусах. Для подсчета значений в м кв. растр следует перепроецировать из географической системы координат (WGS84) в прямоугольную (UTM, для этой территории 35 зона)
Re: Определить площадь растра одного цвета
Добавлено: 30 янв 2019, 22:39
_taras_
Делаю так. Выделяю нужный слой, растр->проекции-> перепроецирование
целевая СК
EPSG:3189-GR96 / UTM zone 35 N
остальное по умолчанию.
Результат
Код: Выделить всё
Значение Pixel count Площадь (m²)
0 77972 0.05416385670209518
1 9411 0.006537424401367386
Почему-то результат слишком маленький 9411×0,006537424401367386=61,523701041 м2?
Где ошибаюсь? Файл исходный приложил.
Re: Определить площадь растра одного цвета
Добавлено: 31 янв 2019, 12:06
darsvid
Потому что вашему начальному растру по ошибке присвоена неверная проекция (на самом деле он в WGS84, но присвоена UTM). Что нужно сделать:
1. Назначить исходному растру правильную проекцию с помощью Raster -> Projections -> Assign Projection
Установить Desired SRS - EPSG:4326, после выполнения растр закрыть и открыть снова и убедиться, что в Properties - Information написано CRS EPSG:4326 - WGS 84 - Geographic
2. Исправленный растр перепроецировать через Raster -> Projections -> Warp (Reproject)
Source CRS - EPSG:4326
Target CRS - EPSG: 32635
Resampling method - Nearest Neighbour
Опять же проверить в свойствах слоя разрешение (у меня получился Pixel Size 78.4702,-78.4702)
3. Посчитать статистику