전체

- np.einsum()이란?numpy 라이브러리에는 각종 연산을 해주는 함수가 있습니다.그중 아인슈타인 표기법을 함수로 만든 np.einsum()이 있습니다.아인슈타인 표기법의 자세한 내용은 다음 링크를 참고해주세요. https://rockt.ai/2018/04/30/einsum numpy 변수 관점에서 설명하면 아인슈타인 표기법은 axis(행과 열)에 따라 특정 규칙으로 계산을 하는 것입니다.좀 더 간단하게 설명하기 위해 3x3 행렬을 생각해봅시다. 이는 3x3 shape의 numpy 배열이라 봐도 되겠죠.행렬 안에 있는 원소 aij라치면 i는 행의 인덱스, j는 열의 인덱스입니다.a31이면 행렬 내에서 3행 1열에 위치한다는 뜻이죠.아인슈타인 표기법에서는 행렬 내의 원소를 ij와 같은 인덱스를 이용..

- 개요저번 포스트에서는 S-DoT 자료를 전처리하고 간단히 그림을 그려보았습니다.2024년 8월 13일에 내린 집중호우는 서울 전역을 강타한 것이 아니라 특정 지역에만 내렸습니다.ASOS/AWS 관측소가 서울에 30개나 있지만 이 정도로는 국지성 호우의 영향을 보기 어렵습니다.하지만 서울에 1000개나 설치한 S-DoT이라면 국지성 호우의 영향을 확인할 수 있겠죠. 집중호우가 내린 곳은 온도가 급격히 떨어지기 때문에 레이더 에코의 위치와 S-DoT의 온도가 감소하는 위치가 거의 같을 겁니다.그러므로 2024년 8월 13일 시간별로 레이더 에코와 S-DoT 온도장을 확인해보겠습니다.우리는 변화를 봐야하기 때문에 정확히는 온도장을 그리는 것이 아니라 1시간 동안 변한 온도의 공간장을 그려야 합니다.만약 진..

- 개요한국에서 기상 관측 자료라 하면 보통 ASOS/AWS 기상 관측 자료를 말합니다.기상청뿐만 아니라 다른 기관이 설치한 ASOS/AWS 관측소의 수는 전국에 500개가 넘는데요.다른 나라와 비교해서 우리나라 정부가 정말 기상 관측소를 정말 촘촘히 설치한 겁니다.심지어 국가가 관리하기 때문에 관측품질도 좋죠.이렇게 촘촘히 설치했지만 서울에 설치되어있는 기상청 ASOS/AWS 관측소 숫자는 약 30개입니다.서울에 25개 행정구가 있으니 1구에 1개 정도 있는 셈이죠.  이외에도 2020년부터 서울에서는 IoT 센서가 설치되어 기상 변수가 측정되고 있습니다.서울 곳곳에 약 1000개의 IoT 센서가 설치되어 있고, 기상 변수 외에 다른 변수도 측정합니다. 기간은 2020년 4월부터 현재까지 시간은 1시간..

"> 자료처리에는 파이썬, 시각화에는 자바스크립트 three.js 라이브러리를 사용했습니다.현재 인공위성 중 절반 이상이 스타링크 위성이라고 합니다.눈에는 잘 보이지 않지만 하늘에는 위성이 참 많네요.

- 개요지금까지 구한 텍스트, 그림으로 html 파일을 만들어봅시다.파이썬으로 html 파일을 직접 작성할 것이고, html/css의 기초는 알고 있으셔야 합니다. - html 파일 만들기파이썬 코드의 정보 중 총 3개가 html 파일로 전달되어야 합니다.1. 시간안개 속보의 2번째 줄을 보면 속보가 나온 시간이 있습니다.그림을 그린 시간보다 10분을 더해줘야 합니다. 2. 안개 현황 텍스트2번째 전 포스트에서 안개 현황에 넣을 텍스트를 만들었습니다.텍스트를 html에 넣기 좋게 태그를 입혀야합니다. 3. 시정 자료 그림저번 포스트에서 vis.png라는 이름으로 저장했습니다. 4. 배치전 아래와 같은 구조로 각 요소를 배치했습니다.점이 찍혀있으면 클래스 이름입니다(e.g., div .container >..

- 개요저번 포스트에서 자료를 다 처리했으니 시정거리 자료로 그림을 그릴 차례입니다.최대한 안개 경보 화면의 그림과 똑같이 그려보죠. - 그림 그리기안개 경보 화면 그림을 보며 어떻게 그릴지 생각해봅시다.위도, 경도로 위치를 잡고 시정거리에 따라 다른 색으로 scatter plot을 그리면 되겠네요.칼라바의 색을 따오고, 칼라바 위치를 조정해야 합니다. 지금 보니 시정 자료뿐만 아니라 한글을 쓰기 위한 font 설정, 도 경계를 그리기 위한 shp파일도 필요하네요.한글 폰트를 넣는 방법은 이 글을 보시고, 한국 행정구역 shp 파일에 대한 것은 이 글을 보시면 됩니다. 먼저 안개 경보 그림의 칼라바 rgb값은 제가 다 땃으니 아래의 값을 쓰세요.import matplotlib.colors as mcolo..

- 개요기상청 API허브에서 필요한 자료를 다운로드 받습니다.먼저 API로 받은 AWS 시정거리 자료를 pandas.DataFrame 형식으로 저장합니다. 다음으로 AWS가 어디 권역에 속하는지 알아야합니다.구역 4를 보시면 "권역: 관측지점 시정거리" 이런 형식의 글이 있습니다.  시정거리가 작은 지역을 권역별로 알려주는 듯합니다.AWS의 위치 정보는 마찬가지로 기상청 API 허브에 있습니다.다만 이 자료를 pandas.DataFrame에 넣기 위해 파싱(parsing)하는 과정이 생각보다 까다롭습니다. 마지막으로 시정거리와 위치 정보를 이용해서 구역4에 넣을 문자열을 생성합니다.  - 시정거리 자료 다운로드 및 처리기상청 API허브 홈페이지로 들어가서  좌측 지상관측탭 > 상단 방재기상관측(AWS)탭..

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

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