Методика

Получение маски

Для получения шейп-файла маски, основы был взяты следующие шейп-файлы:

  1. Для получения береговой линии была взята ледовая обстановка Ber030718_1.shp (нет льда), контуры полигонов были приняты в качестве контура берега.
  2. Для получения шельфа взята  полигональная  тема глубин, все полигоны со значением глубины до 200 метров объединены в один. Получившийся полигон одной стороной прилегает к суше, с другой ограничен глубиной 200 метров.
  3. Полигональная тема  с делением на полосы шириной 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асчет площадей


  1. Средствами ArcView 3.2 выполнялось объединение темы ледовой обстановки и темы mask_ber_1.
  2. При помощи стандартного скрипта CalculateAcreage вычислялись площади полигонов.
  3. Полученные темы получали имя Union_1, Union_2 и т. д. Темы, соответствующие одному ледовому году сохранялись в папке с имененем winYY, например, в папке  win97 содержаться шейп –файлы относящиеся к 1997 календарному году, а в папке win04 — к 2004.
  4. . К таблице полученного объединения применялся запрос  для получения полигонов содержащих  только лёд и не попадающих на сушу или вне Берингова моря.
  5.  
    Шейр-файл Ber060306
     
    Шейп-файл объединения ледовой обстановки и маски
    Шейп-файл ледовой обстановки Шейп-файл объединения ледовой обстановки и маски

     

    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] <> "").

     
    Выбранные полигоны
    Запрос Выбранные в результате запроса полигоны

  1. Выбранные записи  экспортировались в Excel.
  2. В Excel таблицы компоновались в книги. Один ледовый  год — одна книга, например, w97.xls содержит данные ледового года 1997 ода (с 4 ноября 1997 года по 6 июля 1998 г.).
  3. На листе, соответствующем данным одной темы-объединения строилась сводная таблица для подсчета площадей полигонов, занятых льдом.
  4. В каждой книге кроме  есть лист с именем  «Соответствие» указано соответствие между исходными темами ледовой обстановки и новыми, полученными в результате объединения, площадь занимаемая льдом  на широте 63-62º N, на широте 58-57 (всего) и на широте 58-57 (шельф). В случае, если лед на определенной широте в заданную дату отсутствовал, то в таблице появляется сообщение #ССЫЛКА!, в этом случае значение площади льда принималось равным нулю.

    Ледовая обстановка

    Объединение

    63-62

    58-57 всего

    5857 шельф

    ber081020_1

    union_1

    #ССЫЛКА!

    #ССЫЛКА!

    #ССЫЛКА!

    ber081103_1

    union_2

    #ССЫЛКА!

    #ССЫЛКА!

    #ССЫЛКА!

    ber081117_1

    union_3

    225.4

    #ССЫЛКА!

    #ССЫЛКА!

    ber081201_1

    union_4

    9054.5

    #ССЫЛКА!

    #ССЫЛКА!

    ber081208_1

    union_5

    10797.9

    #ССЫЛКА!

    #ССЫЛКА!

    ber081215_1

    union_6

    3429.0

    2050.7

    1883.4

    ber081222_1

    union_7

    2384.3

    563.4

    563.4

    ber081229_1

    union_8

    12411.3

    576.0

    576.0

    ber090608_1

    union_30

    325.3

    2816.0

    2807.8

    ber090615_1

    union_31

    6.1

    2574.6

    2407.4

    ber090622_1

    union_32

    909.8

    2794.7

    1570.3

    ber090629_1

    union_33

    #ССЫЛКА!

    #ССЫЛКА!

    #ССЫЛКА!

    ber090706_1

    union_34

    #ССЫЛКА!

    #ССЫЛКА!

    #ССЫЛКА!

    ber090713_1

    union_35

    #ССЫЛКА!

    #ССЫЛКА!

    #ССЫЛКА!

     

    Данные всех лет сведены в одну таблицу Ice_svod.xls.

    Некоторые темы из набора Ледовая обстановка частично обрезаны, и на часть территории данные о площади льда отсутствуют. В этом случае в таблице площадей в соответствующем столбце появляется пропуск. Для заполнения пропуска использовались  одно предшествующее и одно следующее значения, значения  встроенная excel функция  Предсказ (линейной регрессии). В случае, если отсутствуют подряд два значения, для восстановления значений в этих точка также использовалась функция Предсказ и существующие значения до и после пропущенных.,

    Каждой дате  сопоставлялась одна из 3 декад и данные усреднялись по декадам (при помощи сводных таблиц Excel).
    Уникальное имя декады состоит из 7 цифр (ГГГГММД) — четыре цифры год, две — месяц, последняя — номер декады.

    Среднее по полю площадь 63

     

     

    Среднее по полю площадь 58

     

    год_месяц_декада

    Итог

     

    год_месяц_декада

    Итог

    1997111

    0.0

     

    1997111

    0.0

    1997112

    1198.3

     

    1997112

    826.9

    1997113

    3593.3

     

    1997113

    416.9

    1997121

    6860.6

     

    1997121

    585.2

    1997122

    12219.2

     

    1997122

    3081.5

    1997123

    11344.0

     

    1997123

    3465.3

    1998011

    12087.6

     

    1998011

    4106.4

    1998012

    12208.1

     

    1998012

    4827.8

     

    Для построения графика площадь, занимаемая льдом делилась на площадь шельфа в соответствующей одноградусной полосе. На  широте  58-57 площадь шельфа была принята  4038.738 км2, на широте 63-62 — 12697.481 км2.

  5. Файлы с рассчитанными площадями сохранены в каталоге https://docs.google.com/open?id=0B0AT94539mqXZVNUUWFGVzlXSTg

     

Comments