The new method was applied to sodium airglow data obtained in 2011 at Syowa Station (69°S, 40Â☎), Antarctica. This method can deal with extensive amounts of imaging data obtained on different years and at various observation sites without bias caused by different event extraction criteria for the person processing the data. We have developed a new analysis method for obtaining the power spectrum in the horizontal phase velocity domain from airglow intensity image data to study atmospheric gravity waves. New statistical analysis of the horizontal phase velocity distribution of gravity waves observed by airglow imaging