위성, 레이더로 보는 2024년 11월 폭설사태
·
프로그래밍/시각화
">  GK2, 레이더, ASOS 모두 기상청 API허브에서 받았습니다.자료 처리에는 파이썬, 시각화에는 Javascript의 three.js 라이브러리를 사용했습니다.근데 지금생각해보니 굳이 three.js 써서 만들 필요는 없었던듯
지도로 보는 우리나라 연도별 기상관측장비 (ASOS, AWS) 설치 현황
·
프로그래밍/시각화
예전에 지도에 지도에 ASOS 관측소 위치 표시하기라는 포스트를 쓴 적이 있습니다(https://climate-data-science.tistory.com/5)이 때는 파이썬의 matplotlib 라이브러리를 이용해서 ASOS 관측소의 위치를 한 장의 그림으로 그렸습니다. 이번에는 ASOS뿐만 아니라 AWS의 위치를 표시하되 이를 연도 순으로 차례차례 시각화해보았습니다.자료는 기상청에서 구했고, 관측이 종료된 관측소는 처음부터 제외했습니다.  " data-ke-type="html">HTML 삽입미리보기할 수 없는 소스
3차원으로 보는 서울의 기후 나선
·
프로그래밍/시각화
- 개요예전에 matplotlib을 이용해서 여러 장의 png 파일을 만들고 이를 gif로 합쳐 서울의 기후 나선을 그렸습니다. 원래 NASA Climate Change 유투브에 나온 영상에서는 2차원의 기후 나선을 그려진 뒤 카메라가 3차원으로 움직입니다.하지만 matplotlib으로 3차원 처리를 하고 시간마다 특정 효과를 넣어 동영상으로 저장하기 어렵죠.그렇기 때문에 이러한 작업처리가 용이한 Javascript를 써서 시각화를 했습니다. - 자료 처리matplotlib 기후 나선 포스트에서는 월평균을 기온을 이용했습니다만1년이 12개월이므로 기후 나선은 12각형이 되므로 이쁘지가 않습니다.시간해상도를 올려서 일평균값으로 기후나선을 그리면 최소 365각형이 되어 원처럼 보입니다.Base period는..
47108 서울 ASOS 월평균 기온으로 본 지구온난화
·
프로그래밍/시각화
47108 서울 ASOS는 1907년부터 지금까지 기상 요소를 관측하고 있습니다.무려 100년이 넘게 관측이 되어 100년 동안 정말 지구온난화가 있었는지 확인 가능하겠죠?다만 슬프게도 한국전쟁기간인 1950년부터 1953년까지 결측이 있긴하지만요. 기준 온도를 정의하자.지구온난화를 시각화할 때는 기준 온도를 먼저 정의해야 합니다.어떤 기간보다 지금이 얼마나 뜨거운지 알아야 하거든요. 보통은 이산화탄소 배출이 심해지기 전(산업혁명 전)인 1850년부터의 관측값을 써야하지만 저는 1908년부터 1949년을 기준온도로 쓰겠습니다.1월부터 12월까지 각각의 기준온도를 계산해서 각 월에서 빼줍니다. 예를 들어 아래에 제가 시각화 하는 2024년 7월 평균 기온은 그냥 2024년 7월 평균 기온이 아니라 1908..
2024년 여름 얼마나 더운거지?
·
프로그래밍/시각화
요즘 날씨에 대해 남들에게 물어보면 진짜 너무 덥고 습하다는 말을 자주 듣습니다. 과거와 비교해서 지금은 얼마나 덥고 습한걸까요? 서울(47018) ASOS 관측소의 평균기온, 평균습도 자료를 이용해서6, 7, 8월 상순, 중순, 하순 평균값으로 순위를 매겨봅시다. 1908년부터 2024년 8월 중순까지 시간순서대로 평균습도, 평균기온을 시각화하는 동영상을 만들었습니다." data-ke-type="html">HTML 삽입미리보기할 수 없는 소스  그럼 순위를 표로 살펴보죠(동영상의 마지막 프레임입니다). 평균기온 6위를 보시면 2024년 8월 중순이 있습니다.2024년 8월 상순은 5위에 있네요.또한 평균기온 상위권 중에서 습도도 좀 높은 편이긴합니다.참고로 제가 따로 확인해보니 2024년 6, 7월은 ..
여러 nc 파일을 빠르게 읽으려면?
·
프로그래밍/파이썬
- 개요nc파일은 NetCDF(Network Common Data Form) 파일로 대기과학 분야에서 많이 사용하는 파일 형식입니다.특히 위도, 경도, 연직층의 좌표계를 갖는 자료는 보통 nc 파일로 저장되어 있습니다. 예전에 저는 ERA5 재분석 자료 1달치를 하나의 nc파일로 저장했었는데요.이 nc 파일에는 특정 변수가 [30일 x 24시간 x 위도격자 x 경도격자 x 연직층]의 차원으로 저장되므로 용량이 꽤 큽니다.물론 좁은 영익인 한반도 근처 자료만 받았지만 30년치라 그냥 이 자료를 모두 읽는 것만 해도 시간이 꽤 걸립니다. 그럼 읽는 속도를 비교하며 nc 파일을 가장 빨리 읽는 방법을 알아보겠습니다.결론을 먼저 말하자면 xarray에서 병렬처리로 파일을 읽는 것이 제일 빠릅니다. - 사용할 라..
[Matplotlib] 타이틀(title)에서 한글 사용하기
·
프로그래밍/파이썬
- 개요우리는 matplotlib 라이브러리를 사용해서 다양한 그림을 그립니다.하지만 Matplotlib의 기본 설정에서는 figure의 제목이나 x축, y축 타이틀을 한글로 적을 수 없습니다.그럼 어떻게 해야 matplotlib의 그림에서 한글을 쓸 수 있는지 알아보겠습니다. - Matplotlib에서 한글 사용하기 (개별 설정)한글 폰트를 사용하는 방법은 간단합니다.한글 폰트를 사용하고자 부분에서 fonproperties를 설정해주면 됩니다.import matplotlib.pyplot as pltimport matplotlib.font_manager as fmfn_font = 'NanumGothic.ttf' # 나눔고딕 폰트 사용fig, ax = plt.subplots(1, 2, figsize=(6,..