This paper proposes a method for producing smooth and positive estimates of the intensity function of an inhomogeneous Poisson process based on the shrinkage of wavelet coefficients of the observed counts. The information projection is used in conjunction with the level-dependent thresholds to yield smooth and positive estimates. This work is motivated by and demonstrated within the context of a problem involving gamma-ray burst data in astronomy. Simulation results are also presented in order to show the performance of the information projection estimators.