-This paper presents a novel method for power system harmonic estimation based on the Park transform. The proposed method firstly extends the signal to a group of three-phase signals in a-b-c coordinate. Then, a linear fitting based method is adopted to estimate the fundamental frequency. Afterwards, the Park transform is utilized to convert the three-phase signals from a-b-c coordinate to d-q-0 coordinate. Finally, the amplitude and phase of a harmonic component of interest can be calculated using the d-axis and q-axis components obtained. Simulation studies have been conducted using matrix laboratory (MATLAB) and power system computer aided design/electromagnetic transients including direct current (PSCAD/EMTDC). Simulation studies in MATLAB have considered three scenarios, i.e., no-frequency-deviation scenario, frequency-deviation scenario and the scenario in the presence of inter-harminics. The results have demonstrated that the proposed method achieves very high accuracy in frequency, phase and amplitude estimation under noisy conditions, and suffers little influence of the inter-harmonics. Moreover, comparison studies have proved that the proposed method is superior to FFT and Interpolated FFT with the Hanning Window (IpFFTHW). Finally, a popular case in PSCAD/EMTDC has been employed to further verify the effectiveness of the proposed method.