Hemodialysis data processing method and systemTechnical Field
The invention relates to the technical field of hemodialysis, in particular to a data processing method and system for hemodialysis.
Background
Hemodialysis is one of kidney replacement treatment modes of patients with acute and chronic renal failure, and is characterized in that in-vivo blood is drained to the outside of the body, and is subjected to substance exchange with electrolyte solution (dialysate) with similar concentration of an organism inside and outside the hollow fiber through a dialyzer consisting of innumerable hollow fibers, so that metabolic wastes in the body are removed, balance between electrolyte and acid and alkali is maintained, excessive moisture in the body is removed, and the whole process of back transfusion of the purified blood is called hemodialysis;
In the traditional hemodialysis data processing process, due to the fact that dialysis instruments are operated for a long time, hardware parts of the dialysis instruments often have faults, collected related instrument data are inaccurate, physiological characteristics of different patients are different, collected data are affected by environmental factors of a dialysis room, error sources are inconvenient to identify timely and accurately in the existing hemodialysis data processing process, due to the fact that effective error identification is lacking, all errors are often treated as the same type in the traditional method, targeted correction measures cannot be adopted according to different error sources, and when the situation of complex and changeable errors is faced, deviation of data processing results is easy to occur, and accurate assessment of dialysis sufficiency is affected.
Accordingly, the present invention provides a data processing method and system for hemodialysis to overcome and ameliorate the shortcomings of the prior art.
Disclosure of Invention
In order to solve the technical problems, the invention provides a data processing method and a system for hemodialysis, which are used for solving the corresponding technical problems in the background art.
In order to achieve the above purpose, the invention adopts the technical scheme that the data processing method for hemodialysis comprises the following steps:
Firstly, collecting instrument core parameters in real time through a sensor arranged in a dialysis instrument, collecting physiological data of a patient in real time through a circulating pipeline sensor arranged outside the patient, collecting environmental interference data in real time through a sensor arranged in a dialysis chamber, and performing standardized processing;
Acquiring acquisition frequency of instrument core parameters, acquisition frequency of patient physiological data and acquisition frequency of environmental interference data, performing cubic spline interpolation on low acquisition frequency parameters according to the highest acquisition frequency parameters serving as a reference time axis, generating virtual sampling values at the time points of the highest acquisition frequency parameters, and associating the virtual sampling values of the instrument core parameters, the virtual sampling values of the patient physiological data and the virtual sampling values of the environmental interference data at the same time point to form a time sequence matrix;
Defining an error type as an internal error caused by hardware faults of a dialysis instrument and an external error caused by environmental interference or individual differences of a patient, performing autocorrelation analysis on a virtual sampling value of a core parameter of the instrument based on a time sequence matrix, marking a high-frequency fault point in combination with a history maintenance record of the dialysis instrument, performing cross-correlation analysis on the virtual sampling value of environmental interference data and the virtual sampling value of physiological data of the patient, calculating an interference coupling coefficient, performing longitudinal trend analysis on data of the same patient for multiple treatments, identifying individual specific interference characteristics, and determining an error source according to an autocorrelation analysis result, a cross-correlation analysis result and a longitudinal trend analysis result;
And fourthly, carrying out layered correction on instrument core parameters, patient physiological data and environmental interference data based on error sources, triggering a dialysis instrument self-checking program on internal errors, adopting self-adaptive filtering on external errors, establishing a patient error model on individual specific interference, and outputting calibrated data.
Preferably, the instrument core parameters comprise blood pump rotation speed, transmembrane pressure and ultrafiltration rate, the patient physiological data comprise systolic pressure, venous pressure and activated clotting time, and the environment interference data comprise room temperature, power frequency noise intensity and illumination intensity.
Preferably, the specific process of marking the high-frequency fault point is as follows:
S101, acquiring acquisition frequency of instrument core parameters, acquisition frequency of patient physiological data and acquisition frequency of environmental interference data, taking a time axis corresponding to the highest acquisition frequency parameter as a reference time axis, wherein time points are sequentially as follows;
The actual collection time point and the corresponding value of the low collection frequency parameter are taken as known data points, a cubic spline function is constructed between every two adjacent data points, and based on the constructed cubic spline function, the corresponding function value is calculated at each time point of the reference time axis, so that a virtual sampling value of the low collection frequency parameter on the reference time axis is obtained;
Combining the virtual sampling value of the instrument core parameter, the virtual sampling value of the patient physiological data and the virtual sampling value of the environment interference data at the same time point to form a data row, and sequentially arranging the data rows corresponding to all the time points according to the time sequence of a reference time axis to form a time sequence matrix;
s102, calculating autocorrelation coefficients under different hysteresis orders k (k=0, 1,2,) according to the virtual sampling values of the instrument core parametersThe calculation formula is as follows:
;
Wherein, theRepresenting a virtual sampling value of the instrument core parameter at an ith time point in the time sequence matrix;
representing the average value of the virtual sampling values of the instrument core parameter at all time points;
respectively calculating autocorrelation coefficients of the blood pump rotating speed a, the transmembrane pressure b and the ultrafiltration rate c under different hysteresis orders to obtain a group of autocorrelation coefficient sequences respectively、And;
S103, according to a preset autocorrelation coefficient thresholdFor each instrument core parameter, when the autocorrelation coefficients are the autocorrelation coefficient sequencesIs greater than the absolute value ofAt the time, the abnormal correlation exists between the value of the instrument core parameter and the value of the instrument core parameter at other time points under the hysteresis order k, and the corresponding hysteresis order k and the time point are recorded, wherein,Is the starting position of the point in time involved in the current calculation of the autocorrelation coefficients,Is relative to the starting time pointA time point after k time units;
Will exceedIs preliminarily marked as potential fault points to form a preliminary fault point set;
S104, acquiring a history maintenance record from a maintenance database of the dialysis instrument, and for the primary fault point setEach time point of (3)If inBefore and after the maintenance record related to the hardware fault of the dialysis instrument exists, the time point is thenFrom a set of preliminary fault pointsRemoving, screening to obtain a screened fault point set;
For the filtered fault point setEach time point of (3)Marking is performed in the time series matrix, and the time series matrix marked with the high-frequency fault points is output.
Preferably, the specific process of calculating the interference coupling coefficient is as follows:
S201, obtaining and combining virtual sampling values of environment interference data and virtual sampling values of patient physiological data at the same time point from a time sequence matrix, and respectively calculating cross-correlation coefficients of the virtual sampling values of the environment interference data and the virtual sampling values of the patient physiological data under different hysteresis orders m (m=0, 1, 2.)The calculation formula is as follows:
;
Wherein, theRepresenting a virtual sampling value of the environment interference data at a q-th time point in the time sequence matrix;
representing the average value of the virtual sampling values of the environmental interference data at all n time points;
representing virtual sampling values of the physiological data of the patient at the (q+m) th time point in the time sequence matrix;
representing virtual sample values of the patient physiological data at a q-th time point in the time series matrix;
Representing an average of the virtual sample values of the patient physiological data at all n time points;
Calculating the cross-correlation coefficients of the room temperature and the systolic pressure, the room temperature and the venous pressure, the room temperature and the activated coagulation time, the power frequency noise intensity and the systolic pressure, the power frequency noise intensity and the venous pressure, the power frequency noise intensity and the activated coagulation time, the illumination intensity and the systolic pressure, the illumination intensity and the venous pressure and the activated coagulation time under different hysteresis orders respectively to obtain a plurality of groups of cross-correlation coefficient sequences;
s202, according to a preset cross-correlation coefficient threshold valueFor each group of calculated cross-correlation coefficient sequences, finding out the value with the maximum absolute value of the cross-correlation coefficient and the corresponding hysteresis order, if the maximum absolute value exceeds the cross-correlation coefficient threshold valueIndicating that there is a high correlation between the ambient interference data and the patient physiological data, the maximum absolute value is taken as the interference coupling coefficient of the corresponding ambient interference data and the patient physiological data combination.
Preferably, the specific process of identifying individual specific interference features is as follows:
S301, obtaining virtual sampling values of core parameters of the same patient instrument at all time points, virtual sampling values of physiological data of the patient at all time points and virtual sampling values of environmental interference data at all time points through a time sequence matrix, and sorting data of each treatment into a matrix form according to the time sequence of the treatment, wherein each row represents one treatment, and each column corresponds to the values of the various data at different time points respectively;
S302, calculating the mean value of the physiological data of each patient aiming at the physiological data of each patient in each treatment of the patientThe calculation formula is as follows:
;
wherein u is the virtual sampling point number of the physiological data of the patient in the treatment process;
is the p-th virtual sampling point value;
the contraction pressure, the venous pressure and the activated clotting time of each treatment of the patient are calculated respectively to obtain the average value of physiological data of each patient in each treatment;
s303, calculating the variance of the physiological data of each patient in each treatmentThe calculation formula is as follows:
;
The contraction pressure, the venous pressure and the activated clotting time of each treatment are respectively calculated to obtain the variance of physiological data of each patient in each treatment;
S304, comparing the mean value and the variance of different treatment times of the same patient, and if the mean value and the variance combination of the physiological data of the patient repeatedly appear in multiple treatments, primarily determining the combination as an individual specific interference characteristic;
The treatment times are taken as an abscissa, the mean value of physiological data of a patient is taken as an ordinate, a trend graph of the mean value changing along with the treatment times is drawn, the change trend of the mean value is observed, a linear regression method is used for fitting a change curve of the mean value along with the treatment times, a regression coefficient is calculated, and if the regression coefficient is not zero, the fact that the mean value has the linear change trend is indicated;
taking the treatment times as an abscissa, taking the variance of the physiological data of the patient as an ordinate, drawing a trend graph of variance variation along with the treatment times, observing the variation trend of the variance, observing the periodical variation of the variance along with the increase of the treatment times by adopting a moving average method, and if the physiological data of the patient gradually rises along with the treatment times and the variance of the physiological data of the patient also shows the increasing trend, and the variation mode continuously exists in multiple treatments, further determining that the combination is an individual specific interference characteristic.
Preferably, the specific process of determining the error source is as follows:
S401, acquiring a high-frequency fault point of an autocorrelation analysis mark, combining a dialysis instrument historical maintenance record, and determining an internal error caused by hardware fault of the dialysis instrument if abnormal change exists in a core parameter of the instrument at a certain time point and the relevant maintenance record exists near the time point in the historical maintenance record;
s402, acquiring an interference coupling coefficient obtained through cross-correlation analysis and calculation, selecting a median as a preset standard threshold based on a distribution range of the interference coupling coefficient between the environmental interference data and the physiological data of the patient in a normal state, and determining an external error caused by the environmental interference if the interference coupling coefficient between the environmental interference data and the physiological data of the patient is greater than the preset standard threshold;
S403, acquiring individual specific interference characteristics identified by longitudinal trend analysis, and determining an external error caused by individual differences of patients if the instrument core parameters are not changed abnormally and no relevant maintenance record exists at the time point, and meanwhile, the interference coupling coefficient between the environmental interference data and the physiological data of the patients does not exceed a preset standard threshold value, and the individual specific interference characteristics exist in the physiological data of the patients.
Preferably, the specific process of outputting the calibrated data is as follows:
s501, triggering a dialysis instrument self-checking program when determining that the error source is an internal error caused by the hardware fault of the dialysis instrument, maintaining the hardware of the dialysis instrument and replacing a fault part, and processing environment interference data by adopting an adaptive filtering method when determining that the error source is an external error caused by the environment interference;
S502, when determining that the error source is an external error caused by individual difference of a patient, establishing a patient error model for individual specific interference by a regression analysis method by taking patient age, weight, environmental interference data and instrument core parameters as independent variables and patient physiological data as dependent variables, wherein the formula is as follows:
;
Wherein Z is patient physiological data;
Age of the patient;
is the weight of the patient;
Y is environmental interference data;
x is an instrument core parameter;
are regression coefficients;
is an error term;
Training and optimizing a patient error model through historical patient physiological data, historical environment interference data and historical instrument core parameters, and determining a regression coefficient value to obtain an updated patient error model;
S503, according to the processing results of different error types, correcting the instrument core parameters, the environmental interference data and the physiological data of the patient, recording the corrected data in a data record file, marking corresponding positions in a time sequence matrix, and outputting the corrected data according to a specified format.
A data processing system for hemodialysis comprises a data acquisition unit, a data processing unit, an error analysis unit and an error calibration unit;
the data acquisition unit is used for acquiring instrument core parameters in real time through a sensor arranged in the dialysis instrument, acquiring physiological data of a patient in real time through a circulating pipeline sensor arranged outside the patient, acquiring environmental interference data in real time through a sensor arranged in the dialysis chamber and carrying out standardized processing;
The data processing unit is used for acquiring the acquisition frequency of the instrument core parameters, the patient physiological data and the environment interference data, performing cubic spline interpolation on the low acquisition frequency parameters according to the highest acquisition frequency parameters serving as a reference time axis, generating virtual sampling values at the time points of the highest acquisition frequency parameters, and associating the virtual sampling values of the instrument core parameters, the virtual sampling values of the patient physiological data and the virtual sampling values of the environment interference data at the same time point to form a time sequence matrix;
The error analysis unit is used for defining the error type as an internal error caused by hardware faults of a dialysis instrument and an external error caused by environmental interference or individual difference of a patient, carrying out autocorrelation analysis on a virtual sampling value of a core parameter of the instrument based on a time sequence matrix, detecting periodic deviation, marking a high-frequency fault point in combination with a history maintenance record of the dialysis instrument, carrying out cross-correlation analysis on the virtual sampling value of environmental interference data and physiological data of the patient, calculating an interference coupling coefficient, carrying out longitudinal comparison on data of multiple treatments of the same patient, identifying individual specific interference, and determining an error source according to an autocorrelation analysis result, a cross-correlation analysis result and a longitudinal comparison result;
the error calibration unit is used for carrying out layered correction on instrument core parameters, patient physiological data and environment interference data based on error sources, triggering a dialysis instrument self-checking program on internal errors, adopting self-adaptive filtering on external errors, establishing a patient-specific error model on individual specific interference, and outputting calibrated data.
Compared with the prior art, the invention has the beneficial effects that: the method comprises the steps of acquiring instrument core parameters, patient physiological data and environmental interference data, acquiring acquisition frequencies of the instrument core parameters, the patient physiological data and the environmental interference data, performing tertiary spline interpolation on low acquisition frequency parameters according to the highest acquisition frequency parameters serving as a reference time axis, generating virtual sampling values at the highest acquisition frequency parameter time points, correlating the virtual sampling values of the instrument core parameters, the patient physiological data and the environmental interference data at the same time points to form a time sequence matrix, defining error types as internal errors caused by dialysis instrument hardware faults and external errors caused by environmental interference or patient individual differences, performing autocorrelation analysis on the virtual sampling values of the instrument core parameters, combining dialysis instrument history maintenance records, marking high-frequency fault points, performing cross-correlation analysis on the virtual sampling values of the environmental interference data and the virtual sampling values of the patient physiological data, calculating interference coupling coefficients, performing longitudinal trend analysis on data of the same patient multi-treatment, identifying individual specific interference characteristics, determining error sources according to the autocorrelation analysis results, the error sources, the dialysis instrument core parameters, the patient physiological data and the environmental interference data, and the external error caused by the individual difference, performing self-correlation analysis on the dialysis instrument, setting up an accurate correction on the hemodialysis instrument, and performing self-adaptive processing on the hemodialysis instrument, improving the blood error, correcting the blood error, and improving the blood-dialysis instrument accuracy, and the blood-dialysis data, and the self-correction model, and improving the accuracy, the safety of hemodialysis treatment is improved, and the risk of complications is reduced.
Drawings
FIG. 1 is a schematic diagram of the overall structure of a preferred embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
Embodiments of the invention:
Referring to fig. 1, a data processing method for hemodialysis includes the following steps:
Firstly, collecting instrument core parameters in real time through a sensor arranged in a dialysis instrument, collecting physiological data of a patient in real time through a circulating pipeline sensor arranged outside the patient, collecting environmental interference data in real time through a sensor arranged in a dialysis chamber, and performing standardized processing;
Acquiring acquisition frequency of instrument core parameters, acquisition frequency of patient physiological data and acquisition frequency of environmental interference data, performing cubic spline interpolation on low acquisition frequency parameters according to the highest acquisition frequency parameters serving as a reference time axis, generating virtual sampling values at the time points of the highest acquisition frequency parameters, and associating the virtual sampling values of the instrument core parameters, the virtual sampling values of the patient physiological data and the virtual sampling values of the environmental interference data at the same time point to form a time sequence matrix;
Defining an error type as an internal error caused by hardware faults of a dialysis instrument and an external error caused by environmental interference or individual differences of a patient, performing autocorrelation analysis on a virtual sampling value of a core parameter of the instrument based on a time sequence matrix, marking a high-frequency fault point in combination with a history maintenance record of the dialysis instrument, performing cross-correlation analysis on the virtual sampling value of environmental interference data and the virtual sampling value of physiological data of the patient, calculating an interference coupling coefficient, performing longitudinal trend analysis on data of the same patient for multiple treatments, identifying individual specific interference characteristics, and determining an error source according to an autocorrelation analysis result, a cross-correlation analysis result and a longitudinal trend analysis result;
And fourthly, carrying out layered correction on instrument core parameters, patient physiological data and environmental interference data based on error sources, triggering a dialysis instrument self-checking program on internal errors, adopting self-adaptive filtering on external errors, establishing a patient error model on individual specific interference, and outputting calibrated data.
The instrument core parameters comprise blood pump rotation speed, transmembrane pressure and ultrafiltration rate, the physiological data of patients comprise systolic pressure, venous pressure and activated clotting time, and the environmental interference data comprise room temperature, power frequency noise intensity and illumination intensity.
The specific process of marking the high-frequency fault point is as follows:
S101, acquiring acquisition frequency of instrument core parameters, acquisition frequency of patient physiological data and acquisition frequency of environmental interference data, taking a time axis corresponding to the highest acquisition frequency parameter as a reference time axis, wherein time points are sequentially as follows;
The actual collection time point and the corresponding value of the low collection frequency parameter are taken as known data points, a cubic spline function is constructed between every two adjacent data points, and based on the constructed cubic spline function, the corresponding function value is calculated at each time point of the reference time axis, so that a virtual sampling value of the low collection frequency parameter on the reference time axis is obtained;
Combining the virtual sampling value of the instrument core parameter, the virtual sampling value of the patient physiological data and the virtual sampling value of the environment interference data at the same time point to form a data row, and sequentially arranging the data rows corresponding to all the time points according to the time sequence of a reference time axis to form a time sequence matrix;
s102, calculating autocorrelation coefficients under different hysteresis orders k (k=0, 1,2,) according to the virtual sampling values of the instrument core parametersThe calculation formula is as follows:
;
Wherein, theRepresenting a virtual sampling value of the instrument core parameter at an ith time point in the time sequence matrix;
representing the average value of the virtual sampling values of the instrument core parameter at all time points;
respectively calculating autocorrelation coefficients of the blood pump rotating speed a, the transmembrane pressure b and the ultrafiltration rate c under different hysteresis orders to obtain a group of autocorrelation coefficient sequences respectively、And;
S103, according to a preset autocorrelation coefficient thresholdFor each instrument core parameter, when the autocorrelation coefficients are the autocorrelation coefficient sequencesIs greater than the absolute value ofAt the time, the abnormal correlation exists between the value of the instrument core parameter and the value of the instrument core parameter at other time points under the hysteresis order k, and the corresponding hysteresis order k and the time point are recorded, wherein,Is the starting position of the point in time involved in the current calculation of the autocorrelation coefficients,Is relative to the starting time pointA time point after k time units;
Will exceedIs preliminarily marked as potential fault points to form a preliminary fault point set;
S104, acquiring a history maintenance record from a maintenance database of the dialysis instrument, and for the primary fault point setEach time point of (3)If inBefore and after the maintenance record related to the hardware fault of the dialysis instrument exists, the time point is thenFrom a set of preliminary fault pointsRemoving, screening to obtain a screened fault point set;
For the filtered fault point setEach time point of (3)Marking is performed in the time series matrix, and the time series matrix marked with the high-frequency fault points is output.
The specific process of calculating the interference coupling coefficient is as follows:
S201, obtaining and combining virtual sampling values of environment interference data and virtual sampling values of patient physiological data at the same time point from a time sequence matrix, and respectively calculating cross-correlation coefficients of the virtual sampling values of the environment interference data and the virtual sampling values of the patient physiological data under different hysteresis orders m (m=0, 1, 2.)The calculation formula is as follows:
;
Wherein, theRepresenting a virtual sampling value of the environment interference data at a q-th time point in the time sequence matrix;
representing the average value of the virtual sampling values of the environmental interference data at all n time points;
representing virtual sampling values of the physiological data of the patient at the (q+m) th time point in the time sequence matrix;
representing virtual sample values of the patient physiological data at a q-th time point in the time series matrix;
Representing an average of the virtual sample values of the patient physiological data at all n time points;
Calculating the cross-correlation coefficients of the room temperature and the systolic pressure, the room temperature and the venous pressure, the room temperature and the activated coagulation time, the power frequency noise intensity and the systolic pressure, the power frequency noise intensity and the venous pressure, the power frequency noise intensity and the activated coagulation time, the illumination intensity and the systolic pressure, the illumination intensity and the venous pressure and the activated coagulation time under different hysteresis orders respectively to obtain a plurality of groups of cross-correlation coefficient sequences;
s202, according to a preset cross-correlation coefficient threshold valueFor each group of calculated cross-correlation coefficient sequences, finding out the value with the maximum absolute value of the cross-correlation coefficient and the corresponding hysteresis order, if the maximum absolute value exceeds the cross-correlation coefficient threshold valueIndicating that there is a high correlation between the ambient interference data and the patient physiological data, the maximum absolute value is taken as the interference coupling coefficient of the corresponding ambient interference data and the patient physiological data combination.
The specific process of identifying individual specific interference features is as follows:
S301, obtaining virtual sampling values of core parameters of the same patient instrument at all time points, virtual sampling values of physiological data of the patient at all time points and virtual sampling values of environmental interference data at all time points through a time sequence matrix, and sorting data of each treatment into a matrix form according to the time sequence of the treatment, wherein each row represents one treatment, and each column corresponds to the values of the various data at different time points respectively;
S302, calculating the mean value of the physiological data of each patient aiming at the physiological data of each patient in each treatment of the patientThe calculation formula is as follows:
;
wherein u is the virtual sampling point number of the physiological data of the patient in the treatment process;
is the p-th virtual sampling point value;
the contraction pressure, the venous pressure and the activated clotting time of each treatment of the patient are calculated respectively to obtain the average value of physiological data of each patient in each treatment;
s303, calculating the variance of the physiological data of each patient in each treatmentThe calculation formula is as follows:
;
The contraction pressure, the venous pressure and the activated clotting time of each treatment are respectively calculated to obtain the variance of physiological data of each patient in each treatment;
S304, comparing the mean value and the variance of different treatment times of the same patient, and if the mean value and the variance combination of the physiological data of the patient repeatedly appear in multiple treatments, primarily determining the combination as an individual specific interference characteristic;
The treatment times are taken as an abscissa, the mean value of physiological data of a patient is taken as an ordinate, a trend graph of the mean value changing along with the treatment times is drawn, the change trend of the mean value is observed, a linear regression method is used for fitting a change curve of the mean value along with the treatment times, a regression coefficient is calculated, and if the regression coefficient is not zero, the fact that the mean value has the linear change trend is indicated;
taking the treatment times as an abscissa, taking the variance of the physiological data of the patient as an ordinate, drawing a trend graph of variance variation along with the treatment times, observing the variation trend of the variance, observing the periodical variation of the variance along with the increase of the treatment times by adopting a moving average method, and if the physiological data of the patient gradually rises along with the treatment times and the variance of the physiological data of the patient also shows the increasing trend, and the variation mode continuously exists in multiple treatments, further determining that the combination is an individual specific interference characteristic.
The specific process of determining the error source is as follows:
S401, acquiring a high-frequency fault point of an autocorrelation analysis mark, combining a dialysis instrument historical maintenance record, and determining an internal error caused by hardware fault of the dialysis instrument if abnormal change exists in a core parameter of the instrument at a certain time point and the relevant maintenance record exists near the time point in the historical maintenance record;
s402, acquiring an interference coupling coefficient obtained through cross-correlation analysis and calculation, selecting a median as a preset standard threshold based on a distribution range of the interference coupling coefficient between the environmental interference data and the physiological data of the patient in a normal state, and determining an external error caused by the environmental interference if the interference coupling coefficient between the environmental interference data and the physiological data of the patient is greater than the preset standard threshold;
S403, acquiring individual specific interference characteristics identified by longitudinal trend analysis, and determining an external error caused by individual differences of patients if the instrument core parameters are not changed abnormally and no relevant maintenance record exists at the time point, and meanwhile, the interference coupling coefficient between the environmental interference data and the physiological data of the patients does not exceed a preset standard threshold value, and the individual specific interference characteristics exist in the physiological data of the patients.
The specific process of outputting the calibrated data is as follows:
s501, triggering a dialysis instrument self-checking program when determining that the error source is an internal error caused by the hardware fault of the dialysis instrument, maintaining the hardware of the dialysis instrument and replacing a fault part, and processing environment interference data by adopting an adaptive filtering method when determining that the error source is an external error caused by the environment interference;
S502, when determining that the error source is an external error caused by individual difference of a patient, establishing a patient error model for individual specific interference by a regression analysis method by taking patient age, weight, environmental interference data and instrument core parameters as independent variables and patient physiological data as dependent variables, wherein the formula is as follows:
;
Wherein Z is patient physiological data;
Age of the patient;
is the weight of the patient;
Y is environmental interference data;
x is an instrument core parameter;
are regression coefficients;
is an error term;
Training and optimizing a patient error model through historical patient physiological data, historical environment interference data and historical instrument core parameters, and determining a regression coefficient value to obtain an updated patient error model;
S503, according to the processing results of different error types, correcting the instrument core parameters, the environmental interference data and the physiological data of the patient, recording the corrected data in a data record file, marking corresponding positions in a time sequence matrix, and outputting the corrected data according to a specified format.
A data processing system for hemodialysis comprises a data acquisition unit, a data processing unit, an error analysis unit and an error calibration unit;
The data acquisition unit is used for acquiring instrument core parameters in real time through a sensor arranged in the dialysis instrument, acquiring physiological data of a patient in real time through a circulating pipeline sensor arranged outside the patient, acquiring environmental interference data in real time through a sensor arranged in the dialysis chamber and performing standardized processing;
The data processing unit is used for acquiring the acquisition frequency of the instrument core parameters, the patient physiological data and the environment interference data, carrying out cubic spline interpolation on the low acquisition frequency parameters according to the highest acquisition frequency parameters serving as a reference time axis, generating virtual sampling values at the time points of the highest acquisition frequency parameters, and associating the virtual sampling values of the instrument core parameters, the virtual sampling values of the patient physiological data and the virtual sampling values of the environment interference data at the same time point to form a time sequence matrix;
The error analysis unit is used for defining the error type as an internal error caused by hardware faults of the dialysis instrument and an external error caused by environmental interference or individual difference of a patient, carrying out autocorrelation analysis on a virtual sampling value of a core parameter of the instrument based on a time sequence matrix, detecting periodic deviation, marking a high-frequency fault point in combination with a history maintenance record of the dialysis instrument, carrying out cross-correlation analysis on the virtual sampling value of environmental interference data and physiological data of the patient, calculating an interference coupling coefficient, carrying out longitudinal comparison on data of multiple treatments of the same patient, identifying individual specific interference, and determining an error source according to an autocorrelation analysis result, a cross-correlation analysis result and a longitudinal comparison result;
The error calibration unit is used for carrying out layering correction on instrument core parameters, patient physiological data and environment interference data based on error sources, triggering a dialysis instrument self-checking program on internal errors, adopting self-adaptive filtering on external errors, establishing a patient-specific error model on individual specific interference, and outputting calibrated data.
The method comprises the steps of acquiring instrument core parameters, patient physiological data and environmental interference data, acquiring acquisition frequencies of the instrument core parameters, the patient physiological data and the environmental interference data, performing tertiary spline interpolation on low acquisition frequency parameters according to the highest acquisition frequency parameters serving as a reference time axis, generating virtual sampling values at the highest acquisition frequency parameter time points, correlating the virtual sampling values of the instrument core parameters, the patient physiological data and the environmental interference data at the same time points to form a time sequence matrix, defining error types as internal errors caused by dialysis instrument hardware faults and external errors caused by environmental interference or patient individual differences, performing autocorrelation analysis on the virtual sampling values of the instrument core parameters, combining dialysis instrument history maintenance records, marking high-frequency fault points, performing cross-correlation analysis on the virtual sampling values of the environmental interference data and the virtual sampling values of the patient physiological data, calculating interference coupling coefficients, performing longitudinal trend analysis on data of the same patient multi-treatment, identifying individual specific interference characteristics, determining error sources according to the autocorrelation analysis results, the error sources, the dialysis instrument core parameters, the patient physiological data and the environmental interference data, and the external error caused by the individual difference, performing self-correlation analysis, triggering the hemodialysis instrument, improving the hemodialysis instrument self-adaptive treatment performance, improving the hemodialysis instrument self-adaptive performance, improving the hemodialysis instrument self-correction performance, and the hemodialysis instrument correction performance, and the blood-adaptive treatment model, improving the blood-quality, reducing the risk of complications.
The size of the interval and the threshold is set for the convenience of comparison, and the size of the threshold depends on the number of sample data and the number of cardinalities set for each group of sample data by a person skilled in the art, so long as the proportional relation between the parameter and the quantized value is not affected.
The formulas are all formulas with dimensions removed and numerical calculation, the formulas are formulas with a large amount of data collected for software simulation to obtain the latest real situation, and preset parameters in the formulas are set by a person skilled in the art according to the actual situation;
In the two embodiments provided in the present application, it should be understood that the disclosed apparatus and system may be implemented in other manners, for example, the apparatus embodiments described above are merely illustrative, for example, the modules are divided into only one kind of logic function, and there may be other manners of dividing actually being implemented, for example, a plurality of modules or components may be combined or may be integrated into another system, or some features may be omitted or not implemented;
The foregoing is only a preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any person skilled in the art, who is within the scope of the present invention, should make equivalent substitutions or modifications according to the technical scheme of the present invention and the inventive concept thereof, and should be covered by the scope of the present invention.