1. A method of texture segmentation for a sonar image, S(x,y), comprising using a processor to perform the steps of: obtaining initial estimates of a first set of parameters of said sonar image; obtaining an estimate of an intensity autocorrelation function (ACF) of said image; initializing a second set of parameters, Î˜ i ={Î± i , [Î¼ X i Î¼ Y i ], Î£ i } of said sonar image based on said first set of parameters; initializing a likelihood function based on said first and second sets of parameters; maximizing said likelihood function to obtain updated parameters; calculating an updated intensity ACF based on said updated parameters; and iteratively returning to maximizing said likelihood function until convergence of said updated intensity ACF.