A new simulation method of field emission ion images is proposed. The method can be applied to any crystal shape, size and atomic structure (possibly with point defects and dislocations), without any trial-and-error adjustment of parameters. The algorithm is illustrated by three simple examples: spherical crystals of perfect FCC, BCC and HCP structures.