CN107625519A - ECG treatment method and device - Google Patents
ECG treatment method and device Download PDFInfo
- Publication number
- CN107625519A CN107625519A CN201710870201.0A CN201710870201A CN107625519A CN 107625519 A CN107625519 A CN 107625519A CN 201710870201 A CN201710870201 A CN 201710870201A CN 107625519 A CN107625519 A CN 107625519A
- Authority
- CN
- China
- Prior art keywords
- lead
- image
- waveform
- curve
- papery
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
The embodiments of the invention provide a kind of ECG treatment method and device, it is related to technical field of image processing.Method includes the papery electrocardiogram image read carrying out separating treatment, obtains lead curve regions image；By lead curve regions image using behavior unit separation as multiple lead wave regions images, each lead wave regions image includes multiple lead signals；Each lead wave regions image is divided into the lead waveform sub-district area image of each lead signals of multiple correspondences；Obtain the pixel position of the lead waveform in lead waveform sub-district area image；Based on lead curve regions image and the pixel position, the range value for reducing lead signals is obtained；Based on the sum of all pixels of the lead waveform in the lead waveform sub-district area image got, the phase between reduction lead signals is obtained；Based on pixel position, range value and it is described between the phase, generation reduction lead signals.Ensure the high reduction degree of reduction lead signals, it is simple to operate, more accurate.
Description
Technical field
The present invention relates to technical field of image processing, in particular to a kind of ECG treatment method and device.
Background technology
With digital medical, the development of big data, computer technology, various hospital information management systems widely should
With the patient health information database of foundation, can inquiring about and trace the health and fitness information of patient, realize that resource is total to using internet
Enjoy and expert's remote diagnosis.Electrocardiogram is the important detection methods of Diagnosing Cardiac disease, as a kind of noninvasive, accurately and rapidly
Detection method, substantial amounts of electrocardio report can be all produced in hospital daily, most electrocardio report is recorded by thermal printing paper, by
In the chemical characteristic of thermal printing paper, the reason for storage condition and time remote, printing graph color can thin out, picture quality change
Difference, or even paper sheet break so that typical clinical electrocardiosignal can not preserve for a long time.China human mortality is large number of, possesses
The pathologic electrocardiosignal of magnanimity.If the electrocardio paper shelves of reservation are reported figure is converted to by laser printer or digital camera
Piece, then by image processing techniques, the Wave data in picture is extracted, then achieve with certain form.Electrocardio after digitlization
Signal not only contributes to management of the hospital to patient information, the teaching and scientific research of clinician, and doctor can also utilize what is obtained
Mass data does the quantitative statisticses analysis and diagnosis more to become more meticulous to heart disease, excavates out more, more accurate clinical diagnosis
Standard, improve to patient's condition-inference accuracy, there is provided more preferable therapeutic scheme.
The existing digitized processing method of electrocardiograph paper curve is that the point tracing, window trace method, prediction of curve are searched
Suo Fa, the possible position of next point is mainly speculated according to the known point on electrocardio track, needs to sentence during processing
Whether the gray scale of disconnected next point and current point is close, carries out threshold process；Trace point be in change slow baseline position and
At the raising and lowering acute variation position of crest, to be handled using different model methods, these methods are substantially
All it is a kind of morphologic approximate evaluation method, complex disposal process is computationally intensive, and quantization error is larger.
The content of the invention
It is an object of the invention to provide a kind of ECG treatment method and device, to improve above mentioned problem.In order to realize
Above-mentioned purpose, the technical scheme that the present invention takes are as follows：
In a first aspect, the embodiments of the invention provide a kind of ECG treatment method, methods described includes reading
Papery electrocardiogram image carries out separating treatment, obtains the lead curve regions image of the papery electrocardiogram image；Led described
Join curve regions image using behavior unit separation as multiple lead wave regions images, each lead wave regions image bag
Include multiple lead signals；Each lead wave regions image is divided into the lead of each lead signals of multiple correspondences
Waveform sub-district area image；Obtain the pixel position of the lead waveform in the lead waveform sub-district area image；Led based on described
Join curve regions image and the pixel position, obtain the range value of reduction lead signals；Based on the lead got
The sum of all pixels of lead waveform in waveform sub-district area image, obtain the phase between the reduction lead signals；Based on the pixel
Point position, the range value and it is described between the phase, generate the reduction lead signals.
Second aspect, the embodiments of the invention provide a kind of electrocardiogram processing unit, described device includes first and obtains list
Member, separating element, cutting unit, second acquisition unit, the 3rd acquiring unit, the 4th acquiring unit and generation unit.First obtains
Unit is taken, the papery electrocardiogram image for that will read carries out separating treatment, obtains the lead of the papery electrocardiogram image
Curve regions image.Separating element, for by the lead curve regions image using behavior unit separation as multiple lead waveforms
Area image, each lead wave regions image include multiple lead signals.Cutting unit, for inciting somebody to action each lead
Wave regions image is divided into the lead waveform sub-district area image of each lead signals of multiple correspondences.Second acquisition unit,
For obtaining the pixel position of the lead waveform in the lead waveform sub-district area image.3rd acquiring unit, for based on
The lead curve regions image and the pixel position, obtain the range value of reduction lead signals.4th acquiring unit, use
The sum of all pixels of lead waveform in based on the lead waveform sub-district area image got, obtain the reduction lead letter
Phase between number.Generation unit, for based on the pixel position, the range value and it is described between the phase, generate it is described reduction lead
Join signal.
The embodiments of the invention provide a kind of ECG treatment method and device, and the papery electrocardiogram image read is entered
Row separating treatment, obtain the lead curve regions image of the papery electrocardiogram image；By the lead curve regions image with
Behavior unit separation is multiple lead wave regions images, and each lead wave regions image includes multiple lead signals；
Each lead wave regions image is divided into the lead waveform sub-district area image of each lead signals of multiple correspondences；
Obtain the pixel position of the lead waveform in the lead waveform sub-district area image；Based on the lead curve regions image and
The pixel position, obtain the range value of reduction lead signals；Based in the lead waveform sub-district area image got
Lead waveform sum of all pixels, obtain it is described reduction lead signals between the phase；Based on the pixel position, the range value
And it is described between the phase, generate the reduction lead signals.Ensure the high reduction degree of reduction lead signals, it is simple to operate, more accurate.
Other features and advantages of the present invention will illustrate in subsequent specification, also, partly become from specification
It is clear that or by implementing understanding of the embodiment of the present invention.The purpose of the present invention and other advantages can be by saying what is write
Specifically noted structure is realized and obtained in bright book, claims and accompanying drawing.
Brief description of the drawings
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below by embodiment it is required use it is attached
Figure is briefly described, it will be appreciated that the following drawings illustrate only certain embodiments of the present invention, therefore be not construed as pair
The restriction of scope, for those of ordinary skill in the art, on the premise of not paying creative work, can also be according to this
A little accompanying drawings obtain other related accompanying drawings.
Fig. 1 is a kind of structured flowchart for the electronic equipment that can be applied in the embodiment of the present invention；
Fig. 2 is the flow chart of ECG treatment method provided in an embodiment of the present invention；
Fig. 3 is the pending lead waveform sub-district area image signal in ECG treatment method provided in an embodiment of the present invention
Figure；
Fig. 4 is the lead waveform sub-district area image schematic diagram in ECG treatment method provided in an embodiment of the present invention；
Fig. 5 is the reduction lead signals schematic diagram generated in ECG treatment method provided in an embodiment of the present invention；
Fig. 6 is the papery electrocardiogram I lead schematic diagrames provided in an embodiment of the present invention read；
Fig. 7 is the structured flowchart for the electrocardiogram processing unit that present example provides.
Embodiment
To make the purpose, technical scheme and advantage of the embodiment of the present invention clearer, below in conjunction with the embodiment of the present invention
In accompanying drawing, the technical scheme in the embodiment of the present invention is clearly and completely described, it is clear that described embodiment is
Part of the embodiment of the present invention, rather than whole embodiments.The present invention implementation being generally described and illustrated herein in the accompanying drawings
The component of example can be configured to arrange and design with a variety of.Therefore, the reality of the invention to providing in the accompanying drawings below
The detailed description for applying example is not intended to limit the scope of claimed invention, but is merely representative of the selected implementation of the present invention
Example.Based on the embodiment in the present invention, what those of ordinary skill in the art were obtained under the premise of creative work is not made
Every other embodiment, belongs to the scope of protection of the invention.
It should be noted that：Similar label and letter represents similar terms in following accompanying drawing, therefore, once a certain Xiang Yi
It is defined, then it further need not be defined and explained in subsequent accompanying drawing in individual accompanying drawing.Meanwhile the present invention's
In description, term " first ", " second " etc. are only used for distinguishing description, and it is not intended that instruction or hint relative importance.
Fig. 1 shows a kind of structured flowchart for the electronic equipment 100 that can be applied in the embodiment of the present invention.As shown in figure 1,
Electronic equipment 100 can include memory 102, storage control 104, one or more (one is only shown in Fig. 1) processors
106th, at Peripheral Interface 108, input/output module 110, audio-frequency module 112, display module 114, radio-frequency module 116 and electrocardiogram
Manage device.
Memory 102, storage control 104, processor 106, Peripheral Interface 108, input/output module 110, audio mould
Directly or indirectly electrically connected between block 112, display module 114,116 each element of radio-frequency module, with realize the transmission of data or
Interaction.For example, electrical connection can be realized by one or more communication bus or signal bus between these elements.At electrocardiogram
Reason method includes at least one software that can be stored in the form of software or firmware (firmware) in memory 102 respectively
Functional module, such as the software function module or computer program that the electrocardiogram processing unit includes.
Memory 102 can store various software programs and module, the electrocardiogram processing provided such as the embodiment of the present application
Programmed instruction/module corresponding to method and device.Processor 106 by run storage software program in the memory 102 with
And module, so as to perform various function application and data processing, that is, realize the ECG treatment method in the embodiment of the present application.
Memory 102 can include but is not limited to random access memory (Random Access Memory, RAM), only
Read memory (Read Only Memory, ROM), programmable read only memory (Programmable Read-Only
Memory, PROM), erasable read-only memory (Erasable Programmable Read-Only Memory, EPROM),
Electricallyerasable ROM (EEROM) (Electric Erasable Programmable Read-Only Memory, EEPROM) etc..
Processor 106 can be a kind of IC chip, have signal handling capacity.Above-mentioned processor can be general
Processor, including central processing unit (Central Processing Unit, abbreviation CPU), network processing unit (Network
Processor, abbreviation NP) etc.；It can also be digital signal processor (DSP), application specific integrated circuit (ASIC), ready-made programmable
Gate array (FPGA) either other PLDs, discrete gate or transistor logic, discrete hardware components.It can
To realize or perform disclosed each method, step and the logic diagram in the embodiment of the present application.General processor can be micro-
Processor or the processor can also be any conventional processors etc..
Various input/output devices are coupled to processor 106 and memory 102 by the Peripheral Interface 108.At some
In embodiment, Peripheral Interface 108, processor 106 and storage control 104 can be realized in one single chip.Other one
In a little examples, they can be realized by independent chip respectively.
Input/output module 110 is used to be supplied to user input data to realize interacting for user and electronic equipment 100.It is described
Input/output module 110 may be, but not limited to, mouse and keyboard etc..
Audio-frequency module 112 provides a user COBBAIF, and it may include one or more microphones, one or more raises
Sound device and voicefrequency circuit.
Display module 114 provides an interactive interface (such as user interface) between electronic equipment 100 and user
Or referred to for display image data to user.In the present embodiment, the display module 114 can be liquid crystal display or touch
Control display.If touch control display, it can be that the capacitance type touch control screen or resistance-type for supporting single-point and multi-point touch operation touch
Control screen etc..Single-point and multi-point touch operation is supported to refer to that touch control display can sense on the touch control display one or more
Individual opening position is with caused touch control operation, and the touch control operation that this is sensed transfers to processor 106 to be calculated and handled.
Radio-frequency module 116 is used to receiving and sending electromagnetic wave, realizes the mutual conversion of electromagnetic wave and electric signal, so that with
Communication network or other equipment are communicated.
It is appreciated that structure shown in Fig. 1 is only to illustrate, electronic equipment 100 may also include it is more more than shown in Fig. 1 or
Less component, or there is the configuration different from shown in Fig. 1.Each component shown in Fig. 1 can use hardware, software or its
Combination is realized.
In the embodiment of the present invention, electronic equipment 100 can be used as user terminal, or as server.User terminal
Can be PC (personal computer) computer, tablet personal computer, mobile phone, notebook computer, intelligent television, set top box, vehicle-mounted
The terminal devices such as terminal.
Referring to Fig. 2, the embodiments of the invention provide a kind of ECG treatment method, methods described include step S200,
Step S210, step S220, step S230, step S240 and step S250, step S260.
Step S200：The papery electrocardiogram image read is subjected to separating treatment, obtains the papery electrocardiogram image
Lead curve regions image.
Based on step S200, further, the background grid in the papery electrocardiogram image is removed, obtains the papery
First curve text information image of electrocardiogram image；Processing is filtered to the first curve text information image, obtained
Second curve text information image；The second curve text information image is handled according to projected outline's method, obtains institute
State the lead curve regions image of papery electrocardiogram image.
The thermal printing paper of papery electrocardiogram image is all pink background grid, and the curve and word of printing are black
Color.The picture of printed report is cromogram, according to RGB combined color image principles, the song of pink background grid and black
Line has obvious difference in the threshold value of R, G, B layer picture.Further, the background grid in the papery electrocardiogram image is obtained
With the intensity value ranges of the prospect curve word in the papery electrocardiogram image；According to threshold method, the papery electrocardio is removed
After background grid in figure image, the first curve text information image of the papery electrocardiogram image is obtained.First curve text
Word frame is that P1 is a black and white picture for including curve and word, and curve and word are white, gray value 255, are carried on the back
Scape processing is black, gray value 0.
Further, due to the resolution ratio of image and the resolution differences of laser printer, the first curve text information figure
As there may be scattered noise signal, in order to obtain clean complete curve, the first curve text information image is entered
Row median filter process, obtain the second curve text information image.
Further, according to the particularity of the second curve text information image, in the second curve text information image
Certain interstitial spaces between line of text and curve regions all be present, and be all distributed rectangular, using projected outline's method, by
Two curve text information image projections to Y-axis, calculate projecting direction pixel it is cumulative and, pixel and be 0 after current line projection, under
A line projected pixel and more than 0, illustrates that the line position is set to blank borderline region, so as to find in the second curve text information image
The upper and lower border of lead curve whole region, a lead curve regions image for including lead curve is generated, i.e., described in acquisition
The lead curve regions image of papery electrocardiogram image is P2.
Step S210：By the lead curve regions image using behavior unit separation as multiple lead wave regions images,
Each lead wave regions image includes multiple lead signals.
Based on step S210, further, according to the up-and-down boundary position of the lead curve regions image, will entirely search
Rope region division is multiple deciles, and each decile is as the row field of search for searching for the often maximum amplitude position of capable multiple lead signals
Domain；
Each self-corresponding row region of search of the multiple decile is projected to Y-axis, by each row region of search found
Pixel and maximum position of adding up as the first reference point；
Based on first reference point, the lead end position of the Y-axis of search projection downwards, note in default hunting zone
The lower boundary position of often row lead signals is recorded, until obtaining multiple lead wave regions images.
In the present embodiment, the multiple is 4.The lead curve regions image is 3 × 4+1 printing model, i.e.,：
Often row has 4 lead signals, plus rhythm and pace of moving things lead, is divided into 4 row subregions altogether.It is P2 by the lead curve regions image
It is P2_1, P2_2, P2_3, P2_4 by 4 lead wave regions images of behavior unit separation.
Step S220：Each lead wave regions image is divided into multiple correspondences, and each the lead signals are led
Join waveform sub-district area image.
There are lead separator, each lead between multiple lead signals in each lead wave regions image are adjacent
There is lead symbol above the original position of signal, further, will each lead wave regions image based on step S220
Projected to X-axis；According to the characteristic of the lead symbol and the lead separator, the reference position of the lead symbol is obtained,
Again using the reference position as the second reference point, the point of predetermined threshold value is exceeded toward the default window ranges search peak in left side
Position, the position of the lead separator is recorded as, to be divided into the pending lead of each lead signals of multiple correspondences
Waveform sub-district area image；The pending lead waveform sub-district area image is projected to Y-axis, according to lead symbol and lead waveform
Projected pixel superimposed curves catastrophe characteristicses, obtain lead character pitch beginning and end position；Pending led described
Join waveform sub-district area image to project to X-axis, isolate lead waveform and lead separator, obtain the beginning and end of lead waveform
Position；Beginning and end position based on the lead character pitch, the beginning and end position of the lead waveform, are led
Join waveform sub-district area image, the lead waveform sub-district area image includes the lead waveform.
Specifically, the characteristic of the lead symbol is：After lead symbol projects to X-axis, accumulation and curve it is intensive and
Amplitude is low.The feature of lead separator is：After lead separator projects to X-axis, pixel is cumulative and amplitude is larger, is obvious
Crest.The one lead wave regions image such as P2_1 is divided into the lead ripple of each lead signals of 4 correspondences
Shape sub-district area image is P2_1_1, P2_1_2, P2_1_3, P2_1_4.As shown in figure 3, pending lead waveform sub-district area image
In, contain lead symbol I and part lead separator.As shown in figure 4, lead waveform sub-district area image only includes lead waveform.
Step S230：Obtain the pixel position of the lead waveform in the lead waveform sub-district area image.
Due to printhead resolution ratio, actual print of each sampled point signal is the multiple pixels of a row, is needed
Position by multiple pixel values processing of a row for 1 pixel.In the mutated site of lead waveform raising and lowering, waveform
Pixel-intensive and quantity is more, the amplitude to avoid processing from bringing reduce error, it is necessary to be in the region of peak value in lead waveform, taken
Most upper peripheral pixels, the region of trough is in waveform, takes and most descends peripheral pixels.
Further, the lead waveform in the lead waveform sub-district area image is projected to Y direction, according to lead ripple
The characteristic of shape, finds baseline position, and using the baseline position as the 3rd reference point, the baseline position above processes pixel is
Most upper peripheral pixels position is taken, the following processes pixel of baseline position most descends peripheral pixels position to take, the baseline position
Average location of pixels is nearby taken, to obtain all pixels of lead waveform point position.
Step S240：Based on the lead curve regions image and the pixel position, reduction lead signals are obtained
Range value.
Further, in the lead curve regions image, the position of calibration symbol is obtained；According to the calibration symbol
Position, obtain the number of pixels that is included of height of the calibration symbol；What the height based on the calibration symbol was included
Number of pixels, obtain the pixel resolution of single pixel；Using the pixel position as the 4th reference point, calculate the reduction and lead
Join the relative magnitude of signal, then based on the pixel resolution and the relative magnitude, obtain the range value for reducing lead signals.
For example, the number of pixels that the height of the calibration symbol is included is 120 pixels, 1mv is represented, calibrates symbol
Sensitivity is 10mm/mv, and the pixel resolution of single pixel is：1/120≈8.3×10-3mv；Using the pixel position as
4th reference point, calculate the relative magnitude of the reduction lead signals；The relative magnitude is multiplied by the pixel resolution again
Obtained product is the range value of reduction lead signals, that is, obtains the range value for reducing lead signals.
Step S250：Based on the sum of all pixels of the lead waveform in the lead waveform sub-district area image got, obtain
Obtain the phase between reducing lead signals.
Further, to the actual time-write interval of the papery electrocardiogram image got and the picture of the lead waveform
Plain sum carries out interpolation processing, obtains the phase between the reduction lead signals.
Step S260：Based on the pixel position, the range value and it is described between the phase, generate it is described reduction lead letter
Number.
For example, each lead waveform is 2.5s, sample rate 500, by the sum of all pixels of the lead waveform, enter row interpolation
Handle as 1250 points, you can using full backup as actual samples signal.In the reduction I lead signals and Fig. 6 that are generated such as Fig. 5
The I leads of the original papery electrocardiogram image read, contrast understand that ECG treatment method provided in an embodiment of the present invention can
Completely to restore papery electrocardiographic wave well.
In order to which ECG treatment method beneficial effect provided in an embodiment of the present invention is further illustrated, 1, table 2 is shown in Table
Shown in table 3, case numbering 1-12 represents 12 papery electrocardiograms, is represented with reference to P-wave amplitude on each papery electrocardiogram image
P wave amplitude angle value, reduction P-wave amplitude represents the corresponding P for obtaining reduction lead signals of the ECG treatment method for implementing the present invention
Wave amplitude angle value, P-wave amplitude difference represent with reference to P-wave amplitude and reduced the difference of P-wave amplitude, and each paper is represented with reference to Q wave amplitudes
Q wave amplitude angle value on barycenter electrograph image, reduction Q wave amplitudes represent that the ECG treatment method for implementing the present invention is corresponding and gone back
The Q wave amplitude angle value of former lead signals, Q wave amplitudes difference represents with reference to Q wave amplitudes and reduced the difference of Q wave amplitudes, with reference to R wave amplitudes
Value represents the R wave amplitude angle value on each papery electrocardiogram image, and reduction R-wave amplitude represents to implement the electrocardiogram processing side of the present invention
The corresponding R wave amplitude angle value for obtaining reduction lead signals of method, R-wave amplitude difference represent with reference to R-wave amplitude and reduced the difference of R-wave amplitude
Value, the S wave amplitude angle value on each papery electrocardiogram image is represented with reference to S wave amplitudes, reduction S wave amplitudes represent to implement the present invention's
The corresponding S wave amplitude angle value for obtaining reduction lead signals of ECG treatment method, the expression of S wave amplitudes difference is with reference to S wave amplitudes and also
The difference of former S wave amplitudes, the T wave amplitude angle value on each papery electrocardiogram image is represented with reference to T wave amplitudes, reduces T wave amplitude tables
Show the corresponding T wave amplitude angle value for obtaining reduction lead signals of the ECG treatment method for implementing the present invention, T wave amplitudes difference represents ginseng
Examine T wave amplitudes and reduce the difference of T wave amplitudes.
By the amplitude of P, Q, R, S, T ripple of 1~table of table 3, each amplitude error is below 0.1mv, P, Q, S, T wave amplitude value difference
Away from especially small, less than 0.05mv, except No. 3 cases because P wave amplitudes are too small, therefore diagnosing machinery is measured when calculating P wave amplitudes
Error be present, by the way that phase contrast understands heart rate difference very little, other phase differences between phase, QT between phase, QRS between the heart rate of table 4, PR
Within 30ms, the accuracy of re-computation is very high after restoring data.When weight analysis and processing being carried out so as to doctor or instrument,
It just can guarantee that the accuracy of diagnosis.
The case range value of table 1 contrast 1
The case range value of table 2 contrast 2
The case range value of table 3 contrast 3
| Case is numbered | With reference to T wave amplitudes (mv) | Reduce T wave amplitudes (mv) | T wave amplitudes difference (mv) |
| 1 | 0.2294 | 0.2052 | 0.0242 |
| 2 | 0.4904 | 0.436 | 0.0544 |
| 3 | 0.3788 | 0.3234 | 0.0554 |
| 4 | 0.2168 | 0.1868 | 0.03 |
| 6 | 0.383 | 0.3254 | 0.0576 |
| 7 | 0.2748 | 0.235 | 0.0398 |
| 8 | 0.1096 | 0.094 | 0.0156 |
| 10 | 0.5478 | 0.4674 | 0.0804 |
| 11 | 0.1318 | 0.1046 | 0.0272 |
| 12 | 0.255 | 0.217 | 0.038 |
The phase contrasts between the case of table 4
The embodiments of the invention provide a kind of ECG treatment method and device, and the papery electrocardiogram image read is entered
Row separating treatment, obtain the lead curve regions image of the papery electrocardiogram image；By the lead curve regions image with
Behavior unit separation is multiple lead wave regions images, and each lead wave regions image includes multiple lead signals；
Each lead wave regions image is divided into the lead waveform sub-district area image of each lead signals of multiple correspondences；
Obtain the pixel position of the lead waveform in the lead waveform sub-district area image；Based on the lead curve regions image and
The pixel position, obtain the range value of reduction lead signals；Based in the lead waveform sub-district area image got
Lead waveform sum of all pixels, obtain it is described reduction lead signals between the phase；Based on the pixel position, the range value
And it is described between the phase, generate the reduction lead signals.Ensure the high reduction degree of reduction lead signals, it is simple to operate, more accurate.
Referring to Fig. 7, present example provides a kind of electrocardiogram processing unit 300, described device 300 can include：
First acquisition unit 310, separating element 320, cutting unit 330, second acquisition unit 340, the 3rd acquiring unit the 350, the 4th
Acquiring unit 360 and generation unit 370.
First acquisition unit 310, the papery electrocardiogram image for that will read carry out separating treatment, obtain the papery
The lead curve regions image of electrocardiogram image.
First acquisition unit 310 can include first and obtain subelement 311.
First obtains subelement 311, for removing the background grid in the papery electrocardiogram image, obtains the papery
First curve text information image of electrocardiogram image；Processing is filtered to the first curve text information image, obtained
Second curve text information image；The second curve text information image is handled according to projected outline's method, obtains institute
State the lead curve regions image of papery electrocardiogram image.
First obtains subelement 311, is additionally operable to obtain the background grid in the papery electrocardiogram image and the papery
The intensity value ranges of prospect curve word in electrocardiogram image；According to threshold method, remove in the papery electrocardiogram image
After background grid, the first curve text information image of the papery electrocardiogram image is obtained.
First obtains subelement 311, is additionally operable to carry out median filter process to the first curve text information image, obtains
Obtain the second curve text information image.
Separating element 320, for by the lead curve regions image using behavior unit separation as multiple lead waveform areas
Area image, each lead wave regions image include multiple lead signals.
Separating element 320 can include separating subelement 321.
Separate subelement 321, for the up-and-down boundary position according to the lead curve regions image, by the whole field of search
Domain is divided into multiple deciles, and each decile is as the row region of search for searching for the often maximum amplitude position of capable multiple lead signals；
Each self-corresponding row region of search of the multiple decile is projected to Y-axis, by the pixel in each row region of search found and
Maximum position add up as the first reference point；Based on first reference point, search is thrown downwards in default hunting zone
The lead end position of shadow Y-axis, the lower boundary position of often row lead signals is recorded, until obtaining multiple lead wave regions figures
Picture.
Cutting unit 330, for each lead wave regions image to be divided into multiple correspondences each lead
The lead waveform sub-district area image of signal.
There are lead separator, each lead between multiple lead signals in each lead wave regions image are adjacent
There is lead symbol above the original position of signal.Cutting unit 330 can include segmentation subelement 331.
Split subelement 331, for each lead wave regions image to be projected to X-axis；Accorded with according to the lead
Number and the lead separator characteristic, obtain the reference position of the lead symbol, then second is used as using the reference position
Reference point, the position of the point of predetermined threshold value is exceeded toward the default window ranges search peak in left side, be recorded as the lead and separate
The position of symbol, to be divided into the pending lead waveform sub-district area image of each lead signals of multiple correspondences；Treated described
Handle lead waveform sub-district area image to project to Y-axis, according to the prominent of the projected pixel superimposed curves of lead symbol and lead waveform
Become characteristic, obtain the beginning and end position of lead character pitch；The pending lead waveform sub-district area image is thrown to X-axis
Shadow, lead waveform and lead separator are isolated, obtain the beginning and end position of lead waveform；It is high based on the lead symbol
The beginning and end position of degree, the beginning and end position of the lead waveform, lead waveform sub-district area image is obtained, it is described to lead
Connection waveform sub-district area image includes the lead waveform.
Second acquisition unit 340, for obtaining the pixel position of the lead waveform in the lead waveform sub-district area image
Put.
Second acquisition unit 340 can include second and obtain subelement 341.
Second obtains subelement 341, for the lead waveform in the lead waveform sub-district area image to be projected to Y-axis side
To, according to the characteristic of lead waveform, baseline position is found, using the baseline position as the 3rd reference point, the baseline position
For above processes pixel to take most upper peripheral pixels position, the following processes pixel of baseline position most descends peripheral pixels position to take
Put, the baseline position nearby takes average location of pixels, to obtain all pixels of lead waveform point position.
3rd acquiring unit 350, for based on the lead curve regions image and the pixel position, being reduced
The range value of lead signals.
3rd acquiring unit 350 can include the 3rd and obtain subelement 351.
3rd obtains subelement 351, in the lead curve regions image, obtaining the position of calibration symbol；Root
According to the position of the calibration symbol, the number of pixels that the height of the calibration symbol is included is obtained；Based on the calibration symbol
The number of pixels that is included of height, obtain the pixel resolution of single pixel；Using the pixel position as the 4th reference point,
The relative magnitude of the reduction lead signals is calculated, then based on the pixel resolution and the relative magnitude, obtains reduction and leads
Join the range value of signal.
4th acquiring unit 360, for based on the lead waveform in the lead waveform sub-district area image got
Sum of all pixels, obtain the phase between the reduction lead signals.
4th acquiring unit 360 can include the 4th and obtain subelement 361.
4th obtains subelement 361, for the actual time-write interval by the papery electrocardiogram image to getting
Interpolation processing is carried out with the sum of all pixels of the lead waveform, obtains the phase between the reduction lead signals.
Generation unit 370, for based on the pixel position, the range value and it is described between the phase, generate the reduction
Lead signals.
Above each unit can be that now, above-mentioned each unit can be stored in memory 102 by software code realization.
Above each unit can equally be realized by hardware such as IC chip.
Electrocardiogram processing unit 300 provided in an embodiment of the present invention, its realization principle and caused technique effect and foregoing
Embodiment of the method is identical, and to briefly describe, device embodiment part does not refer to part, refers to corresponding in preceding method embodiment
Content.
In several embodiments provided herein, it should be understood that disclosed apparatus and method, can also pass through
Other modes are realized.Device embodiment described above is only schematical, for example, flow chart and block diagram in accompanying drawing
Show the device of multiple embodiments according to the present invention, method and computer program product architectural framework in the cards,
Function and operation.At this point, each square frame in flow chart or block diagram can represent the one of a module, program segment or code
Part, a part for the module, program segment or code include one or more and are used to realize holding for defined logic function
Row instruction.It should also be noted that at some as in the implementation replaced, the function that is marked in square frame can also with different from
The order marked in accompanying drawing occurs.For example, two continuous square frames can essentially perform substantially in parallel, they are sometimes
It can perform in the opposite order, this is depending on involved function.It is it is also noted that every in block diagram and/or flow chart
The combination of individual square frame and block diagram and/or the square frame in flow chart, function or the special base of action as defined in performing can be used
Realize, or can be realized with the combination of specialized hardware and computer instruction in the system of hardware.
In addition, each functional module in each embodiment of the present invention can integrate to form an independent portion
Point or modules individualism, can also two or more modules be integrated to form an independent part.
If the function is realized in the form of software function module and is used as independent production marketing or in use, can be with
It is stored in a computer read/write memory medium.Based on such understanding, technical scheme is substantially in other words
The part to be contributed to prior art or the part of the technical scheme can be embodied in the form of software product, the meter
Calculation machine software product is stored in a storage medium, including some instructions are causing a computer equipment (can be
People's computer, server, or network equipment etc.) perform all or part of step of each embodiment methods described of the present invention.
And foregoing storage medium includes：USB flash disk, mobile hard disk, read-only storage (ROM, Read-Only Memory), arbitrary access are deposited
Reservoir (RAM, Random Access Memory), magnetic disc or CD etc. are various can be with the medium of store program codes.Need
Illustrate, herein, such as first and second or the like relational terms be used merely to by an entity or operation with
Another entity or operation make a distinction, and not necessarily require or imply between these entities or operation any this reality be present
The relation or order on border.Moreover, term " comprising ", "comprising" or its any other variant are intended to the bag of nonexcludability
Contain, so that process, method, article or equipment including a series of elements not only include those key elements, but also including
The other element being not expressly set out, or also include for this process, method, article or the intrinsic key element of equipment.
In the absence of more restrictions, the key element limited by sentence "including a ...", it is not excluded that including the key element
Process, method, other identical element also be present in article or equipment.
The preferred embodiments of the present invention are the foregoing is only, are not intended to limit the invention, for the skill of this area
For art personnel, the present invention can have various modifications and variations.Within the spirit and principles of the invention, that is made any repaiies
Change, equivalent substitution, improvement etc., should be included in the scope of the protection.It should be noted that：Similar label and letter exists
Similar terms is represented in following accompanying drawing, therefore, once being defined in a certain Xiang Yi accompanying drawing, is then not required in subsequent accompanying drawing
It is further defined and explained.
The foregoing is only a specific embodiment of the invention, but protection scope of the present invention is not limited thereto, any
Those familiar with the art the invention discloses technical scope in, change or replacement can be readily occurred in, should all be contained
Cover within protection scope of the present invention.Therefore, protection scope of the present invention described should be defined by scope of the claims.
It should be noted that herein, such as first and second or the like relational terms are used merely to a reality
Body or operation make a distinction with another entity or operation, and not necessarily require or imply and deposited between these entities or operation
In any this actual relation or order.Moreover, term " comprising ", "comprising" or its any other variant are intended to
Nonexcludability includes, so that process, method, article or equipment including a series of elements not only will including those
Element, but also the other element including being not expressly set out, or it is this process, method, article or equipment also to include
Intrinsic key element.In the absence of more restrictions, the key element limited by sentence "including a ...", it is not excluded that
Other identical element also be present in process, method, article or equipment including the key element.
Claims (10)
1. a kind of ECG treatment method, it is characterised in that methods described includes：
The papery electrocardiogram image read is subjected to separating treatment, obtains the lead curve regions of the papery electrocardiogram image
Image；
By the lead curve regions image using behavior unit separation as multiple lead wave regions images, each lead ripple
Shape area image includes multiple lead signals；
Each lead wave regions image is divided into the lead waveform subregion of each lead signals of multiple correspondences
Image；
Obtain the pixel position of the lead waveform in the lead waveform sub-district area image；
Based on the lead curve regions image and the pixel position, the range value for reducing lead signals is obtained；
Based on the sum of all pixels of the lead waveform in the lead waveform sub-district area image got, the reduction lead is obtained
Phase between signal；
Based on the pixel position, the range value and it is described between the phase, generate the reduction lead signals.
2. according to the method for claim 1, it is characterised in that carry out the papery electrocardiogram image read at separation
Reason, the lead curve regions image of the papery electrocardiogram image is obtained, including：
The background grid in the papery electrocardiogram image is removed, obtains the first curve word letter of the papery electrocardiogram image
Cease image；
Processing is filtered to the first curve text information image, obtains the second curve text information image；
The second curve text information image is handled according to projected outline's method, obtains the papery electrocardiogram image
Lead curve regions image.
3. according to the method for claim 2, it is characterised in that the background grid in the papery electrocardiogram image is removed,
The first curve text information image of the papery electrocardiogram image is obtained, including：
Obtain the prospect curve word in the background grid and the papery electrocardiogram image in the papery electrocardiogram image
Intensity value ranges；
According to threshold method, after removing the background grid in the papery electrocardiogram image, the papery electrocardiogram image is obtained
First curve text information image.
4. according to the method for claim 2, it is characterised in that place is filtered to the first curve text information image
Reason, the second curve text information image is obtained, including：
Median filter process is carried out to the first curve text information image, obtains the second curve text information image.
5. according to the method for claim 1, it is characterised in that by the lead curve regions image with behavior unit separation
For multiple lead wave regions images, including：
According to the up-and-down boundary position of the lead curve regions image, whole region of search is divided into multiple deciles, each
Decile is as the row region of search for searching for the often maximum amplitude position of capable multiple lead signals；
Each self-corresponding row region of search of the multiple decile is projected to Y-axis, by the picture in each row region of search found
Element and maximum position of adding up are as the first reference point；
Based on first reference point, the lead end position of the Y-axis of search projection downwards, record are every in default hunting zone
The lower boundary position of row lead signals, until obtaining multiple lead wave regions images.
6. according to the method for claim 1, it is characterised in that multiple leads in each lead wave regions image
There is lead separator between signal is adjacent, there is lead symbol above the original position of each lead signals, will each lead
Wave regions image is divided into the lead waveform sub-district area image of each lead signals of multiple correspondences, including：
Each lead wave regions image is projected to X-axis；
According to the characteristic of the lead symbol and the lead separator, the reference position of the lead symbol is obtained, then with institute
Reference position is stated as the second reference point, the position toward the default window ranges search peak in left side more than the point of predetermined threshold value,
Be recorded as the position of the lead separator, be divided into multiple correspondences each the lead signals pending lead waveform
Area image；
The pending lead waveform sub-district area image is projected to Y-axis, according to the projected pixel of lead symbol and lead waveform
The catastrophe characteristicses of superimposed curves, obtain the beginning and end position of lead character pitch；
The pending lead waveform sub-district area image is projected to X-axis, lead waveform and lead separator is isolated, is led
Join the beginning and end position of waveform；
Beginning and end position based on the lead character pitch, the beginning and end position of the lead waveform, are led
Join waveform sub-district area image, the lead waveform sub-district area image includes the lead waveform.
7. according to the method for claim 1, it is characterised in that obtain the lead ripple in the lead waveform sub-district area image
The pixel position of shape, including：
Lead waveform in the lead waveform sub-district area image is projected to Y direction, according to the characteristic of lead waveform, searched
To baseline position, using the baseline position as the 3rd reference point, the baseline position above processes pixel is to take most upper peripheral picture
Plain position, the following processes pixel of baseline position most descend peripheral pixels position to take, and the baseline position nearby takes average picture
Plain position, to obtain all pixels of lead waveform point position.
8. according to the method for claim 7, it is characterised in that based on the lead curve regions image and the pixel
Position, the range value of reduction lead signals is obtained, including：
In the lead curve regions image, the position of calibration symbol is obtained；
According to the position of the calibration symbol, the number of pixels that the height of the calibration symbol is included is obtained；
The number of pixels that height based on the calibration symbol is included, obtain the pixel resolution of single pixel；
Using the pixel position as the 4th reference point, the relative magnitude of the reduction lead signals is calculated, then based on the picture
Plain resolution ratio and the relative magnitude, obtain the range value of reduction lead signals.
9. according to the method for claim 1, it is characterised in that based in the lead waveform sub-district area image got
Lead waveform sum of all pixels, obtain it is described reduction lead signals between the phase, including：
Entered by the actual time-write interval of the papery electrocardiogram image to getting and the sum of all pixels of the lead waveform
Row interpolation processing, obtain the phase between the reduction lead signals.
10. a kind of electrocardiogram processing unit, it is characterised in that described device includes：
First acquisition unit, the papery electrocardiogram image for that will read carry out separating treatment, obtain the papery electrocardiogram
The lead curve regions image of image；
Separating element, for by the lead curve regions image using behavior unit separation as multiple lead wave regions images,
Each lead wave regions image includes multiple lead signals；
Cutting unit, for each lead wave regions image to be divided into multiple correspondences, each the lead signals are led
Join waveform sub-district area image；
Second acquisition unit, for obtaining the pixel position of the lead waveform in the lead waveform sub-district area image；
3rd acquiring unit, for based on the lead curve regions image and the pixel position, obtaining reduction lead letter
Number range value；
4th acquiring unit, it is total for the pixel based on the lead waveform in the lead waveform sub-district area image got
Number, obtain the phase between the reduction lead signals；
Generation unit, for based on the pixel position, the range value and it is described between the phase, generate it is described reduction lead letter
Number.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710870201.0A CN107625519B (en) | 2017-09-20 | 2017-09-20 | Electrocardiogram processing method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710870201.0A CN107625519B (en) | 2017-09-20 | 2017-09-20 | Electrocardiogram processing method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107625519A true CN107625519A (en) | 2018-01-26 |
| CN107625519B CN107625519B (en) | 2020-09-08 |
Family
ID=61103824
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710870201.0A Active CN107625519B (en) | 2017-09-20 | 2017-09-20 | Electrocardiogram processing method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107625519B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109480824A (en) * | 2018-12-11 | 2019-03-19 | 武汉中旗生物医疗电子有限公司 | Processing method, device and the server of EGC waveform data |
| CN109480825A (en) * | 2018-12-13 | 2019-03-19 | 武汉中旗生物医疗电子有限公司 | The processing method and processing device of electrocardiogram (ECG) data |
| CN113077472A (en) * | 2021-04-07 | 2021-07-06 | 华南理工大学 | Paper electrocardiogram curve image segmentation method, system, device and medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102038498A (en) * | 2011-01-20 | 2011-05-04 | 天津理工大学 | Method for digitalizing paper electrocardiogram |
| US8473038B2 (en) * | 2008-05-30 | 2013-06-25 | Yuan Ze University | Mobile and web-based 12-lead ECG management |
| WO2014021515A1 (en) * | 2012-07-30 | 2014-02-06 | 가톨릭대학교 산학협력단 | Method and apparatus for evaluating health of fetus |
| CN106202948A (en) * | 2016-07-18 | 2016-12-07 | 上海市第人民医院 | A kind of can the method and system of digitized ultrasoundcardiogram report automatically management |
| CN106326818A (en) * | 2015-06-30 | 2017-01-11 | 东南大学 | Method and device for digitizing paper hydrological data |
-
2017
- 2017-09-20 CN CN201710870201.0A patent/CN107625519B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8473038B2 (en) * | 2008-05-30 | 2013-06-25 | Yuan Ze University | Mobile and web-based 12-lead ECG management |
| CN102038498A (en) * | 2011-01-20 | 2011-05-04 | 天津理工大学 | Method for digitalizing paper electrocardiogram |
| WO2014021515A1 (en) * | 2012-07-30 | 2014-02-06 | 가톨릭대학교 산학협력단 | Method and apparatus for evaluating health of fetus |
| CN106326818A (en) * | 2015-06-30 | 2017-01-11 | 东南大学 | Method and device for digitizing paper hydrological data |
| CN106202948A (en) * | 2016-07-18 | 2016-12-07 | 上海市第人民医院 | A kind of can the method and system of digitized ultrasoundcardiogram report automatically management |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109480824A (en) * | 2018-12-11 | 2019-03-19 | 武汉中旗生物医疗电子有限公司 | Processing method, device and the server of EGC waveform data |
| CN109480825A (en) * | 2018-12-13 | 2019-03-19 | 武汉中旗生物医疗电子有限公司 | The processing method and processing device of electrocardiogram (ECG) data |
| CN113077472A (en) * | 2021-04-07 | 2021-07-06 | 华南理工大学 | Paper electrocardiogram curve image segmentation method, system, device and medium |
| CN113077472B (en) * | 2021-04-07 | 2023-02-10 | 华南理工大学 | Paper electrocardiogram curve image segmentation method, system, device and medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107625519B (en) | 2020-09-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Pilia et al. | ECGdeli-an open source ECG delineation toolbox for MATLAB | |
| US11147500B2 (en) | Electrocardiogram processing system for delineation and classification | |
| Wang et al. | Arrhythmia classification algorithm based on multi-head self-attention mechanism | |
| Qin et al. | An Adaptive and Time‐Efficient ECG R‐Peak Detection Algorithm | |
| US11617528B2 (en) | Systems and methods for reduced lead electrocardiogram diagnosis using deep neural networks and rule-based systems | |
| EP3698709A1 (en) | Electrocardiogram information processing method and electrocardiogram workstation system | |
| Rahman et al. | Robustness of electrocardiogram signal quality indices | |
| EP3698708A1 (en) | Ecg information processing method and ecg workstation | |
| Ravichandran et al. | Novel tool for complete digitization of paper electrocardiography data | |
| Maheshwari et al. | An automated algorithm for online detection of fragmented QRS and identification of its various morphologies | |
| CN102274020A (en) | Low-power consumption portable electrocardiograph monitor and control method thereof | |
| CN107625519A (en) | ECG treatment method and device | |
| Wang et al. | Automatic digital ECG signal extraction and normal QRS recognition from real scene ECG images | |
| Shivashankara et al. | ECG-Image-Kit: a synthetic image generation toolbox to facilitate deep learning-based electrocardiogram digitization | |
| Monedero | A novel ECG diagnostic system for the detection of 13 different diseases | |
| Saini et al. | K-nearest neighbour-based algorithm for P-and T-waves detection and delineation | |
| US20150265175A1 (en) | Non-Invasive Evaluation of Cardiac Repolarisation Instability for Risk Stratification of Sudden Cardiac Death | |
| EP4041073A1 (en) | Systems and methods for electrocardiogram diagnosis using deep neural networks and rule-based systems | |
| CN118021316B (en) | Dynamic electrocardiographic data processing method and device, electronic equipment and storage medium | |
| Ganesh et al. | Combining optical character recognition with paper ECG digitization | |
| CN109949897A (en) | Display methods, device and the terminal device of monitoring curve | |
| CN107622245A (en) | Papery Wave shape extracting method and device | |
| Gupta et al. | Pre-processing based ECG signal analysis using emerging tools | |
| Mitra et al. | Pattern defined heuristic rules and directional histogram based online ECG parameter extraction | |
| Bae et al. | Electrocardiogram fiducial point detector using a bilateral filter and symmetrical point-filter structure |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |