ABSTRACT The importance of Satellite Aperture Radar (SAR) imagery systems is increasing day-by-day in various field such as earth observation, hi-technology war mechanisms, etc. The images captured by SAR imagery systems are mainly used to detect and classify objects captured in the images. Due to the complexity of the image capturing process, SAR images can be highly noisy – often consisting of multiplicative noise, also known as speckle. To detect or classify objects in SAR images, speckle noise must be removed from images. During despeckling process, the preservation of important information in the SAR images while removing noise such as edge or patterns is a crucial task. Despeckling methods are liable to compromise on edge preservation ability while aspiring for good quality denoised images. Many researchers have proposed wavelet transform based despeckling approaches such as Discrete Wavelet Transform (DWT), Undecimated Wavelet Transform (UDWT), Dual Tree Complex Wavelet Transform (DTCWT). In these approaches, finding the best values of the coefficients plays an important role in yielding excellent denoised images with preserved edges. In this paper, we have proposed a novel optimization framework that optimizes thresholding coefficients of DTCWT despeckling method for SAR images. The proposed optimization framework is based on Fruit Fly Optimization (FOA) algorithm. The approach is a multi-objective optimization algorithm that is used to find maximum values for Peak Signal-to-Noise Ratio (PSNR), Mean Structural Similarity Index (MSSIM) and Equivalent number of look (ENL). The maximum value of MSSIM indicates high edge preservation capacity, whereas maximum PSNR value indicates good quality denoising of images. The maximum ENL value represents a good speckle-noise smoothing capability. We have applied our framework over some classical images as well as over SAR images of MSTAR dataset. In our experiments, we found that our proposed framework results in the excellent PSNR 36.87 dB, 35.4 dB and 37.8 dB, MSSIM values 0.92, 0.93 and 0.92, respectively, in the case of lena image, MSTAR dataset image and TerraSAR-X dataset image.