Nothing Special   »   [go: up one dir, main page]

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 PDF

Info

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
Application number
CN202110147521.XA
Other languages
Chinese (zh)
Other versions
CN113010996A (en
Inventor
雷波
雷林
罗才武
康虔
杨蓉
罗润
谢超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of South China
Original Assignee
University of South China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of South China filed Critical University of South China
Priority to CN202110147521.XA priority Critical patent/CN113010996B/en
Publication of CN113010996A publication Critical patent/CN113010996A/en
Application granted granted Critical
Publication of CN113010996B publication Critical patent/CN113010996B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design 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

A method for extracting a soil radon concentration abnormal region based on information entropy-subregion median liner value filtering relates to the technical field of uranium mine resource exploration, combines information entropy and subregion median liner value filtering image processing technology, constructs a functional relation model, determines parameter values of large and small windows according to the functional relation model, uses large and small windows corresponding to the parameter values to perform subregion median liner value filtering processing on a soil radon concentration contour map, and extracts a corresponding soil radon concentration abnormal region according to a subregion median liner value filtering result; the invention also provides a system for extracting the soil radon concentration abnormal region based on the information entropy-median liner value filtering in the subarea. The method avoids the problem that the result of processing the subzone median value filtering data is uncertain due to human experience factors, ensures the accuracy of the subzone median value filtering on identifying abnormal areas, and can provide valuable reference basis for delineating the exploration position and the mineralization range of the uranium mine.

Description

基于熵-子区中位数滤波提取氡浓度异常的方法和系统Method and system for extracting radon concentration anomalies based on entropy-subregion median filtering

技术领域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)

1. The method for extracting the soil radon concentration abnormal region based on the information entropy-median liner value filtering in the sub-region is characterized in that the parameter ratio of a large window to a small window in the median liner value filtering in the sub-region is determined according to the following formula:
S(m/n)|min=∫m/n(-1)×P(i)×lnP(i) (1);
in the formula (1), m is the parameter value of a large window, n is the parameter value of a small window, m/n is the parameter ratio of the large window and the small window, S (m/n) isminP (i) is the minimum information entropy value of a median liner value filtering image in a sub-region corresponding to the condition that the parameter ratio of a large window and a small window is m/n, and P (i) is the probability of the i-th level gray level of the corresponding median liner value filtering image in the sub-region after graying under the condition that the parameter value of the large window is m and the parameter value of the small window is n;
determining parameter values of a large window and a small window according to the parameter ratio determined by the formula (1), using the large window and the small window corresponding to the parameter values to conduct filtering processing on the soil radon concentration contour map to obtain a sub-region median liner value, and extracting a corresponding soil radon concentration abnormal region according to the filtering result of the sub-region median liner value.
2. The method for extracting the radon concentration abnormal region in the soil based on the information entropy-median liner value filtering in the subarea according to claim 1, which is characterized in that: the method also comprises the step of measuring the soil radon concentration and forming a soil radon concentration contour map.
3. The method for extracting the radon concentration abnormal region in the soil based on the information entropy-median liner value filtering in the subarea according to claim 1, which is characterized in that: and the method also comprises the steps of performing median lining value filtering processing on the soil radon concentration contour map in the subarea and performing graying processing on the median lining value filtering result in the subarea by adopting a large window with the parameter value of m and a small window with the parameter value of n.
4. The method for extracting radon concentration abnormal regions in soil based on information entropy-median liner value filtering in subareas according to claim 3, which is characterized by comprising the following steps of: after the corresponding sub-area median substrate value filtering result is grayed, the method also comprises the step of determining the occurrence probability of the i-th level gray level after the sub-area median substrate value filtering image is grayed.
5. The method for extracting the radon concentration abnormal region in the soil based on the information entropy-median liner value filtering in the subarea as claimed in claim 4, wherein: determining the i-th level gray degree occurrence probability of the grayed image of the corresponding median liner value filtering image in the subarea according to the following formula:
P(i)=N(i)/N (2);
in the formula (2), P (i) is the probability of the i-th level gray level after the median-value filtering image in the subarea is grayed; n (i) is the total number of times of the i-th level gray level after the median lining value filtering image graying in the subarea; and N is the total number of pixels of the grayed median contrast value filtering image in the sub-area.
6. System for extracting soil radon concentration abnormal region based on information entropy-median lining value filtering in subregion, its characterized in that includes:
the parameter determining unit calculates the parameter ratio of the large window and the small window in the median liner value filtering in the subarea according to the following formula (1), and determines the parameters of the large window and the small window according to the calculated parameter ratio:
S(m/n)|min=∫m/n(-1)×P(i)×lnP(i) (1);
in the 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 window and the small window, S (m/n) & gtYminP (i) is the minimum information entropy value of a median liner value filtering image in a corresponding sub-area when the parameter ratio of a large window and a small window is m/n, and P (i) is the probability of the i-th level gray level after the corresponding median liner value filtering image in the sub-area is grayed under the condition that the parameter of the large window is m and the parameter of the small window is n;
the image processing unit is used for making a sub-region median liner value filtering treatment on the soil radon concentration contour map by using the large window and the small window corresponding to the parameters according to the large window parameter and the small window parameter determined by the parameter determination unit and outputting a sub-region median liner value filtering result;
and the abnormal region extraction unit is used for extracting a corresponding soil radon concentration abnormal region from the filtering result of the median lining value in the subarea output by the image processing unit.
7. The system for extracting radon concentration abnormal regions in soil based on the information entropy-median liner value filtering in the subareas as claimed in claim 6, wherein: the image processing unit is also used for carrying out sub-region median lining value filtering processing on the soil radon concentration contour map by adopting a large window with a parameter value of m and a small window with a parameter value of n and carrying out graying processing on corresponding sub-region median lining value filtering results.
8. The system for extracting radon concentration abnormal regions in soil based on the information entropy-median liner value filtering in the subareas as claimed in claim 7, wherein: the parameter determining unit is also used for determining the probability of the i-th gray level after the gray level of the filtered image of the median substrate value in the sub-area is grayed after the image processing unit grays the filtered result of the median substrate value in the corresponding sub-area.
9. The system for extracting radon concentration anomaly region in soil based on information entropy-median liner value filtering in subareas according to claim 8, characterized by: the parameter determining unit determines the i-th level gray degree occurrence probability of the grayed median liner value filtering image in the corresponding sub-area according to the following formula:
Figure FDA0002931179970000031
in the formula (2), P (i) is the probability of the i-th level gray level after the median-value filtering image in the subarea is grayed; n (i) is the total times of the ith gray level after the median lining value filtering image graying in the subarea appears; and N is the total number of pixels of the grayed median contrast value filtering image in the sub-area.
CN202110147521.XA 2021-02-03 2021-02-03 Method and system for extracting radon concentration abnormity based on entropy-median filtering in sub-region Active CN113010996B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (3)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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