Non-invasive diagnostic data acquisition methods and systems

Abstract

This invention provides a method and system for acquiring non-invasive diagnostic data, relating to the field of data acquisition technology. The method includes: acquiring raw optical acquisition data and raw contact pressure data, and performing time alignment processing; performing contact pressure change identification processing to generate contact pressure influence domain data; performing optical path disturbance localization processing on the raw optical acquisition data in the aligned data based on the contact pressure influence domain data to generate optical interference fragment data; performing optical compensation construction processing on the corresponding optical interference fragments based on the optical interference fragment data to generate domain-compensated optical data; performing domain-segmentation splicing processing between each contact pressure influence domain based on the domain-compensated optical data to generate continuously corrected optical data; and performing target data extraction processing based on the continuously corrected optical data to generate non-invasive diagnostic data. This invention improves the autonomy and accuracy of non-invasive diagnostic data acquisition.

Description

Technical Field

[0001] This invention relates to the field of data acquisition technology, and in particular to a method and system for acquiring non-invasive diagnostic data. Background Technology

[0002] Existing non-invasive diagnostic data acquisition technologies primarily utilize optical sensing (such as near-infrared spectroscopy NIRS), electrophysiological sensing (such as photoplethysmography, PPG), or millimeter-wave radar-based vital sign detection to perform non-contact or surface-contact measurements of changes in human tissue, hemodynamic parameters, or physiological signals. For example, photoplethysmography-based acquisition systems typically use light-emitting diodes (LEDs) and photodiodes to calculate physiological indicators such as heart rate and blood oxygen saturation by detecting changes in the intensity of transmitted or reflected light. These systems usually require sensors fixed to the skin surface to acquire continuous optical waveforms and combine them with preset optical path models, calibration coefficients, or empirical algorithms to achieve real-time acquisition of non-invasive diagnostic data.

[0003] However, when the aforementioned optical or electrophysiological sensing systems are applied to scenarios involving human movement (such as slight walking or hand tremors), the acquired signals are susceptible to baseline drift and waveform distortion caused by dynamic changes in contact pressure. For example, in wearable photoplethysmography (PPG) devices, minute relative displacement between the sensor and the skin can cause changes in the optical path length, resulting in low-frequency fluctuations in reflected light intensity unrelated to the actual pulse wave. In practical studies, this low-frequency interference may partially overlap with the frequency band of the real physiological signal, making it difficult for traditional filtering algorithms to effectively distinguish between them. This can lead to jumps in heart rate calculations, waveform peak shifts, or pulse waveform distortion, affecting the reliability of subsequent diagnostic models. Therefore, in data acquisition involving dynamic body scenarios, traditional systems often struggle to guarantee the stability and accuracy of non-invasive diagnostic data. Summary of the Invention

[0004] The purpose of this invention is to provide a method and system for acquiring non-invasive diagnostic data, aiming to solve the problems mentioned in the background art.

[0005] To solve the above-mentioned technical problems, the technical solution of the present invention is as follows:

[0006] Firstly, a method for collecting non-invasive diagnostic data, the method comprising:

[0007] The raw optical acquisition data and raw contact pressure data are acquired and time-aligned to form aligned data.

[0008] Based on the original contact pressure data in the alignment data, perform contact pressure change identification processing, and generate contact pressure influence domain data by the contact pressure change amplitude and contact pressure change rate.

[0009] Based on the domain data affected by contact pressure, optical path disturbance localization processing is performed on the original optical acquisition data in the aligned data. By identifying the change in light intensity slope that occurs synchronously with the change in contact pressure, optical interference fragment data is generated.

[0010] Based on the optical interference fragment data, optical compensation construction processing is performed on the corresponding optical interference fragments. By analyzing the change direction and amplitude of the interference fragments and the pressure change trend of the domain affected by the contact pressure of the corresponding fragments, domain-compensated optical data is generated.

[0011] Based on the domain-compensated optical data, domain-segmentation processing is performed between the domains affected by each contact pressure. Continuous correction optical data is generated by using the light intensity difference, phase difference, and waveform period difference between adjacent domains.

[0012] Based on continuously calibrated optical data, target data extraction and processing are performed to generate non-invasive diagnostic data.

[0013] Preferably, based on the original contact pressure data in the aligned data, contact pressure change identification processing is performed to generate contact pressure influence domain data through the contact pressure change amplitude and contact pressure change rate, including:

[0014] The difference between adjacent sampling points in the original contact pressure data in the alignment data is calculated to generate pressure change amplitude data.

[0015] Based on the pressure change amplitude data, the rate of change is calculated according to the change in adjacent time intervals to generate pressure change rate data;

[0016] Based on the pressure change amplitude data and pressure change rate data, compare them with the preset pressure change amplitude threshold and preset pressure change rate threshold respectively, and record the continuous time segment that meets the matching condition of any preset threshold in the comparison results as candidate pressure change segment data.

[0017] Based on the temporal continuity of the candidate pressure change segment data, candidate pressure change segments with time intervals exceeding a preset interval threshold are split into multiple pressure change segment data.

[0018] Based on the start and end points of the pressure change segment data, corresponding segment boundary point data are generated, and contact pressure influence domain data are generated based on the segment boundary point data.

[0019] Preferably, based on the contact pressure-affected domain data, optical path disturbance localization processing is performed on the original optical acquisition data in the alignment data. By identifying the light intensity slope change that occurs synchronously with the contact pressure change, optical interference segment data is generated, including:

[0020] Based on the influence of contact pressure on the domain data, the corresponding original optical acquisition data is extracted in each domain to form domain optical data.

[0021] The light intensity difference between adjacent sampling points of the domain optical data is calculated to generate light intensity change slope data.

[0022] By comparing the slope data of light intensity change with the pressure change trend in the corresponding domain, synchronous data of light intensity change points are generated.

[0023] Based on the continuity of the synchronous change point data of light intensity on the time axis, the continuous synchronous change points are combined into optical disturbance segment data.

[0024] Based on the slope change of the optical disturbance segment data at the boundary, the interference boundary points are identified, and interference segment boundary point data is generated.

[0025] Optical interference segment data is generated based on the optical disturbance segment data and the boundary point data of the interference segment.

[0026] Preferably, based on the optical interference fragment data, optical compensation construction processing is performed on the corresponding optical interference fragments. By analyzing the direction and magnitude of the interference fragment's change, as well as the pressure change trend of the affected region by its contact pressure, regional compensation optical data is generated, including:

[0027] The change direction data of each interference segment is determined based on the optical interference segment data, and compensation trend data opposite to its direction is generated based on the change direction data.

[0028] Based on the variation amplitude data of the interference segments in the optical interference segment data, the target compensation amount data is generated by uniformly distributing the data according to the time length or proportionally distributing it according to the pressure change trend.

[0029] Based on the target compensation amount data, compensation transition zone data is constructed at the start and end positions of the interference segment. The compensation amount is gradually increased or decreased from zero within the compensation transition zone to generate compensation transition data.

[0030] Compensation processing is performed on the optical interference segments based on compensation trend data, target compensation amount data, and compensation transition data to generate compensated interference segment data.

[0031] The compensated interference fragment data is combined with the undisturbed optical acquisition data in a time series to generate domain-compensated optical data.

[0032] Preferably, based on the variation amplitude data of the interference segments in the optical interference segment data, the target compensation amount data is generated by uniformly distributing the data according to the time length or proportionally distributing it according to the pressure change trend, including:

[0033] Extract the variation amplitude data of the interference segments from the optical interference segment data, and generate the time interval data of the interference segments based on the start and end times of the interference segments;

[0034] Based on the time interval data of the interference segment, the change amplitude data is divided into multiple compensation sub-regions in chronological order, so that each compensation sub-region corresponds to a time distribution segment of the change amplitude data, thus generating compensation sub-region time period data;

[0035] Based on the time period data of the compensation sub-region, the change amplitude data is evenly distributed over time and mapped to each compensation sub-region to generate uniformly distributed compensation amount data;

[0036] Based on the pressure change trend of the contact pressure influence domain to which the interference segment belongs, the change amplitude data is distributed according to the relationship of the pressure change trend in each compensation sub-region to generate trend compensation amount data.

[0037] Based on the compatibility conditions between uniformly distributed compensation data and trend compensation data, a compensation distribution method that better matches the waveform recovery characteristics of the interference segment is selected from the two to generate target compensation data.

[0038] Preferably, compensation transition zone data is constructed based on the start and end positions of the interference segment according to the target compensation amount data, so that the compensation amount gradually increases or decreases from zero within the compensation transition zone, generating compensation transition data, including:

[0039] The starting point and ending point data of the interference segment are determined based on the optical interference segment data, and the boundary point data of the interference segment are generated based on the starting point and ending point data.

[0040] Based on the boundary point data of the interference segment, the boundary compensation area is divided into several compensation transition sub-regions, and compensation transition sub-region data is generated.

[0041] Based on the compensation transition sub-region data, the target compensation amount data is gradually processed in a pattern of gradually increasing from zero to the target compensation amount or gradually decreasing from the target compensation amount, to generate compensation transition sequence data.

[0042] Based on the compensation gradient sequence data, continuous adjustment processing is performed between adjacent sampling points to ensure that the compensation amount changes smoothly point by point within the compensation transition zone, thereby generating compensation gradient adjustment data;

[0043] Based on the compensation gradual adjustment data, construct complete compensation transition zone data.

[0044] Preferably, compensation processing is performed on the optical interference segments based on compensation trend data, target compensation amount data, and compensation transition data to generate compensated interference segment data, including:

[0045] Based on the compensation trend data, the optical interference segment data is subjected to orientation correction processing, and the compensation direction is adjusted to be opposite to the change direction of the interference segment to generate orientation correction data.

[0046] Based on the orientation correction data, the target compensation data is used to perform compensation allocation processing on each sampling point inside the interference segment. The compensation amount is allocated according to the time sequence of the sampling point position to generate internal compensation allocation data.

[0047] Based on the internal compensation allocation data and compensation transition data, a compensation gradient superposition process is performed at the boundary of the interference segment to ensure that the compensation amount maintains a continuous transition between the internal part and the boundary part of the interference segment, thereby generating compensation superposition structure data.

[0048] Compensation and fusion processing is performed on the entire interference segment based on the compensated overlay structure data. By fusing the compensated overlay structure data with the original interference segment data in the time series, compensated interference segment data is generated.

[0049] Secondly, a non-invasive diagnostic data acquisition system, the system comprising:

[0050] The data alignment module is used to acquire raw optical acquisition data and raw contact pressure data, and perform time alignment processing to form aligned data;

[0051] The contact pressure change identification module is used to perform contact pressure change identification processing based on the original contact pressure data in the alignment data, and generate contact pressure influence domain data through the contact pressure change amplitude and contact pressure change rate.

[0052] The optical path disturbance positioning module is used to perform optical path disturbance positioning processing on the original optical acquisition data in the aligned data based on the contact pressure influence domain data. By identifying the light intensity slope change that occurs synchronously with the contact pressure change, it generates optical interference fragment data.

[0053] The optical compensation construction module is used to perform optical compensation construction processing on the corresponding optical interference segments based on the optical interference segment data. It generates regional compensation optical data by considering the change direction and amplitude of the interference segments and the pressure change trend of the region affected by the contact pressure of the interference segments.

[0054] The domain stitching module is used to perform domain stitching processing between each contact pressure-affected domain based on the domain compensation optical data. It generates continuous correction optical data by using the light intensity difference, phase difference, and waveform period difference between adjacent domains.

[0055] The target data extraction module is used to perform target data extraction processing based on continuously corrected optical data to generate non-invasive diagnostic data.

[0056] The above-described solution of the present invention has at least the following beneficial effects:

[0057] (a) Improve data stability and accuracy

[0058] This invention generates precisely aligned data by performing time alignment processing on the raw optical acquisition data and raw contact pressure data. This processing method ensures a strict correspondence between the two types of data on the time axis, providing a reliable foundation for subsequent interference identification and compensation. Compared with existing technologies, the time alignment method effectively reduces optical signal distortion and baseline drift caused by relative sensor displacement, significantly improving the stability and accuracy of the data acquisition process.

[0059] (ii) Effectively isolate dynamic interference from real physiological signals

[0060] Traditional non-invasive diagnostic systems are easily affected by dynamic changes in contact pressure during movement. This invention analyzes contact pressure data to identify the amplitude and rate of pressure changes and generates contact pressure influence domain data. This method can accurately divide the influence ranges of different pressure changes, facilitating independent interference identification and compensation within each domain. It effectively avoids the overlap of dynamic interference and physiological signals in traditional systems, thus improving the recognition accuracy of true physiological signals.

[0061] (III) Targeted optical interference compensation

[0062] This invention is based on the synchronicity between changes in contact pressure and changes in optical signals. It uses changes in light intensity slope to identify optical interference segments and performs targeted compensation on these segments. By analyzing the direction, amplitude, and pressure change trends of the interference segments, a compensation amount is generated and smoothed through a compensation transition zone. This allows the recovered optical signal to more accurately reflect the physiological waveform, avoiding waveform distortion or failure caused by simple filtering in traditional systems.

[0063] (iv) Enhance data splicing and processing capabilities

[0064] In this invention, based on domain-compensated optical data, domain-by-domain stitching processing is performed. Adjustments are made using the differences in light intensity, phase, and waveform period between adjacent domains to generate continuously corrected optical data. This stitching method ensures smooth connection of optical data between different domains, avoiding discontinuities or abrupt changes that occur during processing in traditional systems. This results in smooth and consistent data over time, improving the overall usability and diagnostic accuracy of the data.

[0065] (v) Improve the reliability and usability of non-invasive diagnostic data

[0066] This invention generates continuously corrected optical data through precise compensation and processing of interfering segments, and extracts target data for non-invasive diagnosis. This method avoids signal distortion caused by sensor movement, making the extracted non-invasive diagnostic data more stable and reliable, especially suitable for physiological signal acquisition during movement. Compared with existing technologies, this invention can better acquire physiological information such as heart rate and blood oxygen saturation in dynamic scenarios, providing more accurate and effective input for subsequent diagnostic models.

[0067] (vi) Applicable to practical application scenarios

[0068] By implementing the data processing method of this invention, the non-invasive diagnostic data acquisition system can better adapt to dynamically changing real-world application scenarios, such as hand tremors during movement or slight walking. This solution can accurately identify and compensate for optical interference caused by changes in contact pressure, providing more stable and efficient technical support for health monitoring in practical applications. Attached Figure Description

[0069] Figure 1 This is a flowchart of a non-invasive diagnostic data acquisition method provided in an embodiment of the present invention. Detailed Implementation

[0070] Exemplary embodiments of the present disclosure will now be described in more detail with reference to the accompanying drawings. While exemplary embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be implemented in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.

[0071] like Figure 1 As shown, embodiments of the present invention propose a method for acquiring non-invasive diagnostic data, the method comprising:

[0072] The raw optical acquisition data and raw contact pressure data are acquired and time-aligned to form aligned data.

[0073] Based on the original contact pressure data in the alignment data, perform contact pressure change identification processing, and generate contact pressure influence domain data by the contact pressure change amplitude and contact pressure change rate.

[0074] Based on the domain data affected by contact pressure, optical path disturbance localization processing is performed on the original optical acquisition data in the aligned data. By identifying the change in light intensity slope that occurs synchronously with the change in contact pressure, optical interference fragment data is generated.

[0075] Based on the optical interference fragment data, optical compensation construction processing is performed on the corresponding optical interference fragments. By analyzing the change direction and amplitude of the interference fragments and the pressure change trend of the domain affected by the contact pressure of the corresponding fragments, domain-compensated optical data is generated.

[0076] Based on the domain-compensated optical data, domain-segmentation processing is performed between the domains affected by each contact pressure. Continuous correction optical data is generated by using the light intensity difference, phase difference, and waveform period difference between adjacent domains.

[0077] Based on continuously calibrated optical data, target data extraction and processing are performed to generate non-invasive diagnostic data.

[0078] In this embodiment of the invention, by performing time alignment processing on the original optical acquisition data and the original contact pressure data, the two types of data maintain a corresponding relationship on the time axis. This provides a common analytical basis for subsequent identification of the correlation between changes in contact pressure and changes in optical signals, which is beneficial to improving the coherence of the overall data processing flow. Based on this, by identifying the magnitude and rate of change of the contact pressure data, a contact pressure influence domain is further generated. This allows data intervals under different pressure states during the acquisition process to be distinguished, laying the data structure for subsequent independent interference identification and compensation by domain.

[0079] Based on the influence of contact pressure, changes in light intensity slope synchronized with pressure variations are identified from the aligned optical acquisition data. This allows optical anomalies caused by pressure fluctuations in the acquired waveform to be accurately located within their corresponding domains. By comparing synchronicity and judging slope changes, these anomaly locations can be combined into optical interference segments, providing a clear target area for subsequent compensation operations, thus making the compensation operations more targeted.

[0080] In the optical compensation construction process, by calculating the direction and amplitude of change of the interfering segment and the corresponding pressure change trend, a compensation trend that is opposite to the interference change is formed. Based on this, the compensation amount and compensation transition zone are constructed, enabling the restoration processing to take into account the continuity of the segment's interior and boundary portions. Furthermore, by performing splicing operations between the various contact pressure-affected domains and adjusting according to the light intensity difference, phase difference, and period difference, the final restored optical data forms a continuous overall waveform across the various domains, reducing splicing inconsistencies caused by domain compensation. Finally, target diagnostic results are extracted from the continuously corrected optical data, resulting in non-invasive diagnostic data with more stable waveform characteristics and a more reliable parameter basis.

[0081] For example, in actual scenarios involving photoplethysmography (PPG) data acquisition from human fingers, pressure fluctuations may occur due to slight finger movements. When this happens, the contact pressure data first reflects the corresponding changes, thus dividing the data into different pressure-affected regions. Subsequently, within the same region, optical data showing light intensity changes synchronized with the pressure can be identified and marked as interfering segments. Compensation amounts are generated for these interfering segments based on their trends, and smooth transitions are constructed before and after them, ensuring continuity of the corrected data across segments. Pulse waveforms or other diagnostic data are then extracted from the corrected data to help obtain relatively accurate non-invasive physiological information even in the presence of interference.

[0082] In a preferred embodiment of the present invention, the original optical acquisition data and original contact pressure data are acquired and time-aligned to form aligned data, specifically including:

[0083] Optical signals from the examined area are continuously acquired using an optical sensor. The sampled light intensity values ​​are recorded sequentially according to the sampling time, forming raw optical acquisition data. Simultaneously, pressure changes between the sensor and the skin are acquired using a pressure sensor, and the pressure values ​​from each sample are recorded sequentially as raw contact pressure data. The sampling times of both types of data are extracted as time series parameters. By comparing the starting sampling times and sampling intervals of the two types of data, the data with higher sampling density is resampled sequentially to ensure its sampling interval is consistent with the other set of data. Subsequently, the two sets of sampled values ​​are mapped one-to-one according to the adjusted time series, ensuring that the optical sampling points are consistent with their corresponding pressure sampling points. This generates aligned data with a complete overall structure, allowing subsequent steps to perform analysis and processing based on the same time reference.

[0084] In a preferred embodiment of the present invention, based on the sectional compensation optical data, sectional splicing processing is performed between the sectional regions affected by each contact pressure. Continuous correction optical data is generated using the light intensity difference, phase difference, and waveform period difference between adjacent sectional regions. Specifically, this includes:

[0085] All contact pressure-affected domain compensation optical data are rearranged chronologically, recording the start and end positions of each domain. For two adjacent domains, the light intensity values ​​at both ends are extracted at the domain boundary, and the difference between them is calculated to determine if there is a sudden change in the compensated light intensity. The phase change trends of the two waveforms are compared, and the order of the peak positions on the two waveforms is used to determine if there is a phase shift. The period lengths of the two waveforms are further compared to confirm whether the rhythm of the waveforms between the domains is consistent. Based on the magnitude of the light intensity difference, phase difference, and period difference, the time period requiring adjustment is determined, and waveform translation, local stretching, and point-by-point buffering adjustments are performed on the domains with larger differences to gradually bring the two waveforms closer to consistency at the domain boundaries. Finally, all domains are sequentially stitched together to form continuous correction optical data without obvious boundaries for subsequent diagnostic data extraction.

[0086] In a preferred embodiment of the present invention, target data extraction processing is performed based on continuously corrected optical data to generate non-invasive diagnostic data, specifically including:

[0087] The continuously calibrated optical data is preprocessed by first identifying structural features in the overall waveform, such as peaks, troughs, and rising and falling segments. Based on the waveform structure, waveform segments corresponding to continuous cardiac cycles are extracted, and the main morphological parameters of each cycle are analyzed, such as peak height, rising duration, falling duration, and cycle length. Statistical analysis of the feature data from multiple cycles identifies parameters related to periodic changes and physiological information in the acquisition area, such as the periodic frequency or amplitude variation trend of the pulse wave signal. Combining the waveform features within continuous segments, and through the selection and combination of these structural feature parameters, non-invasive diagnostic data for subsequent diagnostic purposes is generated, ensuring that the analysis results reflect the physiological fluctuations of the subject under the current acquisition state.

[0088] In a preferred embodiment of the present invention, based on the original contact pressure data in the alignment data, contact pressure change identification processing is performed to generate contact pressure influence domain data through the contact pressure change amplitude and contact pressure change rate, including:

[0089] The difference between adjacent sampling points in the original contact pressure data in the alignment data is calculated to generate pressure change amplitude data.

[0090] Based on the pressure change amplitude data, the rate of change is calculated according to the change in adjacent time intervals to generate pressure change rate data;

[0091] Based on the pressure change amplitude data and pressure change rate data, compare them with the preset pressure change amplitude threshold and preset pressure change rate threshold respectively, and record the continuous time segment that meets the matching condition of any preset threshold in the comparison results as candidate pressure change segment data.

[0092] Based on the temporal continuity of the candidate pressure change segment data, candidate pressure change segments with time intervals exceeding a preset interval threshold are split into multiple pressure change segment data.

[0093] Based on the start and end points of the pressure change segment data, corresponding segment boundary point data are generated, and contact pressure influence domain data are generated based on the segment boundary point data.

[0094] In this embodiment of the invention, by calculating the difference between adjacent sampling points of the original contact pressure data, the instantaneous amplitude of pressure change can be represented in data form, facilitating the determination of whether pressure fluctuations occur during the acquisition process. Combining the amplitude change with adjacent time intervals, the pressure change rate can be further obtained, expanding the pressure change process from a single amplitude description to a combined amplitude and rate characterization, which helps identify different pressure states of short-term rapid changes and long-term gradual changes. Subsequently, comparing these two types of change data with preset amplitude and rate thresholds allows for the initial differentiation of continuous time segments under different pressure change modes, forming candidate regions. Further segmenting the candidate regions according to temporal continuity avoids misjudgments of small segments due to sampling jitter, resulting in more stable and consistent pressure change segments. By generating segment boundary points, different pressure change modes can be clearly demarcated, forming a complete contact pressure influence domain, enabling subsequent interference identification of optical data to be analyzed within a clearly defined domain.

[0095] In a preferred embodiment of the present invention, the method for setting a preset pressure change amplitude threshold specifically includes:

[0096] Based on the collected raw contact pressure data, statistical analysis is first performed. The fluctuation range of the sampled data is determined by calculating parameters such as the standard deviation, maximum value, and minimum value of multiple samples over a period of time. Then, a pressure change amplitude threshold is set, which can be selected based on actual application requirements. For example, under normal use, the fluctuation range of contact pressure should generally be within an acceptable range; changes exceeding this range are considered abnormal. During implementation, the threshold selection can refer to the equipment's sensitivity or the maximum pressure change fluctuation in typical application scenarios. The threshold is typically set to twice the standard deviation of the contact pressure data or other appropriate multiples to ensure sufficient sensitivity to sudden, large changes. Finally, based on the selected amplitude threshold, if the raw contact pressure change amplitude exceeds this threshold, it is considered a significant pressure change event and subsequent processing is performed.

[0097] In a preferred embodiment of the present invention, the method for setting a preset pressure change rate threshold specifically includes:

[0098] By statistically analyzing the rate of change of contact pressure data, basic statistics such as the average, maximum, and minimum rates of pressure change are calculated. First, the rate of pressure change between every two consecutive sampling points is calculated, i.e., the amount of pressure change per unit time. Then, based on these rate data, an appropriate rate threshold is selected. This threshold is typically set several times the average or standard deviation of the pressure change rate data to ensure it is sufficient to identify drastic pressure changes. For example, if the sampling interval is 1 millisecond, the rate threshold can be set to twice or more the standard rate. The rate threshold can also be adjusted according to the application scenario; if the device requires a response to rapidly changing pressure, the threshold can be set lower. Finally, if the rate of pressure change exceeds the set threshold, the pressure change at that moment is considered a sudden or rapid change, and it is subsequently marked as a potential change segment for further analysis and processing.

[0099] In a preferred embodiment of the present invention, the method for setting the preset interval threshold specifically includes:

[0100] Based on the time interval information of the collected contact pressure data, the time intervals of consecutive sampling points are first statistically analyzed to calculate the average and standard deviation of the time intervals. Next, a preset interval threshold is set, which can be selected according to actual needs. Typically, if the sampling time interval is too long (e.g., the time interval between sampling points exceeds the set threshold), it may indicate loss or anomalies in the sampling process, leading to intermittent pressure changes. To avoid this, the threshold can be set to twice or more the average time interval plus the standard deviation to ensure that the threshold is effective within the normal fluctuation range. For example, if the average sampling time interval is 100 milliseconds and the standard deviation is 20 milliseconds, then the set interval threshold can be 140 milliseconds or higher. If the time interval between consecutive sampling points exceeds the preset threshold, the data segment is considered to have discontinuities or anomalies, and further data splitting and processing are performed.

[0101] In a preferred embodiment of the present invention, the difference between adjacent sampling points is calculated on the original contact pressure data in the alignment data to generate pressure change amplitude data, specifically including:

[0102] Each contact pressure sampling point is sequentially read from the aligned data to obtain a time sequence of all pressure values. For two adjacent pressure values, the numerical difference between them is calculated, and this difference is taken as the pressure change amplitude at the corresponding time. Specifically, the pressure value at the later sampling point is directly subtracted from the pressure value at the earlier sampling point to obtain the magnitude and direction of the pressure change. The change amplitudes of all adjacent sampling points are recorded in chronological order to form a set of pressure change amplitude data describing the magnitude of pressure change over time, which is used for subsequent quantitative analysis of the degree of pressure change.

[0103] In a preferred embodiment of the present invention, pressure change rate data is generated by calculating the rate of change based on the pressure change amplitude data according to the change in adjacent time intervals, specifically including:

[0104] Each pressure change amplitude value is read from the pressure change amplitude data, and the corresponding sampling time is obtained when the amplitude is recorded. The sampling times of two consecutive amplitude changes are interpolated to determine the time interval between them. The amplitude value is then compared with its corresponding time interval. By expressing the change in amplitude within the time interval as the change per unit time, the pressure change rate for each sampling point is generated. All calculated change rates are recorded in chronological order, quantifying the rate of pressure change over time and providing fundamental data for subsequently identifying rapidly changing or gradually changing zones.

[0105] In a preferred embodiment of the present invention, pressure change amplitude data and pressure change rate data are compared with preset pressure change amplitude thresholds and preset pressure change rate thresholds, respectively. Continuous time segments that satisfy any preset threshold matching condition in the comparison results are recorded as candidate pressure change segment data. Specifically, this includes:

[0106] The pressure change amplitude data at each sampling time is evaluated. When the change amplitude exceeds a preset pressure change amplitude threshold, that time is marked as a point where pressure change may occur. Similarly, the pressure change rate data at each sampling time is evaluated. When the change rate exceeds a preset pressure change rate threshold, that time is marked as a point where rapid pressure change may occur. Subsequently, these marked times are analyzed for continuity according to time sequence. Adjacent or very short-interval marked times are combined into continuous time segments, and these segments are recorded as candidate pressure change segment data to initially identify areas where pressure fluctuations may occur.

[0107] In a preferred embodiment of the present invention, based on the temporal continuity of the candidate pressure change segment data, candidate pressure change segments with time intervals exceeding a preset interval threshold are split to generate multiple pressure change segment data, specifically including:

[0108] The process iterates through each time segment in the candidate pressure change data, detecting the sampling time interval between any two adjacent points within that segment. When a certain adjacent time interval exceeds a preset threshold, the segment is split at that interval, forming two independent time periods. This process is repeated for all segments, ultimately dividing excessively long or discontinuous candidate segments into multiple more continuous pressure change segments. This results in more consistent pressure changes within each segment, providing a more accurate structural basis for generating domain boundaries.

[0109] In a preferred embodiment of the present invention, corresponding segment boundary point data is generated based on the start and end points of the pressure change segment data, and contact pressure influence domain data is generated based on the segment boundary point data, specifically including:

[0110] From each pressure change segment, the start and end sampling times of that segment are extracted, and these two times are defined as the corresponding segment boundary points. All segment boundary points are then arranged chronologically to form a complete boundary point sequence. Based on this boundary point sequence, the aligned data is divided into multiple independent pressure influence intervals, each corresponding to a pressure change domain, giving each pressure change pattern a clearly defined range in the data structure. Finally, these domains are labeled and output as contact pressure influence domain data, providing a precise segmentation basis for subsequent interference identification of optical data by domain.

[0111] In a preferred embodiment of the present invention, based on the contact pressure influence domain data, optical path disturbance localization processing is performed on the original optical acquisition data in the alignment data. By identifying the light intensity slope change that occurs synchronously with the contact pressure change, optical interference segment data is generated, including:

[0112] Based on the influence of contact pressure on the domain data, the corresponding original optical acquisition data is extracted in each domain to form domain optical data.

[0113] The light intensity difference between adjacent sampling points of the domain optical data is calculated to generate light intensity change slope data.

[0114] By comparing the slope data of light intensity change with the pressure change trend in the corresponding domain, synchronous data of light intensity change points are generated.

[0115] Based on the continuity of the synchronous change point data of light intensity on the time axis, the continuous synchronous change points are combined into optical disturbance segment data.

[0116] Based on the slope change of the optical disturbance segment data at the boundary, the interference boundary points are identified, and interference segment boundary point data is generated.

[0117] Optical interference segment data is generated based on the optical disturbance segment data and the boundary point data of the interference segment.

[0118] In this embodiment of the invention, by using contact pressure influence domain segmentation to extract optical acquisition data, the optical waveform and pressure changes maintain a range-based correspondence, avoiding mixed analysis across pressure states. Calculating the light intensity difference in the segmented optical data transforms the original waveform's trend into slope data, making sudden increases and decreases in light intensity easier to identify. Comparing the slope of the light intensity change with the corresponding pressure change trend in the segment allows for the selection of optical change points correlated with pressure changes, thus eliminating genuine fluctuations caused by physiological factors and improving the accuracy of interference identification. By judging the temporal continuity of synchronous change points, multiple single-point anomalies can be merged into optical disturbance segments with complete interference characteristics, giving the abnormal segments a continuous structure on the time axis. Furthermore, generating interference segment boundary points based on the slope changes at the boundaries of the disturbance segments allows for precise calibration of the interference start and end positions. The final optical interference segment data enables subsequent compensation steps to directly apply to the abnormal segments, reducing the impact on the normal waveform.

[0119] In a preferred embodiment of the present invention, based on the regional data influenced by contact pressure, corresponding original optical acquisition data is extracted in each regional region to form regional optical data, specifically including:

[0120] Each region in the contact pressure influence regional data is read according to its recorded start and end times, and compared with the optical sampling times recorded in the alignment data to identify all sampling points in the optical data that match the time range of that region. These sampling points are stored in chronological order to form the regional optical data corresponding to that region. The above extraction process is repeated for all pressure change regions, so that the original optical acquisition data is split into multiple regional optical sequences according to the regional structure, providing an independent basis for subsequent interference identification in each region.

[0121] In a preferred embodiment of the present invention, the intensity difference between adjacent sampling points of the domain optical data is calculated to generate intensity change slope data, specifically including:

[0122] The light intensity value of each sampling point in the segmented optical data is read sequentially, and its numerical difference is compared with that of the next adjacent sampling point. The light intensity value of the next sampling point is subtracted from that of the previous sampling point to obtain the change in light intensity within that time interval. This change is recorded according to the sampling time sequence of the optical data, so that each change corresponds to a specific time position, thus forming the light intensity change slope data. By expressing the light intensity difference between adjacent sampling points in this way, the fast-changing and slow-changing segments of the original waveform in the time dimension can be effectively distinguished, providing basic characteristics for subsequent judgment of optical interference.

[0123] In a preferred embodiment of the present invention, synchronous comparison is performed between the light intensity change slope data and the pressure change trend within the corresponding domain to generate light intensity synchronous change point data, specifically including:

[0124] The pressure change trend of a specific region is read from the contact pressure influence data and recorded as a time series. This trend is then compared with the light intensity change slope data on the same time axis to determine if the direction of light intensity change is consistent with the direction of pressure change, and whether the magnitude of the light intensity change shows a positive or negative correlation with the magnitude of pressure change. When the light intensity change slope shows a synchronous relationship with the pressure change trend in terms of change trend, direction, or magnitude, the corresponding sampling location is recorded as a synchronous light intensity change point. Through this comparison process, optical change points that are correlated with pressure changes within that region can be screened, providing a reliable basis for subsequent interference fragment identification.

[0125] In a preferred embodiment of the present invention, based on the continuity of the light intensity synchronous change point data on the time axis, the continuous synchronous change points are combined into optical perturbation segment data, specifically including:

[0126] All light intensity synchronization change points are sorted according to their sampling time, and the sampling time interval between adjacent synchronization change points is determined point by point. When the interval between two adjacent points is less than or equal to a preset continuity threshold, these two points are considered as continuous parts of the same disturbance trend. As time progresses, all synchronization change points that meet the continuity condition are combined into an optical disturbance segment. When the time interval between two synchronization change points exceeds the set threshold, the current disturbance segment is terminated, and a new disturbance segment is recorded starting from the next point. The generated optical disturbance segment data can accurately represent the disturbance segment affected by pressure synchronization.

[0127] In a preferred embodiment of the present invention, identifying interference boundary points based on the slope change of optical perturbation segment data at the boundary, and generating interference segment boundary point data, specifically includes:

[0128] The slope of the light intensity change between the first and second sampling points in each optical perturbation segment is read, and it is determined whether the trend of this change deviates significantly from the overall trend within the perturbation segment. When the deviation exceeds a set boundary deviation threshold, this position is marked as the starting boundary point of the perturbation segment. Similarly, the slope of the light intensity change between the last sampling point and the previous sampling point in the perturbation segment is read. When the trend of this change is inconsistent with the overall trend within the perturbation segment, this position is marked as the ending boundary point of the perturbation segment. This method can determine the accurate start and end positions of the interference segment, thus accurately expressing the range of interference.

[0129] In a preferred embodiment of the present invention, generating optical interference segment data based on optical disturbance segment data and interference segment boundary point data specifically includes:

[0130] The start and end points of the interference segment boundary point data are used as the boundary range of the interference segment. All sampling points within this boundary range are extracted from the optical perturbation segment data, arranging the sampling values ​​within the interference segment in chronological order. The extracted sequence is used as the core content of the interference segment, and the start and end points of the interference segment, along with their corresponding optical values, are combined to form complete optical interference segment data. This optical interference segment data can fully describe the location, duration, and variation pattern of the interference, providing accurate input for subsequent optical compensation construction.

[0131] In a preferred embodiment of the present invention, optical compensation construction processing is performed on the corresponding optical interference segments based on the optical interference segment data. This process generates domain-specific compensated optical data by considering the direction and magnitude of the interference segment's change, as well as the pressure change trend of the contact pressure affecting the domain to which it belongs. The data includes:

[0132] The change direction data of each interference segment is determined based on the optical interference segment data, and compensation trend data opposite to its direction is generated based on the change direction data.

[0133] Based on the variation amplitude data of the interference segments in the optical interference segment data, the target compensation amount data is generated by uniformly distributing the data according to the time length or proportionally distributing it according to the pressure change trend.

[0134] Based on the target compensation amount data, compensation transition zone data is constructed at the start and end positions of the interference segment. The compensation amount is gradually increased or decreased from zero within the compensation transition zone to generate compensation transition data.

[0135] Compensation processing is performed on the optical interference segments based on compensation trend data, target compensation amount data, and compensation transition data to generate compensated interference segment data.

[0136] The compensated interference fragment data is combined with the undisturbed optical acquisition data in a time series to generate domain-compensated optical data.

[0137] In this embodiment of the invention, by analyzing the direction of change in the optical interference segment data, the main trend of the interference waveform can be determined, ensuring that the subsequent compensation direction is consistent with the interference, thus providing a basis for restoring the original waveform. By processing the amplitude of the interference segment and the pressure change trend in its respective domain, a compensation amount reflecting the degree of interference and the law of pressure change can be generated, making the compensation operation both quantitatively based and trend-consistent. Constructing a compensation transition zone allows the compensation amount to gradually change at the start and end positions of the interference segment, avoiding abrupt waveform changes caused by direct compensation. Further compensation processing is performed on the interference segment, comprehensively applying the compensation trend, compensation amount, and the change mode of the transition zone to the original waveform, allowing the interference segment to be adjusted point by point in time, thereby forming a compensation result that is closer to the real physiological waveform. By stitching the compensated interference segment with the undisturbed optical acquisition data, domain-compensated optical data that maintains continuity can be obtained, providing a smoother base signal for subsequent cross-domain stitching.

[0138] In a preferred embodiment of the present invention, the change direction data of each interference segment is determined based on the optical interference segment data, and compensation trend data opposite to its direction is generated based on the change direction data, specifically including:

[0139] All sampling points of the interference segment are sequentially read from the optical interference segment data, and the light intensity values ​​of each sampling point are arranged in chronological order to form an interference waveform sequence. For the light intensity change between each adjacent sampling point, it is determined whether the change is positive or negative, thereby determining the direction of change of the segment at that time position. The direction of change at each position in the entire interference segment is then statistically analyzed to determine whether the interference segment exhibits an overall upward trend, a downward trend, or a composite trend such as an initial upward trend followed by a downward trend. Based on the obtained directional trend, the compensation trend is designed to be opposite to the interference trend. For example, when the overall interference segment shows an upward trend, the compensation trend is made downward; when the interference shows a downward trend, the compensation trend is made upward. The final generated compensation trend data can be used as the directional reference for subsequent compensation amount calculation and compensation operation.

[0140] In a preferred embodiment of the present invention, the compensated interference segment data and the undisturbed optical acquisition data are time-series composited to generate domain-compensated optical data, specifically including:

[0141] First, the compensated interference segment data is sorted according to the original acquisition time, maintaining its order on the time axis. Then, undisturbed optical acquisition data is read from the aligned data and compared with the compensated interference segments according to the sampling time sequence to identify the time interval corresponding to the compensated segment. Within the interval containing the compensated segment, the original interference sampling points are replaced with compensated sampling points, ensuring the interference section is completely covered by the compensated waveform; outside the compensated section, the original normal waveform remains unchanged. Next, an overall coherence check is performed on the stitched data to ensure there are no discontinuities at the boundary between the compensated and normal segments, and necessary fine-tuning is used to keep the intensity changes at the boundary smooth. The final domain-compensated optical data can serve as the basis for subsequent cross-domain stitching, ensuring the entire data segment maintains the continuity of the time sequence and structural consistency after compensation processing.

[0142] In a preferred embodiment of the present invention, target compensation data is generated by uniformly distributing the target compensation amount data according to the time length or proportionally distributing it according to the pressure change trend based on the variation amplitude data of the interference segments in the optical interference segment data, including:

[0143] Extract the variation amplitude data of the interference segments from the optical interference segment data, and generate the time interval data of the interference segments based on the start and end times of the interference segments;

[0144] Based on the time interval data of the interference segment, the change amplitude data is divided into multiple compensation sub-regions in chronological order, so that each compensation sub-region corresponds to a time distribution segment of the change amplitude data, thus generating compensation sub-region time period data;

[0145] Based on the time period data of the compensation sub-region, the change amplitude data is evenly distributed over time and mapped to each compensation sub-region to generate uniformly distributed compensation amount data;

[0146] Based on the pressure change trend of the contact pressure influence domain to which the interference segment belongs, the change amplitude data is distributed according to the relationship of the pressure change trend in each compensation sub-region to generate trend compensation amount data.

[0147] Based on the compatibility conditions between uniformly distributed compensation data and trend compensation data, a compensation distribution method that better matches the waveform recovery characteristics of the interference segment is selected from the two to generate target compensation data.

[0148] In this embodiment of the invention, by extracting the amplitude variation data from the optical interference segment and combining it with the start and end times of the interference segment to generate time interval data, the range of the interference segment on the time axis can be quantitatively expressed, facilitating subsequent processing by splitting it according to time sequence. Dividing this time interval into multiple compensation sub-regions allows for finer-grained distribution of amplitude variation, avoiding the problem of insufficient or excessive compensation in certain areas caused by applying a single compensation amount to the entire segment. By mapping the amplitude variation according to time relationships to generate a uniformly distributed compensation amount, the compensation requirements can be met when the interference is relatively evenly distributed within the segment. Furthermore, distributing the amplitude variation according to the pressure change trend is more adaptable in situations where the correlation between interference and pressure trends is strong, allowing the compensation amount to reflect the direction and strength of the interference change. Combining the adaptability selection of uniform distribution and trend-based methods can improve the accuracy of compensation amount generation, making subsequent compensation processing more consistent with the actual change structure of the interference segment and providing a more reasonable basis for waveform restoration.

[0149] In a preferred embodiment of the present invention, the variation amplitude data of the interference segment is extracted based on the optical interference segment data, and the time interval data of the interference segment is generated based on the start time and end time of the interference segment, specifically including:

[0150] A complete sampling sequence of the interference segment is obtained from the optical interference segment data, and the light intensity value of each sampling point is recorded in chronological order. Then, the light intensity difference between adjacent sampling points is read to determine the light intensity variation amplitude at various times within the interference segment, and these amplitudes are recorded sequentially as amplitude data. Furthermore, the sampling times of the first and last sampling points of the interference segment are read and recorded as the start and end times of the interference segment. Based on these two time points, a time interval corresponding to the interference segment is generated, clearly defining the range of the interference segment in the time dimension and providing a time reference for subsequent compensation calculation and allocation.

[0151] In a preferred embodiment of the present invention, based on the time interval data of the interference segment, the change amplitude data is divided into multiple compensation sub-regions in chronological order, so that each compensation sub-region corresponds to a time distribution segment of the change amplitude data, generating compensation sub-region time interval data, specifically including:

[0152] First, based on the length of the interference segment's time interval, the interval is divided into several consecutive sub-segments of equal duration, ensuring each sub-segment occupies the same length on the time axis. Then, the variation amplitude data is mapped to each sub-segment according to the sampling time sequence, assigning variation amplitude data belonging to a specific time range to the corresponding compensation sub-segment. For each compensation sub-segment, all variation amplitude data contained within it are recorded, and a time range, start and end sampling points, and internal data list are established for each compensation sub-segment, ensuring each compensation sub-segment has a complete structure in both the time and data distribution dimensions.

[0153] In a preferred embodiment of the present invention, based on the time period data of the compensation sub-regions, the variation amplitude data is mapped to each compensation sub-region in a time-uniform distribution to generate uniformly distributed compensation amount data, specifically including:

[0154] Based on the time length of all compensation sub-segments, the compensation allocation ratio that each compensation sub-segment should bear in the entire interference segment is determined. Then, based on the total amount of variation amplitude data, this total amount is evenly distributed to each compensation sub-segment according to the time distribution ratio, so that each compensation sub-segment receives an initial compensation value proportional to its time length. Subsequently, for the sampling points within each sub-segment, this initial compensation value is averaged again according to the number of sampling points, so that each sampling point receives the same compensation amount. The final uniformly distributed compensation amount data is a set of compensation amount sequences continuously distributed over time, providing a reference for subsequent trend compensation comparisons.

[0155] In a preferred embodiment of the present invention, based on the pressure change trend of the contact pressure influence domain to which the interference segment belongs, the change amplitude data is distributed according to the relationship of the pressure change trend in each compensation sub-region to generate trend compensation amount data, specifically including:

[0156] The pressure change trend of the affected segment is extracted from the contact pressure influence domain data, and this trend is unfolded into a trend change sequence in chronological order. Based on the trend change sequence, the pressure value changes within each time sub-segment are determined; for example, sub-segments with faster pressure increases are considered to have a more significant interference effect. Subsequently, the change amplitude data is allocated according to the degree of pressure change in different sub-segments, so that areas with larger pressure changes receive more compensation, while areas with smaller pressure changes receive less compensation. For each sub-segment, the compensation amount is further gradually transitioned from the starting position to the ending position of the sub-segment according to its time position within the sub-segment, so that the trend compensation data presents a distribution structure consistent with the pressure change trend.

[0157] In a preferred embodiment of the present invention, based on the compatibility condition between uniformly distributed compensation data and trend compensation data, a compensation distribution mode that better suits the waveform recovery characteristics of the interference segment is selected from the two to generate target compensation data, specifically including:

[0158] First, structural analysis is performed on both the uniformly distributed compensation data and the trend compensation data to check whether the two distribution patterns match the main direction, location, and intensity of change of the interference segment. For example, when the change of the interference segment is relatively uniform throughout the time interval, the uniformly distributed compensation method is preferred; when the change of the interference segment has a clear synchronous relationship with the pressure trend, the trend compensation method is preferred. Then, the compensation sequence of the selected compensation method is organized chronologically to form complete target compensation data, ensuring that this data is consistent with the time interval and change amplitude of the interference segment, so that it can serve as input data for constructing the compensation transition zone and subsequent compensation processing.

[0159] In a preferred embodiment of the present invention, compensation transition region data is constructed based on the start and end positions of the interference segment according to the target compensation amount data, so that the compensation amount gradually increases or decreases from zero within the compensation transition region, thereby generating compensation transition data, including:

[0160] The starting point and ending point data of the interference segment are determined based on the optical interference segment data, and the boundary point data of the interference segment are generated based on the starting point and ending point data.

[0161] Based on the boundary point data of the interference segment, the boundary compensation area is divided into several compensation transition sub-regions, and compensation transition sub-region data is generated.

[0162] Based on the compensation transition sub-region data, the target compensation amount data is gradually processed in a pattern of gradually increasing from zero to the target compensation amount or gradually decreasing from the target compensation amount, to generate compensation transition sequence data.

[0163] Based on the compensation gradient sequence data, continuous adjustment processing is performed between adjacent sampling points to ensure that the compensation amount changes smoothly point by point within the compensation transition zone, thereby generating compensation gradient adjustment data;

[0164] Based on the compensation gradual adjustment data, construct complete compensation transition zone data.

[0165] In this embodiment of the invention, by determining the start and end positions from the optical interference segment, a clear boundary range can be provided for the construction of the subsequent compensation transition zone, giving the compensation transformation a positional reference. The compensation transition sub-regions are divided according to the waveform changes between the boundary positions, giving the transition zone a structured nature, allowing for differentiated compensation processing based on the variation characteristics within different sub-regions. By generating a compensation gradient sequence in a gradually increasing or decreasing pattern, it is ensured that the compensation amount changes smoothly within the transition zone, avoiding abrupt changes in the waveform. Continuous adjustment of the gradient sequence ensures that the compensation difference between adjacent sampling points remains stable, contributing to a more natural change curve and reducing jump errors. The final constructed compensation transition zone data establishes a smooth connection structure between the interference segment boundary and the compensation amount within the segment, providing a smooth transition for the overall compensation of subsequent segments, and helping to improve the continuity and integrity of the waveform after compensation processing.

[0166] In a preferred embodiment of the present invention, determining the start-point and end-point data of the interference segment based on the optical interference segment data, and generating interference segment boundary point data based on the start-point and end-point data, specifically includes:

[0167] Complete time series and corresponding light intensity values ​​are extracted from optical interference segment data. The light intensity changes at each sampling point are examined sequentially. Points where the light intensity fluctuates significantly and the direction of change differs from or abruptly changes from the previous time period are marked as potential start or end points of interference. Specifically, the first sampling point showing a fluctuation is considered the start point of the interference segment, and the last sampling point showing a fluctuation is considered the end point. If the amplitude of light intensity changes within the interference segment reaches a set threshold and the trend of change changes significantly, these points are further confirmed as boundary points of the interference segment. Using this method, the start and end times of the interference segment are determined based on the characteristics of light intensity value changes, and boundary point data of the interference segment is generated for subsequent processing.

[0168] In a preferred embodiment of the present invention, based on the boundary point data of the interference segment, boundary compensation region division processing is performed on the target compensation amount data to divide the boundary region into several compensation transition sub-regions, generating compensation transition sub-region data, specifically including:

[0169] Based on the previously generated boundary point data of the interference segments, the start and end points of the interference segments are first determined, and then a boundary compensation zone is defined between these two points. Within the compensation zone, it is divided into several consecutive compensation transition sub-zones, each corresponding to a time period. To ensure the smoothness of the compensation transition, the time length of each sub-zone is typically set based on the rate of change of the compensation amount or the characteristics of the pressure change. For areas with drastic changes, the compensation transition sub-zone time is shorter, while for areas with gentler changes, a longer sub-zone time is used. The range, time length, and variation pattern of the compensation amount for each compensation transition sub-zone are clearly recorded and used as the basic data in the compensation process to generate compensation transition sub-zone data.

[0170] In a preferred embodiment of the present invention, based on the compensation transition sub-region data, the target compensation amount data is gradually processed according to a pattern of gradually increasing or gradually decreasing from zero to the target compensation amount to generate compensation gradual sequence data, specifically including:

[0171] First, starting from the beginning of each compensation transition sub-region, the initial value of the target compensation amount is set to zero. Based on the increasing trend of compensation amount within that region, the compensation amount is gradually increased until the target value is reached. Similarly, at the end of the compensation transition region, the compensation amount is gradually reduced to zero, ensuring a smooth transition and avoiding abrupt changes. Based on the time span of each sub-region, the rate of increase in compensation amount, and the pressure change trend, a continuous sequence of compensation amount changes is generated, allowing the compensation amount to gradually change throughout the entire interval, ensuring a smooth compensation effect without abrupt changes. In this way, the generated compensation gradient sequence data naturally reflects the correction process of the optical waveform, avoiding distortion or discontinuities that may occur during the compensation process.

[0172] In a preferred embodiment of the present invention, continuous adjustment processing between adjacent sampling points is performed based on the compensation gradient sequence data, so that the compensation amount changes smoothly point by point within the compensation transition region, generating compensation gradient adjustment data, specifically including:

[0173] The compensation gradient sequence data is read, and the difference between each pair of adjacent sampling points in the sequence is checked to ensure that the change in compensation between any two adjacent sampling points does not exceed a set threshold. If the difference is too large, the compensation amount in that region is adjusted by linear interpolation to make the change smoother. During the adjustment process, if the change in compensation amount is large, the length of the transition region is appropriately increased and the rate of change is slowed down; if the change is small, the transition region is appropriately shortened and the change in compensation amount is accelerated. In this way, the adjusted compensation gradient sequence allows the compensation amount to gradually increase or decrease in the time dimension, ensuring a smooth transition in the compensation process, thereby generating the final compensation gradient adjustment data.

[0174] In a preferred embodiment of the present invention, a complete compensation transition zone data is constructed based on the compensation gradual adjustment data, specifically including:

[0175] The compensation values ​​for each time point are extracted from the compensation gradient adjustment data, and these compensation values ​​are mapped to the time periods of the compensation transition zone according to the time series. The data structure of the entire compensation transition zone is constructed based on the trend and duration of the compensation values, and the changes in compensation values ​​within each time period are recorded. The final compensation transition zone data will contain the complete compensation process from the start to the end of the compensation zone, ensuring smooth changes in compensation values ​​for each time period and meeting the requirements of optical waveform correction. The compensation transition zone data can be used by subsequent compensation processing modules as a smooth transition zone for the compensated waveform, ensuring that there are no unnatural jumps or abrupt changes after waveform correction.

[0176] In a preferred embodiment of the present invention, compensation processing is performed on the optical interference segment based on compensation trend data, target compensation amount data, and compensation transition data to generate compensated interference segment data, including:

[0177] Based on the compensation trend data, the optical interference segment data is subjected to orientation correction processing, and the compensation direction is adjusted to be opposite to the change direction of the interference segment to generate orientation correction data.

[0178] Based on the orientation correction data, the target compensation data is used to perform compensation allocation processing on each sampling point inside the interference segment. The compensation amount is allocated according to the time sequence of the sampling point position to generate internal compensation allocation data.

[0179] Based on the internal compensation allocation data and compensation transition data, a compensation gradient superposition process is performed at the boundary of the interference segment to ensure that the compensation amount maintains a continuous transition between the internal part and the boundary part of the interference segment, thereby generating compensation superposition structure data.

[0180] Compensation and fusion processing is performed on the entire interference segment based on the compensated overlay structure data. By fusing the compensated overlay structure data with the original interference segment data in the time series, compensated interference segment data is generated.

[0181] In this embodiment of the invention, by using compensation trend data to correct the direction of optical interference segments, the compensation direction can be made opposite to the interference direction, ensuring that the compensation produces a cancellation effect rather than enhancing the interference. The compensation amount is sequentially allocated to each sampling point within the interference segment according to the target compensation amount data, ensuring that the compensation magnitude of each sampling point has a clear basis and avoiding local waveform distortion caused by uneven compensation within the segment. Combined with compensation transition zone data, a gradual overlay process is performed at the boundary of the interference segment, enabling a smooth transition in intensity and change trend between the internal compensation area and the boundary compensation area, reducing the possibility of abrupt changes at the boundary before and after compensation. By fusing direction correction data, internal compensation allocation data, and compensation overlay structure data, overall compensation is performed on the interference segment, achieving unified adjustment between different data structures. This ensures that the compensated waveform can both correct the offset caused by interference and maintain the stability of its internal shape. The processed compensated interference segment can form a continuous signal structure with the undisturbed waveform, providing more reliable basic data for subsequent data splicing and diagnostic analysis.

[0182] In a preferred embodiment of the present invention, direction correction processing is performed on the optical interference segment data based on the compensation trend data, adjusting the compensation direction to be opposite to the change direction of the interference segment, and generating direction correction data, specifically including:

[0183] The compensation direction change information is read from the compensation trend data to identify the trend of the interference segment (e.g., rising, falling, etc.). Then, based on the compensation trend data, it is determined that the compensation should be opposite to the direction of the original interference segment. For example, when the interference segment shows an upward trend, the compensation trend should be downward; when the interference segment shows a downward trend, the compensation trend should be upward. Next, each sampling point in the optical interference segment data is read, and the direction of the compensation amount is applied in reverse to each sampling point, modifying the value of the compensation amount to meet the correction requirement of being opposite to the interference direction. In this way, direction correction data is generated, which can correct the optical interference segment to a form consistent with the compensation trend direction, ensuring that the compensated waveform correction meets actual requirements.

[0184] In a preferred embodiment of the present invention, based on the direction correction data, compensation allocation processing is performed on each sampling point within the interference segment using the target compensation amount data. The compensation amount is allocated according to the time sequence of the sampling point positions to generate internal compensation allocation data, specifically including:

[0185] Based on the compensation direction at each time point in the direction correction data, the target compensation amount data is allocated to each sampling point in chronological order. First, the total amount of compensation data and the time range for allocation are determined. Then, according to the location and time distribution of the sampling points, the compensation data is allocated to each sampling point, ensuring that the magnitude of the compensation amount is proportional to its time position. For example, the compensation amount is smaller at the beginning of the time period and gradually increases (or decreases, depending on the compensation trend) as the sampling point approaches the end of the time period. Finally, the compensation amount at each sampling point is recorded, forming a continuous compensation allocation process within the interference segment, ensuring that the compensation effect is consistent with the changing trend of the optical waveform.

[0186] In a preferred embodiment of the present invention, based on internal compensation allocation data and compensation transition data, a compensation gradient superposition process is performed at the boundary of the interference segment to ensure a continuous transition of the compensation amount between the internal part and the boundary part of the interference segment, thereby generating compensation superposition structure data, specifically including:

[0187] First, the internal compensation allocation data is read and the distribution of compensation amounts at each time point is obtained, with particular attention paid to the change process of the compensation amounts. Based on the compensation transition data, the transition range of the compensation amount at the boundary of the interference segment is determined. At the beginning and end of the compensation segment, the compensation amount is gradually adjusted, increasing or decreasing gradually to ensure a smooth transition at the boundary without abrupt changes. Following this gradual approach, the compensation amount is gradually applied, expanding from the interior of the segment towards the boundary, generating a compensation superposition effect at the boundary. This ensures a natural connection between the compensation amount at the boundary and the compensation amount in the interior, without any discontinuities or over-adjustment. Ultimately, the generated compensation superposition structure data smoothly connects the compensation amounts at the beginning and end on the time axis, ensuring the continuous and stable compensation effect throughout the entire interference segment.

[0188] In a preferred embodiment of the present invention, compensation and fusion processing is performed on the entire interference segment based on the compensation superposition structure data. This is achieved by fusing the compensation superposition structure data with the original interference segment data in a time series to generate compensated interference segment data. Specifically, this includes:

[0189] First, the compensated overlay structure data and the original interference segment data are time-aligned to ensure that the two data sequences have the same time base. Next, based on the compensation amount generated in the compensated overlay structure data, the compensation amount is added point-by-point to each corresponding sampling point in the original interference segment data. For sampling points within the compensated section, the original optical signal is replaced with the compensation amount; for unaffected parts, the original data remains unchanged. In particular, attention is paid to the overlay method of the compensation amount to ensure that its introduction does not lead to over-correction or distortion of the data. Finally, the compensated interference segment data is smoothed to form a continuous and smooth waveform on the time axis, consistent with the original data, and the intensity changes within the interference region are corrected, generating the final compensated interference segment data.

[0190] Embodiments of the present invention also provide a non-invasive diagnostic data acquisition system, the system comprising:

[0191] The data alignment module is used to acquire raw optical acquisition data and raw contact pressure data, and perform time alignment processing to form aligned data;

[0192] The contact pressure change identification module is used to perform contact pressure change identification processing based on the original contact pressure data in the alignment data, and generate contact pressure influence domain data through the contact pressure change amplitude and contact pressure change rate.

[0193] The optical path disturbance positioning module is used to perform optical path disturbance positioning processing on the original optical acquisition data in the aligned data based on the contact pressure influence domain data. By identifying the light intensity slope change that occurs synchronously with the contact pressure change, it generates optical interference fragment data.

[0194] The optical compensation construction module is used to perform optical compensation construction processing on the corresponding optical interference segments based on the optical interference segment data. It generates regional compensation optical data by considering the change direction and amplitude of the interference segments and the pressure change trend of the region affected by the contact pressure of the interference segments.

[0195] The domain stitching module is used to perform domain stitching processing between each contact pressure-affected domain based on the domain compensation optical data. It generates continuous correction optical data by using the light intensity difference, phase difference, and waveform period difference between adjacent domains.

[0196] The target data extraction module is used to perform target data extraction processing based on continuously corrected optical data to generate non-invasive diagnostic data.

[0197] It should be noted that this system is a system corresponding to the above method. All implementation methods in the above method embodiments are applicable to this embodiment and can achieve the same technical effect.

[0198] Embodiments of the present invention also provide a computing device, including: a processor and a memory storing a computer program, wherein the computer program, when executed by the processor, performs the method described above. All implementations in the above method embodiments are applicable to this embodiment and can achieve the same technical effects.

[0199] Embodiments of the present invention also provide a computer-readable storage medium storing instructions that, when executed on a computer, cause the computer to perform the method described above. All implementations in the above method embodiments are applicable to this embodiment and can achieve the same technical effects.

[0200] This case study applies to wearable non-invasive diagnostic devices. Addressing users' needs for collecting heart rate and blood oxygen saturation data during light exercise (such as walking or slight hand movements), the "contact pressure-optical signal linkage correction" technology of this invention eliminates motion interference. The corrected diagnostic data is then matched with a standard physiological database in the cloud to output health assessment suggestions, solving the problems of data distortion and low diagnostic reference value during exercise in traditional devices.

[0201] Core equipment and data sources for the case study:

[0202] 1. Data Acquisition Device: A smart bracelet integrating a near-infrared optical sensor (emission wavelength 660nm / 940nm, sampling frequency 50Hz) and a thin-film pressure sensor (sampling frequency 50Hz, measurement range 0-5N), worn at the base of the user's index finger.

[0203] 2. Cloud-based standard database: Stores physiological data of 100,000+ healthy individuals and individuals with common cardiovascular diseases, including: normal heart rate range (60-100 beats / min at rest, 90-130 beats / min during light exercise); normal blood oxygen saturation range (95%-100%); typical pulse waveform templates (grouped by age: 20-30 years, 31-50 years, and over 51 years, including characteristic parameters such as peak / trough ratio, rising slope, and cycle stability); and data quality assessment thresholds after interference correction (e.g., waveform similarity ≥85% is considered a valid match).

[0204] 3. Test subject: A 35-year-old male user who engages in light daily exercise (walking for 30 minutes daily), has no known cardiovascular disease, and is walking during the test (walking speed 4km / h), with a data collection time of 30 seconds.

[0205] Case implementation steps (in conjunction with the technical solution of this invention):

[0206] Step 1: Raw data acquisition and time alignment

[0207] An optical sensor collects the light intensity signal reflected from the finger skin, generating raw optical acquisition data (30 seconds × 50Hz = 1500 sampling points, including light intensity fluctuation signals).

[0208] The pressure sensor collects contact pressure data between the wristband and the skin, generating raw contact pressure data (1500 sampling points, with pressure fluctuating between 0.8N and 2.5N due to walking).

[0209] According to the method of the present invention, the pressure data is resampled (the sampling interval is adjusted to 20ms) based on the sampling clock of the optical sensor to form time-aligned data (each optical sampling point corresponds to one pressure sampling point).

[0210] Step 2: Contact pressure domain division and optical interference localization

[0211] Calculate the difference between adjacent sampling points for the aligned pressure data to generate pressure change amplitude data (maximum change amplitude 0.6N); calculate pressure change rate data (maximum rate 0.3N / ms).

[0212] The preset pressure change amplitude threshold is 0.2N and the rate threshold is 0.1N / ms. Three consecutive candidate pressure change segments are selected (corresponding to pressure fluctuations caused by the user's hand swinging while walking, each segment lasting 1.2-1.8 seconds).

[0213] After being split according to time continuity, three contact pressure influence domains are generated (domain 1: 0-1.5 seconds, domain 2: 8.2-9.8 seconds, domain 3: 22.5-24.0 seconds).

[0214] Optical data is extracted within each domain, the slope of light intensity change is calculated, and the synchronization with the pressure change trend (rising / falling) is compared to locate 3 optical interference segments (completely synchronized with the pressure change segment, with a sudden change in light intensity slope ≥0.05 lux / ms).

[0215] Step 3: Optical compensation and domain stitching

[0216] For each interfering segment:

[0217] 1. Determine the direction of interference (the light intensity in domain 1 shows a decreasing trend, and the compensation trend is set to increase).

[0218] 2. Extract the interference amplitude (maximum decrease of light intensity in domain 1 is 0.8 lux), distribute the compensation amount proportionally according to the pressure change trend (first increase then decrease), and generate target compensation amount data;

[0219] 3. Construct a compensation transition zone with 5 sampling points (0.1 seconds) at the start and end positions of the interference segment, with the compensation amount gradually changing from 0 to the target value;

[0220] By integrating compensation trend, compensation amount and transition zone data, the interference segment is compensated, and then combined with the undisturbed optical data to generate domain-compensated optical data.

[0221] The light intensity difference (maximum difference 0.1 lux), phase difference (≤5 ms), and period difference (≤0.02 seconds) between adjacent subdomains are calculated. Subdomains are stitched together by local waveform translation to generate continuous corrected optical data (without obvious boundary abrupt changes and smooth waveform).

[0222] Step 4: Extraction of target diagnostic data

[0223] Extracting pulse waveform features from continuously corrected optical data:

[0224] Heart rate: The number of pulse cycles within 30 seconds is counted and calculated to be 105 beats / minute;

[0225] Blood oxygen saturation: 97% was calculated based on the 660nm / 940nm light intensity ratio and the corrected waveform peak ratio.

[0226] Pulse waveform characteristic parameters: peak / trough ratio = 3.2, rising slope = 0.08 lux / ms, period stability coefficient of variation = 3.5%;

[0227] The above parameter combination constitutes the final non-invasive diagnostic data, which is then uploaded to the cloud.

[0228] Step 5: Cloud-based data matching and health advice generation

[0229] 1. Matching Dimensions and Algorithms:

[0230] Dimension 1: Heart rate zone matching (105 beats / min for users, and 90-130 beats / min for men aged 31-50 under light exercise conditions in the cloud).

[0231] Dimension 2: Blood oxygen saturation matching (97%, cloud normal range 95%-100%);

[0232] Dimension 3: Pulse waveform similarity matching (using the Dynamic Time Warping (DTW) algorithm, the user's corrected pulse waveform is compared with the standard waveform template of a 31-50 year old male during light exercise in the cloud, with a similarity of 92%).

[0233] Matching degree calculation: weighted score (heart rate 30% + blood oxygen 30% + waveform similarity 40%), final matching degree = 93%.

[0234] 2. Matching degree classification and recommendations:

[0235] Grading criteria: ≥85% is "high match", 60%-84% is "medium match", and <60% is "low match";

[0236] Matching result: High match (93%);

[0237] Health advice:

[0238] "Your current heart rate (105 beats / min), blood oxygen saturation (97%), and pulse waveform during light exercise are all within the range of healthy men aged 31-50, with no obvious abnormalities. It is recommended to maintain the current exercise intensity and continue monitoring. If your heart rate continues to exceed 130 beats / min or your blood oxygen level drops below 95% during exercise, please stop exercising and rest. If necessary, go to the hospital for a cardiovascular function test."

[0239] Case study technical effectiveness verification:

[0240] 1. Interference correction effect: Before correction, optical interference caused by exercise resulted in a heart rate calculation error of ±15 beats / min and blood oxygen saturation fluctuation of ±3%; after correction, the heart rate error was ≤±2 beats / min, the blood oxygen fluctuation was ≤±0.5%, and the data stability was improved by more than 85%.

[0241] 2. Matching accuracy: Because the corrected data eliminates motion interference, the matching degree with the standard data in the cloud increases from 62% before correction to 93%, effectively avoiding "misjudgment caused by interference data";

[0242] 3. Practicality: It is suitable for light exercise scenarios, solving the pain point of "data unavailable during exercise" in traditional wearable devices, and providing users with real-time and reliable health assessment references.

[0243] The above description represents the preferred embodiments of the present invention. It should be noted that those skilled in the art can make various improvements and modifications without departing from the principles of the present invention, and these improvements and modifications should also be considered within the scope of protection of the present invention.

Claims

1. A method for acquiring non-invasive diagnostic data, characterized in that, The method includes: The raw optical acquisition data and raw contact pressure data are acquired and time-aligned to form aligned data. Based on the original contact pressure data in the alignment data, perform contact pressure change identification processing, and generate contact pressure influence domain data by the contact pressure change amplitude and contact pressure change rate. Based on the domain data affected by contact pressure, optical path disturbance localization processing is performed on the original optical acquisition data in the alignment data. By identifying the changes in light intensity slope that occur synchronously with the changes in contact pressure, optical interference fragment data is generated. Based on the optical interference fragment data, optical compensation construction processing is performed on the corresponding optical interference fragments. By analyzing the change direction and amplitude of the interference fragments and the pressure change trend of the domain affected by the contact pressure of the corresponding fragments, domain compensation optical data is generated. Based on the domain-compensated optical data, domain-segmentation processing is performed between the domains affected by each contact pressure. Continuous correction optical data is generated by using the light intensity difference, phase difference, and waveform period difference between adjacent domains. Based on continuously calibrated optical data, target data extraction and processing are performed to generate non-invasive diagnostic data.

2. The method for acquiring non-invasive diagnostic data according to claim 1, characterized in that, Based on the original contact pressure data in the aligned data, contact pressure change identification processing is performed. Contact pressure influence domain data is generated by analyzing the magnitude and rate of contact pressure change, including: The difference between adjacent sampling points in the original contact pressure data in the alignment data is calculated to generate pressure change amplitude data. Based on the pressure change amplitude data, the rate of change is calculated according to the change in adjacent time intervals to generate pressure change rate data; Based on the pressure change amplitude data and pressure change rate data, compare them with the preset pressure change amplitude threshold and preset pressure change rate threshold respectively, and record the continuous time segment that meets the matching condition of any preset threshold in the comparison results as candidate pressure change segment data. Based on the temporal continuity of the candidate pressure change segment data, candidate pressure change segments with time intervals exceeding a preset interval threshold are split into multiple pressure change segment data. Based on the start and end points of the pressure change segment data, corresponding segment boundary point data are generated, and contact pressure influence domain data are generated based on the segment boundary point data.

3. The method for acquiring non-invasive diagnostic data according to claim 1, characterized in that, Based on the contact pressure-affected domain data, optical path disturbance localization processing is performed on the raw optical acquisition data in the aligned data. By identifying changes in light intensity slope that occur synchronously with changes in contact pressure, optical interference fragment data is generated, including: Based on the influence of contact pressure on the domain data, the corresponding original optical acquisition data is extracted in each domain to form domain optical data. The light intensity difference between adjacent sampling points of the domain optical data is calculated to generate light intensity change slope data. By comparing the slope data of light intensity change with the pressure change trend in the corresponding domain, synchronous data of light intensity change points are generated. Based on the continuity of the synchronous change point data of light intensity on the time axis, the continuous synchronous change points are combined into optical disturbance segment data. Based on the slope change of the optical disturbance segment data at the boundary, the interference boundary points are identified, and interference segment boundary point data is generated. Optical interference segment data is generated based on the optical disturbance segment data and the boundary point data of the interference segment.

4. The method for acquiring non-invasive diagnostic data according to claim 1, characterized in that, Based on the optical interference fragment data, optical compensation construction processing is performed on the corresponding optical interference fragments. By analyzing the direction and magnitude of the interference fragment changes and the pressure change trend of the corresponding contact pressure, domain-specific compensated optical data is generated, including: The change direction data of each interference segment is determined based on the optical interference segment data, and compensation trend data opposite to its direction is generated based on the change direction data. Based on the variation amplitude data of the interference segments in the optical interference segment data, the target compensation amount data is generated by uniformly distributing the data according to the time length or proportionally distributing it according to the pressure change trend. Based on the target compensation amount data, compensation transition zone data is constructed at the start and end positions of the interference segment. The compensation amount is gradually increased or decreased from zero within the compensation transition zone to generate compensation transition data. Compensation processing is performed on the optical interference segments based on compensation trend data, target compensation amount data, and compensation transition data to generate compensated interference segment data. The compensated interference fragment data is combined with the undisturbed optical acquisition data in a time series to generate domain-compensated optical data.

5. The method for acquiring non-invasive diagnostic data according to claim 4, characterized in that, Based on the variation amplitude data of the interference segments in the optical interference segment data, the target compensation amount data is generated by uniformly distributing the data according to the time length or proportionally distributing it according to the pressure change trend, including: Extract the variation amplitude data of the interference segments from the optical interference segment data, and generate the time interval data of the interference segments based on the start and end times of the interference segments; Based on the time interval data of the interference segment, the change amplitude data is divided into multiple compensation sub-regions in chronological order, so that each compensation sub-region corresponds to a time distribution segment of the change amplitude data, thus generating compensation sub-region time period data; Based on the time period data of the compensation sub-region, the change amplitude data is evenly distributed over time and mapped to each compensation sub-region to generate uniformly distributed compensation amount data; Based on the pressure change trend of the contact pressure influence domain to which the interference segment belongs, the change amplitude data is distributed according to the relationship of the pressure change trend in each compensation sub-region to generate trend compensation amount data. Based on the compatibility conditions between uniformly distributed compensation data and trend compensation data, a compensation distribution method that better matches the waveform recovery characteristics of the interference segment is selected from the two to generate target compensation data.

6. The method for acquiring non-invasive diagnostic data according to claim 4, characterized in that, Based on the target compensation amount data, compensation transition zone data is constructed at the start and end positions of the interference segment. The compensation amount gradually increases or decreases from zero within the compensation transition zone, generating compensation transition data, including: The starting point and ending point data of the interference segment are determined based on the optical interference segment data, and the boundary point data of the interference segment are generated based on the starting point and ending point data. Based on the boundary point data of the interference segment, the boundary compensation area is divided into several compensation transition sub-regions, and compensation transition sub-region data is generated. Based on the compensation transition sub-region data, the target compensation amount data is gradually processed in a pattern of gradually increasing from zero to the target compensation amount or gradually decreasing from the target compensation amount, to generate compensation transition sequence data. Based on the compensation gradient sequence data, continuous adjustment processing is performed between adjacent sampling points to ensure that the compensation amount changes smoothly point by point within the compensation transition zone, thereby generating compensation gradient adjustment data; Based on the compensation gradual adjustment data, construct complete compensation transition zone data.

7. The method for acquiring non-invasive diagnostic data according to claim 4, characterized in that, Based on the compensation trend data, target compensation amount data, and compensation transition data, compensation processing is performed on the optical interference segments to generate compensated interference segment data, including: Based on the compensation trend data, the optical interference segment data is subjected to orientation correction processing, and the compensation direction is adjusted to be opposite to the change direction of the interference segment to generate orientation correction data. Based on the orientation correction data, the target compensation data is used to perform compensation allocation processing on each sampling point inside the interference segment. The compensation amount is allocated according to the time sequence of the sampling point position to generate internal compensation allocation data. Based on the internal compensation allocation data and compensation transition data, a compensation gradient superposition process is performed at the boundary of the interference segment to ensure that the compensation amount maintains a continuous transition between the internal part and the boundary part of the interference segment, thereby generating compensation superposition structure data. Compensation and fusion processing is performed on the entire interference segment based on the compensated overlay structure data. By fusing the compensated overlay structure data with the original interference segment data in the time series, compensated interference segment data is generated.

8. A non-invasive diagnostic data acquisition system, characterized in that, The system, used in any one of claims 1 to 7, comprises: The data alignment module is used to acquire raw optical acquisition data and raw contact pressure data, and perform time alignment processing to form aligned data; The contact pressure change identification module is used to perform contact pressure change identification processing based on the original contact pressure data in the alignment data, and generate contact pressure influence domain data through the contact pressure change amplitude and contact pressure change rate. The optical path disturbance positioning module is used to perform optical path disturbance positioning processing on the original optical acquisition data in the aligned data based on the contact pressure influence domain data. By identifying the light intensity slope change that occurs synchronously with the contact pressure change, it generates optical interference fragment data. The optical compensation construction module is used to perform optical compensation construction processing on the corresponding optical interference segments based on the optical interference segment data. It generates regional compensation optical data by considering the change direction and amplitude of the interference segments and the pressure change trend of the region affected by the contact pressure of the interference segments. The domain stitching module is used to perform domain stitching processing between each contact pressure-affected domain based on the domain compensation optical data. It generates continuous correction optical data by using the light intensity difference, phase difference, and waveform period difference between adjacent domains. The target data extraction module is used to perform target data extraction processing based on continuously corrected optical data to generate non-invasive diagnostic data.

9. A computing device, characterized in that, include: One or more processors; A storage device for storing one or more programs, which, when executed by one or more processors, cause the one or more processors to implement the method as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a program that, when executed by a processor, implements the method as described in any one of claims 1 to 7.

Images