Curvature-based signal segmentation method for solar-blind ultraviolet photodetectors
12517978 ยท 2026-01-06
Assignee
Inventors
- Jiufei Luo (Chongqing, CN)
- Wei Liu (Chongqing, CN)
- Sheng Lu (Chongqing, CN)
- Song Feng (Chongqing, CN)
- Haiqing Li (Chongqing, CN)
- Hai Zeng (Chongqing, CN)
- Hongcheng Li (Chongqing, CN)
- Xinyu Wang (Chongqing, CN)
- Baojun Yang (Chongqing, CN)
- Denghua Zheng (Chongqing, CN)
Cpc classification
G06F17/16
PHYSICS
International classification
G06F17/16
PHYSICS
Abstract
The present disclosure relates to a curvature-based signal segmentation method for solar-blind ultraviolet (UV) photodetectors, including the steps of: collecting solar-blind UV radiation using solar-blind UV photodetectors according to a preset sampling time and generating the same into sampled signals; performing noise reduction on the sampled signals to generate filtered signals using a fractional calculus digital filter; calculating estimated curvatures of the filtered signals using a method of circle fitting by discrete points, and numerically optimizing the estimated curvatures using an exponential function to obtain curvature vectors of the filtered signals; and calculating weight coefficients and weighting the filtered signals using a moving average filter function and a binarization method according to the curvature vectors.
Claims
1. A curvature-based signal segmentation method for solar-blind ultraviolet (UV) photodetectors, comprising steps performed by a processor: S1: instructing a solar-blind UV photodetector to collect solar-blind UV radiation according to a preset sampling time, and generate sampled signals from the solar-blind UV radiation; S2: instructing a fractional calculus digital filter to receive the sampled signals, and perform noise reduction on the sampled signals to generate filtered signals; S3: in a memory, calculating estimated curvatures of the filtered signals using a method of circle fitting by discrete points, and numerically optimizing the estimated curvatures using an exponential function to obtain curvature vectors of the filtered signals; S4: in the memory, weighting the filtered signals using a moving average filter function and a binarization method according to the curvature vectors of the filtered signals, to obtain a plurality of non-zero discrete time-domain segments; and excluding interference segments from each non-zero discrete time-domain segment based on a mean value and width thereof, to obtain candidate segmentation point signals; and S5: in the memory, setting determination indicators according to the candidate segmentation point signals to determine upper and lower segment boundary points, and segmenting and rectifying the filtered signals by means of the upper and lower segment boundary points to decode signals from solar-blind UV photodetectors.
2. The curvature-based signal segmentation method for solar-blind UV photodetectors according to claim 1, wherein step S3 comprises: S31: defining the filtered signal S={S.sub.0, S.sub.1, . . . , S.sub.n-1}, where n is a length of the filtered signal, and framing the filtered signal according to a preset frame length w and a frame shift r to obtain S={SW.sub.0, SW.sub.1, . . . , SW.sub.i, . . . , SW.sub.(n-w-1)r}, where SW.sub.i represents an i.sub.th frame signal of the filtered signal, and (nw1)r represents the number of frames in the filtered signal; S32: dividing, by means of a midpoint of the length M=<w/2>, the i.sub.th frame signal SW.sub.i={SW.sub.i1, SW.sub.i2, . . . , SW.sub.i<w/2>, . . . , SW.sub.iw} of the filtered signal into left and right segmented frames {SW.sub.i1, SW.sub.i2, . . . , SW.sub.i|<w/2>-1|} and {SW.sub.i|<w/2>+1|, . . . , SW.sub.iw}, where < > represents rounding down; counting mid-values of signals of the left and right segmented frames, denoted as S.sub.iR and S.sub.iL, respectively; and taking (S.sub.iR, M/2)=(x.sub.1, y.sub.1), (S.sub.iM, M)=(x.sub.2, y.sub.2), and (S.sub.iL, 3M/2)=(x.sub.3, y.sub.3) as three discrete points for the i.sub.th frame signal; S33: constructing, based on a circle expansion formula x.sub.2+y.sub.2+ax+by+c=0, an error equation W.sub.1i=.sub.j=1.sup.3(x.sub.j.sup.2+y.sub.j.sup.2+ax.sub.j+by.sub.j+c).sup.2 for a circle of the i.sub.th frame signal and a position error equation W.sub.2i=E.sub.j=1.sup.3(y.sub.j(dx.sub.j+e)).sup.2 for a straight line and a center of the circle using the three discrete points (x.sub.j, y.sub.j)(j=1, 2, 3) of the i.sub.th frame signal of the filtered signal; and calculating partial derivatives 0 for the unknown parameters a, b, c, d, and e to obtain approximate estimated values .sub.l, {circumflex over (b)}.sub.l, .sub.l, {circumflex over (d)}.sub.l, and .sub.l for this frame signal; S34: constructing a fitting circle based on .sub.l, {circumflex over (b)}.sub.l, .sub.l, {circumflex over (d)}.sub.l, and .sub.l, and calculating an estimated curvature for the i.sub.th frame signal of the filtered signal:
3. The curvature-based signal segmentation method for solar-blind UV photodetectors according to claim 1, wherein step S4 comprises: S41: smoothing the curvature vector R using the moving average filter function, setting a threshold =Amax(R) to traverse each element in the smoothed result, and setting weights of elements greater than the threshold as 1 and weights of elements less than or equal to the threshold as 0, to generate a weight coefficient vector =[.sub.0, .sub.1, . . . , .sub.n-1] corresponding to the smoothed result, A being an adjustable parameter; S42: weighing the filtered signal according to the weight coefficient vector , i.e., Q=S, to obtain the plurality of non-zero discrete time-domain segments; and S43: extracting the mean value and width of each non-zero discrete time-domain segment sequentially, discarding the non-zero discrete time-domain segment if the mean value thereof is less than 30% of the threshold and the width is below 150, and otherwise, retaining the non-zero discrete time-domain segment, to obtain the candidate segmentation point signal.
4. The curvature-based signal segmentation method for solar-blind UV photodetectors according to claim 3, wherein A ranges from 0.1-0.3.
5. The curvature-based signal segmentation method for solar-blind UV photodetectors according to claim 1, wherein the setting determination indicators according to the candidate segmentation point signals to determine upper and lower segment boundary points comprises: determining the lower segment boundary point according to the candidate segmentation point signal by setting the determination indicator, segmenting the filtered signal using the lower segment boundary point, estimating the curvature vector for each segment of signal using the method of circle fitting by discrete points, and finally taking a horizontal coordinate corresponding to the maximum curvature as the upper segment boundary point.
6. The curvature-based signal segmentation method for solar-blind UV photodetectors according to claim 5, wherein the determining the lower segment boundary point according to the candidate segmentation point signal by setting the determination indicator comprises: S51: extracting horizontal coordinates corresponding to the maximum values in various non-zero discrete time-domain segments of the candidate segmentation point signal to form a one-dimensional vector ={.sub.0, .sub.1, . . . .sub.k, . . . , .sub.K-1}, where K represents the number of non-zero discrete time-domain segments in the candidate segmentation point signal, performing a first-order forward difference operation on the one-dimensional vector ={.sub.0, .sub.1, . . . .sub.k, . . . , .sub.K-1}, and calculating a mean value to obtain a neighborhood width used for determination of segment boundary points; S52: taking a horizontal coordinate .sub.k corresponding to the maximum value in each non-zero discrete time-domain segment as a center point, and intercepting sample data R.sub..sub.
7. The curvature-based signal segmentation method for solar-blind UV photodetectors according to claim 6, wherein the energy difference is calculated by the formula .sub.k=|R.sub..sub.
8. The curvature-based signal segmentation method for solar-blind UV photodetectors according to claim 6, wherein the parameter t is an adjustable parameter, ranging from 0.1-0.3.
9. The curvature-based signal segmentation method for solar-blind UV photodetectors according to claim 1, wherein the segmenting and rectifying the filtered signal comprises: rectifying the signal through the upper and lower segment boundary points; setting all filtered signal data from the lower segment boundary point to the first upper segment boundary point on the left as 0, and setting all the data from the lower segment boundary point to the first upper segment boundary point on the right as 1 according to the characteristics of the filtered signal, with each lower segment boundary point as a center; and completing the automatic rectification of the filtered signal.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
(1)
(2)
(3)
(4)
(5)
(6)
DETAILED DESCRIPTION
(7) The embodiments of the present disclosure are described by means of particular specific examples below, and other advantages and effects of the present disclosure can be readily understood by those skilled in the art from the contents disclosed in the specification. The present disclosure may also be implemented or applied in different specific embodiments, and the details in the specification may be modified or altered in various ways based on different viewpoints and applications without departing from the spirit of the present disclosure. It is to be noted that the illustrations provided in the following examples are solely intended to schematically explain the basic concept of the present disclosure. The following examples and features thereof may be combined with each other without conflict.
(8) The accompanying drawings are solely for illustrative purposes, representing schematic diagrams rather than physical images, and are not to be construed as limiting the present disclosure. To better illustrate the examples of the present disclosure, some components in the drawings may be omitted, enlarged, or reduced in size, and do not represent the actual dimensions of the product. For those skilled in the art, it is understandable that some well-known structures and descriptions thereof in the drawings may be omitted.
(9) In the accompanying drawings of the examples of the present disclosure, the identical or similar labels correspond to the identical or similar components. In the description of the present disclosure, it is to be understood that, the orientation or state relations indicated by the terms up, down, left, right, front, rear, etc., are based on those shown in the accompanying drawings and merely for the ease of describing the present disclosure and simplifying the description, rather than indicating or implying that the device or element referred to must be in a specific orientation or constructed and operated in a specific orientation, and therefore cannot be interpreted as limiting the present disclosure. For those of ordinary skill in the art, the specific meaning of the above terms may be understood according to specific circumstances.
(10) As shown in
(11) In an example, two sets of data are collected using the solar-blind UV photodetector, with a sampling frequency of fs=200000, a sampling time of 5 s, and a signal length of n=10000000. The time response of the detector under two sets of different wavelengths of light is shown in
(12) In an example, to extract signal segmentation points, a differential order of 0.9 is set, and a fractional calculus digital lowpass filter is employed to perform noise reduction on the sampled signals. The results are shown in Table 3. It can be seen that the signal to noise ratio is improved effectively after filtering. S3: estimated curvatures of the filtered signals are calculated using a method of circle fitting by discrete points, and the estimated curvatures are numerically optimized using an exponential function to obtain curvature vectors of the filtered signals.
(13) In an example, to realize signal segmentation, the method of circle fitting by discrete points is employed to determine the estimated curvatures, which are numerically optimized. The specific calculation steps are as follows. S31: the filtered signal is defined as S={S.sub.0, S.sub.1, . . . , S.sub.n-1}, where n is a length of the filtered signal, and the filtered signal is framed according to a preset frame length w and a frame shift r to obtain S={SW.sub.0, SW.sub.1, . . . , SW.sub.i, . . . , SW.sub.(n-w-1)r}, where SW.sub.i represents an i.sub.th frame signal of the filtered signal, and (nw1)r represents the number of frames in the filtered signal. In some examples, w=1500. S32: by means of a midpoint of the length M=
(14)
for a circle of the i.sub.th frame signal and a position error equation
(15)
for the straight line and a center of the circle are constructed using the three discrete points (x.sub.j, y.sub.j) (j=1, 2, 3) of the i.sub.th frame signal of the filtered signal. Partial derivatives 0 are calculated for the unknown parameters a, b, c, d, and e to obtain approximate estimated values ,
,
,
, and
for this frame signal.
(16) Preferably, the estimation formulas for ,
,
,
, and
are as follows.
(17) ,
,
,
, and
, and an estimated curvature is calculated for the i.sub.th frame signal of the filtered signal:
(18) represents the estimated curvature for the i.sub.th frame signal of the filtered signal, .sub.i is a parameter reflecting the bending direction of a curve of the i.sub.th frame signal, and sign() represents a sign function. A frame shift r=1 is selected, and an estimated curvature sequence ={.sub.0, .sub.1, . . . , .sub.n-1} of the filtered signal is obtained through periodic extension.
(19) Referring to
(20) In an example, by setting a reasonable threshold to generate a weighting coefficient, the filtered signal is weighted to obtain a plurality of non-zero discrete time-domain segments containing interference and segment boundary points. To reduce the burden on the subsequent determination of segment boundary points, the interference can be excluded by extracting the mean value and width of each weighted segment. The specific steps are shown as follows. S41: the curvature vector R is smoothed using the moving average filter function, and a threshold =Amax(R) is set to traverse each element in the smoothed result. Weights of elements greater than the threshold are set as 1 and weights of elements less than or equal to the threshold are set as 0, to generate a weight coefficient vector =[.sub.0, .sub.1, . . . , .sub.n-1] corresponding to the smoothed result.
(21) Preferably, the parameter A is an adjustable parameter, ranging from 0.1-0.3. In this example, A=0.2. S42: the filtered signal is weighted according to the weight coefficient vector , i.e., Q=S, to obtain the plurality of non-zero discrete time-domain segments. S43: the mean value and width of each non-zero discrete time-domain segment are extracted sequentially. The non-zero discrete time-domain segment is discarded if the mean value thereof is less than 30% of the threshold and the width is below 150, and otherwise, the non-zero discrete time-domain segment is retained. S5: determination indicators are set according to the candidate segmentation point signal to determine upper and lower segment boundary points, and the filtered signal is segmented and rectified by means of the upper and lower segment boundary points to obtain signals from solar-blind UV photodetectors.
(22) Preferably, setting determination indicators according to the candidate segmentation point signal to determine upper and lower segment boundary points includes that: the lower segment boundary point is determined according to the candidate segmentation point signal by setting the determination indicator, the filtered signal is segmented using the lower segment boundary point, the curvature vector for each segment of signal is estimated using the method of circle fitting by discrete points, and finally a horizontal coordinate corresponding to the maximum curvature is taken as the upper segment boundary point.
(23) In an example, through the pre-exclusion of the interference segment described in step S43, an inflection point interval containing both upper and lower segment boundary points can be obtained. Since some upper segment boundary points may be lost during the weighting of the filtered signal and larger non-boundary points may still exist in the numerical optimization results, in the algorithm, determination indicators are selected to exclude the upper segment boundary points and interference segments, while simultaneously recording all lower segment boundary points. Based on this, after segmenting the filtered signal using the lower segment boundary point, the discrete curvature estimation method is further utilized to determine the upper segment boundary point, thereby achieving signal rectification. The above content includes the following specific steps. S51: horizontal coordinates corresponding to the maximum values in various non-zero discrete time-domain segments of the candidate segmentation point signals are extracted to form a one-dimensional vector ={.sub.0, .sub.1, . . . , .sub.k, . . . , .sub.K-1}, where K represents the number of non-zero discrete time-domain segments in the candidate segmentation point signal. A first-order forward difference operation is performed on the one-dimensional vector ={.sub.0, .sub.1, . . . .sub.k, . . . , .sub.K-1}, and a mean value is calculated to obtain a neighborhood width used for determination of segment boundary points. S52: according to the waveform characteristics of the signal shown in
(24) Preferably, the energy difference is calculated by the formula E.sub..sub.
(25) Preferably, the parameter t is an adjustable parameter, ranging from 0.1-0.3. In this example, t=0.15.
(26) Preferably, after segmenting the filtered signal using the lower segment boundary point, the method of circle fitting circles by discrete points is employed to estimate the curvature vector of each signal segment. Finally, the horizontal coordinate corresponding to the maximum curvature is taken as the upper segment boundary point, and the upper and lower boundary segmentation points are utilized to achieve signal segmentation. The result is shown in
(27) Preferably, the segmenting and rectifying the filtered signal includes that: the signal is rectified through the upper and lower segment boundary points; all filtered signal data from the lower segment boundary point to the first upper segment boundary point on the left are set as 0, and all the data from the lower segment boundary point to the first upper segment boundary point on the right are set as 1 according to the characteristics of the filtered signal, with each lower segment boundary point as a center, so that the automatic rectification of the filtered signal is completed. The signal is rectified by assigning high and low levels, the result is as shown in
(28) Finally, it is to be noted that the example described above is merely used to illustrate rather than limiting the technical solutions of the present disclosure. Although the present disclosure is described in detail by reference to the preferred example, it is to be understood by those ordinary skilled in the art that modifications or equivalent replacements can be made to the technical solutions of the present disclosure without departing from the principle and scope of the technical solutions of the present disclosure, and these modifications or equivalent replacements are included in the scope of claims of the present disclosure.