CN109801241A - A kind of solar flare image based on modified dark priority algorithm removes cloud method - Google Patents
A kind of solar flare image based on modified dark priority algorithm removes cloud method Download PDFInfo
- Publication number
- CN109801241A CN109801241A CN201910059000.1A CN201910059000A CN109801241A CN 109801241 A CN109801241 A CN 109801241A CN 201910059000 A CN201910059000 A CN 201910059000A CN 109801241 A CN109801241 A CN 109801241A
- Authority
- CN
- China
- Prior art keywords
- image
- value
- dark
- indicate
- cloud
- 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
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000001914 filtration Methods 0.000 claims abstract description 15
- 230000005540 biological transmission Effects 0.000 claims abstract description 13
- 230000002146 bilateral effect Effects 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims abstract description 6
- 238000005457 optimization Methods 0.000 claims abstract description 6
- 238000010586 diagram Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 14
- 239000003595 mist Substances 0.000 abstract description 5
- 125000001475 halogen functional group Chemical group 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000003760 hair shine Effects 0.000 description 1
- 239000005433 ionosphere Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Landscapes
- Image Processing (AREA)
Abstract
A kind of solar flare image based on modified dark priority algorithm removes cloud method, comprising: describes the concept in channel using mathematic(al) representation；Assuming that cloud atlas piece expression formula: I (x)=J (x) t (x)+A (1-t (x)) seeks global atmosphere light value A.Size from dark channel diagram according to brightness takes preceding 0.1% pixel.In these positions, the value of the corresponding point with maximum brightness is found in original foggy image, as A value.Bilateral filtering is used when solving t (x) to improve.Calculate the local mean value and Local standard deviation of dark image D (x, y): according to known airlight vector A, calculation optimization projects graph expression formula.In view of when the value very little of transmission plot t, the value that will lead to J is bigger than normal, so that it is whole excessive to white field to make image, therefore a state value t generally can be set0, when t value is less than t0When, enable t=t0.Therefore the image expression formula finally restored.A kind of solar flare image based on modified dark priority algorithm of the present invention removes cloud method, by modified dark concept to cloud and mist image restoring, goes cloud effect with good.
Description
Technical field
The invention belongs to technical field of image processing, and in particular to a kind of sun based on modified dark priority algorithm
Flare image removes cloud method.
Background technique
The sun is the fixed star most close with human relation, has inseparable pass with the life of the mankind and production activity
System.It is filled with magnetic field in solar atmosphere, stores huge magnetic energy.It, can be by too when the magnetic energy being stored in magnetic field is excessive
Positive outburst activity releases energy, and solar flare is most violent one of solar activity outburst form.And the outburst of solar flare is then by shadow
The ionosphere in earth magnetic field and its overhead is rung, and then influences the mankind's activities such as satellite navigation, radio communication.Therefore it shines to the sun
The understanding of spot phenomenon and observation and people realize space exploration and carry out the real needs of precautionary measures.
Due to can all be influenced by earth atmosphere cloud layer, therefore to flare image when observing solar flare phenomenon
Cloud removing is increasingly important.The dark channel prior that doctor He Mingkai in 2009 proposes is theoretical, can effectively go to image
Cloud processing, but due to having used minimum filtering in dark, obtained transmissivity contains halo effect and blocky effect,
In order to solve this problem, soft-matting algorithm (is chosen using soft-matting and Steerable filter optimization algorithm herein
Compare) Lai Youhua transmissivity, wherein soft-matting can eliminate halo phenomenon and blocky phenomenon well, but at that time
Between complexity greatly increase；Steerable filter Algorithms T-cbmplexity is smaller, but its recovery after image edge region there are still
A degree of cloud, therefore can use bilateral filtering algorithm to optimize soft-matting algorithm and guiding filtering, with optimization
The calculating of transmissivity.
Summary of the invention
In order to solve the technical issues of above-mentioned solar flare image removes cloud, the present invention provides one kind to be helped secretly based on modified
The solar flare image of road priority algorithm removes cloud method, by modified dark concept to cloud and mist image restoring, has good
Go cloud effect.
The technical scheme adopted by the invention is as follows:
A kind of solar flare image based on modified dark priority algorithm removes cloud method, comprising the following steps:
Step 1: the concept in channel being described using mathematic(al) representation, for image J, dark JdarkIt can be expressed as public affairs
Formula:
In above formula, JcThe color channel image of expression image J, Ω (x) one piece of territory of expression, pixel center x,
Y is any point in Ω (x) territory, and c indicates r, tri- channels g, b, then Jc(y) it means that in Ω (x) territory
The image in the channel c at middle y.Its meaning is exactly to find out the minimum value of tri- components of RGB, is then carried out to the width single channel figure minimum
Value filtering.Dark channel prior is pointed out: Jdark→ 0, it is equivalent to Jdark≈ 0, therefore true J can be acquired according to this conditionc。
Step 2: assuming that cloud atlas piece expression formula (2):
I (x)=J (x) t (x)+A (1-t (x)) (2)
Wherein, I (x) represents figure to be processed, and J (x) represents true figure, and t (x) represents transmissivity, and expression can reach
The part light that computer system is not scattered, A indicate global atmosphere light value.
Step 3: seeking global atmosphere light value A.Size from dark channel diagram according to brightness takes preceding 0.1% pixel.?
In these positions, the value of the corresponding point with maximum brightness is found in original foggy image, as A value.
Bilateral filtering is used when solving t (x) to improve.
Step 4: calculating the local mean value and Local standard deviation of dark image D (x, y), then atmosphere light is estimated by the difference of the two
Curtain:
Wherein,Indicate atmosphere light curtain, D (x, y) indicates dark image, B1(x, y) indicates the office of dark image D (x, y)
Portion's mean value, B2(x, y) indicates the Local standard deviation of dark image D (x, y), FB(x, y) indicates to use the algorithmic function of bilateral filtering.
Due toIt is the local mean value of D (x, y) and the difference of Local standard deviation,
Then:
Step 5: according to known global atmosphere light value A, calculation optimization projects graph expression formula (6)
T (x, y) indicates transmission rate matrix in above formula,Indicate atmosphere light curtain.
Step 6: in view of when the value very little of transmission plot t, the value that will lead to J is bigger than normal, to keep image whole to white field
Excessively, therefore generally a state value t can be set0, when t value is less than t0When, enable t=t0.Therefore the image expression formula finally restored
Are as follows:
J (x, y) indicates that treated image, I (x) indicate that original image, t (x, y) indicate transmission rate matrix, A table in formula
Show global atmosphere light value, t0State the threshold value of the transmission plot t chosen.
A kind of solar flare image based on modified dark priority algorithm of the present invention removes cloud method, and technical effect is such as
Under:
1: cloud removing being carried out to solar flare image with this method, and is emulated.Cloud is gone to it according to simulation result
Effect is assessed.The results show that improved dark priority algorithm significant effect in terms of solar flare removes cloud and mist.
2: innovatory algorithm proposed by the present invention substantially reduces and calculates the time, reduces machine calculation complexity.At image
Contrast on effect after reason, improved algorithm have obtained relatively clear image, and improve the block estimated in transmittance figure
Shape phenomenon；While this paper algorithm refines transmissivity, the effect at smoothed image edge is also acted.
3: innovatory algorithm proposed by the present invention compared to primal algorithm the average gray of image, detailed information display,
Picture information quantity and opposite readability show more preferably face to face, it is easier to observe position and the image of solar flare.
Detailed description of the invention
Present invention will be further explained below with reference to the attached drawings and examples:
Fig. 1 is that the solar flare image for the modified dark priority algorithm that the embodiment of the present invention proposes removes the stream of cloud method
Cheng Tu
Fig. 2 is solar flare image before cloud.
Fig. 3 (1) is the solar flare image handled with original dark priority algorithm.
Fig. 3 (2) is the solar flare image handled with modified dark priority algorithm.
Specific embodiment
Dark channel prior algorithm applies to the defogging processing of image first, due to the imaging model of cloud and the imaging model of mist
It is similar, it is all the sum of the energy by the radiation information of target through overdamping and atmosphere light through both overdamping, therefore using dark
Channel prior knowledge carries out cloud removing to image using the algorithm.
Based on improved dark channel prior algorithm.Since primal algorithm is solving transmissivity t (x, y) Shi Yunyong soft-
Matting algorithm, so that computation complexity is greatly increased with the time is calculated, in order to mitigate computation burden, the present invention is with bilateral
Filtering algorithm replaces soft-matting process, ensure that the quality of image procossing while having reached purpose.
A kind of solar flare image based on modified dark priority algorithm removes cloud method, comprising the following steps:
Step 1: the concept in channel being described using mathematic(al) representation, for image J, dark JdarkIt can be expressed as public affairs
Formula:
In above formula, JcThe color channel image of expression image J, Ω (x) one piece of territory of expression, pixel center x,
Y is any point in Ω (x) territory, and c indicates r, tri- channels g, b, then Jc(y) it means that in Ω (x) territory
The image in the channel c at middle y.Its meaning is exactly to find out the minimum value of tri- components of RGB, is then carried out to the width single channel figure minimum
Value filtering.Dark channel prior is pointed out: Jdark→ 0, it is equivalent to Jdark≈ 0, therefore true J can be acquired according to this conditionc。
Step 2: assuming that cloud atlas piece expression formula (2):
I (x)=J (x) t (x)+A (1-t (x)) (2)
Wherein, I (x) represents figure to be processed, and J (x) represents true figure, and t (x) represents transmissivity, and expression can reach
The part light that computer system is not scattered, A indicate global atmosphere light value.
Step 3: seeking global atmosphere light value A.Size from dark channel diagram according to brightness takes preceding 0.1% pixel.?
In these positions, the value of the corresponding point with maximum brightness is found in original foggy image, as A value.
Bilateral filtering is used when solving t (x) to improve.
Step 4: calculating the local mean value and Local standard deviation of dark image D (x, y), then atmosphere light is estimated by the difference of the two
Curtain:
Wherein,Indicate atmosphere light curtain, D (x, y) indicates dark image, B1(x, y) indicates the office of dark image D (x, y)
Portion's mean value, B2(x, y) indicates the Local standard deviation of dark image D (x, y), FB(x, y) indicates to use the algorithmic function of bilateral filtering.
Due toIt is the local mean value of D (x, y) and the difference of Local standard deviation,
Then:
Step 5: according to known global atmosphere light value A, calculation optimization projects graph expression formula (6)
T (x, y) indicates transmission rate matrix in above formula,Indicate atmosphere light curtain.
Step 6: in view of when the value very little of transmission plot t, the value that will lead to J is bigger than normal, to keep image whole to white field
Excessively, therefore generally a state value t can be set0, when t value is less than t0When, enable t=t0.Therefore the image expression formula finally restored
Are as follows:
J (x, y) indicates that treated image, I (x) indicate that original image, t (x, y) indicate transmission rate matrix, A table in formula
Show global atmosphere light value, t0State the threshold value of the transmission plot t chosen.
Fig. 1 is the flow chart that the solar flare image based on modified dark priority algorithm removes cloud method；This method is
Original solar flare image is handled, when solving transmissivity with bilateral filtering algorithm to soft-matting algorithm with
Guiding filtering is replaced, to optimize the calculating of transmissivity.
A solar flare image is now randomly selected, such as Fig. 2 carries out cloud removing to it, programs with MATLAB, algorithm
Implementation process inputs original image first, solves its dark J as shown in Fig. 1 processdark(x)；Again by calculation described in text
Method obtains the value of global atmosphere light, big to estimate then by calculating the local mean value and Local standard deviation of dark image D (x, y)
Gas light curtain.Then according to known airlight vector, the transmissivity t (x) after projection figure is optimized is calculated.Finally according to formula
(7) distribution function for obtaining cloud atlas picture, restored just obtained it is final remove cloud atlas picture, as shown in Fig. 3 (2).It is simultaneously
Compared with the algorithm before improvement, by the solar flare image of selection with original dark priority algorithm at
Reason, obtains cloud atlas picture, as shown in Fig. 3 (1).
In conjunction with gray value characteristic, gray value is higher at solar flare, therefore display is brighter, and it is in not advise that cloud and mist gray value is also higher
Then light tone.Therefore by Fig. 2 and Fig. 3 (1), Fig. 3 (2) no matter comparison is it can be seen that improve front and back, dark priority algorithm removes cloud
Effect is all very significant, but by the comparison of Fig. 3 (1), Fig. 3 (2) it is found that higher using improved algorithm picture clarity, too
Positive solar flare is more obvious, embodies the superiority of innovatory algorithm.
Claims (1)
1. a kind of solar flare image based on modified dark priority algorithm removes cloud method, it is characterised in that including following step
It is rapid:
Step 1: the concept in channel being described using mathematic(al) representation, for image J, dark JdarkIt can be expressed as formula:
In above formula, JcIndicate that the color channel image of image J, Ω (x) indicate one piece of territory, pixel center x, y Ω
(x) any point in territory, c indicate r, tri- channels g, b, then Jc(y) c at y is meant that in Ω (x) territory
The image in channel；Its meaning is exactly to find out the minimum value of tri- components of RGB, then carries out minimum value filter to the width single channel figure
Wave；Dark channel prior is pointed out: Jdark→ 0, it is equivalent to Jdark≈ 0, therefore true J can be acquired according to this conditionc；
Step 2: assuming that cloud atlas piece expression formula (2):
I (x)=J (x) t (x)+A (1-t (x)) (2)
Wherein, I (x) represents figure to be processed, and J (x) represents true figure, and t (x) represents transmissivity, and expression can reach calculating
The part light that machine system is not scattered, A indicate global atmosphere light value；
Step 3: seeking global atmosphere light value A；Size from dark channel diagram according to brightness takes preceding 0.1% pixel；At these
In position, the value of the corresponding point with maximum brightness is found in original foggy image, as A value；
Bilateral filtering is used when solving t (x) to improve；
Step 4: the local mean value and Local standard deviation of dark image D (x, y) are calculated, then atmosphere light curtain is estimated by the difference of the two:
Wherein,Indicate atmosphere light curtain, D (x, y) indicates dark image, B1(x, y) indicates the part of dark image D (x, y)
Value, B2(x, y) indicates the Local standard deviation of dark image D (x, y), FB(x, y) indicates to use the algorithmic function of bilateral filtering；
Due toIt is the local mean value of D (x, y) and the difference of Local standard deviation,Then:
Step 5: according to known global atmosphere light value A, calculation optimization projects graph expression formula (6)
T (x, y) indicates transmission rate matrix in above formula,Indicate atmosphere light curtain；
Step 6: in view of when the value very little of transmission plot t, the value that will lead to J is bigger than normal, thus keep image whole excessive to white field,
Therefore a state value t generally can be set0, when t value is less than t0When, enable t=t0；Therefore the image expression formula finally restored are as follows:
J (x, y) indicates that treated image, I (x) indicate that original image, t (x, y) indicate that transmission rate matrix, A indicate complete in formula
Ball air light value, t0State the threshold value of the transmission plot t chosen.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910059000.1A CN109801241A (en) | 2019-01-22 | 2019-01-22 | A kind of solar flare image based on modified dark priority algorithm removes cloud method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910059000.1A CN109801241A (en) | 2019-01-22 | 2019-01-22 | A kind of solar flare image based on modified dark priority algorithm removes cloud method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109801241A true CN109801241A (en) | 2019-05-24 |
Family
ID=66559968
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910059000.1A Pending CN109801241A (en) | 2019-01-22 | 2019-01-22 | A kind of solar flare image based on modified dark priority algorithm removes cloud method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109801241A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113989142A (en) * | 2021-10-25 | 2022-01-28 | 黑龙江科技大学 | Infrared image processing method |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130071043A1 (en) * | 2011-09-08 | 2013-03-21 | Fujitsu Limited | Image defogging method and system |
| US20160071244A1 (en) * | 2014-09-04 | 2016-03-10 | National Taipei University Of Technology | Method and system for image haze removal based on hybrid dark channel prior |
| CN106157267A (en) * | 2016-07-12 | 2016-11-23 | 中国科学技术大学 | A kind of image mist elimination absorbance optimization method based on dark channel prior |
| CN106204491A (en) * | 2016-07-12 | 2016-12-07 | 中国科学技术大学 | A kind of adapting to image defogging method based on dark channel prior |
| CN106530246A (en) * | 2016-10-28 | 2017-03-22 | 大连理工大学 | Image dehazing method and system based on dark channel and non-local prior |
| CN107133926A (en) * | 2017-04-13 | 2017-09-05 | 中原智慧城市设计研究院有限公司 | A kind of image defogging method based on dark primary priori |
| CN107730472A (en) * | 2017-11-03 | 2018-02-23 | 昆明理工大学 | A kind of image defogging optimized algorithm based on dark primary priori |
| CN108154492A (en) * | 2017-12-25 | 2018-06-12 | 北京航空航天大学 | A kind of image based on non-local mean filtering goes haze method |
-
2019
- 2019-01-22 CN CN201910059000.1A patent/CN109801241A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130071043A1 (en) * | 2011-09-08 | 2013-03-21 | Fujitsu Limited | Image defogging method and system |
| US20160071244A1 (en) * | 2014-09-04 | 2016-03-10 | National Taipei University Of Technology | Method and system for image haze removal based on hybrid dark channel prior |
| CN106157267A (en) * | 2016-07-12 | 2016-11-23 | 中国科学技术大学 | A kind of image mist elimination absorbance optimization method based on dark channel prior |
| CN106204491A (en) * | 2016-07-12 | 2016-12-07 | 中国科学技术大学 | A kind of adapting to image defogging method based on dark channel prior |
| CN106530246A (en) * | 2016-10-28 | 2017-03-22 | 大连理工大学 | Image dehazing method and system based on dark channel and non-local prior |
| CN107133926A (en) * | 2017-04-13 | 2017-09-05 | 中原智慧城市设计研究院有限公司 | A kind of image defogging method based on dark primary priori |
| CN107730472A (en) * | 2017-11-03 | 2018-02-23 | 昆明理工大学 | A kind of image defogging optimized algorithm based on dark primary priori |
| CN108154492A (en) * | 2017-12-25 | 2018-06-12 | 北京航空航天大学 | A kind of image based on non-local mean filtering goes haze method |
Non-Patent Citations (3)
| Title |
|---|
| KAIMING HE等: "Single Image Haze Removal Using Dark Channel Prior", 《IEEE TRANSACTIONS ON PATTERN ANALYSIS AND MACHINE INTELLIGENCE》 * |
| 孙抗等: "基于双边滤波的实时图像去雾技术研究", 《北京理工大学学报》 * |
| 李娟等: "基于局部均值和标准差的图像去雾算法", 《计算机工程与应用》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113989142A (en) * | 2021-10-25 | 2022-01-28 | 黑龙江科技大学 | Infrared image processing method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103955905B (en) | Based on the single image to the fog method that fast wavelet transform and weighted image merge | |
| CN106157267B (en) | Image defogging transmissivity optimization method based on dark channel prior | |
| CN106204491B (en) | A kind of adapting to image defogging method based on dark channel prior | |
| CN102831591B (en) | Gaussian filter-based real-time defogging method for single image | |
| CN107767354A (en) | A kind of image defogging algorithm based on dark primary priori | |
| CN104253930B (en) | A kind of real-time video defogging method | |
| CN104867121B (en) | Image Quick demisting method based on dark primary priori and Retinex theories | |
| CN105046658B (en) | A kind of low-light (level) image processing method and device | |
| CN110232661A (en) | Low illumination colour-image reinforcing method based on Retinex and convolutional neural networks | |
| Wang et al. | Variational single nighttime image haze removal with a gray haze-line prior | |
| CN106846263A (en) | The image defogging method being immunized based on fusion passage and to sky | |
| CN103530848A (en) | Double exposure implementation method for inhomogeneous illumination image | |
| CN102750674A (en) | Video image defogging method based on self-adapting allowance | |
| CN108564549A (en) | A kind of image defogging method based on multiple dimensioned dense connection network | |
| CN110675340A (en) | Single image defogging method and medium based on improved non-local prior | |
| WO2024060576A1 (en) | Image dehazing method based on dark channel prior | |
| CN103578083A (en) | Single image defogging method based on joint mean shift | |
| CN104299198A (en) | Fast image defogging method based on dark channels of pixels | |
| CN108629750A (en) | A kind of night defogging method, terminal device and storage medium | |
| Liang et al. | Learning to remove sandstorm for image enhancement | |
| CN109685735B (en) | Single picture defogging method based on fog layer smoothing prior | |
| CN107977941A (en) | A kind of bright areas color fidelity and the image defogging method of contrast enhancing | |
| CN102968767A (en) | Method for real-time restoration of fog-degraded image with white balance correction | |
| CN108898561B (en) | Defogging method, server and system for foggy image containing sky area | |
| CN104715456B (en) | A kind of defogging method of image |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190524 |
|
| RJ01 | Rejection of invention patent application after publication |