The particle-in-cell code WARP has been enhanced to incorporate both two- and three-dimensional sheath extraction models giving WARP the capability of simulating entire ion beam transport systems including the extraction of beams from plasma sources. In this article, we describe a method of producing initial ion distributions for plasma extraction simulations in electron cyclotron resonance (ECR) ion sources based on experimentally measured sputtering on the source biased disk. Using this initialization method, we present preliminary results for extraction and transport simulations of an oxygen beam and compare them with experimental beam imaging on a quartz viewing plate for the superconducting ECR ion source VENUS.