This paper presents a new wavelet based retrieval approach based on Spherically Invariant Random Vector (SIRV) modeling of wavelet subbands. Under this multivariate model, wavelet coefficients are considered as a realization of a random vector which is a product of the square root of a scalar random variable (called multiplier) with an independent Gaussian vector. We propose to work on the joint distribution of the scalar multiplier and the multivariate Gaussian process. For measuring similarity between two texture images, the geodesic distance is provided for various multiplier priors. A comparative study between the proposed method and conventional models on the VisTex image database is conducted and indicates that SIRV modeling combined with geodesic distance achieves higher recognition rates than classical approaches.