In this study a semi-automated and observer-independent algorithm for quantifying post-stenotic signal loss (PSL) in 3D phase-contrast (PC) magnetic resonance angiography (MRA) of patients with renal artery stenosis is presented. This algorithm was developed on MRA datasets of stenotic phantoms, which were included in a flow circuit with stationary flows. The length and the severity of the PSL (incorporating both length and degree of PSL) in the maximum intensity projections (MIPs) of MRA datasets were proposed for quantifying stenoses. The algorithm was tested in renal arteries of ten patients with renal artery stenosis and seven healthy volunteers. Digital subtraction angiography (DSA) was performed in the patients and served as the gold standard. Stenosis severity showed better correlation with the severity of the PSL than with the length, both for in vitro as in vivo. Spearman correlation coefficients (rS) showed statistically significant correlations between the severity of the PSL and parameters determined by DSA, i.e. percent diameter stenosis (rS = 0.90). The length of the PSL showed no correlation with the diameter stenosis (rS = 0.37).