Patent Application
20250078288
The various aspects and variations thereof relate to a method and device for assessment of functioning of the myocardium tissue for providing pumping function of the heart.
For assessment of the functioning of the heart and subsequent assessment whether this functioning is in accordance with reference data or not to enable a cardiologist to diagnose any pathological issues, various technologies are available. Such techniques are analysis of an electrocardiogram and assessment of information recorded multiple imaging technologies.
It is preferred to provide an improved technique for evaluation of the functioning of the myocardial tissue and pumping function of the heart and to provide data for such evaluation and, optionally, to provide the data to a cardiologist or other medical practitioner, preferably via an electronic display screen.
A first aspect provides, in an electronic computing device, a method of determining a motion of an heart wall. The method comprises receiving a sequence of images acquired consecutively in time at fixed time intervals, the images comprising image data of at least part of a heart having a heart wall, identifying, in each image, at least part of the heart wall and defining, in each image, at least one vector point on the identified heart wall. The method further comprises determining motion of the vector point over time, based on the sequence and receiving a first reference related to the heart. Based on the tracking of the vector point and the first reference, a motion vector is determined indicating motion of the vector point relative to the reference.
By determining one or more motion vectors, motion of one or more parts of the heart wall may be determined, in particular motion of parts of a wall of a specific cavity, like the left ventricle, may be determined relative to the reference. With the reference being provided within the cavity, motion of the wall relative to, for example, a central point in the cavity may be determined. The magnitude of the motion, in particular in the direction to or from the reference, may provide an indication of a quality of a pumping function of the heart and the particular cavity in particular. A part of the wall may move parallel to the tangential plane of the identified wall. In such case, the movement contributes little or not to the pumping function. With the reference provided in the cavity—in the model—and determining the motion vector indicating motion of the vector points towards or away from the reference—a reference point, a reference curve, or a combination thereof—the part of the movement contributing to pumping may be determined. The identified heart wall may be one of the outer wall of a heart, the inner wall of a specific cavity or another reference, for example between the outer wall and the inner wall and is preferably provided as a curve.
The first aspect may comprise identifying a set of vector points on the heart wall, assigning a first sub-set of the set of vector points to a first region of the heart wall, determining movement of the vector points in the first sub-set over time, based on the determined movement of the vector points in the first sub-set, determining a first region motion vector of the first region indicating motion of the region relative to the first reference.
According to cardiological medical standards, an endocardial wall of a particular cavity may be divided in particular section. Providing motion data per standardised region provides improved input to a medical practitioner for setting a diagnosis.
The first aspect may further comprise assigning a second sub-set of the set of vector points to a second region of the heart wall, determining motion of the vector points in the second sub-set over time, and, based on the tracking of the vector points in the second sub-set, determining a second region motion vector of the second region indicating motion of the region relative to the first reference.
By determining motion vectors for multiple regions, more information may be provided to medical practitioners, allowing them to compare motion of various sections.
The first reference may be different from the second reference. This may allow improved accuracy.
The method my comprise determining, based on the determined motion of at least part of the vector points, a motion of the heart over time, and adjusting the first region motion vector based on the determined motion of the heart to compensate for the motion of the heart.
During heart beat cycle, the heart moves in a body of a subject. With the subject being stabilised during image acquisition, this means that the actual motion of the heart is captured by the images. Hence, the position of the heart may not be the same in each acquired image or acquired frame. By determining a common motion in vector points or a motion of the reference—determined relative to at least one position of the identified heart wall contour —, motion of the heart in particular heart beat and motion of the heart from one frame relative to another frame may be determined. This determined heart motion may be used to adjust the determined motion vectors related to vector points on the identified heart wall, so for example on an identified heart wall contour.
The method may comprise determining an inward vector component of the motion vector, the inward vector component being a component of the motion vector parallel to a line between the first reference point and the vector point. As indicated above, in particular inward motion of the endocardial wall as represented by the identified heart wall is relevant to evaluation of a pumping function of a heart.
The motion vector may be based on a maximum displacement of the vector point relative to the first reference point. Displacement of blood and an amount of blood displaced in a heart beat may be determined by the maximum displacement of sections of or points on the heart wall, between the end of the systolic phase and the end of the diastolic phase.
The reference may be provided between the base and the apex, in particular between one third and two third of that trajectory, more in particular between half that trajectory and a point generally at a third from the base on that trajectory. the reference may comprises single reference point, a multitude of reference points, an open curve, a closed curve, or any combination thereof. A closed curve generally defines an area. An open curve may be provided from a first location closer to the base to a second point further away from the base. If the reference comprises multiple points, such points may be provided on a curve, either closes or open. If the reference comprises multiple points, each reference point may correspond to a vector point on the identified heart wall.
The method may further comprise receiving a reference motion value indicating a reference amount of motion of the vector point relative to the reference point, deriving a vector motion value based on the motion vector, calculating a motion deviation value providing an indication of a difference between the vector motion value and the reference motion value; and providing the motion deviation value for displaying on an electronic display. This provides a medical practitioner with data related to a particular subject, compared to a reference group.
A second aspect provides a computer program product comprising computer executable code causing a processing unit of a computer, when the code is loaded in the processing unit or a memory connected to the processing unit, to execute the method of the first aspect.
A third aspect provides a non-transitory computer readable medium comprising the computer program product of the second aspect.
A fourth aspect provides an electronic computing device comprising a data acquisition module arranged to receive a sequence of images acquired consecutively in time at fixed time intervals, the images comprising image data of at least part of a heart having an heart wall, an image processing module arranged to identify, in each image, at least part of the heart wall, define, in each image, at least one vector point on the heart wall, track the vector point over time, based on the sequence, a motion calculation module arranged to receive a first reference point related to the heart, determine, based on the tracking of the vector point and the reference point, a motion vector indicating motion of the vector point relative to the reference point.
The various aspects and variations thereof will be discussed below in further detail in conjunction with drawings. In the drawings:
The imaging controller 116 is connected to an imaging transmitter 132 and an imaging receiver 134. The imaging transmitter 132 sends out energy—electromagnetic waves like X-ray radiation, sounds waves and ultrasound waves in particular, other energy or a combination thereof—to the imaging receiver 134. The imaging transmitter 132 and the imaging receiver 134 are part of a medical imaging system 130. The medical imaging system 130 may be an MRI system, a CT system, a standard X-Ray system, an ultrasonic imaging system, another imaging system, or a combination thereof. The medical imaging system is arranged to obtain imaging data of a heart 190 of a mammal, a human in particular. The imaging data preferably comprises images divided by approximately constant time intervals.
The imaging controller 116 may be arranged to process the imaging data receiver 134 to provide a two-dimensional representation, a three-dimensional reconstruction of the heart 190 or any other part of a body. The imaging controller 116 may further be arranged to provide two-dimensional cross-sections of the three-dimensional reconstruction. Additionally or alternatively, two-dimensional images as acquired are passed through to the central processing unit 112. Some or all steps may also be executed by the central processing unit 112, instead of or in conjunction with the imaging controller 116. The imaging controller and the central processing unit 112 may be included in one and the same physical entity.
The central processing unit 112 comprises an edge detection subunit 122 as a wall detection unit, a segmenting subunit 124 for segmenting data, a motion detection subunit 126 for tracking of data points across frames divided in time and a vector computing subunit 128 for execution vector operations. The various subunits may be hardwired, programmed, in any way, including temporary, permanent, in a volatile manner, in a non-volatile manner, other, or combination thereof. The central processing unit 112 may be further arranged to execute the method according to the first aspect and methods as discussed below, with or without any option indicated and any variation thereof. In any chosen realisation, the unit eventually provides a number of positions on the image, where such positions change during time. The positions are identified by their coordinate values and represent places of interest therein.
The storage module 114 is arranged to have image data stored thereon that may be acquired by means of the medical imaging system 130, as directly obtained as well as in a processed way. Furthermore, the storage module 114 may have data stored thereon comprising code executable by the central processing unit 112, enabling the central processing unit to operate as discussed above.
The procedure starts in a terminator 202 and proceeds to step 204 in which a sequence of images is obtained from the heart 190. The images are preferably obtained with approximately constant intervals between image acquisition moments. The images may be provided as a stream with frames. The images are preferably provided such that in the two-dimensional area of a cardiac cavity under scrutiny is as large as possible. For example, if operation of the left ventricle is to be assessed, the two-dimensional images in the sequence are provided under such angle that the line from base to apex is as long as possible, preferably at the end of the diastolic phase—or any other moment within the heartbeat cycle.
Alternatively or additionally, the area of the left ventricle on images is as large as possible at the end of the diastolic phase. In particular if two-dimensional images are provided based upon a three-dimensional model, this option is feasible. Otherwise, the medical imaging system 130 and the subject under scrutiny are properly positioned relative to one another. In other examples, the various steps discussed below are executed on a three-dimensional model. In the latter case, the images received are three-dimensional images.
In step 206, a wall contour of the heart 190 is identified in images of the stream or other sequence of image data. The wall contour may be detected by means edge detection subunit 122 or another detection subunit, including receiving manual input indicating a contour, that is arranged to detect the applicable wall.
In step 208, the wall thus identified is segmented, which may be done by means of the segmenting subunit 124. The wall may be segmented in accordance with a pre-defined scheme, based on established defined areas of the endocardial wall or the myocardium. Alternatively or additionally, the segments are defined based on size of each segment. In such case, the total number of segments defined may depend on the size of the heart.
In step 210, wall points 312 are defined on the wall contour 310 as vector points. The wall points 312 may be placed on pre-determined locations relative to known segments like, for example, segments of the myocardium or endocardial wall having been specifically defined in the science of cardiology. Alternatively or additionally, the wall points 312 are defined at pre-determined distances relative to one another. Subsequently, if segments have been defined, the wall points 312 are assigned to the segment in which they are defined.
In step 212, a motion reference 320 is defined within the schematic representation 300. The motion reference 320 may be a straight line or a curved line—a curve —, a single point, a finite set of points, other or a combination thereof. Within this context, a curve or a line may be interpreted as a set of points. Points of a set of points, either finite or infinite, may be provided on a curve or a line.
Preferably, the motion reference 320 is provided relative to the apex and/or base of the heart 190 as provided by the schematic representation 300. More in particular, the motion reference 320 is provided between one third and two third of the distance between base and apex, more in particular between half and two third of the distance between base and apex, for example based on the distance between apex and base at the end of the diastolic phase. The motion reference may coincide with a physiological direction, like the line between apex and base or a line perpendicular thereto.
In
Optionally, if the motion reference 320 is not a single point but as a finite set of reference points 324 or a reference curve 322 or reference line 322, particular wall points 312 may be assigned to a particular reference point 324 or to a particular point on the reference curve 322 or reference line 322. In step 214, it is checked whether all applicable images or frames received have been assessed. It may be noted that not every image received has to be assessed—also every second image may be assessed, for example. If not all images have been assessed, the procedure jumps back to step 206.
If all images have been processed as described above, the procedure continues to step 218 in which images may be grouped on per heartbeat basis, by determining cardiac cycle boundaries. The cardiac cycle boundaries may be determined based on size of the left ventricle, based on electrical signals (electrocardiogram), other, or a combination thereof.
In step 220, motion of wall points is obtained over time, for example by means of tracking, the motion detection subunit 126. Such may be done on a per cycle basis, if image data is grouped by cycle. If not, motion of each wall point is tracked over a particular interval in time—over the full acquired stream or at least part thereof. If the motion is tracked on a per-cycle basis, step 222 checks whether all cycles have been tracked. If motion is not per cycle, step 222 may be superfluous.
Based on the motion of the points, the vector computing subunit 128 computes a motion vector per wall point 312 or per segment in step 224. The motion vector may indicate a measure of motion like the maximum displacement or the differential displacement between two phases of the cardiac cycle, in a particular direction, based on a particular reference. In a particular case, an average of measure of motion may be provided as a general measure for a segment, like an average of maxima or an average of minima. The measure of displacement may in particular be determined for one heartbeat cycle. Alternatively, the motion vector may indicate a specific component of a specific vector component of the wall point motion parallel to a line from the wall point 312 to the motion reference 320, from the end of the systolic phase to the end of the diastolic phase. In the latter case, the motion vector or a reference motion vector component is always directed towards the motion reference or away from the motion reference.
The motion vector is thus determined based on a first location of the wall point 312 at the end of the diastolic phase, a first location of the wall point 312 at the end of the systolic phase and the location of the reference 320, optionally on location of reference points 324 provided thereon. In one example, each wall point 312 is associated with a specific, optionally unique, reference point 324. These reference points may be located on the reference curve 322. Furthermore, these reference points may be spaced apart.
Based on the determined motion vectors, a component of the motion vectors parallel to a line between the applicable wall point 312 and the motion reference 320 is determined in step 226. Based on the determined motion vectors determined in step 224, motion vectors per segment may be determined in step 228. Such segment motion vector may be the average or median of the wall point motion vectors for wall points 312 associated with the applicable segment. Either based on the segment motion vectors or based on the vector components determined in step 226 or based on both, a segment motion vector component of the segment motion vector parallel to a line between the applicable segment and the motion reference 320 is determined in step 230. In step 232 is checked whether all segments and/or all wall points 312 have been assessed by the vector computing subunit 128. If not, the process branches back to step 224.
If all segments and/or all wall points 312 have been assessed by the vector computing subunit 128, the overall movement of the heart 190 or of the wall contour 310 during a cycle may be determined in step 234. During a heartbeat cycle, the heart may move from left to right in the body and back. This movement may be detected based on the information provided by the movement of the wall points 312, for example by means of tracking of the wall points 312, and for example by determining common motion components over time of all or at least part of the determined motion vectors. This global movement may applied to obtain the motion of the motion reference 320 or, the other way round, the motion of the motion reference 320 can be directly detected from the images and taken as representing the overall motion. The overall motion as well as the motion of the motion reference 320—one or more curves, parts thereof, points thereon, other points or a combination thereof —, lines can be considered in their entire vector value, or they can be limited to a component of the motion vector parallel to a physiological direction, like the longitudinal or transversal with respect to the orientation of the heart. Based on the determined overall motion of the heart 190 during a heartbeat cycle, a compensation factor may be determined as adjustment data in step 236 to adjust motion of wall points 312 or segments relative to the motion reference 320 over a cycle for motion of the heart 190.
Using the adjustment data, all, some, one—or none—of the motion data values discussed above—vectors, vector components and motion paths—may be adjusted for motion of the heart 190 during a cycle in step 238, yielding adjusted motion data values.
In step 240, reference values for motion values are obtained. The reference values may relate to at least one of motion tracks, wall point motion vectors, wall point motion vector components relative to the motion reference, segment motion vectors, segment motion vector components relative to the motion reference, other motion values or a combination thereof.
In step 242, based on the reference motion values and motion values determined and calculated as discussed above, motion deviation values are calculated. The deviation values may be at least one of a ratio of the motion vector value and the reference motion value or the inverse thereof, a difference between the motion vector value and the reference motion value, an absolute value of the difference between the motion vector value and the reference motion value, a ratio of the difference between the motion vector value and the reference motion value on one hand and the reference motion value on the other hand or the inverse thereof, a ratio of the absolute value of the difference between the motion vector value and the reference motion value on one hand and the reference motion value on the other hand or the inverse thereof, other, or a combination thereof.
With the motion data processed and vector and vector component data calculated, as well as the deviation data having been calculated, at least part of the available data may be displayed on the electronic display screen 144. In step 244, the deviation data is displayed.
In the bullseye structure 510, the deviation data may be plotted by means of different colours, different colour intensities, different shadings, different hedgings, other, or a combination thereof as depicted by
In step 246, vector component data towards the motion reference 320 is displayed in the bullseye plot 500. As shown in
In the bullseye plot 500, the deviation data is displayed such that different parts of the sectors may provide different data. And the vector component data is displayed per sector of the bullseye. If multiple data values for the same entity are available for one sector, further processing of data may be required, for example taking a median or average of the values for all data points within a sector. In an alternative, the deviation data may be displayed on a per-sector basis. In an alternative, the vector data may be displayed with multiple data values per sector.
After at least part of the data thus discussed has been processed and displayed, the procedure ends in step 248.
The various aspects and embodiment relate to the following numbered examples:
| Number | Date | Country | Kind |
|---|---|---|---|
| 2028172 | May 2021 | NL | national |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/NL2022/050248 | 5/6/2022 | WO |
