Method and system for automated center positioning for radial velocity mapping in vector contrast imaging - Patents.com

Abstract

A method 100 for generating a radial velocity map for a target region, the method comprising: (i) receiving 120 a series of ultrasound images of the target region; (ii) generating 130 a grid including a plurality of points representing the target region, each representing a possible radial center in the target region; (iii) generating 140 a radial flux magnitude map for the grid including radial contrast flux for each of the plurality of points, the radial flux magnitude map generating step 140 comprising the following steps repeated until a radial contrast flux for each of the plurality of points has been calculated: selecting 142 a point of the plurality of points as a radial center; calculating 144 a radial velocity field for a first region around the selected point; and determining 146 a radial contrast flux for the first region; and (iv) automatically selecting 150 a radial center for the series of contrast images.

Description

[Technical field]

[0001] The present disclosure is generally directed to methods and systems for generating radial velocity maps in vector contrast imaging. [Background technology]

[0002]

[0002] Contrast-enhanced ultrasound (CEUS) is the ultrasound imaging of tissue of interest following intravenous injection of a contrast agent containing gaseous micro- and / or nanobubbles. CEUS images are collected as microbubbles of contrast agent wash-in and wash-out into and out of a lesion or organ over the duration of the contrast bolus. The dynamic wash-in and wash-out patterns exhibited by the lesion or organ of interest are then utilized to characterize the lesion or disease.

[0003]

[0003] Vector contrast imaging (VCI) has recently emerged as a breakthrough technique for visualizing and quantifying the motion and velocity of microbubbles as vectors (with both magnitude and direction) for the purpose of more confident analysis and diagnosis. There are various features available for mapping and visualization of bubbles using VCI, including velocity magnitude, velocity direction, and radial velocity relative to a given radial center. For easy visualization of vascularity, a temporal accumulation display is often utilized, as the location and / or velocity of contrast microbubbles perfusing a lesion or organ on each frame are accumulated over all selected frames. Summary of the Invention [Problem to be solved by the invention]

[0004]

[0004] However, although VCI can provide a map of the velocity and / or direction of microbubbles, this information is difficult to interpret, especially in the context of complex filling patterns with one or more flow centers. In such a situation, it is physically meaningful and hemodynamically (i.e., ultrasound contrast dynamics) valid to express velocity vectors in the form of radial velocities relative to a given radial center associated with a source or drain of local blood flow during a selected CEUS imaging sequence. Furthermore, for the analysis of VCI, the direction center should be determined before calculating the spatial distribution and temporal accumulation of radial velocities. Selecting the proper direction center can have a significant impact on the interpretation and understanding of VCI, and therefore on the patient's condition and diagnosis. However, selecting the direction center can be quite subjective and inaccurate. [Means for solving the problem]

[0005]

[0005] Therefore, there remains a need for a method and system for automatically and objectively determining a directional center for vector contrast imaging during contrast-enhanced ultrasound. Various embodiments and implementations herein are directed to a method and system for generating a radial velocity map, including automatically determining a directional center for a region being vector contrast imaged during contrast-enhanced ultrasound. The system receives a series of contrast images of a target region of a patient visualized using contrast-enhanced ultrasound. The system then generates a grid that represents some or all of the target region and includes a plurality of points, each of which represents a possible radial center in the target region. The system generates a radial flux magnitude map for the grid, including a radial contrast flux for each of the plurality of points. The radial flux magnitude map is generated by the following steps, which are repeated until the radial contrast flux for each of the plurality of points is calculated: (i) selecting one point of the plurality of points as a radial center; (ii) calculating the radial velocity field for each of the series of contrast images for a first region (e.g., a limited circular region with a predetermined radius) around the selected point; and (iii) determining the radial contrast flux for the first region for a single cumulative image generated from some or all of the series of contrast images. The system then uses the generated radial flux magnitude map to dynamically select a radial center for the series of contrast images, the radial center comprising one point of the plurality of points that has the maximum radial contrast flux relative to the remaining points of the plurality of points. The system then provides a visualization of the selected radial center and / or the generated radial flux magnitude map to the user via a user interface.

[0006] In general, in one aspect, a method is provided for generating a radial velocity map for a target region, the method including the steps of: (i) receiving a series of contrast images of the target region visualized using contrast enhanced ultrasound; (ii) generating a grid representing some or all of the target region, the grid including a plurality of points each representing a possible radial center in the target region; and (iii) generating a radial flux magnitude map for the grid including a radial contrast flux for each of the plurality of points, the step of generating the radial flux magnitude map being repeated until a radial contrast flux for each of the plurality of points has been calculated, the steps being: (1) selecting a point of the plurality of points as a radial center; (2) generating a radial flux magnitude map comprising the steps of: (1) calculating a radial velocity field for a first region around the selected point for each of the series of contrast images; and (3) determining a radial contrast flux for the first region for a single cumulative image generated from some or all of the series of contrast images; and (iv) automatically selecting a radial center for the series of contrast images using the generated radial flux magnitude map, wherein the radial center includes one point of the plurality of points that has the maximum radial contrast flux relative to the remaining points of the plurality of points.

[0007]

[0007] According to one embodiment, the step of determining the radial contrast flux includes determining a radial flux volume by summing the calculated radial velocity fields across all pixels in a first region for a single cumulative image generated from some or all of the series of contrast images, and dividing the determined radial flux volume by time and the area of ​​the first region.

[0008] According to one embodiment, the step of determining the radial flux volume (RFV) is performed according to the formula RFV=Σ(V r ΔS ΔT), where V r is the radial velocity field, ΔS is the area coverage of one pixel in a 2D contrast image, and ΔT is the time interval between two successive images in a sequence of contrast images. Note that ΔS is the volume of one voxel in a 3D sequence of contrast images.

[0009] According to one embodiment, the step of dividing the determined radial flux volume by time and the area of ​​the first region to generate the radial contrast flux (Flux) has the formula Flux=RFV / T / A, where T=ΔT·(N2-N1), where N1 is the start image in the series of contrast images, N2 is the end image in the series of contrast images, and A is the area of ​​the first region for 2D contrast images. Note that A is the volume of the first region for 3D contrast images.

[0010]

[0010] According to one embodiment, the multiple points are equidistantly spaced within the grid.

[0011]

[0011] According to one embodiment, the multiple points include less than all pixels in the target area.

[0012] According to one embodiment, the density of the points in the grid is determined automatically.

[0013] According to one embodiment, the method further comprises the step of providing the selected radial center and / or the generated radial flux magnitude map via a user interface.

[0014] According to a second aspect, there is provided a system for generating a radial velocity map for a target region, the system comprising a series of contrast images of the target region visualized using contrast enhanced ultrasound. The system also includes a processor for: (i) generating a grid including a plurality of points representing some or all of the target area, each representing a possible radial center in the target area; (ii) generating a radial flux magnitude map for the grid including a radial contrast flux for each of the plurality of points, where generating the radial flux magnitude map includes the following steps repeated until a radial contrast flux for each of the plurality of points has been calculated: (a) selecting a point of the plurality of points as a radial center; (b) calculating a radial velocity field for a first region around the selected point for each of the series of contrast images; and (c) determining the radial contrast flux for the first region for a single accumulated image generated from some or all of the series of contrast images; and (iii) selecting a radial center for the series of contrast images using the generated radial flux magnitude map, where the radial center includes a point of the plurality of points having a maximum radial contrast flux relative to the remaining points of the plurality of points. The system also includes a user interface for providing the selected radial center.

[0015] According to one embodiment, the user interface further provides the generated radial flux magnitude map.

[0016]

[0016] It should be appreciated that all combinations of the above concepts and additional concepts described in more detail below (provided that such concepts are not mutually inconsistent) are contemplated as being part of the inventive subject matter disclosed herein. In particular, all combinations of claimed subject matter appearing at the end of this disclosure are contemplated as being part of the inventive subject matter disclosed herein. It should also be appreciated that terms explicitly employed herein that may also appear in the disclosures incorporated by reference should be given the meaning most consistent with the specific concepts disclosed herein.

[0017]

[0017] These and other aspects of the various embodiments will be apparent from and elucidated with reference to the embodiments described hereinafter.

[0018]

[0018] In the drawings, similar reference characters generally refer to the same parts throughout different views. The figures illustrating the features and manner of implementing various embodiments should not be construed as being limited to other possible embodiments falling within the scope of the appended claims. Also, the drawings are not necessarily to scale, with emphasis instead generally being placed on illustrating the principles of various embodiments. [Brief description of the drawings]

[0019] [Figure 1] 1 is a flowchart of a method for generating a radial velocity map according to one embodiment. [Diagram 2]

[0020] FIG. 1 illustrates a schematic representation of a radial velocity mapping system, according to one embodiment. [Diagram 3]

[0021] FIG. 1 illustrates a schematic representation of a series of contrast-enhanced ultrasound images, each including a selected first point ("radial center") and a selected first region ("circular area") surrounding the selected first point, in accordance with one embodiment. [Figure 4A]

[0022] 1 is a radial flux map showing original values ​​according to one embodiment. [Figure 4B]

[0023] 1 is a radial flux map showing absolute values ​​of radial flux according to one embodiment; [Diagram 5]

[0024] 1 is a radial flux magnitude map illustrating locations of maximum flux values ​​in both the inward and outward radial directions, according to one embodiment. [Figure 6]

[0025] 1 is a radial flow velocity map whose radial center is automatically determined by the location of the maximum flux value in the radial flux magnitude map, according to one embodiment. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0020]

[0026] This disclosure describes various embodiments of systems and methods for automatically determining a directional center for a vector contrast imaged region in contrast enhanced ultrasound. More generally, applicants have recognized and appreciated that it would be beneficial to provide a method and system for significantly improving the analysis of vector contrast imaging. Accordingly, a radial velocity mapping system receives a series of contrast images of a target region of a patient visualized using contrast enhanced ultrasound. The system then generates a grid including a plurality of points representing a portion or all of the target region, each representing a possible radial center in the target region. The system generates a radial flux magnitude map for the grid including the radial contrast flux for each of the plurality of points. The radial flux magnitude map is generated by the following steps repeated until the radial contrast flux for each of the plurality of points is calculated: (i) selecting a point of the plurality of points as a radial center; (ii) calculating a radial velocity field for a first region around the selected point for each of the series of contrast images; and (iii) determining the radial contrast flux for the first region for a single accumulated image generated from a portion or all of the series of contrast images. The system then uses the generated radial flux magnitude map to automatically select a radial center for the series of contrast images, the radial center including one point of the plurality of points that has a maximum radial contrast flux relative to the remainder of the plurality of points. The system then provides a visualization of the selected radial center and / or the generated radial flux magnitude map to a user via a user interface.

[0021]

[0027] Thus, according to one embodiment, the method and system described or otherwise contemplated herein generate quantitative flow parameters for objective documentation of directional flow patterns. From the quantitative flow parameter map, the method and system identify a central location for outward flows as a source and a central location for inward flows as a sink. Furthermore, the method and system enable automatic locationing of directional centers for radial velocity imaging based on sources and / or drains in the directional flow pattern.

[0022]

[0028] According to one embodiment, the systems and methods described or otherwise contemplated herein may, in some non-limiting embodiments, be implemented as an element of a product for a contrast-enhanced ultrasound method or system.

[0023]

[0029] With reference to Figure 1, in one embodiment, a flow chart of a method 100 for generating a radial velocity map for a target region using a radial velocity mapping system is provided. The method described with respect to the figure is provided merely as an example and should be understood as not limiting the scope of the present disclosure. The radial velocity mapping system may be any of the systems described or otherwise contemplated herein. The radial velocity mapping system may be a single system or multiple different systems.

[0024]

[0030] In step 110 of the method, a radial velocity mapping system is provided. Referring to an embodiment of a radial velocity mapping system 200 shown in FIG. 2, for example, the system includes one or more of a processor 220, a memory 230, a user interface 240, a communication interface 250, and a storage 260, interconnected via one or more system buses 212. It will be understood that FIG. 2 is abstract in some respects and that the actual organization of the components of the system 200 may be different and more complex than that shown. Furthermore, the radial velocity mapping system 200 may be any of the systems described or otherwise contemplated herein. Other elements and components of the radial velocity mapping system 200 are disclosed and / or contemplated elsewhere herein.

[0025]

[0031] In step 120 of the method, a series of ultrasound images are received. The ultrasound images can be any images that can be utilized for the method. According to one embodiment, a series of ultrasound images is obtained over a period of time via contrast-enhanced ultrasound, where a region of interest in a patient is imaged after the introduction of a contrast agent containing gaseous micro- and / or nano-bubbles. The images are then collected as bubbles from the contrast agent enter and leave one or more components of the region of interest, such as a tumor or an organ, among other possible components. According to one embodiment, the images are such that vector contrast imaging (VCI) can be utilized to visualize and quantify the direction and velocity of the microbubbles. There are various features available for mapping and visualization of bubbles using VCI, including velocity magnitude, velocity direction, and radial velocity relative to a given radial center.

[0026]

[0032] According to one embodiment, the ultrasound device used to obtain the contrast-enhanced ultrasound image is a component of a radial velocity mapping system. However, according to another embodiment, the ultrasound device used to obtain the contrast-enhanced ultrasound image is not a component of the radial velocity mapping system, instead the system receives images from a local or remote ultrasound system. The local or remote ultrasound system, which may or may not be a component of the radial velocity mapping system, obtains images for immediate or future analysis. Thus, the contrast-enhanced ultrasound image is either immediately used or stored in a local or remote storage for use in further steps of the method, before or after image processing. In the case of immediate analysis, the system receives or acquires images from the ultrasound device. In the case of later analysis, the system receives or acquires images from the storage.

[0027]

[0033] In step 130 of the method, the radial velocity mapping system generates a grid that represents a portion or all of the visualized target region and includes a plurality of points in or on the grid. Each of these points is a possible radial center for the gridded target region. To perform vector contrast imaging (VCI) for the target region with one or more blood flow sources and / or drains, the direction and velocity of bubble movement are quantified or visualized with respect to the radial center. The radial center can be manually determined or pre-determined, but this introduces subjectivity in the VCI analysis. Using the methods and systems described or otherwise envisioned herein, the radial velocity mapping system provides an objective method for determining the radial center. According to one embodiment, the grid is generated in the memory of the system. The plurality of points are determined automatically by the system and may be randomly spaced or equidistantly spaced with respect to each other within the grid.

[0028]

[0034] In step 140 of the method, the radial velocity mapping system generates a radial flux magnitude map for the generated grid. The radial flux magnitude map includes radial contrast flux for some or all of the points in the grid. The radial flux magnitude map facilitates selection of the radial center by the system. The radial flux magnitude map can be generated using a variety of mechanisms, including methods described or otherwise contemplated herein.

[0029]

[0035] According to one embodiment, the following method is utilized to generate a radial flux magnitude map, by way of non-limiting example: The method is repeated until each or all of the plurality of points intended to be utilized contain the calculated radial contrast flux.

[0030]

[0036] In step 142, one point of the multiple points is selected as the radial center. This may be based on a random selection or may be based on some input from either the user or the system. For example, the user designates a temporary or permanent point in the grid to be the first radial center. As another example, the system selects the first point as the first radial center by selecting the same point each time or by randomly selecting a point. The next point selected as the second radial center may be based on the selection of the first point or may be randomly or user selected as well.

[0031]

[0037] In step 144, the system calculates the radial velocity field for a first region around the selected point for a series of acquired or received contrast images. According to one embodiment, the system calculates the radial velocity field for a first region around the selected point starting with a start image or frame N1 and ending with an end image or frame N2. nA sequence of acquired or received contrast images may be selected that includes some or all of the images ending with .

[0032]

[0038] According to one embodiment, the region around the first selected point for which the radial velocity field is calculated can be of any shape or size. The shape and / or size of the region can be based on the generated grid, the selected point, the target region, pre-determined or pre-programmed parameters, and / or any other factors. According to one embodiment, the region is a circular area having a first diameter and with the first selected point at the center of the circle.

[0033]

[0039] 3, in one embodiment, a series of contrast-enhanced ultrasound images is provided, each of which includes a selected first point (the "radial center") and a selected first region (the "circular area") surrounding the selected first point. Although shown as a circle, the first region can be any shape and size.

[0034]

[0040] The system performs the selected frames (N1...N n ) the radial velocity field V around the radial center for each R Determine the radial velocity field V R can be calculated using any method or mechanism for determining the radial velocity field. According to one embodiment, the radial velocity field is calculated using an image cross-correlation method (i.e., block matching) for image pixel displacement, or the Hungarian algorithm for microbubble pairing and tracking, although other methods are possible. According to one embodiment, the Hungarian algorithm is more accurate when tracking microbubbles across multiple frames, but is less efficient in computational cost than image cross-correlation based methods. Other methods are possible. The calculated radial velocity field V Rcan be immediately available or can be stored in local or remote memory for later use.

[0035]

[0041] In step 146, the system determines the radial contrast flux for a first region about the selected point for a single cumulative image generated from some or all of the series of contrast images. The radial contrast flux can be determined according to a variety of different methods.

[0036]

[0042] According to one embodiment for determining the radial contrast flux, the system generates a series of contrast images (N1...N n For a single accumulation image generated from some or all of the selected first regions, the calculated radial velocity fields are summed over all pixels within the selected first region, and the determined radial flux volume is divided by the time and the area of ​​the selected first region. For example, determining the radial flux volume (RFV) enclosed by the selected first region may be performed using the formula RFV=Σ(Vr ΔS ΔT) (Equation 1) where V r is the radial velocity field, ΔS is the area coverage of one pixel for a 2D image or the volume coverage of one voxel for a 3D image, and ΔT is the mean velocity field over a series of contrast images (N1...N n ) is the time interval between two consecutive images in

[0037]

[0043] According to one embodiment, dividing the determined radial flux volume by the time and the area of ​​the first region to generate a radial contrast flux (Flux) is performed according to the formula: Flux=RFV / T / A (Formula 2) where T = ΔT (N n −N1), where N1 is the starting image in the series of contrast images, and N nis the end image in the series of contrast images, and A is the area of ​​the first region for a 2D image, or the volume of the first region for a 3D image.

[0038]

[0044] Steps 142, 144, and 146 are repeated for each of the points in the generated grid, after which the method loops out at step 140. For example, at decision point 148 in FIG. 1, the system asks whether steps 142, 144, and 146 are to be repeated for an additional point of the points. If so, the system proceeds to step 142. If not, the system proceeds to step 140. Notably, however, this is only one example of a method for determining radial contrast flux. At the end of step 140, however, the system includes a radial flux magnitude map for the generated grid. The generated radial flux magnitude map may be utilized immediately or may be saved in local or remote storage for later use by the method.

[0039]

[0045] In step 150 of the method, the radial velocity mapping system utilizes the generated radial flux magnitude map to objectively select a radial center. According to one embodiment, the radial velocity mapping system automatically selects the radial center using the generated radial flux magnitude map. For example, the system may select a point within the plurality of points that includes a maximum radial contrast flux relative to the remaining points of the plurality of points.

[0040]

[0046] 4A and 4B, examples of radial flux magnitude maps for 2D contrast images are provided, generated using an image cross-correlation method that can be utilized to select the radial center, for example. FIG. 4A is a radial flux map showing the original values ​​(and thus positive and negative flux) and FIG. 4B is a radial flux map showing the absolute values ​​of the radial flux (mL / s / cm 2 ) in a non-limiting example. For this non-limiting example, the map was constructed over a 2D grid spaced every 10 pixels (0.8 mm) in both dimensions of the ROI in a contrast image of an abnormal human thyroid gland. The raw and absolute values ​​of radial flux were calculated over 200 frames (after initial contrast appearance) over a circular area with a radius of 25 pixels (2 mm). Other parameters are possible.

[0041]

[0047] According to one embodiment, a spatial distribution map of flow parameters (such as radial contrast flux and its magnitude in Fig. 4A and Fig. 4B) can be constructed to quantify the directional flow pattern. For the calculation of radial flux, it is important to select a sufficient diameter of the radial area (as the selected first region) to maintain both low variation in the calculation of flux (lower for larger circle areas) and high sensitivity of spatial distribution (higher for smaller circle areas). It is also important to properly select the timing and length of the image sequence. These parameters can be predetermined or experimentally derived.

[0042]

[0048] According to one embodiment, the radial center in the radial velocity map represents a focal point in a complex flow pattern, and therefore the spatial distribution map of quantitative flow parameters can be utilized for automatic selection of the radial center. As demonstrated in Fig. 5 and Fig. 6 for a two-dimensional (2D) VCI based on image cross-correlation method, the location of the maximum flux magnitude value in the radial flux magnitude map (Fig. 5) can be directly used as the radial center for the radial flow velocity map (Fig. 6). This is an example of automatic selection of the radial center using a quantitative flow map. Thus, Fig. 5 includes a radial flux magnitude map showing the location of the maximum flux value in both the inward radial direction and the outward radial direction. Fig. 6 includes a radial flow velocity map whose radial center is automatically determined by the location of the maximum flux value in the radial flux magnitude map.

[0043]

[0049] According to one non-limiting embodiment, the user can select the grid size and / or area per grid point, and the system automatically calculates the map and automatically locates the highest inflow / outflow points and displays their values. A threshold can be set to limit the algorithm to only the highest points and thus the maximum absolute values. According to one embodiment, the flux can be tied to the divergence theorem (Gauss's law). Thus, the system alternatively calculates the divergence of the velocity vector field and uses this to calculate the flux directly.

[0044]

[0050] According to one embodiment, the radial flux can also be separated into "outward flux", for centrifugal flows only (moving away from the center or axis), and "inward flux", for centripetal flows only (moving towards the center or axis). Outward flux can be calculated using only positive radial velocities (VR>0), and inward flux can be calculated using only negative radial velocities (VR<0). The radial flux can thus be the sum of outward and inward radial fluxes. According to one embodiment, the inward and outward radial fluxes can be particularly useful for Hungarian algorithm-based VCI maps, since closely adjacent opposing flows can be imaged separately via tracking of individual microbubbles across multiple frames on such VCI maps. The radial centers in the inward, outward and normal radial velocity maps can be determined from the location of the maximum absolute value in the inward flux map (using only VR<0), the outward flux map (using only VR>0) and the normal flux map (i.e. the map of the sum of the inward and outward flux), respectively. Many other methods are possible.

[0045]

[0051] In step 160 of the method, the system provides a report including, among other possible information, one or more of the selected radial center and the generated radial flux magnitude map. The system may provide the report to a clinician or other user via a user interface of the system. The display may also include information about the patient, contrast-enhanced ultrasound, and / or other information. Any information may be communicated to another device by wired and / or wireless communication. For example, the system communicates information to a mobile phone, a computer, a laptop, a wearable device, and / or other device configured to enable display and / or other communication of the report. The user interface may be any device or system that allows for communicating and / or receiving information and includes a display, a mouse, and / or a keyboard for receiving user commands.

[0046]

[0052] In an optional step 170 of the method, the information generated by step 160 is provided electronically by the system to a clinician or other decision maker who then utilizes the displayed graphics with the selected radial center and / or the generated radial flux magnitude map for patient care decision making. For example, the clinician or other decision maker can use the displayed graphics, such as a radial flow velocity map with an automatically selected radial center, to interpret blood flow in a lesion, organ, tumor, or other target area. The interpretation of blood flow is, for example, a diagnosis, and thus, this diagnosis can be used to select a treatment or other care option for the patient. Many other implementations are possible, such as automatically providing the information into a clinical decision support (CDS) system for later use.

[0047]

[0053] Referring to Figure 2, there is shown a schematic representation of a radial velocity mapping system 200. System 200 may be any of the systems described or otherwise contemplated herein and may include any of the components described or otherwise contemplated herein. It will be understood that Figure 2 is abstract in some respects and that the actual organization of the components of system 200 may be different and more complex than that shown.

[0048]

[0054] According to one embodiment, the system 200 includes a processor 220 capable of executing instructions or otherwise processing data stored in the memory 230 or storage 260, for example, to perform one or more steps of the method. For example, a computer-readable storage medium embodied in a non-transitory memory includes instructions that, when executed by the computer or processor 220, cause the computer or processor 220 to perform the steps of the method. The processor 220 may be in the form of one or more modules. The processor 220 may take any suitable form, including, but not limited to, a microprocessor, a microcontroller, multiple microcontrollers, a circuit, a field programmable gate array (FPGA), an application specific integrated circuit (ASIC), a single processor, or multiple processors.

[0049]

[0055] Memory 230 may take any suitable form, including non-volatile memory and / or RAM. Memory 230 includes various memories, such as, for example, L1, L2, or L3 cache, or system memory. Thus, memory 230 includes static random access memory (SRAM), dynamic RAM (DRAM), flash memory, read-only memory (ROM), or other similar memory devices. The memory may store, among other things, an operating system. RAM is used by the processor for temporary storage of data. According to one embodiment, the operating system includes code that, when executed by the processor, controls the operation of one or more components of system 200. It will be apparent that in embodiments in which the processor performs one or more of the functions described herein in hardware, software described as corresponding to such functions in other embodiments may be omitted.

[0050]

[0056] User interface 240 includes one or more devices for enabling communication with a user. The user interface may be any device or system that enables communication and / or reception of information and includes a display, mouse, and / or keyboard for receiving user commands. In some embodiments, user interface 240 includes a command line interface or a graphical user interface that is presented to a remote terminal via communication interface 250. The user interface may be located with one or more other components of the system or may be located remotely from the system and in communication via a wired and / or wireless communication network.

[0051]

[0057] The communication interface 250 includes one or more devices for enabling communication with other hardware devices. For example, the communication interface 250 includes a network interface card (NIC) configured to communicate according to an Ethernet protocol. Additionally, the communication interface 250 implements a TCP / IP stack for communication according to a TCP / IP protocol. Various alternative or additional hardware or configurations for the communication interface 250 will be apparent.

[0052]

[0058] Storage 260 includes one or more machine-readable storage media, such as read-only memory (ROM), random access memory (RAM), magnetic disk storage media, optical storage media, flash memory devices, or similar storage media. In various embodiments, storage 260 stores instructions for execution by processor 220 or data on which processor 220 operates. For example, storage 260 stores operating system 261 for controlling various operations of system 200.

[0053]

[0059] It will be apparent that the various information described as being stored in storage 260 may additionally or alternatively be stored in memory 230. In this regard, memory 230 may also be considered to constitute a storage device, and storage 260 may be considered a memory. Various other configurations may be apparent. Moreover, both memory 230 and storage 260 may be considered to be non-transitory machine-readable media. As used herein, the term non-transitory may be understood to exclude transitory signals, but to include all forms of storage, including both volatile and non-volatile memory.

[0054]

[0060] Although system 200 is shown as including one of each described component, various components may be repeated in various embodiments. For example, processor 220 includes multiple microprocessors configured to execute the methods described herein alone or to execute steps or subroutines of the methods described herein such that the multiple processors cooperate to achieve the functionality described herein. Furthermore, when one or more components of system 200 are implemented in a cloud computing system, various hardware components may belong to separate physical systems. For example, processor 220 includes a first processor in a first server and a second processor in a second server. Many other variations and configurations are possible.

[0055]

[0061] According to one embodiment, the electronic medical record system 270 is an electronic medical record database from which information about the patient, including ultrasound images obtained from / for the patient, is obtained or received. The electronic medical record database may be a local or remote database and is in direct and / or indirect communication with the radial velocity mapping system 200. Thus, according to one embodiment, the radial velocity mapping system comprises an electronic medical record database or system 270.

[0056]

[0062] According to one embodiment, the system includes an ultrasound device 280. The ultrasound device may be any ultrasound device capable of capturing contrast-enhanced images as described or otherwise contemplated herein. According to another embodiment, rather than including an ultrasound device, the system is simply in wireless and / or wired communication with a local or remote ultrasound device from which one or more ultrasound images are obtained or otherwise received.

[0057]

[0063] According to one embodiment, storage 260 of system 200 stores one or more algorithms, modules, and / or instructions for performing one or more functions or steps of methods described or otherwise contemplated herein. For example, the system includes radial flux magnitude map instructions 262, radial center selection instructions 263, and / or reporting instructions 264, among other instructions or data.

[0058]

[0064] According to one embodiment, the radial flux magnitude map command 262 instructs the system to generate a radial flux magnitude map. According to one embodiment, the radial flux magnitude map includes radial contrast flux for some or all of the points in the grid. The radial flux magnitude map can be generated using various mechanisms. According to one embodiment, the system selects one point of the points in the grid as a radial center. The system then calculates a radial velocity field for a first region around the selected point for a plurality of series of acquired or received contrast images. The system then determines the radial contrast flux for the first region around the selected point for a single accumulated image generated from some or all of the series of contrast images.

[0059]

[0065] According to one embodiment, a radial center selection command 263 instructs the system to automatically select a radial center using the generated radial flux magnitude map. According to one embodiment, the system can select a point within the plurality of points that includes a maximum radial contrast flux relative to the remaining points of the plurality of points. Methods are possible.

[0060]

[0066] According to one embodiment, the reporting instructions 264 instruct the system to generate and provide information to a user via a user interface, including one or more of the selected radial center, the generated radial flux magnitude map, and the radial flow velocity map with the automatically selected radial center. The display may also include information about the patient, contrast enhanced ultrasound, and / or any other information. Any information may be communicated to another device by wired and / or wireless communication. For example, the system communicates information to a mobile phone, a computer, a laptop, a wearable device, and / or other device configured to enable display and / or other communication of the report. The user interface may be any device or system that allows for communicating and / or receiving information and includes a display, a mouse, and / or a keyboard for receiving user commands.

[0061]

[0067] According to one embodiment, the radial velocity mapping system is configured to process thousands or millions of data points to generate a radial flow velocity map with an automatically selected radial center. For example, generating a radial flux magnitude map for a plurality of points in a generated grid for an image involves millions or billions of calculations. Each generated radial flux magnitude map is a new map that did not exist before and includes millions of data points and millions or billions of calculations. Automatically selecting a radial center and then generating and displaying a radial flow velocity map with an automatically selected radial center further includes millions of additional calculations. Thus, generating and providing a radial flow velocity map with an automatically selected radial center involves a process that involves a large amount of calculations and analysis that the human brain cannot accomplish in one or more lifetimes. By providing improved analysis of contrast-enhanced ultrasound images, the novel radial velocity mapping system and method has a huge positive impact and effect on ultrasound images and analysis compared to prior art systems.

[0062]

[0068] All definitions defined and used herein should be understood to control over dictionary definitions, definitions in documents incorporated by reference, and / or ordinary meanings of the defined terms.

[0063]

[0069] As used in this specification and claims, "a" or "an" should be understood to mean "at least one," unless expressly stated to the contrary.

[0064]

[0070] As used in the specification and claims, the phrase "and / or" should be understood to mean "either or both" of the elements so conjoined, i.e., elements that are conjunctive in some cases and disjunctive in other cases. Multiple elements listed with "and / or" should be construed in the same manner, i.e., "one or more" of the elements so conjoined. Other elements other than the elements specifically identified by the "and / or" clause may optionally be present, whether related or unrelated to those elements specifically identified.

[0065]

[0071] As used herein and in the claims, "or" should be understood to have the same meaning as "and / or" as defined above. For example, when separating items in a list, "or" or "and / or" should be interpreted as being inclusive, i.e., including at least one of an element or list of elements, but also including two or more of them, and optionally including additional unlisted items. Only terms expressly stated to the contrary, such as "only one of" or "exactly one of," or, when used in the claims, "consisting of," refer to the inclusion of exactly one element of an element or list of elements. In general, as used herein, the term "or" should only be interpreted as indicating exclusive alternatives (i.e., "one or the other, but not both") when preceded by terms of exclusivity, such as "either," "one of," "only one of," or "exactly one of."

[0066]

[0072] As used in this specification and claims, the phrase "at least one" in reference to a list of one or more elements should be understood to mean at least one element selected from any one or more of the elements in the list of elements, but not necessarily including at least one of each and every element specifically listed in the list of elements, and not excluding any combination of elements in the list of elements. This definition also allows for the optional presence of elements other than those elements specifically identified in the list of elements to which the phrase "at least one" refers, whether related or unrelated to the specifically identified elements.

[0067]

[0073] Also, unless expressly stated to the contrary, in any method claimed herein that includes two or more steps or actions, it should be understood that the order of the method steps or actions is not necessarily limited to the order in which the method steps or actions are recited.

[0068]

[0074] In the claims, as well as in the specification above, all transitional phrases such as "comprising," "including," "carrying," "having," "including," "involving," "holding," "consisting of," etc., are to be understood to be open-ended, i.e., meaning "including, but not limited to." Only the transitional phrases "consisting of" and "consisting essentially of" shall be closed or semi-closed transitional phrases, respectively.

[0069]

[0075] While several inventive embodiments have been described and illustrated herein, those skilled in the art will readily envision various other means and / or structures for performing the functions described herein and / or obtaining one or more of the results and / or advantages described herein, and each such variation and / or modification is deemed to be within the scope of the inventive embodiments described herein. More generally, those skilled in the art will readily appreciate that all parameters, dimensions, materials, and configurations described herein are exemplary, and that the actual parameters, dimensions, materials, and / or configurations will depend on the particular application or applications for which the inventive teachings are used. Those skilled in the art will recognize, or be able to ascertain using no more than routine experimentation, many equivalents to the specific inventive embodiments described herein. It is therefore to be understood that the above-described embodiments are presented by way of example only, and that within the scope of the appended claims and their equivalents, the inventive embodiments may be practiced otherwise than as specifically described and claimed. The inventive embodiments of the present disclosure are directed to each individual feature, system, article, material, kit, and / or method described herein. Furthermore, any combination of two or more such features, systems, articles, materials, kits, and / or methods is included within the inventive scope of the present disclosure, if such features, systems, articles, materials, kits, and / or methods are not mutually inconsistent.

Claims

1. 1. A method for generating a radial velocity map for a target region, the method comprising: receiving a series of contrast images of the target area visualized using contrast-enhanced ultrasound; generating a grid representing part or all of the target area, the grid including a plurality of points each representing a possible radial center within the target area; generating a radial flux magnitude map for the grid including a radial contrast flux for each of the plurality of points, the step of generating the radial flux magnitude map comprising the following steps repeated until a radial contrast flux for each of the plurality of points has been calculated: selecting one point of the plurality of points as a radial center; calculating, for each of the series of contrast images, a radial velocity field for a first region around the selected point; determining a radial contrast flux for the first region for a single cumulative image generated from some or all of the series of contrast images; generating a radial flux magnitude map having: using the generated radial flux magnitude map to automatically select a radial center for the series of contrast images, the radial center comprising one point of the plurality of points having a maximum radial contrast flux relative to the remainder of the plurality of points; generating a radial velocity map having the automatically selected radial center; A method comprising:

2. determining the radial contrast flux comprises: determining a radial flux volume by summing the calculated radial velocity fields over all pixels within the first region for a single cumulative image generated from some or all of the series of contrast images; dividing the determined radial flux volume by a time and an area of ​​the first region; 2. The method of claim 1, comprising:

3. The step of determining the radial flux volume (RFV) is performed using the formula RFV=Σ(V) ｒ ・ΔS・ΔT) where V ｒ 3. The method of claim 2, wherein σ is the radial velocity field, ΔS is the area coverage of one pixel for a 2D image or the volume coverage of one voxel for a 3D image, and ΔT is the time interval between two consecutive images in the series of contrast images.

4. Dividing the determined radial flux volume by time and the area of ​​the first region to generate a radial contrast flux (Flux) is performed using the formula Flux=RFV / T / A where T = ΔT · (N ２ -N １ ), where N １ is the starting image in the series of contrast images, and N ２ 3. The method of claim 2, wherein ∑ is the end image in the series of contrast images, and A is the area of ​​the first region for a 2D image or the volume of the first region for a 3D image.

5. The method of claim 1 , wherein the plurality of points are equidistantly spaced within the grid.

6. The method of claim 1 , wherein the plurality of points includes fewer than all pixels in the target region.

7. The method of claim 1 , wherein the density of the plurality of points within the grid is determined automatically.

8. The method of claim 1 , further comprising providing the selected radial center and / or the generated radial flux magnitude map via a user interface.

9. 1. A system for generating a radial velocity map for a target region, the system comprising: a series of contrast images of the target area visualized using contrast-enhanced ultrasound; (i) generating a grid representing part or all of the target region, the grid including a plurality of points each representing a possible radial center within the target region; and (ii) generating a radial flux magnitude map for the grid including a radial contrast flux for each of the plurality of points, wherein generating the radial flux magnitude map comprises the following steps repeated until a radial contrast flux for each of the plurality of points has been calculated: (a) selecting one of the plurality of points as a radial center; and (b) calculating, for each of the series of contrast images, a radial velocity field for a first region around the selected point. (c) generating a radial flux magnitude map, the radial flux magnitude map comprising: (i) calculating a radial contrast flux for the first region for a single cumulative image generated from some or all of the series of contrast images; (ii) selecting a radial center for the series of contrast images using the generated radial flux magnitude map, the radial center comprising one point of the plurality of points having a maximum radial contrast flux relative to the remainder of the plurality of points; and (iv) generating a radial velocity map having the automatically selected radial center. a user interface for presenting the generated radial velocity map having the selected radial center and the automatically selected radial center; and A system comprising:

10. determining the radial contrast flux comprises: determining a radial flux volume by summing the calculated radial velocity fields over all pixels within the first region for a single cumulative image generated from some or all of the series of contrast images; dividing the determined radial flux volume by a time and an area of ​​the first region; The system of claim 9 , comprising:

11. The step of determining the radial flux volume (RFV) is performed using the formula RFV=Σ(V) ｒ ・ΔS・ΔT) where Vr is the radial velocity field, ΔS is the area coverage of one pixel for a 2D image or the volume coverage of one voxel for a 3D image, and ΔT is the time interval between two consecutive images in the series of contrast images.

12. The step of dividing the determined radial flux volume by time and the area of ​​the first region to generate a radial contrast flux (Flux) is performed using the formula Flux=RFV / T / A where T = ΔT · (N ２ -N １ ), where N １ is the starting image in the series of contrast images, and N ２ 11. The system of claim 10, wherein A is the area of ​​the first region for a 2D image or the volume of the first region for a 3D image. (i) the plurality of points are equidistantly spaced within the grid; or (ii) the plurality of points includes fewer than all pixels in the target region; The system of claim 9.

14. The system described in claim 9, wherein the user interface further provides the generated radial flux magnitude map.

15. A computer-readable storage medium embodied in non-transitory memory containing instructions that, when executed by a computer, cause the computer to perform the steps of the method of claim 1.

Images