[go: up one dir, main page]

CN109815857B - A Clustering Method for Ionizing Radiation Time Series - Google Patents

A Clustering Method for Ionizing Radiation Time Series Download PDF

Info

Publication number
CN109815857B
CN109815857B CN201910018760.8A CN201910018760A CN109815857B CN 109815857 B CN109815857 B CN 109815857B CN 201910018760 A CN201910018760 A CN 201910018760A CN 109815857 B CN109815857 B CN 109815857B
Authority
CN
China
Prior art keywords
gaussian mixture
mixture model
peak
branch
class
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
CN201910018760.8A
Other languages
Chinese (zh)
Other versions
CN109815857A (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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Priority to CN201910018760.8A priority Critical patent/CN109815857B/en
Publication of CN109815857A publication Critical patent/CN109815857A/en
Application granted granted Critical
Publication of CN109815857B publication Critical patent/CN109815857B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Complex Calculations (AREA)
  • Image Analysis (AREA)
  • Measurement Of Radiation (AREA)

Abstract

一种电离辐射时间序列的聚类方法,包括以下步骤:1)对输入的电离辐射序列数据进行裁剪,截取峰值波段;2)对检测出的波峰段进行Fourier变换;3)根据频域系数计算出频谱;4)这些波峰段序列对应的频谱段都是等长的;对这些频谱段进行k均值聚类,迭代到收敛为止;5)令k均值聚类得到的每类的矢量中心、以及方差矢量作为该k支高斯混合模型的每支的均值和方差;6)使用期望最大算法对高斯混合模型进行迭代,迭代到收敛为止;7)对每一个频谱段使用收敛后的高斯混合模型计算其属于i类别的条件概率,取概率值最大的支属为其类别号,这样就定出了该频谱段对应的输入的辐射时间序列的波峰段的类别。本发明灵活性较好、准确性较高。

Figure 201910018760

A clustering method for ionizing radiation time series, comprising the following steps: 1) trimming input ionizing radiation sequence data, and intercepting peak wave bands; 2) performing Fourier transform on detected wave peak sections; 3) calculating according to frequency domain coefficients 4) The spectrum segments corresponding to these peak segment sequences are all of equal length; perform k-means clustering on these spectrum segments, and iterate until convergence; 5) Make the vector center of each class obtained by k-means clustering, and The variance vector is used as the mean and variance of each branch of the k Gaussian mixture model; 6) Use the expectation-maximization algorithm to iterate the Gaussian mixture model until convergence; 7) Use the converged Gaussian mixture model for each spectrum segment to calculate It belongs to the conditional probability of category i, and the branch with the largest probability value is taken as its category number, so that the category of the peak section of the input radiation time series corresponding to this spectrum section is determined. The invention has better flexibility and higher accuracy.

Figure 201910018760

Description

Clustering method of ionizing radiation time sequence
Technical Field
The invention belongs to the field of data mining and big data processing, and relates to a clustering method of an ionizing radiation time sequence.
Background
With the rise of the application of the internet of things, more and more sensor data are accumulated, and how to analyze and process the data is an important subject faced by people to mine valuable information. These data are basically time series data, i.e. there is a strict chronological order, but the data at each time sampling point have different dimensions, for example, video is multi-dimensional matrix data of each frame, and each time sampling point of the temperature control sensor is continuous one-dimensional data (i.e. scalar quantity), etc. The ionizing radiation time series is the ionizing radiation values of local peripheral areas acquired by ionizing radiation sensors deployed in real environments (water, air and land), and mainly senses radiation intensities of gamma, neutron lines and the like. Ionizing radiation, which causes disturbances of the physiological metabolic functions of the human body by damaging cells to charge them, and electromagnetic radiation, which is harmless and causes at most thermal effects, are very weak and thus negligible in many cases, are different concepts.
The ionizing radiation value is a continuous one-dimensional time sequence, and the value of the ionizing radiation value fluctuates along with the influences of events such as rainfall, snowfall, perspective construction (namely, perspective equipment such as X-ray is used for detecting industrial pipeline facilities), hospital radiotherapy (including cardiac surgery and the like), even nuclear facility leakage and the like. How to automatically cluster and identify the events of the time series fluctuation values so as to distinguish the radiation fluctuation abnormal events is a very important application topic. The traditional time series clustering requires the length of sequence segments participating in clustering to be consistent, which limits the flexibility of application, and the forced use of equal-length preprocessing can cause distortion of a dynamic change mechanism in a sequence.
The time series of ionizing radiation peaks vary in duration, for example, the rise in radiation caused by precipitation must vary with the duration of the precipitation. Therefore, equal-length clustering techniques, such as k-means clustering in the time domain, gaussian mixture model, and other classical techniques, cannot be used for clustering of radiation peaks.
Disclosure of Invention
In order to overcome the defects of poor flexibility and poor accuracy of the conventional ionizing radiation time sequence clustering mode, the invention provides a clustering method of an ionizing radiation time sequence with good flexibility and high accuracy, which does not need to introduce isometric preprocessing and adopts a Gaussian mixture model technology of time sequence segments on a frequency domain to carry out automatic clustering.
The technical scheme adopted by the invention for solving the technical problems is as follows:
a clustering method of ionizing radiation time series comprises the following steps:
1) cutting input ionizing radiation sequence data, and intercepting a peak wave band;
2) carrying out Fourier transformation on the detected peak segment, and replacing the peak segment with filling 0 if the length of the peak segment is not enough for the input length requirement of the Fourier transformation;
3) calculating a frequency spectrum according to the frequency domain coefficients, namely calculating a nonnegative value corresponding to each bin serial number, and obtaining a real part, an imaginary part square and a root opening number corresponding to each bin of the frequency domain coefficients;
4) the spectrum sections corresponding to the wave crest section sequences are all equal in length, namely the same bin number is contained; performing k-means clustering on the frequency spectrum sections until the frequency spectrum sections are converged in an iteration mode;
5) using the vector center and the variance vector of each class obtained by k-means clustering as the mean and the variance of each of the k Gaussian mixture models;
6) iterating the Gaussian mixture model by using an expected maximum algorithm until the Gaussian mixture model is converged;
7) and calculating the conditional probability of each frequency spectrum section belonging to the class i by using the converged Gaussian mixture model, wherein the class i is the ith branch of the Gaussian mixture model, and the branch with the maximum probability value belongs to the class number of the I, so that the class of the wave crest section of the input radiation time sequence corresponding to the frequency spectrum section is determined.
Further, in the step 1), a peak band is detected by using an adaptive threshold method.
The invention has the following beneficial effects: the method does not need to introduce isometric preprocessing, and adopts a Gaussian mixture model technology of time sequence segments on a frequency domain to perform automatic clustering; the flexibility is better, the accuracy is higher.
Drawings
Fig. 1 is a flowchart of a clustering method of an ionizing radiation time series.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
Referring to fig. 1, a clustering method of an ionizing radiation time series includes the following steps:
1) cutting input ionizing radiation sequence data, and intercepting a peak wave band;
2) carrying out Fourier transformation on the detected peak segment, and replacing the peak segment with filling 0 if the length of the peak segment is not enough for the input length requirement of the Fourier transformation;
3) calculating a frequency spectrum according to the frequency domain coefficients, namely calculating a nonnegative value corresponding to each bin serial number, and obtaining a real part, an imaginary part square and a root opening number corresponding to each bin of the frequency domain coefficients;
4) the spectrum sections corresponding to the wave crest section sequences are all equal in length, namely the same bin number is contained; performing k-means clustering on the frequency spectrum sections until the frequency spectrum sections are converged in an iteration mode;
5) using the vector center and the variance vector of each class obtained by k-means clustering as the mean and the variance of each of the k Gaussian mixture models;
6) iterating the Gaussian mixture model by using an expected maximum algorithm until the Gaussian mixture model is converged;
7) and calculating the conditional probability of each frequency spectrum section belonging to the class i by using the converged Gaussian mixture model, wherein the class i is the ith branch of the Gaussian mixture model, and the branch with the maximum probability value belongs to the class number of the I, so that the class of the wave crest section of the input radiation time sequence corresponding to the frequency spectrum section is determined.
Further, in the step 1), a peak band is detected by using an adaptive threshold method.

Claims (2)

1. A method for clustering a time series of ionizing radiations, characterized in that it comprises the following steps:
1) cutting input ionizing radiation sequence data, and intercepting a peak wave band;
2) carrying out Fourier transformation on the detected peak segment, and replacing the peak segment with filling 0 if the length of the peak segment is not enough for the input length requirement of the Fourier transformation;
3) calculating a frequency spectrum according to the frequency domain coefficients, namely calculating a nonnegative value corresponding to each bin serial number, and obtaining a real part, an imaginary part square and a root opening number corresponding to each bin of the frequency domain coefficients;
4) the spectrum sections corresponding to the wave crest section sequences are all equal in length, namely the same bin number is contained; performing k-means clustering on the frequency spectrum sections until the frequency spectrum sections are converged in an iteration mode;
5) enabling the vector center and the variance vector of each class obtained by k-means clustering to serve as the mean and the variance of each branch of the k-branch Gaussian mixture model;
6) iterating the Gaussian mixture model by using an expected maximum algorithm until the Gaussian mixture model is converged;
7) and calculating the conditional probability of each frequency spectrum section belonging to the class i by using the converged Gaussian mixture model, wherein the class i is the ith branch of the Gaussian mixture model, and the branch with the maximum probability value belongs to the class number of the I, so that the class of the wave crest section of the input radiation time sequence corresponding to the frequency spectrum section is determined.
2. The method for clustering temporal sequences of ionizing radiations as set forth in claim 1, wherein in the step 1), the peak sections are detected by an adaptive threshold method.
CN201910018760.8A 2019-01-09 2019-01-09 A Clustering Method for Ionizing Radiation Time Series Active CN109815857B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910018760.8A CN109815857B (en) 2019-01-09 2019-01-09 A Clustering Method for Ionizing Radiation Time Series

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910018760.8A CN109815857B (en) 2019-01-09 2019-01-09 A Clustering Method for Ionizing Radiation Time Series

Publications (2)

Publication Number Publication Date
CN109815857A CN109815857A (en) 2019-05-28
CN109815857B true CN109815857B (en) 2021-05-18

Family

ID=66603293

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910018760.8A Active CN109815857B (en) 2019-01-09 2019-01-09 A Clustering Method for Ionizing Radiation Time Series

Country Status (1)

Country Link
CN (1) CN109815857B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119399606B (en) * 2024-09-19 2025-12-19 上海大学 Method for eliminating off-axis digital holographic zero-order item

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106127220A (en) * 2016-06-01 2016-11-16 苏州大学 A kind of time series classification method and device
CN106384128A (en) * 2016-09-09 2017-02-08 西安交通大学 Method for mining time series data state correlation
CN106446251A (en) * 2016-10-14 2017-02-22 厦门大学 Continuous simulate quantitative time series data mining method
CN107463947A (en) * 2017-07-13 2017-12-12 浙江工业大学 A kind of pattern clustering method of time series

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9417700B2 (en) * 2009-05-21 2016-08-16 Edge3 Technologies Gesture recognition systems and related methods

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106127220A (en) * 2016-06-01 2016-11-16 苏州大学 A kind of time series classification method and device
CN106384128A (en) * 2016-09-09 2017-02-08 西安交通大学 Method for mining time series data state correlation
CN106446251A (en) * 2016-10-14 2017-02-22 厦门大学 Continuous simulate quantitative time series data mining method
CN107463947A (en) * 2017-07-13 2017-12-12 浙江工业大学 A kind of pattern clustering method of time series

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
时空数据挖掘研究进展;刘大有 等;《计算机研究与发展》;20121205;全文 *
时间序列的表示与分类算法综述;原继东,王志海;《计算机科学》;20150315;第42卷(第3期);全文 *

Also Published As

Publication number Publication date
CN109815857A (en) 2019-05-28

Similar Documents

Publication Publication Date Title
Liu et al. Application of deep convolutional neural networks for detecting extreme weather in climate datasets
Lee et al. Applying machine learning methods to detect convection using Geostationary Operational Environmental Satellite-16 (GOES-16) advanced baseline imager (ABI) data
Li et al. De-noising low-frequency magnetotelluric data using mathematical morphology filtering and sparse representation
Hou et al. Deep neural network denoising model based on sparse representation algorithm for ECG signal
Jiang et al. Principal component analysis for extremes and application to US precipitation
CN105205112A (en) System and method for excavating abnormal features of time series data
CN115034254B (en) Nuclide identification method based on HHT (Highway-HHT) band energy characteristics and convolutional neural network
WO2019019565A1 (en) Energy distribution characteristic based mine microseismic signal identification method
CN108742697B (en) Heart sound signal classification method and terminal equipment
CN103268607B (en) A kind of common object detection method under weak supervision condition
CN107657242A (en) A kind of mt noise identification and separation method
CN112380998A (en) Low-voltage transformer area missing data completion method based on matrix completion
CN106446081B (en) A Method of Mining Time-Series Data Correlation Based on Change Consistency
Salmane et al. A method for the automated detection of solar radio bursts in dynamic spectra
CN109815857B (en) A Clustering Method for Ionizing Radiation Time Series
Gaetan et al. A hierarchical model for the analysis of spatial rainfall extremes
CN110269605A (en) A kind of electrocardiosignal noise recognizing method based on deep neural network
CN119580049A (en) A CLIP-based multimodal approach to realize intelligent crop recognition system
Zhang et al. Combining multi-level feature extraction algorithm with residual graph convolutional neural network for partial discharge detection
CN108983301B (en) A kind of K index shows reporting method
Pander A new approach to robust, weighted signal averaging
Wang et al. Extraction of fetal ECG by logarithmic hyperbolic secant adaptive algorithm in alpha-stable noise
CN111616686B (en) Physiological signal imbalance analysis method based on Heaviside kernel function
CN119992263A (en) Ecological environment degradation assessment system, method, storage medium and electronic device
CN118633910A (en) An automatic denoising method for magnetocardiographic signals based on independent component analysis

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