Method for 4D imaging of a heart of a living being, said method including at least the following steps: (a) an acquisition step wherein unfocused ultrasonic waves are transmitted in the heart by a 2D array ultrasonic probe and raw data from backscattered ultrasonic waves are acquired by said 2D array ultrasonic probe; (b) an imaging step wherein a sequence of 3D images is generated from said raw data, said sequence of 3D images forming an animation showing movements of an imaged volume including at least part of the heart; (c) a velocity computing step wherein 3D cartography of at least one parameter related to blood velocity and tissue velocity are automatically computed in said imaged volume and time, based on said sequence of 3D images; (d) a locating step wherein at least one point of interest having a predetermined property is automatically located in said sequence of 3D images, wherein said at least one point of interest is automatically located based solely on said 3D cartography and its temporal profile; (e) a quantification step wherein at least one velocity chosen between blood velocity and tissue velocity is automatically determined at said at least one point of interest and a predetermined quantification parameter is automatically computed, involving said at least one velocity, wherein said at least one velocity is automatically determined at said at least one point of interest based solely on said 3D cartography and its temporal profile; wherein the quantification parameter involves: - either a peak blood velocity in a certain anatomic area and said locating step (d) includes automatically locating said point of interest as a point of maximum blood velocity in said anatomic area and in at least part of the sequence of 3D images; - or a peak tissue velocity in a certain anatomic area and said locating step (d) includes automatically locating said anatomic area in the sequence of 3D images and automatically locating said point of interest as a point of maximum tissue velocity in said anatomic area in the sequence of 3D images. (a) an acquisition step wherein unfocused ultrasonic waves are transmitted in the heart by a 2D array ultrasonic probe and raw data from backscattered ultrasonic waves are acquired by said 2D array ultrasonic probe; (b) an imaging step wherein a sequence of 3D images is generated from said raw data, said sequence of 3D images forming an animation showing movements of an imaged volume including at least part of the heart; (c) a velocity computing step wherein 3D cartography of at least one parameter related to blood velocity and tissue velocity are automatically computed in said imaged volume and time, based on said sequence of 3D images; (d) a locating step wherein at least one point of interest having a predetermined property is automatically located in said sequence of 3D images, wherein said at least one point of interest is automatically located based solely on said 3D cartography and its temporal profile; (e) a quantification step wherein at least one velocity chosen between blood velocity and tissue velocity is automatically determined at said at least one point of interest and a predetermined quantification parameter is automatically computed, involving said at least one velocity, wherein said at least one velocity is automatically determined at said at least one point of interest based solely on said 3D cartography and its temporal profile; - either a peak blood velocity in a certain anatomic area and said locating step (d) includes automatically locating said point of interest as a point of maximum blood velocity in said anatomic area and in at least part of the sequence of 3D images; - or a peak tissue velocity in a certain anatomic area and said locating step (d) includes automatically locating said anatomic area in the sequence of 3D images and automatically locating said point of interest as a point of maximum tissue velocity in said anatomic area in the sequence of 3D images.