- 개요북극 진동 패턴은 다 계산했으니 최대한 NOAA에서 제공하는 그림과 비슷하게 그려봅니다.길게 설명할 건 없고 그냥 코드로 들어가죠.  - NOAA 그림 따라 그리기1. 칼라맵 따오기뭐 어려운 건 없고 NOAA 그림을 캡쳐해서 그림판으로 연 다음색 선택 >>> 색 편집으로 NOAA 그림 칼라맵 색상의 rgb값을 가져온겁니다.colors에 저장하고 mcolors를 이용해서 custom_cmap으로 저장합니다.import matplotlib.colors as mcolorscolors = [ (0.078431373, 0.392156863, 0.823529412), # -45 ~ -40 (0.156862745, 0.509803922, 0.941176471), # -40 ~ -35 (0...

- 개요북극 진동 패턴을 구하기 위해 자료 처리를 했으니 이제 북극 진동 패턴을 구해봐야겠죠.EOF 분석을 한 후 여기서 나온 시계열을 다시 hgtAno에 regression하여 regression map을 구해야 합니다.저번에 구한 hgtAno 변수부터 시작합니다. - EOF 분석 이론EOF는 Empirical Orthgonal Function으로 직역하면 경험직교함수입니다.자세히 설명하기에는 포스트를 하나 새로 쓸 정도로 분량이 많아서 대략적으로만 설명합니다. 대기과학에서 EOF 분석을 사용하는 이유는 특정 지역에서 가장 변동성이 큰 패턴, 다른 말로는 가장 설명을 잘하는 (explainable) 공간 패턴을 찾고, 그 공간 패턴이 어떤 식으로 변하는지(공간 패턴의 시계열) 확인하고자 함입니다.이 때..

- 개요이번 포스트에서는 저번에 정한 자료 수집, 자료 처리 목표를 달성합니다.목표는 자료 수집에서 1개, 자료 처리에서 3개로 총 4가지입니다. - 자료 수집목표 1. NCEP2 재분석 자료 수집: NOAA 사이트에서 월평균 1000 hPa 지위고도 다운로드 (2, 3)먼저 아래 링크로 들어갑니다.https://psl.noaa.gov/data/gridded/data.ncep.reanalysis2.html다음 그림의 설명대로 하시면 NCEP2 재분석 자료인 hgt.mon.mean.nc파일을 다운로드 받을 수 있습니다.  - 자료 처리목표 2. 자료 추출: 북위 20도 이상, 1979~2000년 (4, 8)nc파일을 읽기 위해 xarray 라이브러리를 씁니다.isel과 sel 그리고 slice를 이용하여 ..

- 개요본격적인 코딩을 하기 전에 무엇을 해야하는지 정리해보고자 합니다.NOAA 홈페이지에 북극 진동 패턴을 구하는 과정이 잘 정리되어있습니다. https://www.cpc.ncep.noaa.gov/products/precip/CWlink/daily_ao_index/history/method.shtml - 북극 진동 패턴 구하는 방법론 해석NOAA 홈페이지의 문단을 한 문장씩 읽으면서 북극 진동 패턴을 구하기 위한 정보들을 모아봅시다. To identify the leading teleconnection patterns in the atmospheric circulation, Emperical Orthogonal Function (EOF) was applied to the monthly mean 1000..

- 북극 진동 패턴이란?이제 슬슬 겨울이 다가오고 있습니다.기상, 기후 뉴스에서 매년 엘니뇨 이야기가 나오는 것처럼 겨울철에는, 특히 우리나라 한파가 올 때 "북극 진동"이라는 키워드를 심심치 않게 볼 수 있습니다. 그리고 뉴스에 나오는 논리는 대충 이럴겁니다.1. 북극의 찬 공기를 가두는 소용돌이가 있다.2. 평소에는 이 소용돌이가 북극의 찬 공기를 가두기 때문에 찬 공기는 고위도 지역에 머무른다.3. 이러저러한 이유로 소용돌이가 약해져서 원래는 고위도에 있던 북극의 찬 공기가 중위도로 남하한다.4. 북극의 찬 공기 때문에 우리나라의 온도는 급격히 떨어져서 한파가 발생한다. 소용돌이의 세기를 숫자로 표현한다면 소용돌이가 쌔면 숫자가 크고, 약하면 숫자가 작겠죠?소용돌이의 세기를 y축, 시간을 x축으로 ..

- 개요이제 그림 그리는 법은 다 알았으니 애니메이션을 만들 차례입니다.여기서 애니메이션은 동영상 파일이 아니라여러 사진을 연속으로 보여주는 gif 파일로 저장됩니다.matplotlib의 함수를 쓰면 쉽게 gif 파일을 만들 수 있습니다. - matplotlib으로 그림 그리는 과정에 대한 이해matplotlib의 애니메이션을 만들어주는 함수를 이해하기 위해선matplotlib 라이브러리에서 그림이 어떤 과정으로 그려지는지 알 필요가 있습니다. 쉽게 설명하려고 비유를 했기 때문에 완벽한 설명을 아닙니다. 먼저 fig=plt.figure()는 그림을 그릴 때 쓰는 canvas를 준비하는 것입니다.ax=fig.add_subplot()은 canvas 내에서 그림을 그릴 도화지를 준비하는 과정입니다. 다음으로 ..

- 개요기후 나선에서는 시간에 따라 과거의 기온선부터 현재의 기온선까지 순차적으로 그려져야 합니다.순차적으로 그리기 위해서는 애니메이션 효과를 써야하며, 그 전에 기온선을 그려야 합니다.이번 포스트에서는 기온선을 그려보겠습니다. - 모든 연도의 기온선 그리기1. 선 그리기선을 그리는 것은 어렵지 않습니다만이전 포스트처럼 극좌표계의 위치를 기후나선에 맞게 수정을 해야합니다.극좌표계에서 0도는 가장 오른쪽 방향이고, 각도가 증가할수록 반시계방향으로 돌아갑니다.반면 기후나선에서는 극좌표계의 90도의 위치가 1월이고 2, 3, 4월은 시계방향의 순서로 위치합니다.즉 시작하는 위치와 회전방향을 수정해줘야 합니다.import matplotlib.pyplot as pltfrom matplotlib.animation ..

- 개요저번 포스트에서는 기온 자료의 전처리 작업을 했습니다.이번 포스트에서는 기온 자료를 시각화하기 전에NASA Climate change 기후 나선 그림에서 1, 0, -1도를 나타내는 원과 그림 중앙의 연도, 외곽의 1~12월 문자를 밑그림으로 그립니다.  - 밑그림 그리기원을 그리는 것과 1월부터 12월을 텍스트로 적는 것은 반복 작업이므로 하나하나 적지말고 함수를 만듭니다.저는 극좌표계로 생각하고 이를 x ,y 좌표로 변환해서 그림을 그렸습니다.코드에서 r은 반지름을 의미합니다. 1. 기본 그림matplotlib으로 그림을 그리면 기본 하얀 바탕에 축이 그려집니다.우린 검은 바탕을 써야하고, 축을 지워야 합니다.import matplotlib.pyplot as pltfig = plt.figure..

- 개요이번 포스트에서는 서울 47108 기온 자료의 전처리 작업을 수행합니다.여러분이 서울 470108 일평균 기온 자료를 가지고 있다고 가정하겠습니다. - 일평균 기온 자료를 월평균 기온 자료로 변환주의할 점으로 저는 이미 기온 자료의 기간을 1908년 1월부터 2023년 12월까지,즉 시작은 1월부터 마지막은 12월로 맞추어두었습니다.만약 이렇게 기간을 잡지 않았다면 월평균 자료로 바꾼 다음 시작은 1월, 마지막은 12월로 맞춰주세요.import numpy as npimport pandas as pddf = pd.read_csv('ASOS_108_daily.csv', encoding='cp949')col1 = df['일시']col2 = df['평균기온(℃)']''' pandas data frame으..