Support transformations of coordinates given as numpy arrays about pyproj HOT 6 CLOSED

>>> from shapely.geometry import LineString, Point
>>> import numpy
>>> ll = LineString([Point(ii, ii) for ii in numpy.arange(0, 90, 0.01)])
>>> points = numpy.asarray(ll)
>>> points
array([[0.000e+00, 0.000e+00],
       [1.000e-02, 1.000e-02],
       [2.000e-02, 2.000e-02],
       ...,
       [8.997e+01, 8.997e+01],
       [8.998e+01, 8.998e+01],
       [8.999e+01, 8.999e+01]])
>>> from pyproj import Transformer
>>> transformer = Transformer.from_crs("epsg:4326", "epsg:3857")
>>> points_3857 = numpy.empty(points.shape)
>>> points_3857[:, 0], points_3857[:, 1] = transformer.transform(points[:, 0], points[:, 1])
>>> points_3857
array([[0.00000000e+00, 0.00000000e+00],
       [1.11319491e+03, 1.11319491e+03],
       [2.22638982e+03, 2.22638986e+03],
       ...,
       [1.00154146e+07, 5.26064404e+07],
       [1.00165278e+07, 5.51925525e+07],
       [1.00176410e+07, 5.96135402e+07]])