- 개요최근에 안개가 심하다는 문자를 받았습니다.안개는 기상현상이니 기상청에 들어가서 한 번 확인을 해봐야겠죠? 기상청 홈페이지에 들어가니 새로운 창에 안개와 관한 속보가 뜹니다. 자세히 보면 안개의 현황과 전망이 나오고 전국 관측소의 가시거리 그림이 나오네요. 요즘 무슨 주제할 지 고민이었는데 이거나 한 번 만들어보겠습니다.아마 다른 포스팅에서 쓰인 기술들을 잘 조합하면 만들 수 있을 것 같네요. - 안개 속보 화면 분석 및 계획짜기안개 속보 화면을 보면 크게 그림 부분과 텍스트의 두 부분으로 구성되어있으니 어떻게 작업을 할지 고민해봅시다. 아래의 2단계를 수행하면 되겠습니다.1. 프로그래밍 언어로 그림 파일 만들기2. 그림 파일을 불러와서 html 파일 만들기 제 생각에 기상청에서는 그림을 그리기 위..

GK2A(천리안 2호) 위성 데이터 처리 1. 다운로드 2. NC 파일 읽기 3. 간단한 시각화 지도에 ASOS 관측소 위치 표시하기 바로가기(지도에 ASOS 관측소 위치 표시하기) 관측 자료의 결측을 시각화하기 바로가기(관측 자료의 결측을 시각화하기) [기상청 API] ASOS 시간(hourly) 자료..

">  GK2, 레이더, ASOS 모두 기상청 API허브에서 받았습니다.자료 처리에는 파이썬, 시각화에는 Javascript의 three.js 라이브러리를 사용했습니다.근데 지금생각해보니 굳이 three.js 써서 만들 필요는 없었던듯

기상청은 3개월 전망을 할 때 특히 겨울철 전망에서 엘니뇨/라니냐, 북극진동을 자주 언급합니다.엘니뇨/라니냐는 교과서에서 배울 정도로 대중적인 개념이지만 북극진동은 그렇지 않죠. 여기서 호기심이 더 많은 분이라면 북극진동은 차가운 북극의 공기를 가두는 역할을 하며음의 북극진동일 때 북극 공기를 가두던 제트가 남쪽으로 내려와 차가운 공기도 같이 내려오게 되니우리나라가 춥다, 한파도 발생한다는 설명을 보셨을겁니다. 여러 과학, 기상 뉴스에서 이런 과정으로 그림으로 잘 설명하고 있습니다.하지만 이런 뉴스에서 실제 대기장을 보여주기보다 전문가들의 자문을 받아 모식도 같은 느낌의 그림을 보여줄겁니다. NOAA에서 제공하는 월 북극진동 지수와 NCEP2 재분석 자료의 대기 변수를 이용해서 음의 북극진동일 때 대기장..

- 개요북극 진동 패턴은 다 계산했으니 최대한 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축으로 ..