Using a full vectorial finite element based beam propagation method combined with the imaginary axis propagation technique, the accurate calculation of the full vectorial leaky modes of a curved optical waveguide is addressed. The numerical precision of the proposed modal solution approach is supported by excellent agreement with published results obtained using other rigorous full vectorial approaches.