This paper investigates the mean square error optimal estimation of scale invariant Wigner spectrum for the class of Gaussian locally self-similar processes, by the multitaper method. In this method, the spectrum is estimated as a weighted sum of scale invariant windowed spectrograms. Moreover, it is shown that the optimal multitapers are approximated by the quasi Lamperti transformation of Hermite functions, which is computationally more efficient. Finally, the performance and accuracy of the estimation is studied via simulation. keywords:Locally self-similar processes scale invariant Wigner spectrum multitaper method Hermite functions time-frequency analysis. 60G18 60G99