A non-invasive method for estimating 2-D pressure gradients from ultrasound vector velocity data is presented. The method relies on in-plane vector velocity fields acquired using the Transverse Oscillation method. The pressure gradients are estimated by applying the Navier-Stokes equations for isotropic fluids to the estimated velocity fields. The velocity fields were measured for a steady flow on a carotid bifurcation phantom (Shelley Medical, Canada) with a 70% constriction on the internal branch. Scanning was performed with a BK8670 linear transducer (BK Medical, Denmark) connected to a BK Medical 2202 UltraView Pro Focus scanner. The results are validated through finite element simulations of the carotid flow model where the geometry is determined from MR images. This proof of concept study was conducted at nine ultrasound frames per second. Estimated pressure gradients along the longitudinal direction of the constriction varied from 0 kPa/m to 10 kPa/m with a normalized bias of -9.1% for the axial component and -7.9% for the lateral component. The relative standard deviation of the estimator, given in reference to the peak gradient, was 28.4% in the axial direction and 64.5% in the lateral direction. A study made across the constriction was also conducted. This yielded magnitudes from 0 kPa/m to 7 kPa/m with a normalized bias of -5.7% and 13.9% for the axial and lateral component, respectively. The relative standard deviations of this study were 45.2% and 83.2% in the axial and lateral direction, respectively.