Invariant joint distribution of a stationary random field and its derivatives: Euler characteristic and critical point counts in 2 and 3D