KR101471509B1 - 코로나 질량 방출 분석 장치 및 그 방법 - Google Patents

코로나 질량 방출 분석 장치 및 그 방법 Download PDF

Info

Publication number
KR101471509B1
KR101471509B1 KR1020140099021A KR20140099021A KR101471509B1 KR 101471509 B1 KR101471509 B1 KR 101471509B1 KR 1020140099021 A KR1020140099021 A KR 1020140099021A KR 20140099021 A KR20140099021 A KR 20140099021A KR 101471509 B1 KR101471509 B1 KR 101471509B1
Authority
KR
South Korea
Prior art keywords
corona
image
stereo
corona mass
mass
Prior art date
Application number
KR1020140099021A
Other languages
English (en)
Inventor
오승준
이정덕
이상우
홍순학
김재훈
한진욱
김영규
Original Assignee
(주)에스이랩
대한민국(국립전파연구원 우주전파센터장)
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by (주)에스이랩, 대한민국(국립전파연구원 우주전파센터장) filed Critical (주)에스이랩
Priority to KR1020140099021A priority Critical patent/KR101471509B1/ko
Application granted granted Critical
Publication of KR101471509B1 publication Critical patent/KR101471509B1/ko

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/23—Clustering techniques
    • G06F18/231—Hierarchical techniques, i.e. dividing or merging pattern sets so as to obtain a dendrogram
    • G06T3/16—
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10032—Satellite or aerial image; Remote sensing

Abstract

본 발명은 코로나 질량 방출 분석 장치 및 그 방법을 개시한다. 즉, 본 발명은 소호 위성과 스테레오 위성의 코로나 그래프의 실시간 수집 체계를 구축하고, 상기 실시간 수집된 코로나 그래프를 분석하여 코로나 질량 방출 발생을 자동으로 인식함으로써 신뢰성 있는 데이터를 제공하고 코로나 질량 방출의 방향, 방출폭, 초기 속도를 자동 인식하고, 3차원 CME 물리량을 제공하여 사용자의 편의성을 향상시킬 수 있다.

Description

코로나 질량 방출 분석 장치 및 그 방법{Apparatus for analyzing coronal mass ejection and method thereof}
본 발명은 코로나 질량 방출 분석 장치 및 그 방법에 관한 것으로, 특히 소호 위성과 스테레오 위성의 코로나 그래프의 실시간 수집 체계를 구축하고, 상기 실시간 수집된 코로나 그래프를 분석하여 코로나 질량 방출 발생을 자동으로 인식하는 코로나 질량 방출 분석 장치 및 그 방법에 관한 것이다.
코로나 질량 방출(Coronal Mass Ejection: CME)은 수십억톤의 태양 물질을 품은 거대한 가스구름이 수백만 Km의 우주 공간으로 뻗어나가는 현상이다.
상기 코로나 질량 방출을 감지하기 위해서 SIDC(Solar Influence Data Center)에서 수행한 CACTus(Computer Aided CME Tracking)가 있다.
이러한 CACTus에서 사용하는 SOHO LASCO 자료는 생성이 주기적이 않아 Hough 변환을 통해 코로나 질량 방출을 찾아내기 위해서는 SOHO LASCO 영상의 결손 부분에 대한 영상의 내삽과정이 필요하고, SOHO LASCO 자료만 이용하므로 3차원 분석이 어렵다.
한국등록특허 제10-1021530호 [명칭: 코로나 물질 방출현상을 이용한 지자기 폭풍 예보 시스템 및 그 방법]
본 발명의 목적은 소호 위성과 스테레오 위성의 코로나 그래프의 실시간 수집 체계를 구축하는 코로나 질량 방출 분석 장치 및 그 방법을 제공하는 데 있다.
본 발명의 다른 목적은 실시간 수집된 코로나 그래프를 분석하여 코로나 질량 방출 발생을 자동으로 인식하고, 산출된 CME 물리량을 근거로 3차원 CME 물리량을 산출하는 코로나 질량 방출 분석 장치 및 그 방법을 제공하는 데 있다.
본 발명의 실시예에 따른 코로나 질량 방출 분석 장치는 소호 위성 및 스테레오 A 위성과 스테레오 B 위성과 통신하는 코로나 질량 방출 분석 장치에 있어서, 상기 소호 위성으로부터 제공되는 소호 라스코 C3 이미지를 수신하고, 상기 스테레오 A 위성 및 상기 스테레오 B 위성으로부터 각각 제공되는 스테레오 COR2 이미지를 수신하는 통신부; 및 상기 수신된 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지에 대한 전처리 과정을 수행하고, 상기 전처리된 이미지를 근거로 코로나 질량 방출을 검출하고, 상기 검출된 코로나 질량 방출의 물리량을 산출하고, 상기 산출된 코로나 질량 방출의 물리량을 근거로 3차원 코로나 질량 방출의 물리량을 산출하는 제어부를 포함할 수 있다.
본 발명과 관련된 일 예로서 상기 제어부는, 상기 수신된 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지에 대한 좌표를 변환하고, 상기 좌표 변환된 복수의 이미지를 근거로 3차원 데이터 큐브를 생성하고, 상기 생성된 3차원 데이터 큐브로부터 배경장을 제거하는 상기 전처리 과정을 수행할 수 있다.
본 발명과 관련된 일 예로서 상기 제어부는, 상기 전처리된 이미지를 근거로 코로나 질량 방출 프론트 후보를 검출하고, 상기 검출된 코로나 질량 방출 프론트 후보에서 별 또는 노이즈에 의해 생긴 후보를 제거하고, 상기 노이즈가 제거된 코로나 질량 방출 프론트 후보들에 대해 계층 군집화(hierarchical clustering) 방법을 적용하여 그룹핑 수행하고, 상기 그룹핑된 각 그룹에 대하여 주어진 시각에 가장 큰 거리를 프로시딩 프론트(proceeding front)로 인식할 수 있다.
본 발명과 관련된 일 예로서 상기 제어부는, 코로나 질량 방출 프론트의 [r, θ, t] 좌표를 산출하고, 상기 산출된 정보를 [x, y, t] 공간에서 그룹핑하고, 상기 각 추출된 개별 그룹에 대하여 2차원 코로나 질량 방출 물리량을 산출할 수 있다.
본 발명과 관련된 일 예로서 상기 코로나 질량 방출의 물리량은, 상기 코로나 질량 방출의 방출 방향, 상기 코로나 질량 방출의 방출 폭, 상기 코로나 질량 방출의 방출 속도 및 상기 코로나 질량 방출의 방출 시각 중 적어도 하나를 포함할 수 있다.
본 발명과 관련된 일 예로서 상기 3차원 코로나 질량 방출의 물리량은, 태양 중심에서의 시각, 속도, 표준 편차, 상기 방출된 코로나 질량 방출의 폭, 방위각 및 천정각 중 적어도 하나를 포함할 수 있다.
본 발명의 실시예에 따른 코로나 질량 방출 분석 방법은 소호 위성 및 스테레오 A 위성과 스테레오 B 위성과 통신하는 코로나 질량 방출 분석 장치의 코로나 질량 방출 분석 방법에 있어서, 통신부를 통해, 상기 소호 위성으로부터 제공되는 소호 라스코 C3 이미지를 수신하고, 상기 스테레오 A 위성 및 상기 스테레오 B 위성으로부터 각각 제공되는 스테레오 COR2 이미지를 수신하는 단계; 제어부를 통해, 상기 수신된 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지에 대한 전처리 과정을 수행하는 단계; 상기 제어부를 통해, 상기 전처리된 이미지를 근거로 코로나 질량 방출을 검출하는 단계; 상기 제어부를 통해, 상기 검출된 코로나 질량 방출의 물리량을 산출하는 단계; 및 상기 제어부를 통해, 상기 산출된 코로나 질량 방출의 물리량을 근거로 3차원 코로나 질량 방출의 물리량을 산출하는 단계를 포함할 수 있다.
본 발명과 관련된 일 예로서 상기 전처리 과정을 수행하는 단계는, 상기 수신된 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지에 대해 직교 좌표계를 극좌표계로 변환하는 과정; 상기 극좌표계로 변환된 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지를 묶어 3차원 데이터 큐브를 생성하는 과정; 및 상기 생성된 3차원 데이터 큐브로부터 배경장을 제거하는 과정;을 포함할 수 있다.
본 발명과 관련된 일 예로서 상기 코로나 질량 방출을 검출하는 단계는, 상기 전처리된 이미지를 근거로 코로나 질량 방출 프론트 후보를 검출하는 과정; 상기 검출된 코로나 질량 방출 프론트 후보에서 별 또는 노이즈에 의해 생긴 후보를 제거하는 과정; 상기 노이즈가 제거된 코로나 질량 방출 프론트 후보들에 대해 계층 군집화 방법을 적용하여 그룹핑 수행하는 과정; 및 상기 그룹핑된 각 그룹에 대하여 주어진 시각에 가장 큰 거리를 프로시딩 프론트로 인식하는 과정;을 포함할 수 있다.
본 발명과 관련된 일 예로서 상기 코로나 질량 방출의 물리량을 산출하는 단계는, 프로시딩 프론트로 인식된 코로나 질량 방출 프론트의 [r, θ, t] 좌표를 산출하는 과정; 상기 산출된 정보를 [x, y, t] 공간에서 그룹핑하는 과정; 및 상기 각 추출된 개별 그룹에 대하여 2차원 코로나 질량 방출 물리량을 산출하는 과정;을 포함할 수 있다.
본 발명은 소호 위성과 스테레오 위성의 코로나 그래프의 실시간 수집 체계를 구축함으로써, 신뢰성 있는 데이터를 제공할 수 있는 효과가 있다.
또한, 본 발명은 실시간 수집된 코로나 그래프를 분석하여 코로나 질량 방출 발생을 자동으로 인식하고, 산출된 CME 물리량을 근거로 3차원 CME 물리량을 산출함으로써, 코로나 질량 방출의 방향, 방출폭, 초기 속도를 자동 인식하고, 3차원 CME 물리량을 제공하여 사용자의 편의성을 향상시키는 효과가 있다.
도 1은 본 발명의 실시예에 따른 코로나 질량 방출 분석 장치의 구성을 나타낸 블록도이다.
도 2는 본 발명의 실시예에 따른 폭 및 문턱값 산출을 의한 예를 나타낸 도이다.
도 3은 본 발명의 실시예에 따른 3차원 CME 물리량을 위한 좌표계를 나타낸 도이다.
도 4는 본 발명의 실시예에 따른 3차원 방출폭 산출을 위한 좌표계 설정을 나타낸 도이다.
도 5는 본 발명의 실시예에 따른 코로나 질량 방출 분석 방법을 나타낸 흐름도이다.
도 6은 본 발명의 실시예에 따른 소호 라스코 C3 직교좌표 이미지를 극좌표계로 변환한 결과를 나타낸 도이다.
도 7은 본 발명의 실시예에 따른 스테레오-A COR2 직교좌표 이미지를 극좌표계로 변환한 결과를 나타낸 도이다.
도 8은 본 발명의 실시예에 따른 스테레오-B COR2 직교좌표 이미지를 극좌표계로 변환한 결과를 나타낸 도이다.
도 9는 본 발명의 실시예에 따른 데이터 큐브의 구조를 나타낸 도이다.
도 10은 본 발명의 실시예에 따른 배경장 추출의 예를 나타낸 도이다.
도 11은 본 발명의 실시예에 따른 배경장 제거 전 및 배경장 제거 후의 비교를 나타낸 도이다.
도 12는 본 발명의 실시예에 따른 CME 프론트 후보의 예를 나타낸 도이다.
도 13은 본 발명의 실시예에 따른 노이즈 제거(또는 데이터 클리닝) 전/후 비교의 예를 나타낸 도이다.
도 14는 본 발명의 실시예에 따른 계층 군집화를 통한 그룹의 구별 예를 나타낸 도이다.
도 15는 본 발명의 실시예에 따른 3차원 공간상에서 구별된 CME 그룹의 예를 나타낸 도이다.
도 16은 본 발명의 실시예에 따른 2차원 CME 물리량 도출 결과를 나타낸 도이다.
도 17 및 도 18은 본 발명의 실시예에 따른 CME 검출 결과를 나타낸 도이다.
본 발명에서 사용되는 기술적 용어는 단지 특정한 실시예를 설명하기 위해 사용된 것으로, 본 발명을 한정하려는 의도가 아님을 유의해야 한다. 또한, 본 발명에서 사용되는 기술적 용어는 본 발명에서 특별히 다른 의미로 정의되지 않는 한, 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자에 의해 일반적으로 이해되는 의미로 해석되어야 하며, 과도하게 포괄적인 의미로 해석되거나, 과도하게 축소된 의미로 해석되지 않아야 한다. 또한, 본 발명에서 사용되는 기술적인 용어가 본 발명의 사상을 정확하게 표현하지 못하는 잘못된 기술적 용어일 때에는 당업자가 올바르게 이해할 수 있는 기술적 용어로 대체되어 이해되어야 할 것이다. 또한, 본 발명에서 사용되는 일반적인 용어는 사전에 정의되어 있는 바에 따라, 또는 전후 문맥상에 따라 해석되어야 하며, 과도하게 축소된 의미로 해석되지 않아야 한다.
또한, 본 발명에서 사용되는 단수의 표현은 문맥상 명백하게 다르게 뜻하지 않는 한 복수의 표현을 포함한다. 본 발명에서 "구성된다" 또는 "포함한다" 등의 용어는 발명에 기재된 여러 구성 요소들 또는 여러 단계를 반드시 모두 포함하는 것으로 해석되지 않아야 하며, 그 중 일부 구성 요소들 또는 일부 단계들은 포함되지 않을 수도 있고, 또는 추가적인 구성 요소 또는 단계들을 더 포함할 수 있는 것으로 해석되어야 한다.
또한, 본 발명에서 사용되는 제 1, 제 2 등과 같이 서수를 포함하는 용어는 구성 요소들을 설명하는데 사용될 수 있지만, 구성 요소들은 용어들에 의해 한정되어서는 안 된다. 용어들은 하나의 구성 요소를 다른 구성 요소로부터 구별하는 목적으로만 사용된다. 예를 들어, 본 발명의 권리 범위를 벗어나지 않으면서 제 1 구성 요소는 제 2 구성 요소로 명명될 수 있고, 유사하게 제 2 구성 요소도 제 1 구성 요소로 명명될 수 있다.
이하, 첨부된 도면을 참조하여 본 발명에 따른 바람직한 실시예를 상세히 설명하되, 도면 부호에 관계없이 동일하거나 유사한 구성 요소는 동일한 참조 번호를 부여하고 이에 대한 중복되는 설명은 생략하기로 한다.
또한, 본 발명을 설명함에 있어서 관련된 공지 기술에 대한 구체적인 설명이 본 발명의 요지를 흐릴 수 있다고 판단되는 경우 그 상세한 설명을 생략한다. 또한, 첨부된 도면은 본 발명의 사상을 쉽게 이해할 수 있도록 하기 위한 것일 뿐, 첨부된 도면에 의해 본 발명의 사상이 제한되는 것으로 해석되어서는 아니 됨을 유의해야 한다.
도 1은 본 발명의 실시예에 따른 코로나 질량 방출 분석 장치(10)의 구성을 나타낸 블록도이다.
도 1에 도시한 바와 같이, 코로나 질량 방출 분석 장치(10)는 통신부(100), 저장부(200), 표시부(300) 및 제어부(400)로 구성된다. 도 1에 도시된 코로나 질량 방출 분석 장치(10)의 구성 요소 모두가 필수 구성 요소인 것은 아니며, 도 1에 도시된 구성 요소보다 많은 구성 요소에 의해 코로나 질량 방출 분석 장치(10)가 구현될 수도 있고, 그보다 적은 구성 요소에 의해서도 코로나 질량 방출 분석 장치(100)가 구현될 수도 있다.
상기 통신부(100)는 유/무선 통신망을 통해 내부의 임의의 구성 요소 또는 외부의 임의의 적어도 하나의 단말기와 통신 연결한다. 이때, 상기 외부의 임의의 단말기는 위성(예를 들어 소호 위성, 스테레오 위성 등 포함), 서버 등을 포함할 수 있다. 여기서, 무선 인터넷 기술로는 무선랜(Wireless LAN: WLAN), 와이브로(Wireless Broadband: Wibro), 와이맥스(World Interoperability for Microwave Access: Wimax), HSDPA(High Speed Downlink Packet Access), IEEE 802.16, 롱 텀 에볼루션(Long Term Evolution: LTE), 광대역 무선 이동 통신 서비스(Wireless Mobile Broadband Service: WMBS) 등이 포함될 수 있다. 또한, 근거리 통신 기술로는, 블루투스(Bluetooth), 와이 파이(Wi-Fi), RFID(Radio Frequency Identification), 적외선 통신(Infrared Data Association: IrDA), UWB(Ultra Wideband), 지그비(ZigBee), 인접 자장 통신(Near Field Communication: NFC), 초음파 통신(Ultra Sound Communication: USC), 가시광 통신(Visible Light Communication: VLC), BLE(Bluetooth Low Energy) 등이 포함될 수 있다. 또한, 유선 통신 기술로는, 전력선 통신(Power Line Communication: PLC), USB 통신, 이더넷(Ethernet), 시리얼 통신(serial communication), 광/동축 케이블 등이 포함될 수 있다.
또한, 상기 통신부(100)는 유니버설 시리얼 버스(Universal Serial Bus: USB)를 통해 상기 단말기와 정보를 상호 전송할 수 있다.
또한, 상기 통신부(100)는 소호 위성, 스테레오 위성 등을 포함하는 상기 위성으로부터 제공되는 이미지(또는 코로나 그래프 이미지)를 수신(또는 수집)한다. 여기서, 상기 이미지(또는 영상)는 소호 위성(또는 소호 라스코 위성)에서 생성되는 소호 라스코 C3(또는 SOHO LASCO C3) 이미지, 스테레오 A 위성 및 스테레오 B 위성에서 생성되는 스테레오 COR2(STEREO COR2) 이미지 등을 포함한다.
상기 저장부(200)는 다양한 사용자 인터페이스(User Interface: UI), 그래픽 사용자 인터페이스(Graphic User Interface: GUI) 등을 저장한다.
또한, 상기 저장부(200)는 상기 코로나 질량 방출 분석 장치(10)가 동작하는데 필요한 데이터와 프로그램 등을 저장한다.
또한, 상기 저장부(200)는 플래시 메모리 타입(Flash Memory Type), 하드 디스크 타입(Hard Disk Type), 멀티미디어 카드 마이크로 타입(Multimedia Card Micro Type), 카드 타입의 메모리(예를 들면, SD 또는 XD 메모리 등), 자기 메모리, 자기 디스크, 광디스크, 램(Random Access Memory: RAM), SRAM(Static Random Access Memory), 롬(Read-Only Memory: ROM), EEPROM(Electrically Erasable Programmable Read-Only Memory), PROM(Programmable Read-Only Memory) 중 적어도 하나의 저장매체를 포함할 수 있다. 또한, 코로나 질량 방출 분석 장치(10)는 인터넷(internet)상에서 저장부(200)의 저장 기능을 수행하는 웹 스토리지(web storage)를 운영하거나, 또는 상기 웹 스토리지와 관련되어 동작할 수도 있다.
또한, 상기 저장부(200)는 상기 제어부(400)의 제어에 의해 상기 통신부(100)를 통해 수신된 이미지(또는 코로나 그래프 이미지)를 저장한다.
상기 표시부(300)는 상기 제어부(400)의 제어에 의해 상기 저장부(200)에 저장된 사용자 인터페이스 및/또는 그래픽 사용자 인터페이스를 이용하여 다양한 메뉴 화면 등과 같은 다양한 콘텐츠를 표시할 수 있다. 여기서, 상기 표시부(300)에 표시되는 콘텐츠는 다양한 텍스트 또는 이미지 데이터(각종 정보 데이터 포함)와 아이콘, 리스트 메뉴, 콤보 박스 등의 데이터를 포함하는 메뉴 화면 등을 포함한다. 또한, 상기 표시부(300)는 터치 스크린 일 수 있다.
또한, 상기 표시부(300)는 액정 디스플레이(Liquid Crystal Display: LCD), 박막 트랜지스터 액정 디스플레이(Thin Film Transistor-Liquid Crystal Display: TFT LCD), 유기 발광 다이오드(Organic Light-Emitting Diode: OLED), 플렉시블 디스플레이(Flexible Display), 3차원 디스플레이(3D Display), 전자잉크 디스플레이(e-ink display), LED(Light Emitting Diode) 중에서 적어도 하나를 포함할 수 있다.
또한, 상기 표시부(300)는 상기 제어부(400)의 제어에 의해 상기 통신부(100)를 통해 수신된 이미지(또는 코로나 그래프 이미지)를 표시한다.
상기 제어부(400)는 상기 코로나 질량 방출 분석 장치(10)의 전반적인 제어 기능을 실행한다.
또한, 상기 제어부(400)는 상기 통신부(100)를 통해 상기 소호 위성으로부터 제공되는 소호 라스코 C3(또는 SOHO LASCO C3) 이미지, 상기 스테레오 A 위성 및 스테레오 B 위성으로부터 제공되는 스테레오 COR2(STEREO COR2) 이미지 등을 수집(또는 아카이브(archive))한다. 여기서, 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지는 FITS 파일 형식, JPEG 파일 형식 등으로 제공되며, 2048×2048, 1024×1024, 512×512, 256×2356, 128×128 해상도 등의 형태로 제공된다.
본 발명의 실시예에서는 JPEG 형태의 이미지를 근거로 설명하고 있으나, 이에 한정되는 것은 아니며 FITS 형태의 이미지도 본 발명의 기술적 특징에 적용할 수 있다.
즉, 상기 제어부(400)는 상기 위성 또는 상기 서버로부터 현재 날짜 및 이전 날짜의 파일 리스트를 상기 통신부(100)를 통해 다운로드한다.
또한, 상기 제어부(400)는 상기 다운로드된 현재 날짜 및 이전 날짜의 파일 리스트에 대해서 HTML 코드 파싱을 수행하여, 이전 날짜로부터 현재 날짜까지의 이미지 리스트를 추출(또는 확인)한다.
또한, 상기 제어부(400)는 상기 추출된 이미지 리스트로부터 가장 최근에 다운로드 받은 이미지 이후의 이미지를 상기 위성 또는 상기 서버로부터 상기 통신부(100)를 통해 다운로드한다.
이때, 상기 소호 라스코 C3 이미지의 경우, 영상 처리 과정에서 생기는 오류로 인하여 실제 영상이 제대로 나오지 않는 경우가 존재한다. 또한, 해당 영상들은 영상 처리 과정에서 오류를 일으키게 되므로 미리 제거하는 과정이 필요하다.
따라서, 상기 제어부(400)는 상기 다운로드한 이미지의 용량을 확인(예를 들어, 해당 이미지에 오류가 있는 경우, 정상적인 용량에 비해서 매우 작은 용량을 지니는 특징이 있음)하고, 상기 확인된 이미지의 용량이 미리 설정된 값(예를 들어 120kB) 이하인 경우, 오류가 있는 이미지로 분류하여 해당 이미지를 제외(또는 삭제)한다.
또한, 상기 제어부(400)는 상기 다운로드한 이미지 파일 정보를 업데이트하고, 상기 업데이트된 이미지 파일 정보를 상기 저장부(200)에 저장한다.
이와 같이, 상기 제어부(400)는 미리 설정된 시간 간격(예를 들어 10분)으로 상기 위성 또는 상기 서버로부터 상기 소호 라스코 C3 이미지, 상기 스테레오 COR2 이미지 등을 수집한다.
또한, 상기 제어부(400)는 상기 수신된(또는 수집된) 이미지에 대한 전처리 과정을 수행한다. 여기서, 상기 전처리 과정은 상기 이미지에 대한 좌표 변환 과정, 상기 좌표 변환된 이미지를 근거로 3차원 데이터 큐브를 생성하는 과정, 상기 생성된 3차원 데이터 큐브로부터 배경장을 제거하는 과정 등을 포함한다.
즉, 상기 제어부(400)는 [x, y] 직교 좌표계의 이미지를 [r, θ]의 극좌표계로 변환하기 위해서 각 이미지에 대해서 이미지의 중심 좌표를 생성한다. 이때, 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지에는 태양의 디스크가 표시(또는 포함)되어 있다. 따라서, 상기 제어부(400)는 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지에 표시된(또는 포함된) 태양의 디스크를 근거로 태양 디스크의 영상 좌표를 추출한다.
또한, 상기 제어부(400)는 상기 추출된 영상 좌표를 근거로 원형 회귀를 통해 각 이미지별로 태양 중심 좌표 및 각 이미지에서의 태양 디스크 반지름(r_sun)을 추출한다.
또한, 상기 제어부(400)는 상기 추출된 이미지 중심 좌표(xc, yc)를 근거로 [x, y] 직교 좌표계로부터 [r, θ] 극좌표계로 변환한다. 이때, 상기 좌표 변환에 사용되는 이미지의 크기가 512×512인 경우, 상기 제어부(400)는 이미지를 충분히 커버할 수 있도록 반지름의 범위(r_sun-범위)는 0에서 255 픽셀의 범위를 갖도록 설정하고, 방위각(θ, 또는 θ_2d)은 1도 단위로 0도에서 359도의 범위를 갖도록 설정한다.
또한, 상기 제어부(400)는 상기 반지름(r_sun) 및 상기 방위각에 대응하는 직교 좌표계를 아래의 [수학식 1]과 같이 나타낼 수 있다.
Figure 112014073165270-pat00001
따라서, 상기 제어부(400)에 의해 극좌표계로 변환된 영상은 512×512에서 250×360 크기의 이미지로 변환된다.
여기서, 상기 [수학식 1]에서는 코로나 질량 방출(Coronal Mass Ejection: CME)의 방향을 나타내는 방출 방향(Position Angle: PA, 북쪽 방향으로부터 반시계 방향의 각도)을 사용하지 않고, 상기 방위각(θ)은 x축으로부터 반시계 방향의 각을 나타낸다. 따라서, 상기 방출 방향은 PA=θ+90°으로부터 산출(또는 변환)할 수 있다.
또한, 상기 제어부(400)는 상기 극좌표계로 변환된 이미지들(예를 들어 상기 극좌표계로 변환된 소호 라스코 C3 이미지, 상기 극좌표계로 변환된 스테레오 COR2 이미지 등 포함)을 묶어 새로운 3차원 데이터 큐브를 생성한다.
즉, 상기 제어부(400)는 상기 극좌표계로 변환된 복수의 이미지를 미리 설정된 특정 각도를 기준으로 펼치고, 상기 펼쳐진 복수의 이미지를 시간별로 정리하여(또는 묶어서) 상기 새로운 3차원 데이터 큐브를 생성한다. 이때, 상기 3차원 데이터 큐브의 시간축의 범위는 미리 설정된 시각인 현재 시각에서 24시간 이전의 자료로부터 현재 시각까지의 이미지에 해당한다. 또한, 결측이나 이미지의 오류 등으로 데이터 사이의 시간 간격이 일정하지 않아 시간축 방향 데이터 큐브의 크기가 매번 달라지므로, 상기 제어부(400)는 데이터 큐브의 시간 정보를 상기 저장부(200)에 별도로 기록하여 보관한다.
또한, 상기 제어부(400)는 상기 생성된 3차원 데이터 큐브로부터 상기 배경장을 제거한다. 여기서, 상기 배경장은 코로나 질량 방출과 달리 시간에 따라 크게 변하지 않는 부분을 의미한다. 이러한 상기 배경장은 코로나 질량 방출 검출에 방해가 되므로, 최대한 제거해야 한다.
즉, 상기 제어부(400)는 상기 생성된 3차원 데이터 큐브에서 [r, θ]에 대해 시계열 이미지 값을 추출한다.
또한, 상기 제어부(400)는 아래의 [수학식 2]와 같은 다항식 회귀를 근거로 상기 추출된 시계열 이미지 값에서 상기 배경장을 추정한다. 이때, 다항식은 2차 및 3차를 적용한다.
Figure 112014073165270-pat00002
여기서, y는 각 시각에 따른 이미지 화소 값을 나타내고, A0, A1, A2, A3는 각각 상기 다항식을 피팅할 때 사용하는 계수로, 설계자의 설계에 따라 피팅 시 설정할 수 있다.
본 발명의 실시예에서는 상기 다항식에 대해서 최대 3차로 설명하고 있으나, 이에 한정되는 것은 아니며, 설계자의 설계에 따라 상기 다항식의 차수를 다양하게 설계할 수 있다.
이와 같이, 차수가 높을수록 배경장의 시간에 따른 변화에 민감하게 반응하나, 너무 높은 차수를 적용하면, 규모가 큰 코로나 질량 방출의 경우 코로나 질량 방출 자체를 제거할 가능성이 있으므로, 차수는 최대 3차가 적당할 수 있다. 또한, 본 발명의 실시예에서는 2차가 가장 적절하다. 또한, 코로나 질량 방출에 해당하는 정보 제거 방지를 위해 1회 차로 구해진 배경장에 비해 편차가 큰 값을 제외한 후, 다시 다항식을 적용하여 배경장을 산출할 수 있다.
또한, 상기 제어부(400)는 상기 생성된 3차원 데이터 큐브로부터 상기 추정된 배경장을 제거한다.
또한, 상기 제어부(400)는 상기 전처리된 이미지(또는 상기 배경장이 제거된 3차원 데이터 큐브)를 근거로 코로나 질량 방출을 검출한다. 여기서, 상기 코로나 질량 방출 검출은 코로나 질량 방출 프론트(front) 후보 검출 과정, 상기 검출된 코로나 질량 방출 프론트 후보에서 별 또는 노이즈에 의해 생긴 후보를 제거하는 과정, 상기 노이즈가 제거된 코로나 질량 방출 프론트 후보들에 대하여 시공간상의 클러스터링을 통해 그룹핑을 한 이후에 가장 앞서 진행하는 프로시딩 프론트(proceeding front)를 검출하여 상기 코로나 질량 방출을 검출하는 과정 등을 포함한다.
즉, 상기 제어부(400)는 상기 전처리된 이미지(또는 상기 배경장이 제거된 3차원 데이터 큐브)의 [θ, t]에서 이미지 r-방향으로의 편차를 정의하는 에지 검출 방법을 사용하여 상기 코로나 질량 방출 프론트 후보를 검출한다.
즉, 상기 제어부(400)는 상기 코로나 질량 방출 프론트의 후보를 검출하기 위해서 아래의 [수학식 3]와 같이 이미지 r-방향의 편차를 정의한다. 상기 [수학식 3]은 이미지 프로세싱 중 커널(kernerl) 프로세싱과 유사하다.
Figure 112014073165270-pat00003
여기서,
Figure 112014073165270-pat00004
는 r-방향 i번째의 화소에서의 편차값이고, i는 r-방향에서의 화소 위치이고, k는 r-방향에서의 화소 위치이고, vk는 k번째 픽셀의 화소값이고, 폭(width)은 r-방향의 편차 계산 시 사용하는 폭(또는 화소의 개수)이다.
상기 제어부(400)는 상기 정의한 이미지 편차가 미리 설정된 문턱값(threshold) 이상의 값을 가지는 경우, 코로나 질량 방출 프론트의 후보로 선정한다. 여기서, 상기 제어부(400)는 상기 코로나 질량 방출 프론트의 경우, 미리 설정된 기준 픽셀에 대해서 태양에 가까운 쪽의 값의 평균값과 태양에서 먼 쪽의 평균값을 비교하며, 상기 비교 결과(또는 상기 이미지 편차값)가 양의 값인 경우에는 코로나 질량 방출의 바깥쪽에 해당하고, 상기 비교 결과가 음의 값인 경우에는 코로나 질량 방출의 안쪽에 해당한다.
이때, 상기 제어부(400)는 코로나 질량 방출을 잘 결정할 수 있도록 폭(width) 및 문턱값을 설정한다. 즉, 폭을 작게 잡을수록 세밀한 형태에 민감하게 반응하나, 너무 낮게 잡을 경우 노이즈가 많아진다. 또한, 문턱값이 낮을수록 희미한 형태를 많이 잡아낼 수 있으나, 너무 낮을 경우 역시 노이즈가 많아진다. 이때, 문턱값의 경우는 양의 값(또는 코로나 질량 방출의 바깥쪽 부분)을 취하도록 설정한다.
예를 들어, 도 2에 도시된 바와 같이, 상기 제어부(400)는 폭 및 문턱값을 산출한다. 여기서, 상기 도 2의 좌측에 도시된 그림(210)은 주어진 시각 및 방향각에 대한 r-방향의 이미지 값을 나타낸다. 또한, 상기 도 2에 표시된 화살표(211)는 수동으로 사람이 찾아낸 코로나 질량 방출 프론트에 해당한다. 또한, 상기 도 2의 우측에 도시된 그림(220)은 주어진 시각 및 방향각에 대한 이미지의 r-방향 편차를 나타낸다. 상기 도 2의 우측에 도시된 그림(220)에서의 황색선, 녹색선 및 청색선은 각각 폭이 1, 2, 3으로 주어졌을 때의 이미지 편차 값에 해당한다.
또한, 상기 도 2에 도시된 바와 같이, 상기 이미지(또는 상기 코로나 그래프 이미지)에서 반경 방향의 크기가 작은 형태에서는 지정한 폭이 작을수록 r-방향의 편차가 크게 나타나지만 큰 차이는 없다. 따라서, 본 발명의 실시예에 따라, 설계자는 폭이 작은 경우 별 또는 우주선(cosmic ray)의 영향을 받을 수 있으므로, 상기 별 또는 우주선의 영향을 고려하여 폭을 결정할 수 있다.
또한, 상기 제어부(400)는 상기 검출된 코로나 질량 방출 프론트 후보에 대하여 노이즈에 해당하는 점을 제거하는 데이터 클리닝(data cleaning) 단계를 수행한다.
즉, 상기 제어부(400)는 상기 검출된 코로나 질량 방출 프론트 후보 중에서 미리 설정된 노이즈 판별 기준에 해당하는 상기 코로나 질량 방출 프론트 후보를 제거한다. 여기서, 상기 미리 설정된 노이즈 판별 기준은 미리 설정된 r-방향 10 픽셀, 미리 설정된 θ-방향 10도의 범위 내에서 검출되는 코로나 질량 방출 프론트의 개수가 미리 설정된 3개 이하인 경우 노이즈로 구별한다.
또한, 상기 제어부(400)는 상기 노이즈가 제거된 코로나 질량 방출 프론트 후보들에 대하여 시공간상의 클러스터링인 계층 군집화(hierarchical clustering) 방법을 통해 그룹핑을 하여, 각 그룹을 구별한다. 이때, 상기 클러스터 사이의 거리는 최근이웃법(nearest neighbor)을 적용한다. 여기서, 상기 최근이웃법은 각 클러스터의 소속 멤버 중에서 가장 거리가 가까운 것을 두 클러스터 사이의 거리로 정의하는 것이다. 또한, 상기 제어부(400)는 만약 두 클러스터 간의 거리가 주어진 기준보다 작으면, 두 클러스터는 하나의 클러스터로 통합한다.
또한, 상기 계층 군집화는 [r, t] 평면에서 수행하며, r과 t의 차원이 다르므로 단위 차이를 고려하여 아래의 [수학식 4]에 의해 2개의 코로나 질량 방출 프론트 후보점 사이의 거리를 정의한다. 여기서, 가중치(ω)(예를 들어 50)를 설정함으로써, 클러스터의 통합이 실제 거리(r)에 중점을 둘 것인지 시간(t)에 중점을 둘 것인지를 결정할 수 있다.
Figure 112014073165270-pat00005
여기서, d는 2개의 코로나 질량 방출 프론트 후보점 사이의 거리를 나타내고, r1 및 r2는 상기 코로나 질량 방출 프론트의 실제 거리를 각각 나타내고, t1 및 t2는 상기 코로나 질량 방출 프론트의 생성 시각(예를 들어 hour 단위)을 각각 나타낸다.
또한, 상기 제어부(400)는 각 클러스터 그룹 사이의 거리가 미리 설정된 기준값(예를 들어 30) 이상인 경우, 해당 클러스터 그룹들을 서로 다른 그룹으로 구별한다.
또한, 상기 제어부(400)는 각 그룹에 대하여 주어진 시각에 가장 큰 거리를 프로시딩 프론트로 인식(또는 검출)한다.
또한, 상기 제어부(400)는 상기 검출된 코로나 질량 방출(또는 상기 프로시딩 프론트로 인식된 상기 코로나 질량 방출 프론트)의 물리량(또는 2차원 코로나 질량 방출의 물리량)을 산출(또는 도출)한다. 여기서, 상기 코로나 질량 방출(또는 상기 프로시딩 프론트로 인식된 상기 코로나 질량 방출 프론트)의 물리량 산출은 상기 코로나 질량 방출 프론트의 [r, θ, t] 좌표를 산출(또는 추출)하고, 상기 산출된 정보를 [x, y, t] 공간에서 그룹핑하는 과정, 상기 각 추출된 개별 그룹에 대하여 2차원 코로나 질량 방출 물리량을 산출하는 과정 등을 포함한다. 이때, 상기 산출되는 2차원 코로나 질량 방출의 물리량은 코로나 질량 방출의 방출 방향(또는 방위각), 방출 폭, 방출 속도, 방출 시각 등을 포함하며, 시선 방향에 따른 투영된 값을 의미한다. 또한, 상기 방출 시각(또는 발생 시각)은 코로나 질량 방출이 일정한 속도로 방출한다는 가정에서 코로나 질량 방출이 태양 중심에 위치할 때의 시각을 의미한다.
또한, 상기 제어부(400)는 상기 소호 라스코 C3 이미지 및 2개의 상기 스테레오 COR2 이미지에 대해서 각각 수행하고, 각각의 이미지에 대한 결과를 상기 저장부(200)에 저장한다.
즉, 상기 제어부(400)는 최종적으로 검출된 상기 코로나 질량 방출 프론트의 [r, θ, t] 좌표를 각 그룹별로 추출한다. 또한, 상기 제어부(400)는 상기 추출된 [r, θ, t] 좌표를 다시 직교 좌표 공간 및 시간 공간에 해당하는 [x, y, t] 좌표로 변환한다. 또한, 상기 제어부(400)는 상기 변환된 [x, y, t] 3차원 공간에서 다시 계층 군집화를 수행한다. 여기서, 상기 계층 군집화에서 그룹간 거리는 최근접이웃 방법을 적용한다.
이와 같이, 상기 제어부(400)는 상기 최종적으로 검출된 상기 코로나 질량 방출 프론트에 대해서 시공간을 통합한 코로나 질량 방출 그룹을 생성한다.
또한, 상기 제어부(400)는 상기 생성된 코로나 질량 방출 그룹을 근거로 2차원 코로나 질량 방출 물리량을 산출(또는 도출)한다. 여기서, 상기 2차원 코로나 질량 방출 물리량은 코로나 질량 방출의 방출 방향(Position Angle: PA), 코로나 질량 방출의 방출 폭(width), 코로나 질량 방출의 방출 속도, 코로나 질량 방출의 방출 시각 등을 포함한다. 이때, 상기 코로나 질량 방출의 방출 방향의 기준은 북쪽 방향을 기준으로 반시계 방향 각도를 의미한다. 또한, 상기 코로나 질량 방출의 방출 폭은 전체 방출 폭(full width)을 의미한다. 또한, 상기 코로나 질량 방출의 방출 속도(단위 km/s)는 시간에 따른 상기 코로나 질량 방출 프론트의 위치에 대해 선형회귀 방법을 적용한 속도를 의미한다. 또한, 상기 코로나 질량 방출의 방출 시각은 상기 코로나 질량 방출의 방출 속도에 의해 상기 코로나 질량 방출이 태양 중심에 있을 때의 시각을 의미한다.
또한, 상기 2차원 코로나 질량 방출 물리량은 상기 코로나 질량 방출의 실제 방출 방향을 고려하지 않고, 투영된 이미지상의 속도, 방향, 방출폭을 의미한다. 또한, 상기 제어부(400)는 상기 검출된 각 코로나 질량 방출 그룹별로 방출 방향에 따른 방출 속도, 방출 시각을 계산하고, 고사분위(upper quartile) 값을 코로나 질량 방출 그룹에 대한 방출 속도로 정의한다.
또한, 실제 각 방위각 별로 계산된 방출 속도 및 방출 시각은 코로나 질량 방출 이미지의 상태에 의해 에러가 많은 값을 나타내므로 튀는 값(또는 노이즈)을 제거하기 위해서 고사분위 값을 선택한다.
또한, 상기 제어부(400)는 상기 산출된 2차원 코로나 질량 방출의 물리량을 저장부(200)에 저장하거나 상기 표시부(300)에 표시한다. 이때, 상기 저장부(200)에 저장되는 2차원 코로나 질량 방출의 물리량은 상기 제어부(400)의 제어에 의해, 상기 소호 라스코 C3 이미지인 경우 [YYYY][MM][DD]_[hh]_la_cme.txt의 형태로 상기 저장부(200)에 저장되고, 상기 스테레오 COR2 이미지인 경우 [YYYY][MM][DD]_[hh]_sa_cme.txt의 형태로 상기 저장부(200)에 저장될 수 있다.
또한, 상기 제어부(400)는 상기 산출된 2차원 코로나 질량 방출의 물리량을 근거로 3차원 코로나 질량 방출의 물리량을 산출(또는 도출)한다. 여기서, 상기 3차원 코로나 질량 방출의 물리량은 태양 중심에서의 시각, 속도(또는 상기 산출된 3차원 공간상의 속도), 표준 편차(또는 상기 산출된3차원 공간상의 표준 편차), 상기 방출된 코로나 질량 방출의 폭, 방위각, 천정각 등을 포함한다.
또한, 상기 제어부(400)는 상기 산출된 3차원 코로나 질량 방출의 물리량을 상기 저장부(200)에 저장하거나 상기 표시부(300)에 표시한다.
이때, 도 3에 도시한 바와 같이, 상기 3차원 코로나 질량 방출의 물리량을 산출하기 위해서 설정한 좌표계는 [x, y, z] 직교 좌표계에서 고려할 때, x-축은 지구 방향, y-축은 태양의 서쪽 방향, z-축은 태양의 북극 방향으로 정의한다. 다만, 실제 소호 라스코 코로나 그래프에서 보는 [x', y'] 좌표는 상기 정의된 3차원 좌표계에서는 [y,z] 평면에 투영된 좌표계에 해당한다.
따라서, 상기 코로나 질량 방출이 방위각(φ) 및 천정각(θ, 또는 θ_3d)으로 방출되는 경우, 상기 제어부(400)는 상기 소호 라스코의 코로나 그래프에서 투영된 코로나 질량 방출 프론트의 중심을 아래의 [수학식 5]와 같이 정의할 수 있다.
Figure 112014073165270-pat00006
여기서, 상기 소호 라스코 코로나 그래프에서 상기 코로나 질량 방출 프론트의 투영된 태양 중심으로부터의 거리(r')는 아래의 [수학식 6]과 같이 나타낼 수 있다.
Figure 112014073165270-pat00007
또한, 상기 제어부(400)는 상기 투영된 코로나 그래프에서 상기 2차원 코로나 질량 방출의 물리량 산출에서 산출한 투영된 속도(V') 및 상기 3차원 코로나 질량 방출의 물리량 산출에서 산출한 방출 속도(V)를 근거로 실제 3차원 공간상의 속도는 아래의 [수학식 7]과 같이 나타낼 수 있다.
Figure 112014073165270-pat00008
또한, 도 4에 도시된 바와 같이, 상기 제어부(400)는 상기 투영된 방출폭으로부터 실제 3차원 상의 방출폭을 계산하기 위해서 상기 투영된 평면과 실제 방출된 원뿔 형태의 코로나 질량 방출 빔에 대한 근사를 하여, 상기 투영된 평면과 상기 방출된 원뿔 사이의 각을 라고 ω할 때, 투영되어 관측되는 방출폭(α')과 실제 방출폭(α) 사이의 관계를 아래의 [수학식 8]과 같이 나타낸다.
Figure 112014073165270-pat00009
이때, 상기 cosω는 아래의 [수학식 9]와 같이 나타낸다.
Figure 112014073165270-pat00010
따라서, 실제 3차원 공간상에서 방출된 코로나 질량 방출의 방출폭은 아래의 [수학식 10]과 같이 나타낸다.
Figure 112014073165270-pat00011
따라서, 태양에서 코로나 질량 방출이 방위각(φ) 및 천정각(θ) 방향으로 방출되었을 경우, 소호 라스코 위성에서 관측되는 경우, 방출 속도는 아래의 [수학식 11]과 같이 나타낼 수 있다.
Figure 112014073165270-pat00012
또한, 스테레오-A 위성의 위치가
Figure 112014073165270-pat00013
로 주어지는 경우, 상기 스테레오-A 위성에서 관측되는 방출 속도는 아래의 [수학식 12]와 같이 나타낼 수 있다.
Figure 112014073165270-pat00014
또한, 스테레오-B 위성의 위치가
Figure 112014073165270-pat00015
로 주어지는 경우, 상기 스테레오-B 위성에서 관측되는 방출 속도는 아래의 [수학식 13]과 같이 나타낼 수 있다.
Figure 112014073165270-pat00016
또한, 상기 방출된 코로나 질량 방출이 세 위성에서 동시에 관측되는 경우, 실제 방출된 코로나 질량 방출의 속도는 같아야 하므로, 방위각(φ) 및 천정각(θ)이 설정되는 경우,
Figure 112014073165270-pat00017
인 관계가 성립해야 한다.
따라서, 상기 제어부(400)는 상기 방위각(φ) 및 천정각(θ)에 대해서 각 위성에서 관측된 코로나 질량 방출의 투영된 방출 속도로부터 실제 속도를 산출한다.
또한, 상기 제어부(400)는 상기 3개의 위성에서 산출한 실제 3차원 공간상에서의 코로나 질량 방출 속도(
Figure 112014073165270-pat00018
)의 표준편차가 가장 작은 점을 3차원 공간상에서의 코로나 질량 방출의 방출 속도로 결정한다.
또한, 코로나 질량 방출이 2개의 위성에서만 관측되는 경우, 상기 제어부(400)는 2개의 위성에서 산출된 방출 속도의 표준편차가 가장 작은 값을 3차원 공간상에서의 코로나 질량 방출의 방출 속도로 결정한다.
또한, 실제 방출된 방향 방위각(φ) 및 천정각(θ)은 하나지만, 수학적으로 위상이 반대인 경우에 의해서 실제 최소값을 가지는 방위각(φ) 및 천정각(θ)의 조합은 4개가 존재한다.
또한, 상기 제어부(400)는 상기 4개의 경우 중에서 상기 관측된 코로나 질량 방출의 방출 방향으로 실제 방출 방향을 결정한다.
또한, 상기 제어부(400)는 상기 코로나 질량 방출의 방출 방향과 실제 방출 속도가 정해지면(또는 결정되면/산출되면), 투영 효과에 의한 오차를 최소화하기 위해서 상기 실제 3차원 공간상의 방출 속도와 상기 투영된 방출 속도의 비가 가장 작은 위성 관측값을 기준으로 상기 코로나 질량 방출의 방출폭을 결정(또는 산출)한다.
또한, 상기 제어부(400)는 상기 산출된 3차원 코로나 질량 방출의 물리량을 저장부(200)에 저장하거나 상기 표시부(300)에 표시한다. 이때, 상기 저장부(200)에 저장되는 3차원 코로나 질량 방출의 물리량은 상기 제어부(400)의 제어에 의해, cme_3d_[YYYY][MM][DD]_[hh].dat의 아스키 텍스트 파일 형태로 저장될 수 있다.
이와 같이, 소호 위성과 스테레오 위성의 코로나 그래프의 실시간 수집 체계를 구축할 수 있다.할 수 있다.
또한, 이와 같이, 실시간 수집된 코로나 그래프를 분석하여 코로나 질량 방출 발생을 자동으로 인식하고, 산출된 CME 물리량을 근거로 3차원 CME 물리량을 산출할 수 있다.
이하에서는, 본 발명에 따른 코로나 질량 방출 분석 방법을 도 1 내지 도 18을 참조하여 상세히 설명한다.
도 5는 본 발명의 실시예에 따른 코로나 질량 방출 분석 방법을 나타낸 흐름도이다.
먼저, 통신부(100)는 위성(예를 들어 소호 위성, 스테레오 위성 등 포함)(미도시) 또는 서버(미도시)로부터 제공되는 이미지(또는 코로나 그래프 이미지)를 수신(또는 수집)한다. 여기서, 상기 이미지(또는 영상)는 소호 위성에서 생성되는 소호 라스코 C3(또는 SOHO LASCO C3) 이미지, 스테레오 A 위성 및 스테레오 B 위성에서 생성되는 스테레오 COR2(STEREO COR2) 이미지 등을 포함한다(S510).
이후, 제어부(400)는 수신된(또는 수집된) 이미지(예를 들어 상기 소호 라스코 C3 이미지, 상기 스테레오 COR2 이미지 등 포함)에 대한 전처리 과정을 수행한다. 여기서, 상기 전처리 과정은 상기 이미지에 대한 좌표 변환 과정, 상기 좌표 변환된 이미지를 근거로 3차원 데이터 큐브를 생성하는 과정, 상기 생성된 3차원 데이터 큐브로부터 배경장을 제거하는 과정 등을 포함한다.
즉, 상기 제어부(400)는 [x, y] 직교 좌표계의 이미지를 [r, θ]의 극좌표계로 변환하기 위해서 각 이미지에 대해서 이미지의 중심 좌표를 생성한다.
이때, 상기 제어부(400)는 각 이미지에 대하여 이미지의 중심 좌표를 생성한다. 즉, 소호 라스코 C3이미지 및 스테레오 COR2 이미지에 표시된(또는 포함된) 태양의 디스크를 근거로 태양 디스크의 영상 좌표를 추출한다. 또한, 상기 제어부(400)는 상기 추출된 영상 좌표를 근거로 원형 회귀를 통해 각 이미지별로 태양 중심 좌표 및 각 이미지에서의 태양 디스크 반지름을 추출한다. 또한, 상기 제어부(400)는 상기 추출된 이미지 중심 좌표(xc, yc)를 근거로 [x, y] 직교 좌표계로부터 [r, θ] 극좌표계로 변환한다.
일 예로, 도 6에 도시된 바와 같이, 상기 제어부(400)는 상기 소호 위성에서 생성되는 소호 라스코 C3 이미지인 직교 좌표 이미지(610)를 극좌표계(620)로 변환한다.
다른 일 예로, 도 7 및 도 8에 도시된 바와 같이, 상기 제어부(400)는 상기 스테레오 A 위성 및 스테레오 B 위성에서 생성되는 스테레오 COR2 이미지인 직교 좌표 이미지(710, 810)를 각각 극좌표계(720, 820)로 변환한다.
또한, 상기 제어부(400)는 상기 극좌표계로 변환된 이미지들(예를 들어 상기 극좌표계로 변환된 소호 라스코 C3 이미지, 상기 극좌표계로 변환된 스테레오 COR2 이미지 등 포함)을 묶어 새로운 3차원 데이터 큐브를 생성한다.
즉, 상기 제어부(400)는 상기 극좌표계로 변환된 복수의 이미지를 미리 설정된 특정 각도를 기준으로 펼치고, 상기 펼쳐진 복수의 이미지를 시간별로 정리하여(또는 묶어서) 상기 새로운 3차원 데이터 큐브를 생성한다.
일 예로, 도 9에 도시된 바와 같이, 상기 제어부(400)는 상기 극좌표계로 변환된 소호 라스코 C3 이미지 및 복수의 스테레오 COR2 이미지를 묶어서 상기 새로운 3차원 데이터 큐브(910)를 생성한다.
또한, 상기 제어부(400)는 상기 생성된 3차원 데이터 큐브로부터 상기 배경장을 제거한다.
즉, 상기 제어부(400)는 상기 생성된 3차원 데이터 큐브에서 [r, θ]에 대해 시계열 이미지 값을 추출한다. 또한, 상기 제어부(400)는 [수학식 2]를 근거로 상기 추출된 시계열 이미지 값에서 상기 배경장을 추정한다. 또한, 상기 제어부(400)는 상기 추정된 배경장을 제거한다.
일 예로, 도 10에 도시된 바와 같이, 상기 제어부(400)는 상기 [r, θ]에 대해 시계열 이미지 값을 추출하고, 상기 추출된 시계열 이미지 값에서 상기 배경장(상기 도 10에서 적색 부분)(1010)을 추정한다.
또한, 도 11에 도시된 바와 같이, 상기 제어부(400)는 상기 추정된 배경장을 제거한다. 여기서, 상기 도 11에 도시된 좌측 그림은 배경장 제거 전을 나타내고, 우측 그림은 배경장을 제거한 후의 그림을 나타낸다. 또한, 상기 도 11에서 상단의 그림은 주어진 방위각에 대하여 시각 및 반경에 따른 영상의 값을 나타내고, 상기 도 11에서 하단의 그림은 주어진 시각에 대하여 방위각 및 반경에 따른 영상의 값을 나타낸다. 또한, 상기 도 11의 그림에서 코로나 질량 방출은 우측 상단 방향의 직선 형태로 나타나게 되는데, 배경장을 제거하기 전에는 잘 보이지 않는다. 또한, 상기 도 11의 하단의 그림에서는 배경장 제거 전에는 강하게 나타났던 부분이 배경장 제거 후에는 거의 사라진 것을 볼 수 있으며, 실제 코로나 질량 방출 신호에는 거의 영향을 미치지 않는다(S520).
이후, 상기 제어부(400)는 상기 전처리된 이미지(또는 상기 배경장이 제거된 3차원 데이터 큐브)를 근거로 코로나 질량 방출을 검출한다. 여기서, 상기 코로나 질량 방출 검출은 코로나 질량 방출 프론트 후보 검출 과정, 상기 검출된 코로나 질량 방출 프론트 후보에서 별 또는 노이즈에 의해 생긴 후보를 제거하는 과정, 상기 노이즈가 제거된 코로나 질량 방출 프론트 후보들에 대하여 시공간상의 클러스터링을 통해 그룹핑을 한 이후에 가장 앞서 진행하는 프로시딩 프론트를 검출하여 상기 코로나 질량 방출을 검출하는 과정 등을 포함한다.
즉, 상기 제어부(400)는 상기 전처리된 이미지(또는 상기 배경장이 제거된 3차원 데이터 큐브)의 [θ, t]에서 이미지 r-방향으로의 편차를 정의하는 에지 검출 방법을 사용하여 상기 코로나 질량 방출 프론트 후보를 검출한다.
일 예로, 도 12에 도시된 바와 같이, 상기 제어부(400)는 상기 [θ, t]에서 이미지 r 방향으로의 편차를 정의하는 에지 검출 방법을 사용하여 상기 코로나 질량 방출 프론트 후보를 검출한다. 여기서, 상기 도 12의 상단의 그림들은 폭이 3 픽셀이고 문턱값이 30 인 경우의 상기 코로나 질량 방출 프론트 후보이고, 상기 도 12의 하단의 그림들은 폭이 2 픽셀이고 문턱값이 25인 경우의 상기 코로나 질량 방출 프론트 후보이다.
또한, 상기 제어부(400)는 상기 검출된 코로나 질량 방출 프론트 후보에 대하여 노이즈에 해당하는 점을 제거한다.
즉, 상기 제어부(400)는 상기 검출된 코로나 질량 방출 프론트 후보 중에서 미리 설정된 노이즈 판별 기준에 해당하는 상기 코로나 질량 방출 프론트 후보를 제거한다. 여기서, 상기 미리 설정된 노이즈 판별 기준은 미리 설정된 r-방향 10 픽셀, 미리 설정된 θ-방향 10도의 범위 내에서 검출되는 코로나 질량 방출 프론트의 개수가 미리 설정된 3개 이하인 경우 노이즈로 구별한다.
일 예로, 도 13에 도시된 바와 같이, 상기 제어부(400)는 상기 검출된 코로나 질량 방출 프론트 후보에 대하여 노이즈에 해당하는 점을 제거한다. 여기서, 상기 도 13의 상단에 나타낸 그림들은 노이즈 제거 전 상태이고, 상기 도 13의 하단에 나타낸 그림들은 노이즈 제거 후 상태이다. 또한, 상기 도 13의 하단에 나타낸 그림들에서 볼 수 있듯이, 노이즈에 해당하는 점들이 잘 제거되고, 코로나 질량 방출은 잘 남아 있는 것을 확인할 수 있다.
또한, 상기 제어부(400)는 상기 노이즈가 제거된 코로나 질량 방출 프론트 후보들에 대하여 시공간상의 클러스터링인 계층 군집화(hierarchical clustering) 방법을 통해 그룹핑을 하여, 각 그룹을 구별한다. 또한, 상기 제어부(400)는 각 그룹에 대하여 주어진 시각에 가장 큰 거리(r)를 프로시딩 프론트로 인식(또는 검출)한다.
일 예로, 도 14에 도시된 바와 같이, 상기 제어부(400)는 계층 군집화를 통해 상기 노이즈가 제거된 코로나 질량 방출 프론트 후보들에 대한 그룹핑을 수행한다(S530).
이후, 상기 제어부(400)는 상기 검출된 코로나 질량 방출(또는 상기 프로시딩 프론트로 인식된 상기 코로나 질량 방출 프론트)의 물리량(또는 2차원 코로나 질량 방출의 물리량)을 산출(또는 도출)한다. 여기서, 상기 코로나 질량 방출(또는 상기 프로시딩 프론트로 인식된 상기 코로나 질량 방출 프론트)의 물리량 산출은 상기 코로나 질량 방출 프론트의 [r, θ, t] 좌표를 산출(또는 추출)하고, 상기 산출된 정보를 [x, y, t] 공간에서 그룹핑하는 과정, 상기 각 추출된 개별 그룹에 대하여 2차원 코로나 질량 방출 물리량을 산출하는 과정 등을 포함한다. 이때, 상기 산출되는 2차원 코로나 질량 방출의 물리량은 코로나 질량 방출의 방출 방향(또는 방위각), 방출 폭, 방출 속도, 방출 시각 등을 포함한다.
즉, 상기 제어부(400)는 최종적으로 검출된 상기 코로나 질량 방출 프론트의 [r, θ, t] 좌표를 각 그룹별로 추출한다. 또한, 상기 제어부(400)는 상기 추출된 [r, θ, t] 좌표를 다시 직교 좌표 공간 및 시간 공간에 해당하는 [x, y, t] 좌표로 변환한다. 또한, 상기 제어부(400)는 상기 변환된 [x, y, t] 3차원 공간에서 다시 계층 군집화를 수행한다.
이와 같이, 상기 제어부(400)는 상기 최종적으로 검출된 상기 코로나 질량 방출 프론트에 대해서 시공간을 통합한 코로나 질량 방출 그룹을 생성한다.
일 예로, 도 15에 도시된 바와 같이, 상기 제어부(400)는 상기 최종적으로 검출된 상기 코로나 질량 방출 프론트의 [r, θ, t] 좌표를 각 그룹별로 추출하고, 상기 추출된 [r, θ, t] 좌표를 [x, y, t] 좌표로 변환하고, 상기 변환된 [x, y, t] 3차원 공간에서 다시 계층 군집화를 수행하여 3차원 공간상에서 구별된 코로나 질량 방출 그룹을 생성한다. 여기서, 상기 도 15에 도시된 색상은 각각의 그룹을 나타낸다.
또한, 상기 제어부(400)는 상기 생성된 코로나 질량 방출 그룹을 근거로 2차원 코로나 질량 방출 물리량을 산출(또는 도출)한다. 여기서, 상기 2차원 코로나 질량 방출 물리량은 코로나 질량 방출의 방출 방향, 코로나 질량 방출의 방출 폭, 코로나 질량 방출의 방출 속도, 코로나 질량 방출의 방출 시각 등을 포함한다. 이때, 상기 코로나 질량 방출의 방출 방향의 기준은 북쪽 방향을 기준으로 반시계 방향 각도를 의미한다. 또한, 상기 코로나 질량 방출의 방출 폭은 전체 방출 폭을 의미한다. 또한, 상기 코로나 질량 방출의 방출 속도는 시간에 따른 상기 코로나 질량 방출 프론트의 위치에 대해 선형회귀 방법을 적용한 속도를 의미한다. 또한, 상기 코로나 질량 방출의 방출 시각은 상기 코로나 질량 방출의 방출 속도에 의해 상기 코로나 질량 방출이 태양 중심에 있을 때의 시각을 의미한다.
일 예로, 도 16에 도시된 바와 같이, 상기 제어부(400)는 상기 생성된 코로나 질량 방출 그룹을 근거로 상기 도 16의 우측 상단에 나타낸 상기 도 16의 좌측 3차원 그림에 표시한 청색 그룹에 대한 상기 코로나 질량 방출의 물리량을 산출하고, 상기 도 16의 우측 하단에 나타낸 상기 도 16의 좌측 3차원 그림에 표시한 적색 그룹에 대한 상기 코로나 질량 방출의 물리량을 산출한다(S540).
이후, 상기 제어부(400)는 상기 산출된 2차원 코로나 질량 방출의 물리량을 근거로 3차원 코로나 질량 방출의 물리량을 산출(또는 도출)한다. 여기서, 상기 3차원 코로나 질량 방출의 물리량은 태양 중심에서의 시각, 속도(또는 상기 산출된 3차원 공간상의 속도), 표준 편차(또는 상기 산출된3차원 공간상의 표준 편차), 상기 방출된 코로나 질량 방출의 폭, 방위각, 천정각 등을 포함한다.
일 예로, 도 17 및 도 18에 도시된 바와 같이, 상기 제어부(400)는 상기 산출된 2차원 코로나 질량 방출의 물리량을 근거로 3차원 코로나 질량 방출의 물리량을 산출(또는 도출)한다(S550).
본 발명의 실시예는 앞서 설명된 바와 같이, 소호 위성과 스테레오 위성의 코로나 그래프의 실시간 수집 체계를 구축하여, 신뢰성 있는 데이터를 제공할 수 있다.
또한, 본 발명의 실시예는 앞서 설명된 바와 같이, 실시간 수집된 코로나 그래프를 분석하여 코로나 질량 방출 발생을 자동으로 인식하고, 산출된 CME 물리량을 근거로 3차원 CME 물리량을 산출하여, 코로나 질량 방출의 방향, 방출폭, 초기 속도를 자동 인식하고, 3차원 CME 물리량을 제공하여 사용자의 편의성을 향상시킬 수 있다.
전술된 내용은 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자라면 본 발명의 본질적인 특성에서 벗어나지 않는 범위에서 수정 및 변형이 가능할 것이다. 따라서, 본 발명에 개시된 실시예들은 본 발명의 기술 사상을 한정하기 위한 것이 아니라 설명하기 위한 것이고, 이러한 실시예에 의하여 본 발명의 기술 사상의 범위가 한정되는 것은 아니다. 본 발명의 보호 범위는 아래의 청구범위에 의하여 해석되어야 하며, 그와 동등한 범위 내에 있는 모든 기술 사상은 본 발명의 권리범위에 포함되는 것으로 해석되어야 할 것이다.
본 발명은 소호 위성과 스테레오 위성의 코로나 그래프의 실시간 수집 체계를 구축하고, 상기 실시간 수집된 코로나 그래프를 분석하여 코로나 질량 방출 발생을 자동으로 인식함으로써 신뢰성 있는 데이터를 제공하고 코로나 질량 방출의 방향, 방출폭, 초기 속도를 자동 인식하고, 3차원 CME 물리량을 제공하여 사용자의 편의성을 향상시키는 것으로, 태양의 코로나 질량 방출 관찰 분야, 시계열 분석 분야 등에서 광범위하게 이용될 수 있다.
10: 코로나 질량 방출 분석 장치 100: 통신부
200: 저장부 300: 표시부
400: 제어부

Claims (10)

  1. 소호 위성 및 스테레오 A 위성과 스테레오 B 위성과 통신하는 코로나 질량 방출 분석 장치에 있어서,
    상기 소호 위성으로부터 제공되는 소호 라스코 C3 이미지를 수신하고, 상기 스테레오 A 위성 및 상기 스테레오 B 위성으로부터 각각 제공되는 스테레오 COR2 이미지를 수신하는 통신부; 및
    상기 수신된 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지에 대한 전처리 과정을 수행하고, 상기 전처리된 이미지를 근거로 코로나 질량 방출을 검출하고, 상기 검출된 코로나 질량 방출의 물리량을 산출하고, 상기 산출된 코로나 질량 방출의 물리량을 근거로 3차원 코로나 질량 방출의 물리량을 산출하는 제어부를 포함하는 코로나 질량 방출 분석 장치.
  2. 제1항에 있어서,
    상기 제어부는,
    상기 수신된 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지에 대한 좌표를 변환하고, 상기 좌표 변환된 복수의 이미지를 근거로 3차원 데이터 큐브를 생성하고, 상기 생성된 3차원 데이터 큐브로부터 배경장을 제거하는 상기 전처리 과정을 수행하는 것을 특징으로 하는 코로나 질량 방출 분석 장치.
  3. 제1항에 있어서,
    상기 제어부는,
    상기 전처리된 이미지를 근거로 코로나 질량 방출 프론트 후보를 검출하고, 상기 검출된 코로나 질량 방출 프론트 후보에서 별 또는 노이즈에 의해 생긴 후보를 제거하고, 상기 노이즈가 제거된 코로나 질량 방출 프론트 후보들에 대해 계층 군집화(hierarchical clustering) 방법을 적용하여 그룹핑 수행하고, 상기 그룹핑된 각 그룹에 대하여 주어진 시각에 가장 큰 거리를 프로시딩 프론트(proceeding front)로 인식하는 것을 특징으로 하는 코로나 질량 방출 분석 장치.
  4. 삭제
  5. 제1항에 있어서,
    상기 코로나 질량 방출의 물리량은,
    상기 코로나 질량 방출의 방출 방향, 상기 코로나 질량 방출의 방출 폭, 상기 코로나 질량 방출의 방출 속도 및 상기 코로나 질량 방출의 방출 시각 중 적어도 하나를 포함하는 것을 특징으로 하는 코로나 질량 방출 분석 장치.
  6. 제1항에 있어서,
    상기 3차원 코로나 질량 방출의 물리량은,
    태양 중심에서의 시각, 속도, 표준 편차, 상기 방출된 코로나 질량 방출의 폭, 방위각 및 천정각 중 적어도 하나를 포함하는 것을 특징으로 하는 코로나 질량 방출 분석 장치.
  7. 소호 위성 및 스테레오 A 위성과 스테레오 B 위성과 통신하는 코로나 질량 방출 분석 장치의 코로나 질량 방출 분석 방법에 있어서,
    통신부를 통해, 상기 소호 위성으로부터 제공되는 소호 라스코 C3 이미지를 수신하고, 상기 스테레오 A 위성 및 상기 스테레오 B 위성으로부터 각각 제공되는 스테레오 COR2 이미지를 수신하는 단계;
    제어부를 통해, 상기 수신된 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지에 대한 전처리 과정을 수행하는 단계;
    상기 제어부를 통해, 상기 전처리된 이미지를 근거로 코로나 질량 방출을 검출하는 단계;
    상기 제어부를 통해, 상기 검출된 코로나 질량 방출의 물리량을 산출하는 단계; 및
    상기 제어부를 통해, 상기 산출된 코로나 질량 방출의 물리량을 근거로 3차원 코로나 질량 방출의 물리량을 산출하는 단계를 포함하는 코로나 질량 방출 분석 방법.
  8. 제7항에 있어서,
    상기 전처리 과정을 수행하는 단계는,
    상기 수신된 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지에 대해 직교 좌표계를 극좌표계로 변환하는 과정;
    상기 극좌표계로 변환된 상기 소호 라스코 C3 이미지 및 상기 스테레오 COR2 이미지를 묶어 3차원 데이터 큐브를 생성하는 과정; 및
    상기 생성된 3차원 데이터 큐브로부터 배경장을 제거하는 과정;을 포함하는 것을 특징으로 하는 코로나 질량 방출 분석 방법.
  9. 제7항에 있어서,
    상기 코로나 질량 방출을 검출하는 단계는,
    상기 전처리된 이미지를 근거로 코로나 질량 방출 프론트 후보를 검출하는 과정;
    상기 검출된 코로나 질량 방출 프론트 후보에서 별 또는 노이즈에 의해 생긴 후보를 제거하는 과정;
    상기 노이즈가 제거된 코로나 질량 방출 프론트 후보들에 대해 계층 군집화 방법을 적용하여 그룹핑 수행하는 과정; 및
    상기 그룹핑된 각 그룹에 대하여 주어진 시각에 가장 큰 거리를 프로시딩 프론트로 인식하는 과정;을 포함하는 것을 특징으로 하는 코로나 질량 방출 분석 방법.
  10. 삭제
KR1020140099021A 2014-08-01 2014-08-01 코로나 질량 방출 분석 장치 및 그 방법 KR101471509B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020140099021A KR101471509B1 (ko) 2014-08-01 2014-08-01 코로나 질량 방출 분석 장치 및 그 방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020140099021A KR101471509B1 (ko) 2014-08-01 2014-08-01 코로나 질량 방출 분석 장치 및 그 방법

Publications (1)

Publication Number Publication Date
KR101471509B1 true KR101471509B1 (ko) 2014-12-11

Family

ID=52678490

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020140099021A KR101471509B1 (ko) 2014-08-01 2014-08-01 코로나 질량 방출 분석 장치 및 그 방법

Country Status (1)

Country Link
KR (1) KR101471509B1 (ko)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101575834B1 (ko) 2015-10-06 2015-12-09 대한민국(국립전파연구원 우주전파센터장) 태양 흑점군 추적 장치 및 그 방법
KR101786361B1 (ko) * 2016-06-29 2017-10-17 경희대학교 산학협력단 헤일로 cme의 3차원 인자 산출 시스템 및 그 방법
KR101914270B1 (ko) 2017-09-20 2018-11-01 대한민국(국립전파연구원 우주전파센터장) Halo cme의 특성 인자 산출 장치 및 그 방법
CN110533100A (zh) * 2019-07-22 2019-12-03 南京大学 一种基于机器学习进行cme检测和跟踪的方法
KR20220059258A (ko) * 2020-11-02 2022-05-10 경희대학교 산학협력단 태양 코로나그래프 관측 영상의 롤 각도 보정 방법 및 그 장치

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20100060705A (ko) * 2008-11-28 2010-06-07 한국 천문 연구원 코로나 물질 방출현상을 이용한 지자기 폭풍 예보 시스템 및 그 방법

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20100060705A (ko) * 2008-11-28 2010-06-07 한국 천문 연구원 코로나 물질 방출현상을 이용한 지자기 폭풍 예보 시스템 및 그 방법

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
CME 및 고속태양풍의 근지구 영향과정 규명연구(한국천문연구원) *
태양활동 및 우주환경연구(한국천문연구원 보현산천문대 2000년 발표자료) *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101575834B1 (ko) 2015-10-06 2015-12-09 대한민국(국립전파연구원 우주전파센터장) 태양 흑점군 추적 장치 및 그 방법
KR101786361B1 (ko) * 2016-06-29 2017-10-17 경희대학교 산학협력단 헤일로 cme의 3차원 인자 산출 시스템 및 그 방법
KR101914270B1 (ko) 2017-09-20 2018-11-01 대한민국(국립전파연구원 우주전파센터장) Halo cme의 특성 인자 산출 장치 및 그 방법
CN110533100A (zh) * 2019-07-22 2019-12-03 南京大学 一种基于机器学习进行cme检测和跟踪的方法
KR20220059258A (ko) * 2020-11-02 2022-05-10 경희대학교 산학협력단 태양 코로나그래프 관측 영상의 롤 각도 보정 방법 및 그 장치
KR102450395B1 (ko) 2020-11-02 2022-09-30 경희대학교 산학협력단 태양 코로나그래프 관측 영상의 롤 각도 보정 방법 및 그 장치

Similar Documents

Publication Publication Date Title
KR101471509B1 (ko) 코로나 질량 방출 분석 장치 및 그 방법
US11709058B2 (en) Path planning method and device and mobile device
US11636289B2 (en) Method, apparatus, and device for classifying LiDAR point cloud data, and storage medium
EP2849117B1 (en) Methods, apparatuses and computer program products for automatic, non-parametric, non-iterative three dimensional geographic modeling
CN103901884B (zh) 信息处理方法和信息处理设备
US9576373B2 (en) Geospatial imaging system providing segmentation and classification features and related methods
KR101531246B1 (ko) 기상 레이더 영상 내 대류 세포와 낙뢰 매칭 시스템 및 그의 제어 방법
CN113359782B (zh) 一种融合lidar点云与图像数据的无人机自主选址降落方法
CN112953383A (zh) 检测光伏板的方法、设备及无人机
CN105139434B (zh) 露天矿山道路网自动构建方法及系统
CN112305559A (zh) 基于地面定点激光雷达扫描的输电线距离测量方法、装置、系统和电子设备
WO2023138163A1 (zh) 一种基于深度图像修复的室内移动机器人玻璃检测与地图更新方法
AU2019349986A1 (en) Apparatus and method for three-dimensional object recognition
Nardinocchi et al. Fully automatic point cloud analysis for powerline corridor mapping
CN111337898A (zh) 激光点云的处理方法、装置、设备及存储介质
CN109685835B (zh) 一种海量卫星影像数据输电线路火点提取的方法和装置
KR101575834B1 (ko) 태양 흑점군 추적 장치 및 그 방법
Kurdi et al. Full series algorithm of automatic building extraction and modelling from LiDAR data
CN115909253A (zh) 一种目标检测、模型训练方法、装置、设备及存储介质
Zheng et al. Forest farm fire drone monitoring system based on deep learning and unmanned aerial vehicle imagery
KR101464718B1 (ko) 태양 흑점군 추적 장치 및 그 방법
Hujebri et al. Automatic building extraction from lidar point cloud data in the fusion of orthoimage
KR101859672B1 (ko) Halo cme 추출 장치 및 그 방법
Qayyum et al. Measuring height of high-voltage transmission poles using unmanned aerial vehicle (UAV) imagery
KR101581836B1 (ko) Euv 디밍 현상 감지 장치 및 그 방법

Legal Events

Date Code Title Description
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20190903

Year of fee payment: 5

R401 Registration of restoration
FPAY Annual fee payment

Payment date: 20191202

Year of fee payment: 6