The problem of systematically determining groups of canonical transformations that leave invariant Hamilton's equations of motion for bound Keplerian system in the plane is formulated and examined in detail. Previously unknown invariance groups are found, and a procedure for finding further ones is outlined. The methods used are directly generalizable to other systems in spaces of arbitrary dimension.