I think Heppe's formula is a reasonable approach to approximate the maximum voltage gradient (V/m) in any direction where measurements are calculated one meter apart only in an even grid. It is not the voltage gradient in a diagonal direction necessarily, but the maximum voltage gradient that is needed. It is an approximation unless a potential is calculated one meter away from x in all directions.
In the case of concentric circular equipotential lines around a single rod, this gives a gradient that is too high, as GordS points out. In this case, Vstep=dv/dx=dv/dy.
If there were a single conductor of infinite length parallel to the x-axis, then you would have straight equipotential lines parallel to the x-axis. Va=Vx=Vb, and dv/dx=0 and sqrt[(dv/dx)²+(dv/dy)²]=dv/dy. In this case, the maximum voltage gradient is in the y direction. The voltage gradient in a diagonal direction would be less than dv/dy, and thus would not be Vstep.
In practical terms, I think if you calculate surface potentials every meter, then it is adequate to take Vstep at x to be equal to the maximum of |Vx-Va|, |Vx-Vb|, |Vx-Vc|, or |Vx-Vd|.