Методика
Получение маски
Для получения шейп-файла маски, основы был взяты следующие шейп-файлы:
Шейп-файл Маска Для получения береговой линии была взята ледовая обстановка Ber030718_1.shp (нет льда), контуры полигонов были приняты в качестве контура берега.
Для получения шельфа взята полигональная тема глубин, все полигоны со значением глубины до 200 метров объединены в один. Получившийся полигон одной стороной прилегает к суше, с другой ограничен глубиной 200 метров.
Полигональная тема с делением на полосы шириной 1 градус.
В результате комбинации этих трёх тем был получен шейп-файл mask_ber_1(https://docs.google.com/open?id=0B0AT94539mqXaXpneURQUEZ3ZTg).
Таблица этого файла содержит поля:
Newfield1 — может принимать значения:
shelf — часть Берингова моря прилегающая к суше (до глубины 200 м.)
out — остальная часть моря
land — суша..
Name —широта, может принимать значения 65-64, 64-63… 51-50, указывает, в каких широтах расположен полигон.
Аrеа — площадь полигона.
Newfield2 — составной идентификатор полигона, позволяющий определить к какой широте и к какому типу относится полигон (позволяет раскрасить тему в уникальные цвета), имя состоит из 3 символов TLL, где T (тип)может принимать значения 0 (море), 1 (шельф), 3 (океан), 4 (суша), LL широта южной границы, может принимать значения от 52 до 66.
Файлы маски сохранены в каталоге https://drive.google.com/#folders/0B0AT94539mqXNWlLajExS2Roc1k
Pасчет площадей
Средствами ArcView 3.2 выполнялось объединение темы ледовой обстановки и темы mask_ber_1.
При помощи стандартного скрипта CalculateAcreage вычислялись площади полигонов.
Полученные темы получали имя Union_1, Union_2 и т. д. Темы, соответствующие одному ледовому году сохранялись в папке с имененем winYY, например, в папке win97 содержаться шейп –файлы относящиеся к 1997 календарному году, а в папке win04 — к 2004.
. К таблице полученного объединения применялся запрос для получения полигонов содержащих только лёд и не попадающих на сушу или вне Берингова моря.
Шейп-файл ледовой обстановки
Шейп-файл объединения ледовой обстановки и маски
Select ([aобщ_спл] <> "НЕТ ЛЬДА" ) and
([aобщ_спл] <> "") and
([aобщ_спл] <> "НЕТ ДАННЫХ") and
([aобщ_спл] <> "ОТКРЫТАЯ ВОДА")and
([aобщ_спл] <> "НЕТ ЛЬДА" ) and
([aобщ_спл] <> "") and
([aобщ_спл] <> "НЕТ ДАННЫХ") and
([aобщ_спл] <> "ОТКРЫТАЯ ВОДА" ) and
([Newfield1] <> "out_be")and
([Newfield1] <> "land") and
([Newfield1] <> "") and
([Newfield1] <> "out_be")and
([Newfield1] <> "land") and
([Newfield1] <> "").
Запрос
Выбранные в результате запроса полигоны
Выбранные записи экспортировались в Excel.
В Excel таблицы компоновались в книги. Один ледовый год — одна книга, например, w97.xls содержит данные ледового года 1997 ода (с 4 ноября 1997 года по 6 июля 1998 г.).
На листе, соответствующем данным одной темы-объединения строилась сводная таблица для подсчета площадей полигонов, занятых льдом.
В каждой книге кроме есть лист с именем «Соответствие» указано соответствие между исходными темами ледовой обстановки и новыми, полученными в результате объединения, площадь занимаемая льдом на широте 63-62º N, на широте 58-57 (всего) и на широте 58-57 (шельф). В случае, если лед на определенной широте в заданную дату отсутствовал, то в таблице появляется сообщение #ССЫЛКА!, в этом случае значение площади льда принималось равным нулю.
Данные всех лет сведены в одну таблицу Ice_svod.xls.
Некоторые темы из набора Ледовая обстановка частично обрезаны, и на часть территории данные о площади льда отсутствуют. В этом случае в таблице площадей в соответствующем столбце появляется пропуск. Для заполнения пропуска использовались одно предшествующее и одно следующее значения, значения встроенная excelфункция Предсказ (линейной регрессии). В случае, если отсутствуют подряд два значения, для восстановления значений в этих точка также использовалась функция Предсказ и существующие значения до и после пропущенных.,
Каждой дате сопоставлялась одна из 3 декад и данные усреднялись по декадам (при помощи сводных таблиц Excel).
Уникальное имя декады состоит из 7 цифр (ГГГГММД) — четыре цифры год, две — месяц, последняя — номер декады.
Для построения графика площадь, занимаемая льдом делилась на площадь шельфа в соответствующей одноградусной полосе. На широте 58-57 площадь шельфа была принята 4038.738 км2, на широте 63-62 — 12697.481 км2.
Файлы с рассчитанными площадями сохранены в каталоге https://docs.google.com/open?id=0B0AT94539mqXZVNUUWFGVzlXSTg