Very, very, minor, but I think
```
x = 2.*np.pi/( np.product(image.array.shape) * image.scale )
```
is possibly a shade quicker than using a sqrt, so I'll make that edit unless you have a strong objection!