Beginning with a review of the existing literature on the computation of spectral distances on noncommutative spaces like Moyal plane and fuzzy sphere, adaptable to Hilbert–Schmidt operatorial formulation, we carry out a correction, revision and extension of the algorithm provided in [1] i.e. [F. G. Scholtz and B. Chakraborty, J. Phys. A, Math. Theor. 46 (2013) 085204] to compute the finite Connes’ distance between normal states. The revised expression, which we provide here, involves the computation of the infimum of an expression which involves the “transverse” [Formula: see text] component of the algebra element in addition to the “longitudinal” component [Formula: see text] of [1], identified with the difference of density matrices representing the states, whereas the expression given in [1] involves only [Formula: see text] and corresponds to the lower bound of the distance. This renders the revised formula less user-friendly, as the determination of the exact transverse component for which the infimum is reached remains a nontrivial task, but under rather generic conditions it turns out that the Connes’ distance is proportional to the Hilbert-Schmidt norm of [Formula: see text], leading to considerable simplification. In addition, we can determine an upper bound of the distance by emulating and adapting the approach of [P. Martinetti and L. Tomassini, Commun. Math. Phys. 323 (2013) 107–141]. We then look for an optimal element for which the upper bound is reached. We are able to find one for the Moyal plane through the limit of a sequence obtained by finite-dimensional projections of the representative of an element belonging to a multiplier algebra, onto the subspaces of the total Hilbert space, occurring in the spectral triple and spanned by the eigen-spinors of the respective Dirac operator. This is in contrast with the fuzzy sphere, where the upper bound, which is given by the geodesic of a commutative sphere, is never reached for any finite [Formula: see text]-representation of [Formula: see text]. Indeed, for the case of maximal noncommutativity ([Formula: see text]), the finite distance is shown to coincide exactly with the above-mentioned lower bound, with the transverse component playing no role. This, however, starts changing from [Formula: see text] onwards and we try to improve the estimate of the finite distance and provide an almost exact result, using our revised algorithm. The contrasting features of these types of noncommutative spaces becomes quite transparent through the analysis, carried out in the eigen-spinor bases of the respective Dirac operators.