CN113010996B - Method and system for extracting radon concentration abnormity based on entropy-median filtering in sub-region - Google Patents
Method and system for extracting radon concentration abnormity based on entropy-median filtering in sub-region Download PDFInfo
- Publication number
- CN113010996B CN113010996B CN202110147521.XA CN202110147521A CN113010996B CN 113010996 B CN113010996 B CN 113010996B CN 202110147521 A CN202110147521 A CN 202110147521A CN 113010996 B CN113010996 B CN 113010996B
- Authority
- CN
- China
- Prior art keywords
- median
- sub
- value
- parameter
- value filtering
- 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.)
- Active
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 94
- 229910052704 radon Inorganic materials 0.000 title claims abstract description 62
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000002689 soil Substances 0.000 claims abstract description 61
- 230000002159 abnormal effect Effects 0.000 claims abstract description 46
- 238000000605 extraction Methods 0.000 claims description 5
- 239000000758 substrate Substances 0.000 claims 4
- 229910052770 Uranium Inorganic materials 0.000 abstract description 8
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 abstract description 8
- 239000000284 extract Substances 0.000 abstract description 4
- 230000033558 biomineral tissue development Effects 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 description 17
- 238000010586 diagram Methods 0.000 description 3
- 230000005856 abnormality Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 101100032908 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) RAD7 gene Proteins 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Image Processing (AREA)
- Measurement Of Radiation (AREA)
Abstract
Description
技术领域technical field
本发明涉及铀矿勘探技术领域,特别涉及一种基于信息熵-子区中位数衬值滤波提取土壤氡浓度异常区域的方法和系统。The invention relates to the technical field of uranium ore exploration, in particular to a method and a system for extracting abnormal areas of soil radon concentration based on information entropy-sub-area median contrast value filtering.
背景技术Background technique
随着核能在能源结构的比重增加,对铀资源的需求越来越大,促进了铀矿资源的勘探与开发。With the increasing proportion of nuclear energy in the energy structure, the demand for uranium resources is increasing, which promotes the exploration and development of uranium resources.
子区中位数衬值滤波法是在EDA(勘查数据分析)技术和滤波技术的基础上建立的无需对原始数据作处理而提取异常信息的一种化学勘探数据处理方法,常用于1:20万、1:5万等金属矿床地球化学弱小异常提取。需要指出的是,子区中位数衬值滤波法(SAMCF)数据处理对工作经验要求较高,例如在目前的工作实践中,子区中位数衬值滤波大、小窗口的参数选择就在很大程度上是依据工作人员的经验,不同大、小窗口条件下子区中位数衬值滤波法获取的正、负异常分布范围不同,异常信息提取区域存在不确定性。The sub-area median contrast filtering method is a chemical exploration data processing method established on the basis of EDA (exploration data analysis) technology and filtering technology without processing the original data and extracting abnormal information. It is commonly used in 1:20 10,000, 1:50,000 and other metal deposits are extracted from weak geochemical anomalies. It should be pointed out that the data processing of the sub-region median contrast filter method (SAMCF) requires high work experience. To a large extent, based on the experience of the staff, the distribution range of positive and negative anomalies obtained by the median contrast filtering method of sub-regions under different large and small window conditions is different, and there is uncertainty in the area of abnormal information extraction.
发明内容SUMMARY OF THE INVENTION
本发明的目的之一是提供一种基于信息熵-子区中位数衬值滤波提取土壤氡浓度异常区域的方法,以解决依赖工作经验选择大、小窗口参数而导致异常信息提取区域存在不确定性的问题。One of the purposes of the present invention is to provide a method for extracting abnormal areas of soil radon concentration based on information entropy-sub-area median contrast value filtering, so as to solve the problem of inconsistencies in abnormal information extraction areas caused by relying on work experience to select large and small window parameters. Certainty issues.
为解决上述技术问题,本发明采用如下技术方案:基于信息熵- 子区中位数衬值滤波提取土壤氡浓度异常区域的方法,按下式确定子区中位数衬值滤波中大、小窗口参数比值:In order to solve the above-mentioned technical problems, the present invention adopts the following technical scheme: based on the information entropy-sub-area median contrast value filtering to extract the method for the abnormal area of soil radon concentration, the following formula is used to determine the medium and small values of the sub-area median contrast value filtering. Window parameter ratio:
S(m/n)|min=∫m/n(-1)×P(i)×lnP(i) (1);S(m/n)| min =∫ m/n (-1)×P(i)×lnP(i) (1);
式(1)中,m为大窗口的参数值,n为小窗口的参数值,m/n为大、小窗口的参数比,S(m/n)|min为大、小窗口的参数比为m/n时所对应的子区中位数衬值滤波图像的最小信息熵值,P(i)为大窗口参数值为m、小窗口参数值为n的条件下所对应的子区中位数衬值滤波图像灰度化后第i级灰度出现的概率;In formula (1), m is the parameter value of the large window, n is the parameter value of the small window, m/n is the parameter ratio of the large and small windows, and S(m/n)| min is the parameter ratio of the large and small windows. is the minimum information entropy value of the median-contrast filter image of the corresponding sub-area when m/n, P(i) is the corresponding sub-area under the condition that the large window parameter value is m and the small window parameter value is n The probability of the occurrence of the i-th grayscale after the grayscale of the bit-contrast filter image;
根据式(1)所确定的参数比值确定大、小窗口的参数值,用与所述参数值相对应的大、小窗口对土壤氡浓度等值线图作子区中位数衬值滤波处理,根据子区中位数衬值滤波结果,提取对应的土壤氡浓度异常区域。Determine the parameter values of the large and small windows according to the parameter ratio determined by the formula (1), and use the large and small windows corresponding to the parameter values to filter the median line value of the soil radon concentration contour map in the sub-region. , and extract the corresponding abnormal area of soil radon concentration according to the filtering result of the median lining value of the sub-area.
进一步来说,还包括测量土壤氡浓度,形成土壤氡浓度等值线图的步骤。Further, it also includes the steps of measuring soil radon concentration and forming a soil radon concentration contour map.
进一步来说,还包括采用参数值为m的大窗口、参数值为n的小窗口对土壤氡浓度等值线图作子区中位数衬值滤波处理并对子区中位数衬值滤波结果作灰度化处理的步骤。Further, it also includes the use of a large window with a parameter value of m and a small window with a parameter value of n to filter the sub-region median contrast value of the soil radon concentration contour map and filter the sub-region median contrast value. The result is a grayscale processing step.
进一步来说,在将相应的子区中位数衬值滤波结果灰度化处理后,还包括确定该子区中位数衬值滤波图像灰度化后第i级灰度出现概率的步骤。Further, after the corresponding sub-area median-contrast filtering result is grayscaled, the step of determining the occurrence probability of the i-th grayscale after the sub-area median-contrast filtering image is grayscaled.
进一步来说,按下式确定相应的子区中位数衬值滤波图像灰度化后第i级灰度出现概率:Further, the probability of occurrence of the i-th grayscale after grayscale of the corresponding sub-region median-contrast filter image is determined by the following formula:
P(i)=N(i)/N (2);P(i)=N(i)/N(2);
式(2)中,P(i)为该子区中位数衬值滤波图像灰度化后第i级灰度出现的概率;N(i)为该子区中位数衬值滤波图像灰度化后第i级灰度出现的总次数;N为该子区中位数衬值滤波图像灰度化后像素总数目。In formula (2), P(i) is the probability of the occurrence of the i-th grayscale of the median-contrast filter image of the sub-region after graying; N(i) is the gray-scale of the median-contrast filter image of the sub-region. The total number of occurrences of the i-th grayscale after the grayscale; N is the total number of pixels after grayscale of the median-contrast filter image in this sub-region.
基于上述提取土壤氡浓度异常区域的方法,本发明还提出一种基于信息熵-子区中位数衬值滤波提取土壤氡浓度异常区域的系统,采用的技术方案如下,该系统包括:Based on the above-mentioned method for extracting abnormal areas of soil radon concentration, the present invention also proposes a system for extracting abnormal areas of soil radon concentration based on information entropy-sub-area median contrast filtering. The adopted technical scheme is as follows, and the system includes:
参数确定单元,按下式(1)计算子区中位数衬值滤波中大、小窗口的参数比值,并根据计算出的参数比值确定大、小窗口的参数:The parameter determination unit calculates the parameter ratio of the large and small windows in the median contrast filter of the sub-region as follows, and determines the parameters of the large and small windows according to the calculated parameter ratio:
S(m/n)|min=∫m/n(-1)×P(i)×lnP(i) (1);S(m/n)| min =∫ m/n (-1)×P(i)×lnP(i) (1);
式(1)中,m为大窗口的参数值,n为小窗口的参数值,m/n为大、小窗口的参数比,S(m/n)|min为大、小窗口的参数比为m/n时所对应的子区中位数衬值滤波图像的最小信息熵值,P(i)为大窗口参数值为m、小窗口参数值为n的条件下所对应的子区中位数衬值滤波图像灰度化后第i级灰度出现的概率。In formula (1), m is the parameter value of the large window, n is the parameter value of the small window, m/n is the parameter ratio of the large and small windows, and S(m/n)| min is the parameter ratio of the large and small windows. is the minimum information entropy value of the median-contrast filter image of the corresponding sub-area when m/n, P(i) is the corresponding sub-area under the condition that the large window parameter value is m and the small window parameter value is n The probability of the occurrence of the i-th grayscale after the grayscale of the bit-contrast filter image.
图像处理单元,根据所述参数确定单元所确定的大、小窗口参数,用与参数相对应的大、小窗口对土壤氡浓度等值线图作子区中位数衬值滤波处理并输出子区中位数衬值滤波结果。The image processing unit, according to the large and small window parameters determined by the parameter determination unit, uses the large and small windows corresponding to the parameters to filter the subregion median contrast value of the soil radon concentration contour map and output the subregion. District median contrast filter results.
异常区域提取单元,根据图像处理单元输出的子区中位数衬值滤波结果,从中提取对应的土壤氡浓度异常区域。The abnormal area extraction unit extracts the corresponding abnormal area of soil radon concentration according to the filtering result of the median contrast value of the sub-area output by the image processing unit.
图像处理单元还用于采用参数值为m的大窗口、参数值为n的小窗口对土壤氡浓度等值线图作子区中位数衬值滤波处理并对相应的子区中位数衬值滤波结果作灰度化处理。The image processing unit is also used to filter the sub-area median contrast value of the soil radon concentration contour map by using a large window with a parameter value of m and a small window with a parameter value of n, and perform the corresponding sub-area median contrast value filtering processing. The value filtering result is grayscaled.
参数确定单元还用于在图像处理单元对子区中位数衬值滤波结果作灰度化处理后确定该子区中位数衬值滤波图像灰度化后第i级灰度出现的概率。The parameter determining unit is further configured to determine the probability of the occurrence of the i-th grayscale of the sub-region median contrast value filtered image after grayscale processing of the subregion median contrast value filtering result by the image processing unit.
参数确定单元按下式确定对应的子区中位数衬值滤波图像灰度化后第i级灰度出现概率:The parameter determination unit determines the probability of occurrence of the i-th grayscale after the grayscale of the corresponding sub-area median-contrast filter image is as follows:
P(i)=N(i)/N (2);P(i)=N(i)/N(2);
式(2)中,P(i)为该子区中位数衬值滤波图像灰度化后第i级灰度出现的概率;N(i)为该子区中位数衬值滤波图像灰度化后第i级灰度出现的总次数;N为该子区中位数衬值滤波图像灰度化后像素总数目。In formula (2), P(i) is the probability of the occurrence of the i-th grayscale of the median-contrast filter image in the sub-region after grayscale; N(i) is the median-contrast-filtered image gray in the sub-region. is the total number of occurrences of the i-th grayscale after the grayscale; N is the total number of pixels after grayscale of the median-contrast filter image in this sub-region.
本发明所取得的有益效果在于:本发明将整个测量区域看作一个提供多元信息的信息源,将信息熵与子区中位数衬值滤波图像处理技术相结合,构建了子区中位数衬值滤波中大、小窗口比值与图像熵的函数关系模型,通过函数关系模型能够计算得出子区中位数衬值滤波中最佳大、小窗口的参数比值,并由此确定后续子区中位数衬值滤波处理时的大、小窗口参数值,该方法避免了人为经验因素导致子区中位数衬值滤波数据处理结果存在不确定性的问题,保证了子区中位数衬值滤波对异常区域识别的准确度,能为圈定铀矿勘探位置和矿化范围提供有价值的参考依据。The beneficial effects obtained by the present invention are as follows: the present invention regards the entire measurement area as an information source that provides multiple information, combines the information entropy with the sub-area median contrast filter image processing technology, and constructs the sub-area median The functional relationship model between the ratio of large and small windows and the image entropy in the contrast filter, the parameter ratio of the optimal large and small windows in the median contrast filter of the subregion can be calculated through the functional relationship model, and the subsequent subregions can be determined accordingly. This method avoids the problem of uncertainty in the data processing results of the median contrast value filtering in the sub-area caused by human experience factors, and ensures the median of the sub-area. The accuracy of the contrast filter to identify abnormal areas can provide a valuable reference for delineating the exploration location and mineralization range of uranium deposits.
附图说明Description of drawings
图1为测量区域土壤氡浓度等值线图;Figure 1 is a contour map of soil radon concentration in the measurement area;
图2为实施例中不同窗口参数条件下测量区土壤氡浓度子区中位数衬值滤波结果;Fig. 2 is the filtering result of median contrast value of soil radon concentration sub-area in measurement area under different window parameter conditions in the embodiment;
图3为实施例中不同差值方法与图像熵值关系图;Fig. 3 is the relation diagram of different difference value method and image entropy value in the embodiment;
图4为实施例中不同差值方法下土壤氡浓度栅格图;Fig. 4 is a grid diagram of soil radon concentration under different difference methods in the embodiment;
图5为实施例中大、小窗口参数比与图像信息熵关系图。FIG. 5 is a graph showing the relationship between the parameter ratio of large and small windows and the image information entropy in the embodiment.
具体实施方式Detailed ways
为了便于本领域技术人员的理解,下面结合实施例与附图对本发明作进一步的说明,实施方式提及的内容并非对本发明的限定。In order to facilitate the understanding of those skilled in the art, the present invention will be further described below with reference to the embodiments and the accompanying drawings, and the contents mentioned in the embodiments are not intended to limit the present invention.
本发明相对于现有技术的改进之处主要在于将整个测量区域看作一个提供多元信息的信息源,将信息熵与子区中位数衬值滤波图像处理技术相结合,构建了子区中位数衬值滤波中大、小窗口比值与图像熵值的函数关系模型,通过函数关系模型计算得出子区中位数衬值滤波中最佳大、小窗口的参数比值,并由此确定后续子区中位数衬值滤波处理时的大、小窗口参数值,从而避免人为经验因素导致子区中位数衬值滤波数据处理结果存在不确定性的问题。Compared with the prior art, the improvement of the present invention mainly lies in that the whole measurement area is regarded as an information source that provides multivariate information, and the information entropy is combined with the sub-area median contrast filtering image processing technology to construct a sub-area The functional relationship model between the ratio of large and small windows and the image entropy value in the median contrast filter, and the parameter ratio of the optimal large and small windows in the median contrast filter of the sub-region is calculated through the functional relationship model, and is determined from this. The large and small window parameter values in the subsequent sub-region median contrast value filtering processing, so as to avoid the problem of uncertainty in the sub-region median contrast value filtering data processing results caused by human experience factors.
下面就上述改进进行详细描述。The above improvements will be described in detail below.
与现有技术不同,本发明所涉基于信息熵-子区中位数衬值滤波提取土壤氡浓度异常区域的方法是按下式确定子区中位数衬值滤波中大、小窗口参数比值:Different from the prior art, the method for extracting abnormal areas of soil radon concentration based on information entropy-sub-area median contrast value filtering according to the present invention is to determine the ratio of large and small window parameters in the sub-area median contrast value filtering according to the following formula: :
S(m/n)|min=∫m/n(-1)×P(i)×lnP(i) (1);S(m/n)| min =∫ m/n (-1)×P(i)×lnP(i) (1);
上述式(1)中,m为大窗口的参数值,n为小窗口的参数值,m/n为大、小窗口的参数比,S(m/n)|min为大、小窗口的参数比为m/n时所对应的子区中位数衬值滤波图像的最小信息熵值,P(i)为大窗口参数值m、小窗口参数值为n的条件下所对应的子区中位数衬值滤波图像灰度化后第i级灰度出现的概率;之后根据式(1)所确定的参数比值确定大、小窗口的参数值,再用与前述参数值相对应的大、小窗口对土壤氡浓度等值线图作子区中位数衬值滤波处理,最后根据子区中位数衬值滤波结果,提取对应的土壤氡浓度异常区域。In the above formula (1), m is the parameter value of the large window, n is the parameter value of the small window, m/n is the parameter ratio of the large and small windows, and S(m/n)| min is the parameter of the large and small windows When the ratio is m/n, the minimum information entropy value of the median-contrast filter image of the corresponding sub-area, P(i) is the parameter value m of the large window and the parameter value of the small window is the corresponding sub-area under the condition of n The probability of occurrence of the i-th grayscale after the grayscale of the bit-contrast filter image; then determine the parameter values of the large and small windows according to the parameter ratio determined by the formula (1), and then use the large and small windows corresponding to the aforementioned parameter values. In the small window, the contour map of soil radon concentration is filtered by the median contrast value of the sub-region, and finally the corresponding abnormal area of soil radon concentration is extracted according to the filtering result of the median contrast value of the sub-region.
上述基于信息熵-子区中位数衬值滤波提取土壤氡浓度异常区域的方案的总体思路为:与现有的子区中位数衬值滤波方式一样,首先,需要划定测量区域并测量该区域的土壤氡浓度,并形成测量区的土壤氡浓度等值线图。然后,采用多组参数比值不同的大、小窗口(大窗口参数为m×m,小窗口参数为n×n)为对土壤氡浓度等值线图作子区中位数衬值滤波处理,并对相应的子区中位数衬值滤波结果作灰度化处理,进而确定对应的子区中位数衬值滤波图像灰度化后第i级灰度出现的概率。The general idea of the above scheme for extracting abnormal areas of soil radon concentration based on information entropy-sub-area median contrast value filtering is as follows: the same as the existing sub-area median contrast value filtering method, first of all, it is necessary to delineate the measurement area and measure Soil radon concentration in this area, and form a contour map of soil radon concentration in the measurement area. Then, multiple groups of large and small windows with different parameter ratios (the parameter of the large window is m×m, and the parameter of the small window is n×n) are used to filter the median line value of the soil radon concentration contour map in the sub-region. The corresponding sub-area median-contrast filtering results are processed to grayscale, and then the probability of the i-th grayscale appearing after the corresponding sub-area median-contrast filtering image is grayscaled is determined.
前述第i级灰度根据实际勘探要求确定,具体来说,子区中位数衬值滤波图像灰度化后第i级灰度出现概率由下式计算得出, P(i)=N(i)/N;式中,P(i)为该子区中位数衬值滤波图像灰度化后第 i级灰度出现的概率;N(i)为该子区中位数衬值滤波图像灰度化后第i 级灰度出现的总次数;N为该子区中位数衬值滤波图像灰度化后像素总数目。The above-mentioned i-th grayscale is determined according to the actual exploration requirements. Specifically, the occurrence probability of the i-th grayscale after the grayscale of the median contrast filter image of the sub-region is calculated by the following formula, P(i)=N( i)/N; in the formula, P(i) is the probability of the occurrence of the i-th grayscale after the grayscale of the median-contrast filter image in the sub-area; N(i) is the median-contrast filter of the sub-area The total number of occurrences of the i-th grayscale after the image is grayed; N is the total number of pixels after the grayscale of the median-contrast filter image in this sub-area.
然后参考式(1)的函数关系模型计算当大、小窗口参数比值为m/n、对应子区中位数衬值滤波图像灰度化后第i级灰度出现的概率为P(i)时的图像信息熵值,根据图像信息熵值与大、小窗口参数比值m/n之间的关系确定最优(图像信息熵值最小时)的大、小窗口的参数比,进而确定大、小窗口参数值。接着,用与前述参数值相对应的大、小窗口对土壤氡浓度等值线图再次作子区中位数衬值滤波处理并输出最终的子区中位数衬值滤波结果。Then, with reference to the functional relationship model of formula (1), when the ratio of the parameters of the large and small windows is m/n, and the corresponding sub-area median contrast filter image is grayscaled, the probability of the ith level grayscale appearing is P(i) According to the relationship between the image information entropy value and the ratio m/n of the large and small window parameters, determine the optimal (when the image information entropy value is the smallest) the parameter ratio of the large and small windows, and then determine the large and small windows. Small window parameter value. Next, use the large and small windows corresponding to the aforementioned parameter values to filter the sub-region median contrast value again on the soil radon concentration contour map and output the final sub-region median contrast value filtering result.
最后,根据最终的子区中位数衬值滤波结果,从中提取对应的土壤氡浓度异常区域。需要指出的是,在本发明中,从子区中位数衬值滤波结果中提取土壤氡浓度异常区域的方式与现有技术相比没有实质性的改变,为简化表述,对该部分内容不再赘述。Finally, according to the filtering result of the final sub-area median lining value, the corresponding abnormal area of soil radon concentration is extracted from it. It should be pointed out that, in the present invention, the method of extracting the abnormal area of soil radon concentration from the filtering result of the median value of the sub-area has no substantial change compared with the prior art. Repeat.
下面结合具体的实例对本发明作进一步说明,先划定长8km、宽 7.5km的区域作为测量区,然后采用静电收集法RAD7型α能谱氡气检测仪对测量区进行野外土壤氡浓度测量,测线距离为500m、测量点距离为100m、土壤深度70cm,形成如图1所示的土壤氡浓度等值线图。The present invention will be further described below in conjunction with specific examples, first delineate an area with a length of 8km and a width of 7.5km as a measurement area, and then adopt the electrostatic collection method RAD7 type alpha energy spectrum radon gas detector to measure the field soil radon concentration in the measurement area, The distance of the survey line is 500m, the distance of the measurement point is 100m, and the soil depth is 70cm, forming a contour map of soil radon concentration as shown in Figure 1.
子区中位数衬值滤波法是在EDA(勘查数据分析)技术和滤波技术基础上开发出来的,以稳健统计学为基础,无需对原始数据作任何处理而提取异常的化探数据处理方法,公式如下:The sub-area median contrast filtering method is developed on the basis of EDA (exploration data analysis) technology and filtering technology. It is based on robust statistics and does not require any processing of the original data to extract abnormal geochemical exploration data processing methods. , the formula is as follows:
Fu=Qu+1.5Sh F u =Q u +1.5S h
F1=Q1-1.5Sh F 1 =Q 1 -1.5S h
CF,P=MWc/Fu C F, P = M Wc /F u
CF,N=MWc/F1 C F,N =M Wc /F 1
式中:Fu表示异常点下限;F1表示异常点上限;Qu表示上4分点;Q1表示下4分点,Sh表示内散度。当CF,P>1,即为正异常;当 CF,N<1,即为负异常。In the formula: F u represents the lower limit of abnormal points; F 1 represents the upper limit of abnormal points; Q u represents the upper 4 points; Q 1 represents the lower 4 points, and Sh represents the inner divergence. When CF, P > 1, it is a positive abnormality; when CF, N <1, it is a negative abnormality.
首先,采用多组参数比值不同的大、小窗口(大窗口参数为6× 6、9×9,小窗口参数为2×2、3×3)对图1所示的土壤氡浓度等值线图作子区中位数衬值滤波处理,并对相应的子区中位数衬值滤波结果做灰度化处理,得到如图2所示的子区中位数衬值滤波结果,当大窗口小于9×9时,子区中位数衬值滤波结果图幅边界效益较强;当小窗口增大时,子区中位数衬值滤波结辨识的异常区域普遍减少。当采用小窗口3×3、大窗口9×9,子区中位数衬值滤波辨识测量区西部和西南部正异常7处,而采用小窗口2×2、大窗口9×9时,子区中位数衬值滤波辨识测量区西部和西南部正异常8处。根据上述采用不同的参数比值的大、小窗口分析发现子区中位数衬值滤波中正异常信号范围与矿床在空间上的分布存在联系。First, use multiple sets of large and small windows with different parameter ratios (large window parameters are 6 × 6, 9 × 9, and small window parameters are 2 × 2, 3 × 3) to compare the soil radon concentration contours shown in Figure 1. Figure 2 shows the sub-area median contrast filter processing, and performs grayscale processing on the corresponding sub-area median contrast value filtering results to obtain the sub-area median contrast value filtering results shown in Figure 2. When the window is smaller than 9×9, the effect of the sub-area median contrast filter results is stronger; when the small window increases, the abnormal areas identified by the sub-area median contrast filter generally decrease. When a small window of 3×3 and a large window of 9×9 are used, and the median contrast filter of the sub-area is used to identify 7 positive anomalies in the west and southwest of the measurement area, and when a small window of 2×2 and a large window of 9×9 are used, the 8 positive anomalies in the west and southwest of the measurement area were identified by filtering the median contrast value of the area. According to the above analysis of large and small windows using different parameter ratios, it is found that there is a relationship between the range of positive anomalous signals in the median contrast filter of sub-regions and the distribution of deposits in space.
接下来,确定相应的子区中位数衬值滤波图像灰度化后第i级灰度出现概率,由下式计算得出:P(i)=N(i)/N;式中,P(i)为该子区中位数衬值滤波图像灰度化后第i级灰度出现的概率;N(i)为该子区中位数衬值滤波图像灰度化后第i级灰度出现的总次数;N为该子区中位数衬值滤波图像灰度化后像素总数目。当所有图像灰度化后像素点的概率分布都相同时,图像距离焦点最远,成像最模糊,图像熵值最大;反之,图像熵值最小,图像最为清晰[19]。Next, determine the occurrence probability of the i-th grayscale after grayscale of the corresponding sub-area median-contrast filter image, which is calculated by the following formula: P(i)=N(i)/N; in the formula, P (i) is the probability of the occurrence of the i-th grayscale of the median-contrast filter image of this sub-region after grayscale; N(i) is the i-th grayscale of the sub-region’s median-contrast filtered image after grayscale The total number of occurrences of the degree; N is the total number of pixels after grayscale of the median contrast filter image in this sub-area. When the probability distribution of all pixels after grayscale is the same, the image is the farthest from the focus, the image is the most blurred, and the image entropy value is the largest; otherwise, the image entropy value is the smallest and the image is the clearest [19] .
然后通过不同的差值方法建立大、小窗口参数与信息熵值的关系图,在Gridata中常用Linear、Cubic、Neareast和V4四种插值方法,如图3所示,由图可看出当大、小窗口参数比值较小时,信息熵值较为离散,存在局部差异;随着大、小窗口参数比值增大,信息熵偏差不大;接着对土壤氡浓度等值线图在Cubic、Linear、Nearest、V4的差值方法下做栅格化处理,获得如图4所示的土壤氡浓度栅格图。Then, the relationship diagram between large and small window parameters and information entropy value is established through different difference methods. Linear, Cubic, Neareast and V4 interpolation methods are commonly used in Gridata, as shown in Figure 3. , When the ratio of small window parameters is small, the information entropy value is relatively discrete, and there are local differences; as the ratio of large and small window parameters increases, the information entropy deviation is not large; , V4 difference method to do rasterization, to obtain the soil radon concentration grid map shown in Figure 4.
再构建大、小窗口参数比值与图像信息熵值的函数关系模型,如式(1)所示;Then build a functional relationship model between the ratio of large and small window parameters and the image information entropy value, as shown in formula (1);
S(m/n)|min=∫m/n(-1)×P(i)×lnP(i) (1);S(m/n)| min =∫ m/n (-1)×P(i)×lnP(i) (1);
通过上述函数关系模型计算当大、小窗口参数比值为m/n、对应子区中位数衬值滤波图像灰度化后第i级灰度出现的概率为P(i)时的图像信息熵值,根据图像信息熵值与大、小窗口参数比值m/n之间的关系,确定最优(图像信息熵值最小时)的大、小窗口参数比值,进而确定大、小窗口参数值;然后用与前述参数值相对应的大、小窗口对土壤氡浓度等值线图再次作子区中位数衬值滤波处理并输出最终的子区中位数衬值滤波结果;最后根据最终的子区中位数衬值滤波结果,从中提取对应的土壤氡浓度异常区域。Calculate the image information entropy when the ratio of the parameters of the large and small windows is m/n, and the probability of the occurrence of the i-th grayscale after the grayscale of the corresponding sub-region median contrast filter image is P(i) by the above functional relationship model According to the relationship between the image information entropy value and the ratio m/n of the large and small window parameters, determine the optimal (when the image information entropy value is the smallest) the ratio of the large and small window parameters, and then determine the large and small window parameter values; Then use the large and small windows corresponding to the aforementioned parameter values to filter the sub-area median contrast value on the soil radon concentration contour map again, and output the final sub-area median contrast value filtering result; The filtering result of the median lining value of the sub-region is used to extract the corresponding abnormal area of soil radon concentration.
通过结合图2和图5分析发现,位于图5中的大小窗口参数比值起始段,图像信息熵值存在低值,但是信息熵值变化幅度大,呈现多期次周期性变化,稳定性较差。但是图2表明,小窗口2×2和3×3、大窗口6×6时,子区中位数衬值滤波结果在测量区东北区域存在明显的边界效益,即“伪正异常”区域。总体来说,测量区信息熵—子区中位数衬值滤波最佳大、小窗口参数比值为15~20。Through the analysis of Figure 2 and Figure 5, it is found that the image information entropy value has a low value in the initial segment of the size window parameter ratio in Figure 5, but the information entropy value changes greatly, showing multiple periodic changes, and the stability is relatively low. Difference. However, Figure 2 shows that when the small window is 2×2 and 3×3, and the large window is 6×6, the median contrast filtering result of the sub-area has obvious boundary benefits in the northeast area of the measurement area, that is, the “pseudo-positive abnormal” area. In general, the optimal ratio of large and small window parameters for the information entropy of the measurement area to the median contrast value of the sub-area filter is 15-20.
综上所述,测量区信息熵-子区中位数衬值滤波模型中信息熵值为2.5~3.6,图像信息熵值随大、小窗口参数周期性波动,随大小窗口参数比值总体呈现“周期波动—下凹型”多项式变化规律。大小窗口参数比值较小或者过大时,信息熵-子区中位数衬值滤波图像信息熵值大,子区中位数衬值滤波结果可信度降低,测量区最佳区中位数衬值滤波大小窗口参数比值为15~20,图像信息熵值最小;最后通过最小信息熵-子区中位数衬值滤波图像结果表明测量区西部异常区域为本区深部铀资源勘探的有利区域。To sum up, the information entropy value in the measurement area information entropy-sub-area median contrast filter model is 2.5 to 3.6, and the image information entropy value fluctuates periodically with the parameters of large and small windows, and the ratio of the parameters of the large and small windows generally presents " Periodic fluctuation - concave type" polynomial variation law. When the ratio of the size window parameter is small or too large, the information entropy of the filtered image with the median contrast value of the sub-area is large, the reliability of the filtering result of the median-contrast value of the sub-area is reduced, and the median of the optimal area of the measurement area is reduced. The contrast filter size window parameter ratio is 15 to 20, and the image information entropy value is the smallest. Finally, the minimum information entropy-sub-area median contrast filter image results show that the abnormal area in the west of the measurement area is a favorable area for deep uranium resource exploration in this area. .
需要指出的是,在本发明中,从子区中位数衬值滤波结果中提取土壤氡浓度异常区域的方式与现有技术相比没有实质性的改变,为简化表述,对该部分内容不再赘述。It should be pointed out that, in the present invention, the method of extracting the abnormal area of soil radon concentration from the filtering result of the median value of the sub-area has no substantial change compared with the prior art. Repeat.
基于上述提取土壤氡浓度异常区域的方法还可以提出一种基于信息熵-子区中位数衬值滤波提取土壤氡浓度异常区域的系统,该系统包括:Based on the above-mentioned method for extracting abnormal areas of soil radon concentration, a system for extracting abnormal areas of soil radon concentration based on information entropy-sub-area median contrast value filtering can also be proposed, and the system includes:
参数确定单元,按下式(1)计算子区中位数衬值滤波中大、小窗口的参数比值,并根据计算出的参数比值确定大、小窗口的参数:The parameter determination unit calculates the parameter ratio of the large and small windows in the median contrast filter of the sub-region as follows, and determines the parameters of the large and small windows according to the calculated parameter ratio:
S(m/n)|min=∫m/n(-1)×P(i)×lnP(i) (1);S(m/n)| min =∫ m/n (-1)×P(i)×lnP(i) (1);
式(1)中,m为大窗口的参数值,n为小窗口的参数值,m/n为大、小窗口的参数比,S(m/n)|min为大、小窗口的参数比为m/n时所对应的子区中位数衬值滤波图像的最小信息熵值,P(i)为大窗口参数值为m、小窗口参数值为n的条件下所对应的子区中位数衬值滤波图像灰度化后第i级灰度出现的概率。In formula (1), m is the parameter value of the large window, n is the parameter value of the small window, m/n is the parameter ratio of the large and small windows, and S(m/n)| min is the parameter ratio of the large and small windows. is the minimum information entropy value of the median-contrast filter image of the corresponding sub-area when m/n, P(i) is the corresponding sub-area under the condition that the large window parameter value is m and the small window parameter value is n The probability of the occurrence of the i-th grayscale after the grayscale of the bit-contrast filter image.
图像处理单元,根据所述参数确定单元所确定的大、小窗口参数,用与参数相对应的大、小窗口对土壤氡浓度等值线图作子区中位数衬值滤波处理并输出子区中位数衬值滤波结果。The image processing unit, according to the large and small window parameters determined by the parameter determination unit, uses the large and small windows corresponding to the parameters to filter the subregion median contrast value of the soil radon concentration contour map and output the subregion. District median contrast filter results.
异常区域提取单元,根据图像处理单元输出的子区中位数衬值滤波结果,从中提取对应的土壤氡浓度异常区域。The abnormal area extraction unit extracts the corresponding abnormal area of soil radon concentration according to the filtering result of the median contrast value of the sub-area output by the image processing unit.
图像处理单元还用于采用参数值为m的大窗口、参数值为n的小窗口对土壤氡浓度等值线图作子区中位数衬值滤波处理并对相应的子区中位数衬值滤波结果作灰度化处理。The image processing unit is also used to filter the sub-area median contrast value of the soil radon concentration contour map by using a large window with a parameter value of m and a small window with a parameter value of n, and perform the corresponding sub-area median contrast value filtering processing. The value filtering result is grayscaled.
参数确定单元还用于在图像处理单元对子区中位数衬值滤波结果作灰度化处理后确定该子区中位数衬值滤波图像灰度化后第i级灰度出现的概率。The parameter determining unit is further configured to determine the probability of the occurrence of the i-th grayscale of the sub-region median contrast value filtered image after grayscale processing of the subregion median contrast value filtering result by the image processing unit.
参数确定单元按下式确定对应的子区中位数衬值滤波图像灰度化后第i级灰度出现概率:The parameter determination unit determines the probability of occurrence of the i-th grayscale after the grayscale of the corresponding sub-area median-contrast filter image is as follows:
P(i)=N(i)/N (2);P(i)=N(i)/N(2);
式(2)中,P(i)为该子区中位数衬值滤波图像灰度化后第i级灰度出现的概率;N(i)为该子区中位数衬值滤波图像灰度化后第i级灰度出现的总次数;N为该子区中位数衬值滤波图像灰度化后像素总数目。In formula (2), P(i) is the probability of the occurrence of the i-th grayscale of the median-contrast filter image of the sub-region after graying; N(i) is the gray-scale of the median-contrast filter image of the sub-region. The total number of occurrences of the i-th grayscale after the grayscale; N is the total number of pixels after grayscale of the median-contrast filter image in this sub-region.
本领域技术人员应当明白,上述基于信息熵-子区中位数衬值滤波用于实现提取土壤氡浓度异常区域的系统可以封装在一个计算机软件系统中并储存于存储装置中由计算装置来执行,本发明并不限于特定的硬件和软件的结合。Those skilled in the art should understand that the above-mentioned system for extracting abnormal areas of soil radon concentration based on information entropy-sub-area median contrast value filtering can be packaged in a computer software system and stored in a storage device for execution by a computing device , the present invention is not limited to a specific combination of hardware and software.
本发明将整个测量区域看作一个提供多元信息的信息源,将信息熵与子区中位数衬值滤波图像处理技术相结合,构建了子区中位数衬值滤波中大、小窗口比值与图像熵的函数关系模型,通过函数关系模型能够计算得出子区中位数衬值滤波中最佳大、小窗口的参数比值,并由此确定后续子区中位数衬值滤波处理时的大、小窗口参数值,该方法避免了人为经验因素导致子区中位数衬值滤波数据处理结果存在不确定性的问题,保证了子区中位数衬值滤波对异常区域识别的准确度,能为圈定铀矿勘探位置和矿化范围提供有价值的参考依据。The invention regards the entire measurement area as an information source providing multiple information, and combines the information entropy with the sub-area median contrast filtering image processing technology to construct a large- and small-window ratio in the sub-area median contrast filtering. The functional relationship model with the image entropy, the parameter ratio of the optimal large and small windows in the median contrast value filtering of the sub-region can be calculated through the functional relationship model, and the subsequent sub-region median contrast value filtering processing can be determined. The method avoids the problem of uncertainty in the data processing results of the median contrast value filtering in the subregion caused by human experience factors, and ensures the accuracy of the identification of abnormal regions by the median contrast value filtering in the subregion. It can provide a valuable reference for delineating the exploration location and mineralization range of uranium deposits.
上述实施例为本发明较佳的实现方案,除此之外,本发明还可以其它方式实现,在不脱离本技术方案构思的前提下任何显而易见的替换均在本发明的保护范围之内。The above-mentioned embodiment is a preferred implementation scheme of the present invention. In addition, the present invention can also be implemented in other ways, and any obvious replacements are within the protection scope of the present invention without departing from the concept of the technical solution.
为了让本领域普通技术人员更方便地理解本发明相对于现有技术的改进之处,本发明的一些附图和描述已经被简化,并且为了清楚起见,本申请文件还省略了一些其它元素,本领域普通技术人员应该意识到这些省略的元素也可构成本发明的内容。In order to make it easier for those skilled in the art to understand the improvements of the present invention relative to the prior art, some drawings and descriptions of the present invention have been simplified, and for the sake of clarity, some other elements are also omitted in this application document, One of ordinary skill in the art would realize that these omitted elements may also constitute the subject matter of the present invention.
Claims (9)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110147521.XA CN113010996B (en) | 2021-02-03 | 2021-02-03 | Method and system for extracting radon concentration abnormity based on entropy-median filtering in sub-region |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110147521.XA CN113010996B (en) | 2021-02-03 | 2021-02-03 | Method and system for extracting radon concentration abnormity based on entropy-median filtering in sub-region |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113010996A CN113010996A (en) | 2021-06-22 |
CN113010996B true CN113010996B (en) | 2022-07-22 |
Family
ID=76385219
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110147521.XA Active CN113010996B (en) | 2021-02-03 | 2021-02-03 | Method and system for extracting radon concentration abnormity based on entropy-median filtering in sub-region |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113010996B (en) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8879813B1 (en) * | 2013-10-22 | 2014-11-04 | Eyenuk, Inc. | Systems and methods for automated interest region detection in retinal images |
CN111681124A (en) * | 2020-07-06 | 2020-09-18 | 核工业北京地质研究院 | A three-dimensional radon anomaly identification method and system for deep sandstone-type uranium mineralization information |
CN111931132A (en) * | 2020-09-27 | 2020-11-13 | 浙江浙能技术研究院有限公司 | Flue gas NOX concentration abnormal value processing method based on standard deviation and median absolute deviation |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10311571B2 (en) * | 2014-10-17 | 2019-06-04 | Stichting Maastricht Radiation Oncology “Maastro-Clinic” | Image analysis method supporting illness development prediction for a neoplasm in a human or animal body |
-
2021
- 2021-02-03 CN CN202110147521.XA patent/CN113010996B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8879813B1 (en) * | 2013-10-22 | 2014-11-04 | Eyenuk, Inc. | Systems and methods for automated interest region detection in retinal images |
CN111681124A (en) * | 2020-07-06 | 2020-09-18 | 核工业北京地质研究院 | A three-dimensional radon anomaly identification method and system for deep sandstone-type uranium mineralization information |
CN111931132A (en) * | 2020-09-27 | 2020-11-13 | 浙江浙能技术研究院有限公司 | Flue gas NOX concentration abnormal value processing method based on standard deviation and median absolute deviation |
Non-Patent Citations (7)
Title |
---|
A system and analytical method for testing the electrical parameters of soil;Yue Han 等;《Transportation Safety and Environment》;20200930;57-71 * |
地形对砂岩型铀矿氡气测量的干扰作用及其修正方法;赵宁博等;《物探与化探》;20170815(第04期);667-671 * |
地球化学异常提取的自适应衬值滤波法;金俊杰等;《物探与化探》;20110815(第04期);526-531 * |
子区中位数衬值滤波法及弱小异常识别;史长义等;《物探与化探》;19990820(第04期);250-257 * |
子区中位数衬值滤波法在川西斑岩型铜矿区地球化学异常的筛选与查证中的应用;费光春等;《物探与化探》;20080215(第01期);66-69 * |
子区中位数衬值滤波法对弱、小异常识别的应用:以玉镇山林场锌异常提取为例;陈馥等;《吉林地质》;20170930;51-54 * |
甘肃省花海盆地土壤氡分布特征及控制因素研究;雷波等;《西北地质》;20200905(第03期);106-112 * |
Also Published As
Publication number | Publication date |
---|---|
CN113010996A (en) | 2021-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107239659B (en) | Method for calculating soil erosion based on improved K-RUSLE model and soil forming rate | |
Pollard et al. | Harmonising topographic & remotely sensed datasets, a reference dataset for shoreline and beach change analysis | |
CN104915640B (en) | A Method for Automatic Recognition and Parameter Extraction of Structural Surfaces in Borehole Images | |
CN108982522B (en) | Method and apparatus for detecting pipe defects | |
CN108520070B (en) | Method and device for screening interest points of electronic map | |
CN107729293B (en) | A Geospatial Anomaly Detection Method Based on Multivariate Adaptive Regression | |
Kienholz et al. | A new semi-automatic approach for dividing glacier complexes into individual glaciers | |
CN108230365B (en) | SAR image change detection method based on multi-source difference image content fusion | |
CN107392887B (en) | A Change Detection Method of Heterogeneous Remote Sensing Image Based on Homogeneous Pixel Transformation | |
CN117057499B (en) | Landslide sensitivity evaluation method, system and equipment in mountain gorge valley region | |
CN111783784A (en) | Method and device for detecting building cavity, electronic equipment and storage medium | |
US20100254219A1 (en) | Detection of Features in Seismic Images | |
CN116895019A (en) | Remote sensing image change detection method and system based on dynamic weighted cross entropy loss | |
Leonard et al. | Measuring channel planform change from image time series: A generalizable, spatially distributed, probabilistic method for quantifying uncertainty | |
CN113010996B (en) | Method and system for extracting radon concentration abnormity based on entropy-median filtering in sub-region | |
CN108986079B (en) | Crack extraction and optimization method for complex fracture network | |
CN109254320A (en) | Seismic properties optimization and sand body Overlay District prediction technique based on image processing method | |
CN118410631A (en) | Spatial prediction of element content in black rock series high background soil and division of impact area | |
CN116561563B (en) | Slope displacement prediction method and related device based on residual prediction model | |
Wei et al. | Susceptibility assessment of earth fissure related to groundwater extraction using machine learning methods combined with weights of evidence | |
CN115424131B (en) | Cloud detection optimal threshold selection method, cloud detection method and cloud detection system | |
CN114152302B (en) | Lake water storage variation estimation method, system, electronic device and medium | |
Liu et al. | Landslide displacement prediction model based on multisource monitoring data fusion | |
CN106778515B (en) | Automatic identification method for axial magnetic flux leakage array signal of flange | |
CN109800625A (en) | A kind of form frame-line detection method based on region growing |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |