Line artifact reduction is provided in multi-beam scanning for Doppler imaging. An analytic solution uses estimated velocities for collocated receive scan lines from different transmit beams to reduce or remove the biases caused by multi-beam receive. For each location, a line or curve is fit to the estimated velocities as a function of distance of the transmit scan lines to the location. The velocity at an intercept of this line with the distance of zero (intercept with the origin) indicates the unbiased velocity. This approach allows a solution even where all of the velocity estimates have a bias with the same sign, such as due to the location being on a same side of the different transmit scan lines.