Accurate arterial stiffness measurement would improve diagnosis and monitoring for many diseases. Atherosclerotic plaques and aneurysms are expected to involve focal changes in vessel wall properties; therefore, a method to image the stiffness variation would be a valuable clinical tool. The pulse wave inverse problem (PWIP) fits unknown parameters from a computational model of arterial pulse wave propagation to ultrasound-based measurements of vessel wall displacements by minimizing the difference between the model and measured displacements. The PWIP has been validated in phantoms, and this study presents the first in vivo demonstration. The common carotid arteries of five healthy volunteers were imaged five times in a single session with repositioning of the probe and subject between each scan. The 1D finite difference computational model used in the PWIP spanned from the start of the transducer to the carotid bifurcation, where a resistance outlet boundary condition was applied to approximately model the downstream reflection of the pulse wave. Unknown parameters that were estimated by the PWIP included a 10-segment linear piecewise compliance distribution and 16 discrete cosine transformation coefficients for each of the inlet boundary conditions. Input data was selected to include pulse waves resulting from the primary pulse and dicrotic notch. The recovered compliance maps indicate that the compliance increases close to the bifurcation, and the variability of the average pulse wave velocity estimated through the PWIP is on the order of 11%, which is similar to that of the conventional processing technique which tracks the wavefront arrival time (13%).