Use astropy for calculating pixel scale
Fix estimate_angle to use Y axis for estimating the angle

Needs test cases
I'm not too sure how well the approximation in estimate_angle works for
non-orthagonal axes, but it is identical to what was used before in pull
request #34

Also add a test for loading some imagecoords with WCS instead of Header