Strain evaluation is of major interest in clinical cardiology as it can quantify the cardiac function. Myocardial elastography, a radio-frequency (RF)-based cross-correlation method, has been developed to evaluate the local strain distribution in the heart in vivo. However, inhomogeneities such as RF ablation lesions or infarction require a three-dimensional approach to be measured accurately. In addition, acquisitions at high volume rate are essential to evaluate the cardiac strain in three dimensions. Conventional focused transmit schemes using 2D matrix arrays, trade off sufficient volume rate for beam density or sector size to image rapid moving structure such as the heart, which lowers accuracy and precision in the strain estimation. In this study, we developed 3D myocardial elastography at high volume rates using diverging wave transmits to evaluate the local axial strain distribution in three dimensions in three open-chest canines before and after radio-frequency ablation. Acquisitions were performed with a 2.5 MHz 2D matrix array fully programmable used to emit 2000 diverging waves at 2000 volumes/s. Incremental displacements and strains enabled the visualization of rapid events during the QRS complex along with the different phases of the cardiac cycle in entire volumes. Cumulative displacement and strain volumes depict high contrast between non-ablated and ablated myocardium at the lesion location, mapping the tissue coagulation. 3D myocardial strain elastography could thus become an important technique to measure the regional strain distribution in three dimensions in humans.