We introduce an algorithm based on a method of snapshots for computing approximate balanced truncations for discrete-time, linear time-periodic systems. The algorithm is applicable to stable (including neutrally stable) systems with very high dimension, allowing either the output or the input to be very high dimensional, as well. This algorithm is illustrated with a numerical example.