Color Doppler imaging with line artifact reduction is provided in multi-beam scanning. Doppler estimates representing the same spatial location but calculated from spatially distinctive transmit beam groups are combined through weighted linear interpolation. Methods of calculating the linear interpolation weights are provided based on geometric relationships and optimization functions. Complete overlapping and superposition among receive beams in the interpolation region are not required. Partial interpolation among the receive beams, where only the estimates of the outer receive scan lines may overlap and be interpolated while estimates for scan lines closer to the transmit scan line are not interpolated, allowing for more rapid frame rate.
