Understanding the properties of interstellar turbulence is a great intellectual challenge and the urge to solve this problem is partially motivated by a necessity to explain the star formation mystery. This review deals with a recently suggested inversion technique as applied to atomic hydrogen. This technique allows to determine 3D turbulence statistics through the variations of 21 cm intensity. We claim that a radio interferometer is an ideal tool for such a study as its visibility function is directly related to the statistics of galactic HI. Next, we show how galactic rotation curve can be used to study the turbulence slice by slice and relate the statistics given in galactic coordinates and in the velocity space. The application of the technique to HI data reveals a shallow spectrum of the underlying HI density that is not compatible with a naive Kolmogorov picture. We show that the random density corresponding to the found spectrum tends to form low contrast filaments that are elongated towards the observer.