Tachycardia diagnosis and treatment equipment

Abstract

The invention discloses tachycardia diagnosis and treatment equipment which comprises an ultrasonic transduction module, a control module, a strain rate module, a three-dimensional imaging module, a spiral wave end point positioning module and a diagnosis module. The ultrasonic transduction module emits ultrasonic waves and receives echo signals, the control module generates three-dimensional displacement field data through the echo signals, and after the strain rate is calculated through the strain rate module, the three-dimensional imaging module stores and maps the three-dimensional displacement field data into three-dimensional data. The spiral wave endpoint positioning module obtains an effective endpoint set through subimage detection and clustering screening, the diagnosis module judges tachycardia based on a preset threshold value, the phased array module is triggered to generate a treatment instruction containing sound intensity, phase and position information, and finally accurate focusing treatment is achieved through the ultrasonic transduction module. Through the non-invasive ultrasonic technology, multi-module cooperative processing and real-time data feedback, non-invasive diagnosis of tachycardia is achieved, and accuracy and safety are both achieved.

Description

Technical Field

[0001] This application relates to the field of ultrasonic diagnosis and treatment technology, and particularly to a tachycardia diagnosis and treatment device. Background Art

[0002] Spiral wave electrical signals are the main generation mechanism of tachycardia. Clinically, if spiral wave electrical signals in the target area can be detected and eliminated in time, tachycardia can be effectively treated.

[0003] Currently, implantable cardioverter defibrillators (ICDs) are generally used to treat tachycardia. The ICD effectively eliminates spiral waves and turbulence by stimulating with a large voltage to depolarize all myocardial cells simultaneously, achieving the effect of treating tachycardia.

[0004] However, the ICD needs to be implanted into the patient's body, which is likely to cause rejection in the patient's body, has a high cost, and has relatively large potential safety hazards to the body due to electrical shock stimulation.

[0005] Therefore, how to ensure the accuracy and safety of treatment is an urgent problem to be solved. Summary of the Invention

[0006] In view of this, a tachycardia diagnosis and treatment device provided by an embodiment of this application can reduce the treatment cost and ensure the safety of treatment. The tachycardia diagnosis and treatment device provided by an embodiment of this application is implemented as follows: A tachycardia diagnosis and treatment device provided by an embodiment of this application includes an ultrasonic transducer module, a control module, a strain rate module, a three-dimensional imaging module, a spiral wave endpoint positioning module, and a diagnosis module that are connected in sequence. The diagnosis module is also connected to the control module, where: The ultrasonic transducer module is configured to emit ultrasonic waves after receiving the transmission signal sent by the control module, and receive the echo signal after the ultrasonic waves are reflected by the target area, and send the echo signal to the control module; The control module is configured to send a transmission signal to the ultrasonic transducer module and receive the echo signal, generate three-dimensional displacement field data corresponding to the target area according to the echo signal, and send the three-dimensional displacement field data to the strain rate module; The strain rate module is configured to perform deformation tensor calculation on the received three-dimensional displacement field data to obtain a strain rate, and send the strain rate to the three-dimensional imaging module; The three-dimensional imaging module is configured to obtain a strain rate data set according to the received strain rate, perform three-dimensional mapping on the strain rate data set to obtain three-dimensional imaging data, and send the three-dimensional imaging data to the spiral wave endpoint positioning module; The spiral wave endpoint positioning module is configured to generate a plurality of sub-images based on the received three-dimensional imaging data, perform spiral wave endpoint detection on the plurality of sub-images to obtain the spiral wave endpoints of the plurality of sub-images, perform clustering and screening on the spiral wave endpoints of the plurality of sub-images to obtain a filtered set of spiral wave endpoints, and send the filtered set of spiral wave endpoints to the diagnosis module; The diagnosis module is configured to determine whether tachycardia exists in the target area according to a preset frequency threshold and the received filtered set of spiral wave endpoints; when tachycardia exists in the target area, the diagnosis module sends a treatment instruction and the filtered set of spiral wave endpoints to the control module; The control module is further configured to generate control information according to the received filtered set of spiral wave endpoints after receiving the treatment instruction from the control module, and send the control information to the ultrasonic transducer module, where the control information includes the sound intensity, phase of the ultrasonic wave, and the position of the spiral wave endpoint; The ultrasonic transducer module is further configured to adjust the sound intensity and phase of the ultrasonic wave according to the received control information and then emit the ultrasonic wave to the target area to perform ultrasonic focusing on the position of the spiral wave endpoint.

[0007] In some embodiments, the strain rate module includes a first receiving unit, a first calculating unit, a conversion unit, and a first sending unit that are connected in sequence. The strain rate module and the three-dimensional imaging module, where; The first receiving unit is configured to receive the three-dimensional displacement field data of the control module; The first calculating unit is configured to perform deformation tensor calculation on the three-dimensional displacement field data to obtain a gradient deformation tensor and a Green-Lagrange strain tensor; perform eigenvalue decomposition and invariant combination operations on the gradient deformation tensor and the Green-Lagrange strain tensor respectively to obtain the tensor invariant of the target area; The conversion unit is configured to perform strain rate conversion on the tensor invariant to obtain the strain rate of the target area; The first sending unit is configured to send the strain rate of the target area to the three-dimensional imaging module.

[0008] In some embodiments, the spiral wave endpoint positioning module includes a second receiving unit, a establishing unit, a cutting unit, a second calculating unit, a screening unit, and a second sending unit and the diagnosis module, where; The second receiving unit is configured to receive the three-dimensional imaging data of the three-dimensional imaging module; The establishing unit is configured to establish a three-dimensional rectangular coordinate system according to the three-dimensional imaging data; The cutting unit is used to cut the three-dimensional imaging data perpendicular to the coordinate axes along the x, y, and z axes of the three-dimensional rectangular coordinate system to generate a plurality of sub-images; The second calculation unit is used to adopt the Moore field tracking algorithm for each of the plurality of sub-images to obtain the intersection line between the foreground and the background in each sub-image, where the foreground is the area containing spiral waves in the sub-image, and the background is the area without spiral waves in the sub-image; when the intersection line meets the preset curvature threshold, the spiral wave endpoints of the sub-image corresponding to the intersection line are obtained by fitting the curvature circle algorithm; The screening unit is used to perform clustering screening on the spiral wave endpoints of the plurality of sub-images to obtain a set of screened spiral wave endpoints; The second sending unit is used to send the set of screened spiral wave endpoints to the diagnosis module.

[0009] In some embodiments, the screening unit includes an acquisition element and a screening element connected in sequence, where; The acquisition element is used to acquire the corresponding plurality of spiral wave endpoints in the plurality of sub-images to obtain an initial set of spiral wave endpoints; The screening element is used to perform a screening process on the spiral wave endpoints in the initial set of spiral wave endpoints whose standard deviation of spatial distribution exceeds a preset distance threshold to obtain a set of screened spiral wave endpoints.

[0010] In some embodiments, the preset distance threshold is 0.5 mm.

[0011] In some embodiments, the ultrasonic transducer module includes at least two groups of ultrasonic transducers.

[0012] In some embodiments, the operating frequency range of the ultrasonic transducers is 1 MHz to 10 MHz.

[0013] In some embodiments, the time of ultrasonic focusing is 20 ms.

[0014] In some embodiments, the preset frequency threshold is the occurrence frequency of 85% of the spiral wave endpoints.

[0015] In some embodiments, the imaging interval time of the three-dimensional mapping is 10 ms.

[0016] A tachycardia diagnosis and treatment device provided by an embodiment of the present application includes an ultrasonic transducer module, a control module, a strain rate module, a three-dimensional imaging module, a spiral wave endpoint positioning module, and a diagnosis module. The ultrasonic transducer module emits ultrasonic waves and receives echo signals. The control module generates three-dimensional displacement field data through the echo signals. After calculating the strain rate by the strain rate module, it is stored by the three-dimensional imaging module and mapped into three-dimensional data. The spiral wave endpoint positioning module obtains an effective endpoint set through sub-image detection and clustering screening. The diagnosis module determines tachycardia based on a preset threshold, triggers the phased array module to generate a treatment instruction containing sound intensity, phase, and position information, and finally the ultrasonic transducer module realizes precise focusing treatment. Through non-invasive ultrasonic technology, multi-module collaborative processing, and real-time data feedback, non-invasive diagnosis of tachycardia is achieved, with both accuracy and safety. The technical problems proposed in the background art are solved. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following will briefly introduce the drawings required for describing the embodiments of the present application or the prior art. Obviously, the following drawings are some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings.

[0018] Figure 1 It is a device structure diagram of a tachycardia diagnosis and treatment device provided by an embodiment of the present application; Figure 2 It is a device structure diagram of another tachycardia diagnosis and treatment device provided by an embodiment of the present application; Figure 3 It is a device structure diagram of another tachycardia diagnosis and treatment device provided by an embodiment of the present application; Figure 4 It is a device structure diagram of another tachycardia diagnosis and treatment device provided by an embodiment of the present application. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0019] To make the objectives, technical solutions, and advantages of the embodiments of the present application clearer, the following will further describe the specific technical solutions of the present application in detail with reference to the drawings in the embodiments of the present application. The following embodiments are used to illustrate the present application but are not intended to limit the scope of the present application.

[0020] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by those of ordinary skill in the technical field to which this application belongs. The terms used herein are only for the purpose of describing the embodiments of the present application and are not intended to limit the present application.

[0021] In the following description, reference is made to "some embodiments", which describe a subset of all possible embodiments. However, it is understood that "some embodiments" may be the same subset or different subsets of all possible embodiments and may be combined with each other without conflict.

[0022] It should be noted that the terms "first / second / third" involved in the embodiments of the present application are used to distinguish similar or different objects and do not represent a specific order for the objects. It is understood that "first / second / third" can be interchanged with a specific order or sequence when permitted, so that the embodiments of the present application described herein can be implemented in an order other than that illustrated or described herein.

[0023] In view of this, the embodiments of the present application provide a device structure diagram of a tachycardia diagnosis and treatment device, as Figure 1 shown, including an ultrasonic transducer module 101, a control module 102, a strain rate module 103, a three-dimensional imaging module 104, a spiral wave endpoint positioning module 105, and a diagnosis module 106 connected in sequence. The diagnosis module 106 is also connected to the control module 102, where: The ultrasonic transducer module 101 is configured to transmit ultrasonic waves after receiving the transmission signal sent by the control module, and receive the echo signal after the ultrasonic waves are reflected by the target area, and send the echo signal to the control module.

[0024] In the embodiments of the present application, the ultrasonic transducer module 101 may be a transducer array composed of a plurality of piezoelectric ceramic elements. Each element can convert an electrical signal into an ultrasonic wave signal and transmit it, receive the echo signal after the ultrasonic wave signal is reflected by the target area, and can also convert the received ultrasonic wave into an electrical signal. The control module 102 generates a transmission signal, and this signal determines key parameters such as the frequency, amplitude, and phase of the ultrasonic wave transmission. In the present application, the operating frequency range of the ultrasonic transducer can be set to 1 MHz to 10 MHz. For example, the transmission frequency is set to 5 MHz. This frequency can not only ensure the effective penetration of ultrasonic waves through human tissues but also obtain a relatively clear image of the target area; the amplitude is set to an appropriate value, such as 10 Vpp, after calibration according to the patient's physical condition and imaging depth requirements, to ensure that the transmitted ultrasonic wave has sufficient energy to reach the target area. In the present application, there are at least two groups of ultrasonic transducers.

[0025] The control module 102 is connected to the ultrasonic transducer module 101 through a high-speed data line, and accurately transmits the transmission signal to the ultrasonic transducer module 101. The drive circuit in the ultrasonic transducer module 101 amplifies the received transmission signal in power to ensure that the signal can effectively drive the piezoelectric ceramic element. Then, according to the timing control information in the transmission signal, the drive circuit sequentially triggers each piezoelectric ceramic element in the transducer array. For example, for imaging of a target area, a linear scanning mode can be adopted, that is, the piezoelectric ceramic elements arranged linearly are sequentially triggered in order, so that the ultrasonic waves generated by them are superimposed to form a beam with a specific directivity, and the beam is focused on a specific area of the target area to ensure the accuracy and clarity of imaging.

[0026] The control module 102 is used to send a transmission signal to the ultrasonic transducer module 101 and receive an echo signal, generate three-dimensional displacement field data of the target area according to the echo signal, and send the three-dimensional displacement field data to the strain rate module 103.

[0027] In the embodiment of the present application, the control module 102 determines the parameters of the transmission signal according to the needs of tachycardia diagnosis and treatment and a preset scanning scheme, so as to balance the penetration depth and imaging resolution of ultrasonic waves.

[0028] Through digital beamforming technology, the control module 102 precisely controls the transmission signals of each array element. By calculating the delay time and phase of different array elements, the transmitted ultrasonic waves are focused on a specific area of the target area in space. For example, using a phased array algorithm, according to the approximate position and depth of the target area, the transmission time of each array element is adjusted, so that the ultrasonic waves form a focal point at the target position, improving the clarity of imaging.

[0029] Through a high-speed data transmission interface, the control module 102 sends the generated transmission signal to the ultrasonic transducer module 101. After receiving the echo signal reflected from the target area, the control module 102 preprocesses the echo signal and then performs acquisition.

[0030] The control module 102 analyzes the preprocessed echo signal and reconstructs a three-dimensional image of the target area using ultrasonic imaging algorithms (such as synthetic aperture imaging algorithm, coherent plane wave imaging algorithm, etc.). In this process, by comparing the phase and amplitude changes of the echo signals at different times, the displacement information of the target area tissue is determined.

[0031] According to the displacement information obtained by analyzing the echo signal, the control module 102 calculates the displacement field of the target area in three-dimensional space. The cross-correlation algorithm or phase difference algorithm can be used to calculate the displacement of the tissue. For example, by calculating the cross-correlation coefficient of the corresponding points in two adjacent frames of images, the maximum correlation point is found, thereby determining the displacement of the tissue. Integrating the displacement information at different positions, the three-dimensional displacement field data of the target area is obtained.

[0032] The control module 102 sends the generated three-dimensional displacement field data of the target area to the strain rate module 103 through the data transmission interface.

[0033] The strain rate module 103 is used to calculate the deformation tensor for the received three-dimensional displacement field data to obtain the strain rate, and send the strain rate to the three-dimensional imaging module 104.

[0034] In the embodiment of the present application, the strain rate module 103 receives the three-dimensional displacement field data through the data transmission interface with the control module 102.

[0035] In order to make the three-dimensional displacement field data more continuous in space and facilitate subsequent calculations, interpolation processing can also be performed on the three-dimensional displacement field data in the present application. Common interpolation methods include linear interpolation, cubic spline interpolation, etc. Through interpolation processing, reasonable displacement values can be estimated between data points, improving the resolution of the data.

[0036] The strain rate module 103 calculates the deformation tensor of the target area using the deformation tensor calculation formula according to the relevant theories in continuum mechanics. For an object in three-dimensional space, the deformation tensor is a second-order tensor, usually represented by a 3×3 matrix. Its calculation is based on the spatial derivatives of the displacement field data. For example, the partial derivatives of the displacement field in the x, y, and z directions constitute the elements of the deformation tensor.

[0037] In actual calculations, since the displacement field data is discrete, numerical difference methods need to be used to approximately calculate the derivatives of the displacement field. For example, the central difference method is used to calculate the first-order partial derivatives, and the derivative values are estimated through the displacement differences between adjacent data points.

[0038] According to the relationship between the deformation tensor and the strain rate, the strain rate module 103 calculates the strain rate of the target area. The strain rate describes the deformation speed of an object per unit time, and it is related to the time derivative of the deformation tensor. In actual calculations, the strain rate can be obtained by performing difference calculations on the deformation tensors at different times.

[0039] The strain rate usually includes multiple components, such as longitudinal strain rate, transverse strain rate, and shear strain rate, etc. The strain rate module 103 calculates these components respectively to comprehensively describe the deformation of the target area. For example, the longitudinal strain rate reflects the stretching or compression speed of the target area tissue in the long axis direction, the transverse strain rate reflects the deformation speed in the short axis direction, and the shear strain rate reflects the degree of shear deformation of the tissue.

[0040] The strain rate module 103 sends the calculated strain rate data of the target area to the three-dimensional imaging module 104 through the data transmission interface.

[0041] A three-dimensional imaging module 104 is configured to obtain a strain rate data set based on the received strain rate, perform three-dimensional mapping on the strain rate data set to obtain three-dimensional imaging data, and send the three-dimensional imaging data to a spiral wave endpoint localization module 105.

[0042] In an embodiment of the present application, the three-dimensional imaging module 104 may receive strain rate data of a target area at different times from the strain rate module 103 according to a preset time interval. Specifically, the preset time interval in the present application is 10 ms. It may also be to obtain all the strain rate data, which is not limited in the present application. The three-dimensional imaging module 104 is connected to the strain rate module 103 through a high-speed and stable data transmission interface.

[0043] The three-dimensional imaging module 104 may set up a dedicated database to store the collected strain rate data. The database adopts an efficient data structure, such as a relational database (such as MySQL) or a non-relational database (such as MongoDB), to optimize the storage and query efficiency of the data.

[0044] Organize the data according to time sequence and spatial position, and each strain rate data point will be associated with its corresponding timestamp and three-dimensional spatial coordinates.

[0045] Divide the three-dimensional space where the target area is located into regular grids, such as cubic grids. Each grid cell has a unique coordinate identifier, corresponding to a specific area within the target area.

[0046] For each grid cell, calculate the strain rate value of the grid cell using an interpolation algorithm (such as trilinear interpolation) based on the existing strain rate data points around it. In this way, the discrete strain rate data points are extended to a continuous three-dimensional strain rate distribution.

[0047] Map the calculated strain rate value of each grid cell to the corresponding position in the three-dimensional space to form three-dimensional imaging data of the strain rate data set. This three-dimensional imaging data can intuitively display the strain rate distribution of the target area at different times.

[0048] After completing the three-dimensional mapping, the three-dimensional imaging module 104 will send the generated three-dimensional imaging data to the spiral wave endpoint localization module 105.

[0049] A spiral wave endpoint localization module 105 is configured to generate multiple sub-images based on the received three-dimensional imaging data, perform spiral wave endpoint detection on the multiple sub-images to obtain spiral wave endpoints of the multiple sub-images, perform clustering and screening on the spiral wave endpoints of the multiple sub-images to obtain a filtered set of spiral wave endpoints, and send the filtered set of spiral wave endpoints to a diagnosis module 106.

[0050] In the embodiment of the present application, the spiral wave endpoint localization module 105 is connected to the three-dimensional imaging module 104 through a high-speed and stable data interface to receive three-dimensional imaging data.

[0051] According to the anatomical structure and imaging characteristics of the target area, the three-dimensional imaging data is segmented according to certain rules. For example, it can be divided into cube sub-blocks of the same size according to spatial position, and each sub-block corresponds to a sub-image. It can also be segmented according to the functional characteristics of different regions of the target area, such as the atrium, ventricle, etc.

[0052] Reasonably set the size and overlap rate of the sub-images. If the sub-images are too large, it may lead to inaccurate detection of spiral wave endpoints; if they are too small, the computational complexity will increase. Generally, the size of the sub-images can be set to a side length of 10 - 20 voxels, and the overlap rate can be set to 10% - 30% to ensure the continuity of information between adjacent sub-images.

[0053] Extract features related to spiral wave endpoints, such as phase gradient, curvature, etc., for each sub-image. For example, by calculating the phase values of each point in the image, the phase gradient can be obtained, and there are usually obvious changes in the phase gradient at the spiral wave endpoints.

[0054] Adopt a special detection algorithm to identify spiral wave endpoints, such as an algorithm based on phase singularity detection. This algorithm analyzes the continuity of the phase in the sub-image and finds the points where the phase undergoes a 2π jump, and these points are the possible spiral wave endpoints.

[0055] Use a clustering algorithm (such as a density-based clustering algorithm) to cluster the spiral wave endpoints detected in all sub-images. The density-based clustering algorithm can divide them into different clusters according to the distance and density between the endpoints.

[0056] Set screening criteria to remove noise points and misdetected points. For example, for clusters with too few points, they can be considered as clusters composed of noise points and be removed; at the same time, remove the endpoints located at the image edge or in positions that do not conform to physiological characteristics.

[0057] Send the screened set of spiral wave endpoints to the diagnostic module in a specific data format.

[0058] The diagnostic module 106 is used to determine whether there is tachycardia in the target area according to a preset frequency threshold and the received screened set of spiral wave endpoints; when there is tachycardia in the target area, the diagnostic module sends a treatment instruction and the screened set of spiral wave endpoints to the control module.

[0059] In the embodiment of the present application, the diagnostic module is connected to the spiral wave endpoint localization module 105 through a reliable data transmission channel to receive the screened set of spiral wave endpoints.

[0060] Preset a period for counting the occurrence frequency of spiral wave endpoints, for example, set it to 1 second. This period duration can not only effectively capture the dynamic changes of spiral wave endpoints, reflecting the instantaneous electrophysiological state of the target area, but also avoid misjudgment due to being too long or too short.

[0061] During the preset period, the diagnostic module monitors and counts the received set of spiral wave endpoints in real time. Whenever a new spiral wave endpoint data point enters the set, the corresponding counter is incremented by 1. As time progresses, based on the accumulated value of the counter and the preset period duration, according to the formula "frequency = number of endpoint occurrences ÷ preset period", the occurrence frequency of spiral wave endpoints is accurately calculated.

[0062] Compare the calculated occurrence frequency of spiral wave endpoints with the preset frequency threshold. In this application, if the occurrence frequency of spiral wave endpoints is greater than 85% of the occurrence frequency within three consecutive heartbeat cycles, it means that the electrical activity disorder in the target area is relatively frequent, and it is sufficient to determine that there is tachycardia in the target area; conversely, if it is less than or equal to the threshold, it is determined that there is no tachycardia in the target area and it is in a relatively normal rhythm state, and there is no need to send a treatment instruction and the filtered set of spiral wave endpoints to the control module 102.

[0063] Once the judgment is completed, the diagnostic module immediately outputs the judgment result in a standardized data format (such as a boolean value, "true" indicates the presence of tachycardia, "false" indicates the absence). This result can be directly called by other modules (such as the control module 102) in the subsequent diagnosis and treatment process to quickly initiate corresponding treatment or monitoring strategies.

[0064] After the diagnostic module determines that there is tachycardia in the target area, it will perform an internal verification of the judgment result again to ensure the accuracy of the judgment. At the same time, organize and screen the data of the set of spiral wave endpoints, check the integrity and accuracy of the data, and perform necessary format conversion and encoding on the data so that it can be stably transmitted between modules.

[0065] The diagnostic module generates a treatment instruction according to the preset treatment plan and the specific situation of the current tachycardia. The treatment instruction contains various key information, such as the start time of treatment, duration, approximate energy intensity range of ultrasonic waves, etc. For example, for mild tachycardia, the treatment duration may be set to 10 seconds and the ultrasonic energy intensity is at a medium level; for a more serious situation, the treatment duration may be extended to 20 seconds and the energy intensity is correspondingly increased.

[0066] Package the generated treatment instruction and the filtered set of spiral wave endpoints to form a complete data packet. Send the data packet to the control module 102 through a stable and reliable data transmission channel.

[0067] The control module 102 is further configured to generate control information according to the received filtered set of spiral wave endpoints after receiving the treatment instruction of the control module, and send the control information to the ultrasonic transducer module 101. The control information includes the sound intensity, phase of the ultrasonic wave, and the position of the spiral wave endpoints.

[0068] In the embodiment of the present application, the control module 102 receives the data packet from the diagnostic module through the internal high-speed data bus. The data packet contains the treatment instruction and the filtered set of spiral wave endpoints.

[0069] The control module 102 matches the most suitable treatment plan from the preset treatment plan database according to the parameters in the received treatment instruction (such as the treatment start time, duration, and approximate energy intensity range). For example, if the treatment instruction indicates mild tachycardia and the set treatment duration is 10 seconds and the energy intensity is at a medium level, the module will select the corresponding detailed treatment parameter combination.

[0070] According to the positions of the endpoints in the filtered set of spiral wave endpoints and the severity of the tachycardia, combined with the matched treatment plan, accurately calculate the ultrasonic sound intensity required for each ultrasonic transducer element. For the endpoints close to the key parts of the target area or in the more severe tachycardia area, the sound intensity can be appropriately increased.

[0071] In order to enable the ultrasonic wave to accurately focus on the spiral wave endpoint position, the control module 102 uses the phased array principle to calculate the phase of the ultrasonic wave emitted by each element according to the three-dimensional spatial position of the endpoint and the layout of the ultrasonic transducer elements. By precisely adjusting the phase, the ultrasonic waves emitted by each element are in-phase superimposed at the spiral wave endpoint, enhancing the focusing effect.

[0072] Associate the calculated sound intensity, phase with the specific position information of the spiral wave endpoints to form a complete set of treatment control information. Each spiral wave endpoint corresponds to a specific set of sound intensity, phase, and position parameters.

[0073] Organize and encapsulate the generated treatment control information in a format recognizable by the ultrasonic transducer module 101. Then send the data packet to the ultrasonic transducer module 101 through the high-speed data interface.

[0074] The ultrasonic transducer module 101 is configured to adjust the sound intensity and phase of the ultrasonic wave according to the received control information and then emit the ultrasonic wave towards the target area to perform ultrasonic focusing on the position of the spiral wave endpoints. In the present application, the ultrasonic focusing time can be set to 20 ms.

[0075] In the embodiment of the present application, the ultrasonic transducer module 101 receives the treatment control information data packet sent by the control module 102 through the high-speed data interface and parses the treatment control information data packet.

[0076] The power adjustment circuit inside the ultrasonic transducer module 101 precisely adjusts the amplitude of the drive signal for each ultrasonic transducer element according to the parsed sound intensity information. For example, by changing the gain of the power amplifier, the sound intensity of the ultrasonic waves emitted by each element meets the requirements of the treatment control information. At the same time, a real-time monitoring and feedback mechanism is adopted to ensure the stability and accuracy of the sound intensity.

[0077] Using the phase control circuit, according to the parsed phase information, the phase of the drive signal for each ultrasonic transducer element is adjusted. By controlling the delay time of the signal, precise phase adjustment is achieved to ensure that the ultrasonic waves emitted by each element are in-phase superimposed at the spiral wave endpoint position.

[0078] Based on the position information of the spiral wave endpoint, the ultrasonic transducer module 101 adopts a phased array focusing algorithm to coordinate the emission time and parameters of each ultrasonic transducer element. So that the ultrasonic waves emitted by all elements are focused in space to the specified spiral wave endpoint position, forming a high-intensity ultrasonic beam.

[0079] During the process of emitting ultrasonic waves, the ultrasonic transducer module 101 monitors the emission state in real time, including parameters such as sound intensity, phase, and emission frequency. If it is found that the parameters deviate from the preset values, adjustments are made in a timely manner to ensure the ultrasonic focusing effect on the spiral wave endpoint. At the same time, according to the parameters such as the treatment start time and duration specified in the treatment instruction, the emission duration of the ultrasonic waves is precisely controlled.

[0080] The embodiment of the present application realizes non-invasive diagnosis and treatment through the whole-process ultrasonic technology, avoiding the risks of invasive methods. Through three-dimensional displacement field, strain rate analysis, and spiral wave endpoint clustering and screening, the target area is accurately located and tachycardia is quantitatively diagnosed. Multi-module collaborative processing and real-time data feedback achieve non-invasive diagnosis of tachycardia, with both accuracy and safety.

[0081] In the above Figure 1 shown situation, the embodiment of the present application also provides a device structure diagram of a tachycardia diagnosis and treatment device. As Figure 2 shown, the strain rate module 103 includes a first receiving unit 1031, a first calculating unit 1032, a conversion unit 1033, and a first sending unit 1034 connected in sequence, where: The first receiving unit 1031 is used to receive the three-dimensional displacement field data sent by the control module 102.

[0082] In the embodiment of the present application, the first receiving unit 1031 receives the three-dimensional displacement field data through the high-speed data interface with the control module 102.

[0083] The first calculation unit 1032 is configured to calculate the deformation tensor for the three-dimensional displacement field data to obtain the gradient deformation tensor and the Green-Lagrange strain tensor; perform eigenvalue decomposition and invariant combination operations on the gradient deformation tensor and the Green-Lagrange strain tensor respectively to obtain the tensor invariants of the target region.

[0084] In the embodiment of the present application, the first calculation unit 1032 performs spatial derivative calculation on the three-dimensional displacement field data to obtain the gradient deformation tensor. In actual calculation, since the displacement field data is discrete, numerical difference methods are usually used to approximately calculate the derivative. For example, the central difference method is used to calculate the first-order partial derivatives of the displacement field in each direction, and these partial derivatives constitute the elements of the gradient deformation tensor. For a three-dimensional space, the gradient deformation tensor is a 3×3 matrix, and its elements reflect the change rate of the displacement field in different directions.

[0085] According to the gradient deformation tensor, the calculation is performed using the calculation formula of the Green-Lagrange strain tensor. The Green-Lagrange strain tensor is a symmetric second-order tensor that describes the strain state of an object during the deformation process and takes into account the geometric relationship before and after the object's deformation. The calculation formula involves the product and summation operations of the gradient deformation tensor. By performing specific combinations and operations on the elements of the gradient deformation tensor, the elements of the Green-Lagrange strain tensor are obtained.

[0086] The first calculation unit 1032 performs eigenvalue decomposition on the calculated gradient deformation tensor and Green-Lagrange strain tensor respectively. The process of decomposing a tensor into eigenvalues and eigenvectors, where the eigenvalues reflect the degree of stretching or deformation of the tensor in different directions, and the eigenvectors represent these directions. Through eigenvalue decomposition, the principal strains and principal strain directions of the tensor can be obtained, providing important information for subsequent analysis.

[0087] According to the results obtained from the eigenvalue decomposition, invariant combination operations are performed. The invariants of a tensor refer to the quantities that remain unchanged under coordinate system transformation and reflect the inherent characteristics of the tensor. For a second-order tensor, there are three basic invariants. By performing specific combination operations (such as summation, multiplication, etc.) on the eigenvalues, the invariants can be obtained. These invariants can comprehensively reflect the deformation characteristics of the target region tissue and provide a basis for subsequent strain rate calculation.

[0088] The conversion unit 1033 is configured to perform strain rate conversion on the tensor invariants to obtain the strain rate of the target region.

[0089] In an embodiment of the present application, the conversion unit 1033 performs strain rate conversion according to the obtained tensor invariants by using a pre-established strain rate conversion model or formula. For example, a generalized Maxwell model or a Kelvin-Voigt model can be adopted to relate stress, strain, and strain rate, and the strain rate is solved through known tensor invariants.

[0090] Taking the generalized Maxwell model as an example, this model is composed of multiple Maxwell units connected in series, and each Maxwell unit contains a spring and a damper. Through mechanical analysis of the model, the relational expression of stress-strain-strain rate can be obtained. Given the material parameters of the target region tissue and the tensor invariants calculated through the previous steps, substituting them into the constitutive equation to obtain the strain rate.

[0091] The first sending unit 1034 is configured to send the strain rate of the target region to the three-dimensional imaging module 104.

[0092] In an embodiment of the present application, the first sending unit 1034 formats the calculated strain rate data of the target region to meet the receiving requirements of the three-dimensional imaging module 104. Then, the strain rate data is sent to the three-dimensional imaging module 104 through the data transmission interface.

[0093] In the embodiment of the present application, through the calculation of the deformation tensor, the deformation of the target region tissue is accurately described, providing an accurate quantitative basis for diagnosis and treatment. The eigenvalue decomposition and invariant combination operation yield tensor invariants, mining the deep essential features of the target region deformation, stably and reliably assisting in diagnosis, and being unaffected by measurement conditions.

[0094] In the above Figure 2 In the case shown, the embodiment of the present application further provides a tachycardia diagnosis and treatment device diagram. As Figure 3 shown, the spiral wave endpoint positioning module 105 includes a second receiving unit 1051, a establishing unit 1052, a cutting unit 1053, a second calculating unit 1054, a screening unit 1055, and a second sending unit 1056 that are connected in sequence, where: The second receiving unit 1051 is configured to receive the three-dimensional imaging data sent by the three-dimensional imaging module 104.

[0095] In an embodiment of the present application, the second receiving unit 1051 is connected to the three-dimensional imaging module 104 through a high-speed and stable data transmission interface to receive the three-dimensional imaging data sent by it.

[0096] The establishing unit 1052 is configured to establish a three-dimensional rectangular coordinate system according to the three-dimensional imaging data.

[0097] In the embodiment of the present application, the establishment unit 1052 establishes a three-dimensional rectangular coordinate system inside the module according to the characteristics of the three-dimensional imaging data and the general anatomical structure of the target area. Taking a fixed point in the target area (such as the center of the target area or a specific anatomical landmark point) as the origin of the coordinate system, the directions of the x, y, and z axes are determined. For example, the long axis direction of the target area is set as the x axis, the short axis direction is set as the y axis, and the direction perpendicular to the plane of the target area is set as the z axis.

[0098] The cutting unit 1053 is used to cut the three-dimensional imaging data perpendicular to the coordinate axes along the x, y, and z axes of the three-dimensional rectangular coordinate system, generating a plurality of sub-images.

[0099] In the embodiment of the present application, the cutting unit 1053 cuts the three-dimensional imaging data perpendicular to the coordinate axes along the x, y, and z axes of the established three-dimensional rectangular coordinate system. When cutting, a suitable cutting interval is set according to the actual requirements and the limitations of computing resources. For example, cutting is performed every 5 voxels (or adjusted according to the imaging resolution and data volume), so as to divide the three-dimensional imaging data into a plurality of sub-images with the same or similar sizes. Each sub-image contains the imaging information of a part of the target area, which is convenient for subsequent detailed detection of the spiral wave endpoints.

[0100] The second calculation unit 1054 is used to adopt the Moore neighborhood tracking algorithm for each sub-image among the plurality of sub-images to obtain the intersection line between the foreground and the background in each sub-image, where the foreground is the area containing the spiral wave in the sub-image and the background is the area without the spiral wave in the sub-image; when the intersection line meets the preset curvature threshold, the spiral wave endpoints of the sub-image corresponding to the intersection line are obtained through the fitting curvature circle algorithm.

[0101] In the embodiment of the present application, for each generated sub-image, the second calculation unit 1054 applies the Moore neighborhood tracking algorithm to determine the intersection line between the foreground (the area containing the spiral wave) and the background (the area without the spiral wave) therein. Specifically, first, a starting point is selected in the sub-image, usually an edge point or a point with obvious features is selected as the starting point. Then, according to the rules of Moore neighborhood tracking, the neighborhood points (such as 8-neighborhood or 4-neighborhood) around the starting point are sequentially checked, and it is determined whether the neighborhood point belongs to the foreground or the background according to the preset gray value or other feature thresholds. The boundary between the foreground and the background is traced point by point until the starting point is returned, thereby obtaining the complete intersection line.

[0102] After the second calculation unit 1054 obtains the intersection lines between the foreground and the background in each sub-image, it calculates the curvatures of the points on the intersection lines. By setting a preset curvature threshold, the points on the intersection lines whose curvatures satisfy the threshold are screened out. For these qualified points, the fitting curvature circle algorithm is used, that is, with this point as the center, by fitting the points within a certain range around it, an approximate curvature circle is obtained. The center of this curvature circle is the corresponding spiral wave endpoint. In this way, the spiral wave endpoints can be accurately extracted from the intersection lines.

[0103] The screening unit 1055 is used to perform clustering screening on the spiral wave endpoints of multiple sub-images to obtain a set of screened spiral wave endpoints.

[0104] In the embodiment of the present application, the screening unit 1055 performs clustering screening on all the spiral wave endpoints obtained from multiple sub-images. Using a clustering algorithm, according to the spatial distance and density distribution between the spiral wave endpoints, they are divided into different clusters. Appropriate clustering parameters (such as the neighborhood radius and the minimum number of points) are set to remove noise points and isolated points (that is, those points that do not belong to any obvious cluster), and the representative spiral wave endpoints are retained, so as to obtain a set of screened spiral wave endpoints.

[0105] The second sending unit 1056 is used to send the set of screened spiral wave endpoints to the diagnosis module.

[0106] In the embodiment of the present application, the second sending unit 1056 formats the set of screened spiral wave endpoints so that it conforms to the data format that the diagnosis module can receive and process. Then, the set is sent to the diagnosis module through the data transmission interface.

[0107] The embodiment of the present application combines the segmentation of sub-images with specific algorithms to reduce data complexity and interference, and improve the detection accuracy of spiral wave endpoints. With the help of the three-dimensional coordinate system and cutting operations, the spatial positions of the endpoints are clarified for accurate positioning of treatment. Clustering screening removes noise and misdetected endpoints, ensures the quality of the endpoint set, and reduces interference.

[0108] In the above Figure 3 shown case, the embodiment of the present application also provides a tachycardia diagnosis and treatment device diagram. As Figure 4 shown, the screening unit 1055 includes an acquisition element and a screening element 10552 connected in sequence, where: The acquisition element 10551 is used to acquire the corresponding multiple spiral wave endpoints in multiple sub-images to obtain an initial set of spiral wave endpoints.

[0109] In an embodiment of the present application, after the component 10551 completes the detection of the spiral wave endpoints of multiple sub-images, the detected spiral wave endpoints in each sub-image are summarized. Each spiral wave endpoint has corresponding coordinate values (x, y, z) in a three-dimensional rectangular coordinate system, and these coordinate values clarify its position in the three-dimensional space of the target area. By integrating the coordinate information of the spiral wave endpoints of all sub-images, an initial spiral wave endpoint set is obtained.

[0110] The screening component 10552 is used to screen out the spiral wave endpoints in the initial spiral wave endpoint set whose standard deviation of spatial distribution exceeds a preset distance threshold, and a screened spiral wave endpoint set is obtained.

[0111] In an embodiment of the present application, for all endpoints in the initial spiral wave endpoint set, the screening component 10552 needs to calculate their standard deviation of spatial distribution. The specific steps are as follows: Calculate the average value of the coordinate values of all endpoints on the x, y, and z axes respectively to obtain the spatial mean point . The calculation formula is: (1).

[0112] where n is the number of endpoints in the initial spiral wave endpoint set, is the coordinate of the i-th endpoint.

[0113] For each endpoint , calculate the sum of the squares of its distances to the mean point , and then calculate the average value to obtain the variance .

[0114] The calculation formula is: (2).

[0115] Take the square root of the variance to obtain the standard deviation of spatial distribution .

[0116] The preset distance threshold is a preset value. In the present application, the preset distance threshold is 0.5 mm. This threshold is used to determine whether the spatial distribution of the spiral wave endpoints is too scattered. If the standard deviation of the spatial distribution of the endpoints exceeds this threshold, it indicates that there may be some abnormal endpoints far from the main distribution area, and these endpoints are likely to be caused by noise, misdetection, etc.

[0117] Compare the calculated standard deviation of spatial distribution with the preset distance threshold. For each endpoint in the initial spiral wave endpoint set, calculate its distance to the mean point.

[0118] (3).

[0119] If the distance corresponding to a certain endpoint causes the standard deviation of the spatial distribution of the entire set to exceed a preset distance threshold, then that endpoint is screened out from the initial set of spiral wave endpoints. Repeat this process until the standard deviation of the spatial distribution of the remaining set of endpoints does not exceed the preset distance threshold. At this time, the obtained set is the screened set of spiral wave endpoints.

[0120] The screened set of spiral wave endpoints contains the spiral wave endpoints that are retained after screening, with relatively concentrated spatial distribution and representativeness. Output this set in a suitable data format (such as a list form, where each element is the three-dimensional coordinates of an endpoint), and it can be sent to the diagnostic module later for further tachycardia diagnosis and treatment plan formulation.

[0121] In the embodiments of the present application, by screening out abnormal endpoints, removing noise and false detection interference, the data quality of the set of spiral wave endpoints is improved, and the actual distribution of spiral waves is accurately presented. The abnormal area is accurately located, making the treatment plan more targeted and personalized, and improving the treatment effect.

[0122] It should be noted that in each embodiment of the present application, each functional unit can be integrated in one processing unit, can exist physically alone, or two or more units can be integrated in one unit.

[0123] It should be understood that the "one embodiment" or "an embodiment" or "some embodiments" mentioned throughout the specification means that the specific features, structures or characteristics related to the embodiment are included in at least one embodiment of the present application. Therefore, the "in one embodiment" or "in an embodiment" or "in some embodiments" that appear throughout the specification do not necessarily refer to the same embodiment. In addition, these specific features, structures or characteristics can be combined in one or more embodiments in any suitable manner. It should be understood that in various embodiments of the present application, the order numbers of the above processes do not mean the order of execution. The order of execution of each process should be determined by its function and internal logic, and should not constitute any limitation to the implementation process of the embodiments of the present application. The serial numbers of the embodiments of the present application above are only for description and do not represent the advantages and disadvantages of the embodiments. The above descriptions of the various embodiments tend to emphasize the differences between the various embodiments, and the same or similar parts can be referred to each other. For the sake of brevity, they will not be repeated herein.

[0124] The term "and / or" in this article is only a description of the association relationship of the associated objects, indicating that there can be three relationships. For example, object A and / or object B can represent: object A exists alone, object A and object B exist simultaneously, and object B exists alone.

[0125] It should be noted that in this article, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, such that a process, method, article or device comprising a series of elements not only includes those elements but also includes other elements not expressly listed, or further includes elements inherent to such process, method, article or device. Without further limitation, an element defined by the statement "comprising one..." does not exclude the existence of additional identical elements in the process, method, article or device comprising such element.

[0126] In several embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. The above-described embodiments are merely illustrative. For example, the division of the modules is only a logical function division. In actual implementation, there may be other division methods, such as multiple modules or components can be combined, or can be integrated into another system, or some features can be ignored or not executed. Additionally, the couplings, direct couplings, or communication connections between the various components shown or discussed can be through some interfaces. The indirect couplings or communication connections of devices or modules can be electrical, mechanical, or in other forms.

[0127] The modules described above as separate components may or may not be physically separated, and the components shown as modules may or may not be physical modules; they can be located in one place or distributed to multiple network units; some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0128] Furthermore, in each embodiment of this application, the various functional modules can all be integrated in one processing unit, or each module can be separately a unit alone, or two or more modules can be integrated in one unit; the above-mentioned integrated modules can be implemented in the form of hardware, or in the form of a combination of hardware and software functional units.

[0129] Those of ordinary skill in the art can understand that all or part of the steps to implement the above method embodiments can be completed by hardware related to program instructions. The foregoing program can be stored in a computer-readable storage medium. When the program is executed, it performs the steps including the above method embodiments; and the foregoing storage medium includes various media that can store program codes, such as removable storage devices, read-only memory (ROM), magnetic disks, or optical discs.

[0130] Alternatively, if the above integrated units of the present application are implemented in the form of software functional modules and sold or used as independent products, they can also be stored in a computer-readable storage medium. Based on such an understanding, the technical solution of the embodiments of the present application essentially or the part that contributes to the related art can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for causing an electronic device to execute all or part of the methods described in the various embodiments of the present application. The foregoing storage medium includes: various media such as removable storage devices, ROMs, magnetic disks, or optical discs that can store program codes.

[0131] The features disclosed in the several method or device embodiments provided by the present application can be combined arbitrarily without conflict to obtain new method embodiments or device embodiments.

[0132] The above is only the implementation manner of the present application, but the protection scope of the present application is not limited thereto. Any person skilled in the art within the technical scope disclosed by the present application can easily think of changes or substitutions, which should all be covered by the protection scope of the present application. Therefore, the protection scope of the present application shall be subject to the protection scope of the claims.

Claims

1. A tachycardia diagnosis and treatment device, characterized in that: It includes an ultrasonic transducer module, a control module, a strain rate module, a three-dimensional imaging module, a spiral wave endpoint positioning module and a diagnosis module connected in sequence, wherein the diagnosis module is also connected to the control module, wherein: The ultrasonic transducer module is used to transmit ultrasonic waves after receiving the transmission signal sent by the control module, and receive the echo signal after the ultrasonic waves are reflected by the target area, and send the echo signal to the control module; The control module is used to send a transmission signal to the ultrasonic transducer module and receive the echo signal, generate three-dimensional displacement field data corresponding to the target area according to the echo signal, and send the three-dimensional displacement field data to the strain rate module; The strain rate module is used to perform deformation tensor calculation on the received three-dimensional displacement field data to obtain the strain rate, and send the strain rate to the three-dimensional imaging module; The three-dimensional imaging module is used to obtain a strain rate data set according to the received strain rate, perform three-dimensional mapping on the strain rate data set to obtain three-dimensional imaging data, and send the three-dimensional imaging data to the spiral wave endpoint positioning module; The spiral wave endpoint positioning module is used to generate a plurality of sub-images according to the received three-dimensional imaging data, perform spiral wave endpoint detection on the plurality of sub-images to obtain spiral wave endpoints of the plurality of sub-images, perform clustering screening on the spiral wave endpoints of the plurality of sub-images to obtain a screened spiral wave endpoint set, and send the screened spiral wave endpoint set to the diagnosis module; The diagnostic module is used to determine whether tachycardia exists in the target area according to a preset frequency threshold and the received filtered spiral wave endpoint set; when tachycardia exists in the target area, the diagnostic module sends a treatment instruction and the filtered spiral wave endpoint set to the control module; The control module is further configured to generate control information according to the received filtered spiral wave endpoint set after receiving the treatment instruction of the control module, and send the control information to the ultrasonic transducer module, wherein the control information includes the sound intensity and phase of the ultrasonic wave and the position of the spiral wave endpoint; The ultrasonic transducer module is also used to adjust the sound intensity and phase of the ultrasonic wave according to the received control information and then transmit the ultrasonic wave to the target area to perform ultrasonic focusing on the position of the spiral wave endpoint.

2. A tachycardia diagnosis and treatment device according to claim 1, characterized in that: The strain rate module comprises a first receiving unit, a first calculating unit, a converting unit and a first sending unit which are connected in sequence, wherein; The first receiving unit is used to receive the three-dimensional displacement field data of the control module; The first calculation unit is used to perform deformation tensor calculation on the three-dimensional displacement field data to obtain a gradient deformation tensor and a Green-Lagrange strain tensor; perform eigenvalue decomposition and invariant combination operations on the gradient deformation tensor and the Green-Lagrange strain tensor respectively to obtain tensor invariants of the target area; The conversion unit is used to perform strain rate conversion on the tensor invariant to obtain the strain rate of the target area; The first sending unit is used to send the strain rate of the target area to the three-dimensional imaging module.

3. The tachycardia diagnosis and treatment device according to claim 1, characterized in that: The spiral wave endpoint positioning module includes a second receiving unit, a establishing unit, a cutting unit, a second calculating unit, a screening unit and a second sending unit which are connected in sequence, wherein; The second receiving unit is used to receive the three-dimensional imaging data of the three-dimensional imaging module; The establishing unit is used to establish a three-dimensional rectangular coordinate system according to the three-dimensional imaging data; The cutting unit is used to cut the three-dimensional imaging data along the x, y, and z axis directions of the three-dimensional rectangular coordinate system respectively perpendicular to the coordinate axis to generate a plurality of sub-images; The second calculation unit is used to use the Moore field tracking algorithm for each sub-image in the multiple sub-images to obtain the intersection line between the foreground and the background in each sub-image, wherein the foreground is the area in the sub-image containing the spiral wave, and the background is the area in the sub-image without the spiral wave; when the intersection line meets the preset curvature threshold, the spiral wave endpoint of the sub-image corresponding to the intersection line is obtained by fitting the curvature circle algorithm; The screening unit is used to perform cluster screening on the spiral wave endpoints of the plurality of sub-images to obtain a screened spiral wave endpoint set; The second sending unit is used to send the screened spiral wave endpoint set to the diagnosis module.

4. A tachycardia diagnosis and treatment device according to claim 3, characterized in that: The screening unit comprises a sequentially connected acquisition element and a screening element, wherein; The acquisition element is used to acquire a plurality of spiral wave endpoints corresponding to the plurality of sub-images to obtain an initial spiral wave endpoint set; The screening element is used to screen out the spiral wave endpoints whose standard deviation of spatial distribution in the initial spiral wave endpoint set exceeds a preset distance threshold, so as to obtain a screened spiral wave endpoint set.

5. The tachycardia diagnosis and treatment device according to claim 4, characterized in that: The preset distance threshold is 0.5 mm.

6. The tachycardia diagnosis and treatment device according to claim 1, characterized in that: The ultrasonic transducer module includes at least two groups of ultrasonic transducers.

7. The tachycardia diagnosis and treatment device according to claim 6, characterized in that: The operating frequency range of the ultrasonic transducer is 1 MHz to 10 MHz.

8. The tachycardia diagnosis and treatment device according to claim 1, characterized in that: The ultrasound focusing time is 20 ms.

9. The tachycardia diagnosis and treatment device according to claim 1, characterized in that: The preset frequency threshold is 85% of the frequency of occurrence of the spiral wave endpoints.

10. The tachycardia diagnosis and treatment device according to claim 6, characterized in that: The imaging interval of the three-dimensional mapping is 10 ms.

Images