A method to expedite data acquisition for multiple spatial-temporal analyses of tissue perfusion by contrast-enhanced ultrasound

For semiquantitative analyses of tissue perfusion using contrast-enhanced ultrasound the acquisition and processing of time intensity curves (TIC) is required. These TICs can be computed for each pixel of an image plane, yielding parametric images of classification numbers like "blood volume" and "flow rate." The expenditure of time for data acquisition and analysis typically limits semiquantitative perfusion imaging to a single image plane in 2-D. 3-D techniques, however, provide a higher diagnostic value since more information (e.g., of an entire lesion) is obtained. Moreover, spatial compounding, being a 2-D-technique where an object is imaged from different viewing angles, is known to improve image quality by reducing artifacts and speckle noise. Both techniques, 3-D and compounding, call for optimized acquisition and processing of TICs in several image planes (3-D) or in several (overlapping) sections of the same image plane (compounding) to decrease the time needed for data acquisition. Here, an approach of interleaved imaging is presented which is applicable, among others, to contrast perfusion imaging using the replenishment method. The total acquisition time is decreased by sequentially scanning image planes twice for short time spans-first, immediately after microbubble destruction to record the initial rise of the TICs, and second, a sufficient time thereafter to assess final values of the TIC. Data from both periods are combined to fit a model function from which parameters are extracted such as perfusion rate and blood volume. This approach was evaluated by in vitro measurements on a perfusion-mimicking phantom for both, individual images such as would be used for volume reconstruction in 3-D and compound images obtained from full angle spatial compounding (FASC, 360deg). An error analysis is conducted to derive the deviation of the extracted parameters of the proposed method compared with the conventional one. These deviations are entailed by a reduct- on in acquisition time of the proposed method, which can be adjusted by several parameters, depending on the prevailing flow. Optimization strategies are proposed to find optimal values for those settings.