We propose a new superresolution imaging technique for interferometry using sparse modeling, utilizing two regularization terms: the 1 -norm and a new function named Total Squared Variation (TSV) of the brightness distribution. TSV is an edge-smoothing variant of Total Variation (TV), leading to reducing the sum of squared gradients. First, we demonstrate that our technique may achieve super-resolution of ∼ 30% compared to the traditional CLEAN beam size using synthetic observations of two point sources. Second, we present simulated observations of three physically motivated static models of Sgr A* with the Event Horizon Telescope (EHT) to show the performance of proposed techniques in greater detail. We find that 1 +TSV regularization outperforms 1 +TV regularization with the popular isotropic TV term and the Cotton-Schwab CLEAN algorithm, demonstrating that TSV is well-matched to the expected physical properties of the astronomical images, which are often nebulous. Remarkably, in both the image and gradient domains, the optimal beam size minimizing root-mean-squared errors is 10 % of the traditional CLEAN beam size for 1 +TSV regularization, and non-convolved reconstructed images have smaller errors than beam-convolved reconstructed images. This indicates that the traditional post-processing technique of Gaussian convolution in interferometric imaging may not be required for the 1 +TSV regularization. We also propose a feature extraction method to detect circular features from the image of a black hole shadow with the circle Hough transform (CHT) and use it to evaluate the performance of the image reconstruction. With our imaging technique and the CHT, the EHT can constrain the radius of the black hole shadow with an accuracy of ∼ 10 − 20 % in present simulations for Sgr A*, suggesting that the EHT would be able to provide useful independent measurements of the mass of the supermassive black holes in Sgr A* and also another primary target, M87.