Curvelet transform is a multidirectional multiscale transform that enables sparse representations for signals. Curvelet-based feature extraction for Synthetic Aperture Radar (SAR) naturally enables utilizing spatial locality; the use of curvelet-based feature extraction is a novel method for SAR clustering. The implemented method is based on curvelet subband Gaussian distribution parameter estimation and cascading these estimated values. The implemented method is compared against original data, polarimetric decomposition features and speckle noise reduced data with use of k-means, fuzzy c-means, spatial fuzzy c-means and self-organizing maps clustering methods. Experimental results show that the curvelet subband Gaussian distribution parameter estimation method with use of self-organizing maps has the best results among other feature extraction-clustering performances, with up to 94.94% overall clustering accuracies. The results also suggest that the implemented method is robust against speckle noise.