Use consistent ordering of x and y when transforming coordinates

[charparr]

For our collective sanity, consider using always_xy=True everywhere we call a pyproj transformer.

https://pyproj4.github.io/pyproj/stable/api/transformer.html

always_xy ([bool](https://docs.python.org/3/library/functions.html#bool), default=False) – If true, the transform method will accept as input and return as output coordinates using the traditional GIS order, that is longitude, latitude for geographic CRS and easting, northing for most projected CRS.

To implement, we'd just run a find and replace on those pyproj API calls and watch to see what breaks.