This work presents aspects of volume visualization of seismic horizons in 3-D seismic data. We consider both the direct and indirect approaches of volume visualization. In the direct approach we investigate the problem of selecting horizons using transfer functions. We present the 2-D opacity technique, which seeks to increase the ability to select horizons for visualization. We compare the use of instantaneous phase, adjusted phase and unwrapped phase as the second dimension, while seismic amplitude is the first dimension. Also in the direct approach, we show that the seismic amplitude gradient is not a good approximation for the normal vectors in seismic horizons. We suggest the gradient of instantaneous phase as a solution to this problem. In the indirect volume visualization approach we introduce a new optimization model to overcome the seismic horizon tracking problem. We present a heuristic method based on a greedy strategy to find solutions that are good approximations of the horizon of interest, even for complex geological structures.