как подсчитать частоту изменения посевов

[li_sasha]

Привет всем.
есть задача, нужно вычислить частоту изменения посевов в период с 2001-2009 годов.
может кто то уже делал подобное?
- из софта имею ERDAS 2010 и ArcGIS 10.
- начальные данные - растровые, разделенные на 9 цветов. Каждый цвет соответствует какому либо виду посева.

В ERDAS вычислил "выявление изменений" (использовал инструмент "Change detection"). То есть получил растр с разницей между 2мя годами.
Что нужно сделать дальше? может я не с той стороны начал?

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

- какие данные собираетесь использовать?
- что такое "посев"?
- что такое "частота изменения посевов"?
- что из себя представляют "начальные данные", какова их роль?

[li_sasha]

- какие данные собираетесь использовать?
- что такое "посев"?
- что такое "частота изменения посевов"?
- что из себя представляют "начальные данные", какова их роль?

Ответ:
- данные в виде растра. То есть tif, в котором имеются около 9 цветов. каждый цвет это какой либо посев, например красный цвет это кукуруза, синий вода и т.д. Имеются растры с 2001 по 2009 год. А точнее каждый растр показывает что в этом году было посеяно.
- посев это вид сельскохозяйственной культуры который был посеян на какой либо территории. Например пискель красного цвета означает что там был посеян хлопок. (В принцыпе, не важно что означают каждые цвета, мне нужно подсчитать частоту изменения этих пикселей)
- под "частотой изменения посевов" я подразумеваю на сколько часто изменился цвет пикселя на каждом участке. Например, возьмем один пиксель с 2001 года и один пиксель того же участка за 2002 год и так далее, если за 2001 год он был красным, то за 2002 год он стал синим, значит он изменился. мне нужно подсчитать частоту изменения, то есть на сколько часто он именился. Например, был хлопок, стал рисом, потом стал пшеницей и так далее.
- начальные данный представляют из себя TIF файлы, каждый файл это показания за год. Данные за 2001 по 2009 год.

заранее извеняюсь, если слишком много написал, просто хотелось объяснить как можно проще, так сказать нулями и единицами

[Voltron]

Смотрите в сторону этой статьи

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

не стоит экономить на постановке задачи, уже понятнее, но еще вопросы:

- в каком виде вам нужно получить результат? вам нужно посчитать частоту по каждому пикселю? Или некий интегральный показатель по всей территории?
- ваше описание того, как вы собираетесь считать "частоту изменений" плохо формализовано, из него непонятно, что именно вы хотите получить в результате? Целочисленное число показывающее количество изменений? Т.е. максимальное число 8 будет обозначать что пиксель меняется с культуры на культуру ежегодно, так? Или если он возвращается к начальной то это уже не изменение и вам нужно количество уникальных культур?

[li_sasha]

не стоит экономить на постановке задачи, уже понятнее, но еще вопросы:

- в каком виде вам нужно получить результат? вам нужно посчитать частоту по каждому пикселю? Или некий интегральный показатель по всей территории?
- ваше описание того, как вы собираетесь считать "частоту изменений" плохо формализовано, из него непонятно, что именно вы хотите получить в результате? Целочисленное число показывающее количество изменений? Т.е. максимальное число 8 будет обозначать что пиксель меняется с культуры на культуру ежегодно, так? Или если он возвращается к начальной то это уже не изменение и вам нужно количество уникальных культур?

ответ:
во первых огромное спасибо за ответ. из ваших вопросов уже вытекают ответы.
итак:
- мне нужно подсчитать результат по каждому пикселю.
- На английском это звучит так:
- "Frequency of crop rotations for a pixel during 2001-2009"
- "Predominant crop or land use for a pixel during 2001-2009"
как я понял, мне нужна частота изменения культуры и преобладающая культура в пикселе. То есть, сколько раз изменились культы и какая культура больше всего была посажена на данном пикселе.

[li_sasha]

я подсчитал статистику преобладающих культур. Но с подсчетом количества изменений всё еще проблемы. В ERDAS, используя инструмент MATRIX был получен результат в котором показывается какая культура изменилась на другую. Другими словами, если был цвет 1 то стал цвет 2 и т.д. Как мне подсчитать количество изменений на каждом пикселе?

Заранее благодарю всех ответивших

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

алгоритм может быть такой

последовательно вычитаете след.год из предыдущего d1=1-2, d2=2-3, d3=3-4 и т.д.
пересчитываете разности: if d1 <> 0 then 1 и т.д.
складываете все разности вместе result = d1+d2+...d8

result - количество изменений на пиксел, минимум 0, максимум 8

[li_sasha]

алгоритм может быть такой

последовательно вычитаете след.год из предыдущего d1=1-2, d2=2-3, d3=3-4 и т.д.
пересчитываете разности: if d1 <> 0 then 1 и т.д.
складываете все разности вместе result = d1+d2+...d8

result - количество изменений на пиксел, минимум 0, максимум 8

спасибо за ответ.
а какими инструментами в ERADS или ArcGIS это можно сделать?

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

в аркгисе в spatial analyst есть калькулятор, должно быть достаточно

[li_sasha]

в аркгисе в spatial analyst есть калькулятор, должно быть достаточно

спасибо огромное за помощь. Вы мне очень помогли. я сделал это, но немного другим способом.

У меня появилась следующая задача. Мне нужно рассчитать наибольшее число повторений.
Например, возьмем один пиксель. Он изменился 9 раз с такой последовательностью 1,1,1,1,1,0,0,1,1. Мне нужно найти наибольшее количество повторений . В данном примере наибольшее количество повторений рана 5-ти.
Не знаю как это сделать.

заранее благодарю за ответ.