CN104615990A - Method for automatically recognizing macula based on Huairou full-disk single-color image - Google Patents
Method for automatically recognizing macula based on Huairou full-disk single-color image Download PDFInfo
- Publication number
- CN104615990A CN104615990A CN201510067024.3A CN201510067024A CN104615990A CN 104615990 A CN104615990 A CN 104615990A CN 201510067024 A CN201510067024 A CN 201510067024A CN 104615990 A CN104615990 A CN 104615990A
- Authority
- CN
- China
- Prior art keywords
- sunspot
- huairou
- full
- macula
- solar
- Prior art date
- Legal status (The legal status 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 status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/10—Terrestrial scenes
- G06V20/13—Satellite images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Astronomy & Astrophysics (AREA)
- Remote Sensing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Analysis (AREA)
Abstract
The invention relates to a method for automatically recognizing macula based on a Huairou full-disk single-color image, and aims at solving the problems of extremely low efficiency caused by manually recognizing the profile of macula on disk and that mass data cannot be met in the solar physics in the domestic astronomy field. The method comprises the steps of extracting gradient information of the disk by the morphology bot-hat transformation according to the full-disk single-color image of the Huairou solar observing station of National Astronomical Observatories of China; removing the influence of inherent limb darkening rule of the disk by the fragmentation threshold method, so as to automatically recognize macula. Compared with other algorithms for recognizing the macula based on high-quality solar space data, the method has the characteristic that the macula of the ground images with relatively low quality and suffering from the interference of instrument noise can be precisely recognized. The method is applicable to automatic recognition of macula in the single-color images or white light images at other wave bands and can provide directly service to the fields of generator theory, solar action predication and space environment monitoring.
Description
Technical field
The invention belongs to the image processing techniques in Solar Physics field, particularly a kind of automatic identifying method being applicable to sunspot in sun image.
Background technology
Sunspot is sun feature the most eye-catching on day face, is the most important index characterizing solar activity level.Analyze sunspot to studying the long-term mechanics of the sun, forecasting the harm that the power of solar flare and monitoring solar activity bring the mankind, all play very important effect.The attribute studying sunspot accurately must identify it.Early stage people are identified by naked eyes.But, in recent years along with observing capacity constantly promotes, new large-scale sun observation equipment constantly comes into operation, a large amount of sunspot observation data is had every day to produce, based on manual type because work efficiency is very low, and accuracy cannot ensure, cannot the demand of satisfying magnanimity data, utilize Computer Automatic Recognition sunspot to become active demand and trend.
Sun observation base, Huairou (HSOS) is one of important base of State Astronomical Observatory, CAS's astronomical sight and research, full-time Vector Magnetic Field telescope of base equipment 10cm bore, full-time the filtergram recording sunspot can be obtained every day, carry out the automatic identification of sunspot based on these images, can be relevant investigation and application and service is provided.
Automatically being identified in of sunspot is domesticly still belonged to blank, proposes certain methods abroad, such as: Zharkov et al. (2005) etc. utilizes frontier probe method and threshold method to extract the profile of sunspot; Curto et al. (2008) utilizes top cap conversion to extract the candidate region of sunspot, then screens these regions in conjunction with the parameter restriction of iteration; Suruchi et al. (2014) then make use of the method for level set to identify sunspot.On the whole, these methods mainly through extracting the gradient of sunspot, then carry out a definite limitation to identify sunspot to gradient.
Above all methods are all carry out sunspot identification based on the spatial data that quality is higher.But full-time photosphere filtergram of Huairou belongs to ground data, is subject to atmospheric interference, limited resolution during observation; The more important thing is, partial data is subject to noise of instrument interference, there is larger noise of instrument, and the position of noise of instrument in different pieces of information is not fixed near center, day face, and existing algorithm does not process for these problems.The present invention is directed to full-time the filtergram of China, propose a kind of new sunspot automatic identifying method, the interference of noise of instrument can be eliminated, improve the precision of sunspot identification in ground data, improve the utilization rate of China's sun observation data, for further investigation solar cycle, sun forecast and space weather monitoring provide convenient.
Leading reference annex
Summary of the invention
The object of the invention is for the some shortcomings existed in the filtergram image of full-time of Huairou (compared to external spatial data, resolution is lower, and be subject to the interference of noise of instrument), a kind of sunspot automatic identifying method proposed, has filled up domestic technological gap in this respect.
For achieving the above object, the present invention takes following design proposal:
(1) first, the heterogeneity of day face existence itself and the interference of cloud layer, the way based on global threshold is made to carry out sunspot identification infeasible, in order to overcome this problem, we have employed the morphology bot-hat transformation that can extract dark object under bright background, and (formula is as (1), (2), (3), (4), A represents original image, B and g representative structure element), remove a day face background, extract the gradient information on day face;
Dilation operation:
Erosion operation:
Closed operation:
Bot-hat transformation: h=A-(A.g) (4)
(2) secondly, owing to being subject to the impact of the intrinsic limb darkening rule of the sun, the gradient information extracted exists near Mian center large, the problem that solar limb is little, to this, we adopt fragmentation threshold to process: utilizing higher thresholds near Mian center, adopt comparatively Low threshold to split respectively at solar limb place, obtain the candidate region of sunspot; Meanwhile, although noise of instrument is near Mian center, its Grad is lower relative to sunspot gradient, therefore, adopts higher thresholds to split, also substantially can eliminate noise of instrument in center, day face.Threshold value H chooses formula as (5), and (x, y) represents the coordinate of each pixel, and R represents day radius surface;
(3) last, because the center gray-scale value in each sunspot region is high, edge gray-scale value is low, and the difference limiting the gray scale of maximum, the minimum pixel of each candidate region at this must be greater than 5 pixels, and remove the candidate region not meeting this condition, finally obtain the sunspot on day face.
In sum, the present invention's advantage compared with existing several method is:
1. process of the present invention is picture quality lower ground face observation data, and the present invention is sensitiveer for the extraction at sunspot edge.
2. the present invention easily causes the noise of instrument of interference to eliminate to sunspot identification by day face, and this is that other technologies are not considered.
3. the sunspot that the present invention also can be used for other solar space data identifies automatically.
Accompanying drawing explanation
Fig. 1 is full-time filtergram original graph.
Fig. 2 carries out the result after closed operation process to original graph.
Fig. 3 is the gradient information on day face.
Fig. 4 is the binary map after adopting fragmentation threshold.
Fig. 5 is the recognition result of sunspot
Embodiment
Embodiments of the invention accompanying drawings is as follows:
The concrete implementation step of sunspot automatic identifying method of the present invention is as follows:
First, to raw data (see Fig. 1), closed operation is utilized to process, obtain a clean day face (see Fig. 2), the structural element be wherein set using to be radius be 15 circle, this, radius of hairy nevus was slightly larger radius ratio day, can eliminate black moles all on day face.
The second, Fig. 2 and Fig. 1 is subtracted each other, obtains the gradient information (see Fig. 3) on day face;
3rd, other method is utilized to extract the profile (developing relative program) of solar limb, calculate the radius in day face, and then the sun face in Fig. 3 is divided into two regions: the region in day 0.8, face radius is set to I1, the region of a day face 0.8-1 solar radius is set to I2.Threshold value H1 is adopted respectively to I1 and I2, H2 splits Fig. 3, obtain binary map 4.According to statistics, we choose H1=22 and H2=10.
4th, extract each region on Fig. 4, the difference meeting pixel maxima and minima is greater than to the region of 5 pixels, thinks sunspot, and with red-label, the region not meeting this condition is then abandoned.The sunspot of extraction is superimposed in Fig. 1 by we, obtains Fig. 5, and the red area in this figure represents the sunspot of identification of the present invention.
The various embodiments described above can in addition some changes under not departing from the scope of the present invention, therefore above explanation comprises and should be considered as exemplary, and is not used to the protection domain limiting the present patent application patent.
Claims (4)
1. the sunspot automatic identifying method based on full-time the filtergram in Huairou, it is characterized in that: compare other algorithms based on high-quality solar space data identification sunspot, the present invention is directed to quality lower and by noise of instrument interference ground image, higher to the accuracy of identification of sunspot.The step of the method is as follows:
(1) first, the heterogeneity of day face itself existence and be subject to the interference of cloud layer, the way based on global threshold is made to carry out sunspot identification infeasible, in order to overcome this problem, we have employed the morphology bot-hat transformation that can extract dark object under bright background, and (formula is as (1), (2), (3), (4), A represents original image, B and g representative structure element), remove a day face background, extract the gradient information on day face;
(2) secondly, owing to being subject to the impact of the intrinsic limb darkening rule of the sun, the gradient information extracted exists near Mian center large, the problem that solar limb is little, to this, we adopt fragmentation threshold to process: utilizing higher thresholds near Mian center, adopt comparatively Low threshold to split respectively at solar limb place, obtain the candidate region of sunspot; Meanwhile, although noise of instrument is near Mian center, its Grad is lower relative to sunspot gradient, therefore, adopts higher thresholds to split, also substantially can eliminate noise of instrument in center, day face;
(3) last, because the center gray-scale value in each sunspot region is high, edge gray-scale value is low, and the difference limiting the gray scale of maximum, the minimum pixel of each candidate region at this must be greater than 5 pixels, and remove the candidate region not meeting this condition, finally obtain the sunspot on day face.
2. a kind of sunspot automatic identifying method based on full-time the filtergram in Huairou according to claim 1, it is characterized in that: step 1) in, propose to utilize morphology bot-hat transformation to extract the gradient information in day face, morphologic formula is as shown in (1), (2), (3), (4).
Dilation operation:
Erosion operation:
Closed operation:
Bot-hat transformation: h=A-(A.g) (4)
3. a kind of sunspot automatic identifying method based on full-time the filtergram in Huairou according to claim 1, it is characterized in that: step 2) in, for the interference of the intrinsic limb darkening rule in day face, and day face noise of instrument place to go, propose to utilize the method for fragmentation threshold to split, threshold value H chooses formula as (5).
4. a kind of sunspot automatic identifying method based on full-time the filtergram in Huairou according to claim 1, it is characterized in that: step 3) in, to each sunspot region, the difference defining the gray scale of maximum, the minimum pixel in region must be greater than 5 pixels.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510067024.3A CN104615990A (en) | 2015-02-10 | 2015-02-10 | Method for automatically recognizing macula based on Huairou full-disk single-color image |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510067024.3A CN104615990A (en) | 2015-02-10 | 2015-02-10 | Method for automatically recognizing macula based on Huairou full-disk single-color image |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104615990A true CN104615990A (en) | 2015-05-13 |
Family
ID=53150427
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510067024.3A Pending CN104615990A (en) | 2015-02-10 | 2015-02-10 | Method for automatically recognizing macula based on Huairou full-disk single-color image |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104615990A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106293663A (en) * | 2015-05-26 | 2017-01-04 | 中国科学院国家天文台 | Sun observation data in Huairou processes software |
CN106570506A (en) * | 2016-10-26 | 2017-04-19 | 昆明理工大学 | Solar activity recognition method based on scale transformation model |
CN107580160A (en) * | 2016-06-30 | 2018-01-12 | 比亚迪股份有限公司 | Remove method, image processor and the camera arrangement in sunspot region in image |
CN108009471A (en) * | 2017-10-25 | 2018-05-08 | 昆明理工大学 | It is a kind of that method for distinguishing is known based on genetic algorithm and the sunspot of simulated annealing |
CN110851627A (en) * | 2019-09-24 | 2020-02-28 | 昆明理工大学 | Method for describing sun black subgroup in full-sun image |
CN113052202A (en) * | 2021-01-29 | 2021-06-29 | 昆明理工大学 | Method for classifying sun black subgroup in full-sun image |
CN113207902A (en) * | 2021-04-06 | 2021-08-06 | 东莞市生命伞生物科技有限公司 | CAZ physical black cluster family harmful microorganism and virus resisting technology |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101173984A (en) * | 2007-11-21 | 2008-05-07 | 中国科学院上海技术物理研究所 | Spaceborne target detection tracing camera in sun viewing blind zone |
CN103605171A (en) * | 2013-09-10 | 2014-02-26 | 国家电网公司 | All-sky imaging instrument and cloud layer characteristic analysis method based on all-sky imaging instrument |
-
2015
- 2015-02-10 CN CN201510067024.3A patent/CN104615990A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101173984A (en) * | 2007-11-21 | 2008-05-07 | 中国科学院上海技术物理研究所 | Spaceborne target detection tracing camera in sun viewing blind zone |
CN103605171A (en) * | 2013-09-10 | 2014-02-26 | 国家电网公司 | All-sky imaging instrument and cloud layer characteristic analysis method based on all-sky imaging instrument |
Non-Patent Citations (5)
Title |
---|
J.J.CURTO等: "Automatic Sunspots Detection on Full-Disk Solar Images using Mathematical Morphology", 《SPRINGER》 * |
S.I.ZHARKOV等: "The latitudinal distributions of sunspot areas and magnetic", 《ADVANCES IN SPACE RESEARCH》 * |
朱士虎: "形态学高帽变换与低帽变换功能扩展及应用", 《计算机工程与应用》 * |
李真真等: "亮背景下形态学星点目标提取算法", 《光电工程》 * |
郑泽忠: "基于高分辨率航空影像高速公路汽车目标检测算法研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106293663A (en) * | 2015-05-26 | 2017-01-04 | 中国科学院国家天文台 | Sun observation data in Huairou processes software |
CN107580160A (en) * | 2016-06-30 | 2018-01-12 | 比亚迪股份有限公司 | Remove method, image processor and the camera arrangement in sunspot region in image |
CN106570506A (en) * | 2016-10-26 | 2017-04-19 | 昆明理工大学 | Solar activity recognition method based on scale transformation model |
CN106570506B (en) * | 2016-10-26 | 2020-10-27 | 昆明理工大学 | Solar activity recognition method based on scale transformation model |
CN108009471A (en) * | 2017-10-25 | 2018-05-08 | 昆明理工大学 | It is a kind of that method for distinguishing is known based on genetic algorithm and the sunspot of simulated annealing |
CN110851627A (en) * | 2019-09-24 | 2020-02-28 | 昆明理工大学 | Method for describing sun black subgroup in full-sun image |
CN113052202A (en) * | 2021-01-29 | 2021-06-29 | 昆明理工大学 | Method for classifying sun black subgroup in full-sun image |
CN113207902A (en) * | 2021-04-06 | 2021-08-06 | 东莞市生命伞生物科技有限公司 | CAZ physical black cluster family harmful microorganism and virus resisting technology |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104615990A (en) | Method for automatically recognizing macula based on Huairou full-disk single-color image | |
CN105956557B (en) | A kind of sequential remote sensing image cloud covered areas domain automatic testing method of object-oriented | |
CN103761717B (en) | A kind of urban water-body extracting method based on panchromatic remote sensing image | |
CN102063727B (en) | Covariance matching-based active contour tracking method | |
CN104197933B (en) | Method for enhancing and extracting star sliding in telescope field of view | |
CN107341809B (en) | Multi-source remote sensing image-oriented coastline automatic extraction and deviation correction method | |
CN102509151A (en) | Video-processing-based crowd density and distribution estimation method | |
CN102750683B (en) | Filtering method for sea surface stripe noise and sea surface stripe cloud in moderate resolution imaging spectroradiometer (MODIS) remote sensing image | |
CN104966274B (en) | A kind of On Local Fuzzy restored method using image detection and extracted region | |
CN105138992A (en) | Coastline detection method based on regional active outline model | |
CN104966291A (en) | Cloud cluster automatic detection method based on foundation cloud atlas | |
CN109801304B (en) | Building index construction and building area rapid extraction method for inhibiting bare land | |
CN116503470A (en) | Method for segmenting, positioning and analyzing states of insulator string disc in infrared image | |
CN103345742B (en) | Based on a kind of detection method improving the Remote Sensing Imagery Change of Markov random field model | |
CN104463822A (en) | Multi-focus image fusing method and device based on multi-scale overall filtering | |
CN103971369A (en) | Optic disc positioning method for retina image | |
CN109472220A (en) | A kind of substation's worker safety helmet detection method and its system based on Faster R-CNN | |
CN104778469B (en) | Automatic extraction method of sun outline | |
CN103487844B (en) | Intelligent coronal mass emission (CME) event observation procedure | |
Lu et al. | An image recognition algorithm based on thickness of ice cover of transmission line | |
CN104182938B (en) | Solar facula repairing method of all-sky nephogram | |
He et al. | Glacier identification from Landsat8 OLI imagery using deep U-NET | |
CN104616265A (en) | Time domain reconstruction method of remote sensing sequence data | |
Zhen-Yan | Chinese character recognition method based on image processing and hidden markov model | |
CN104268836A (en) | Watershed segmentation mark point extraction method based on local area homogeneity indexes |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150513 |