For the purpose of this discussion we define stretching of a single-dimension space as any transformation from one space to another that has the following property: at least an arbitrarily chosen sequence of two consecutive, equal in length, intervals in the input space is transformed into a sequence of two consecutive, not equal in length, intervals in the output space. Stretching an x-space to a y-space will be denoted as
|y = f(x)||(6)|
where is a positive real number whose value does not matter for the purpose of this discussion. As it can be seen in Fig. 2, if we keep the same sampling rate (), aliasing can occur when doing the reverse transformation, from x to y. In order to avoid aliasing, we need to compute , the largest accceptable sampling rate in the y domain. This can sometimes lead to a larger number of samples in the y domain, and thus to larger computational expense. This can be limited to some extent if the signal in the x-space has been bandpassed, as is often the case with seismic data, with the largest frequency present in the data () smaller than the Nyquist frequency given by the sampling rate (fNy). Thus, we can replace in our calculations with
which will result in a larger than that computed using , the sampling rate in the x space.
In order to compute , we will consider two points in the x space, as seen in Fig. 2, such as
The largest sampling rate in the y-space that will not result in aliasing is , the minimum possible value of . Suppose there is a value xm that minimizes . Then,
In particular, in the case of log-stretch, given by equation (1), if tm plays the role of xm from the equation above, then
because tm is the equivalent of xa from eq. (7) and Fig. 2. Thus, we get