.all-contributorsrc

 {
   "projectName": "pyproj",
-  "projectOwner": "snowman2",
+  "projectOwner": "pyproj4",
   "repoType": "github",
   "repoHost": "https://github.com",
   "files": [
@@ -308,6 +308,15 @@
       "contributions": [
         "doc"
       ]
+    },
+    {
+      "login": "edpop",
+      "name": "Eduard Popov",
+      "avatar_url": "https://avatars3.githubusercontent.com/u/13479292?v=4",
+      "profile": "https://github.com/edpop",
+      "contributions": [
+        "doc"
+      ]
     }
   ],
   "contributorsPerLine": 7

.gitignore

 pyproj/proj_dir/
 pyproj/*.c
+pyproj/*/*.c
+pyproj/*/*.html
 pyproj/*.html
 proj-*/
 

.travis.yml

@@ -11,7 +11,7 @@ env:
   global: 
     - PROJ_BASE_DIR=$HOME/proj_install
     - CYTHON_COVERAGE=True
-    - PROJSOURCE=6.2.1
+    - PROJSOURCE=6.3.1
     # Following generated with
     - WHEELHOUSE_UPLOADER_USERNAME=travis-worker
     # Following generated by
@@ -26,29 +26,33 @@ matrix:
       os: osx
       env:
         - PYTHON=3.6.8
-        - PROJSOURCE=6.2.1
     - python: 3.5
-      env:
-        - PROJSOURCE=6.2.1
+    - python: 3.6
     - python: 3.6
       env:
         - PROJSOURCE=6.2.0
     - python: 3.6
       env:
         - PROJSOURCE=6.2.1
+    - python: 3.6
+      env:
+        - PROJSOURCE=6.3.0
     - python: 3.7
       env:
-        - PROJSOURCE=6.2.1
         - DOC=true
     - python: 3.8
-      env:
-        - PROJSOURCE=6.2.1
-    - python: "nightly"
+    - python: 3.8
       env:
         - PROJSOURCE=git
+    # - python: "nightly"
+    #   env:
+    #     - PROJSOURCE=git
     
   allow_failures:
-    - python: "nightly"
+    # - python: "nightly"
+    #   env:
+    #     - PROJSOURCE=git
+    - python: 3.8
       env:
         - PROJSOURCE=git
 

MANIFEST.in

@@ -8,6 +8,7 @@ include pyproj/*.pyx
 include pyproj/*.pxd
 include pyproj/*.pxi
 include test/sample.out
+include test/conftest.py
 recursive-include docs *
 prune docs/_build
 prune pyproj/proj_dir

Makefile

@@ -56,6 +56,8 @@ clean-setup: ## run python setup.py clean
 
 clean-cython: ## clean the cython files
 	rm -f pyproj/*.so
+	rm -f pyproj/*/*.so
+	rm -f pyproj/*/*.c
 	rm -f pyproj/*.c
 
 lint: ## check style with flake8

README.md

@@ -5,7 +5,7 @@ Python interface to [PROJ](http://proj.org) (cartographic projections and coordi
 <p align="center">
 <a href="https://gitter.im/pyproj4-pyproj/community?utm_source=badge&utm_medium=badge&utm_campaign=pr-badge&utm_content=badge"><img alt="Join the chat at https://gitter.im/pyproj4-pyproj/community" src="https://badges.gitter.im/pyproj4-pyproj/community.svg"></a>
 <a href="#contributors"><img alt="All Contributors" src="https://img.shields.io/badge/all_contributors-28-orange.svg?style=flat-square"></a>
-<a href="https://travis-ci.org/pyproj4/pyproj"><img alt="Build Status" src="https://travis-ci.org/pyproj4/pyproj.svg"></a>
+<a href="https://travis-ci.org/pyproj4/pyproj"><img alt="Build Status" src="https://travis-ci.org/pyproj4/pyproj.svg?branch=master"></a>
 <a href="https://ci.appveyor.com/project/jswhit/pyproj"><img alt="Build Status" src="https://ci.appveyor.com/api/projects/status/8xkka4s97uwhkc64/branch/master?svg=true"></a>
 <a href="https://coveralls.io/github/pyproj4/pyproj?branch=master"><img alt="Coverage Status" src="https://coveralls.io/repos/github/pyproj4/pyproj/badge.svg?branch=master"></a>
 <a href="https://badge.fury.io/py/pyproj"><img alt="PyPI" src="https://badge.fury.io/py/pyproj.svg"></a>
@@ -73,6 +73,7 @@ Thanks goes to these wonderful people ([emoji key](https://allcontributors.org/d
   </tr>
   <tr>
     <td align="center"><a href="https://github.com/glostis"><img src="https://avatars0.githubusercontent.com/u/25295717?v=4" width="100px;" alt="Guillaume Lostis"/><br /><sub><b>Guillaume Lostis</b></sub></a><br /><a href="https://github.com/pyproj4/pyproj/commits?author=glostis" title="Documentation">📖</a></td>
+    <td align="center"><a href="https://github.com/edpop"><img src="https://avatars3.githubusercontent.com/u/13479292?v=4" width="100px;" alt="Eduard Popov"/><br /><sub><b>Eduard Popov</b></sub></a><br /><a href="https://github.com/pyproj4/pyproj/commits?author=edpop" title="Documentation">📖</a></td>
   </tr>
 </table>
 

appveyor.yml

@@ -21,19 +21,19 @@ environment:
     #   PYTHON_ARCH: "64"
     #   VS_VERSION: Visual Studio 14
     #   APPVEYOR_BUILD_WORKER_IMAGE: Visual Studio 2015
-    #   PROJSOURCE: 6.2.1
+    #   PROJSOURCE: 6.3.0
     - PYTHON: "C:\\Python36-x64"
       PYTHON_VERSION: "3.6"
       PYTHON_ARCH: "64"
       VS_VERSION: Visual Studio 14
       APPVEYOR_BUILD_WORKER_IMAGE: Visual Studio 2015
-      PROJSOURCE: 6.2.1
+      PROJSOURCE: 6.3.1
   # - PYTHON: "C:\\Python37-x64"
     #   PYTHON_VERSION: "3.7"
     #   PYTHON_ARCH: "64"
     #   VS_VERSION: Visual Studio 14
     #   APPVEYOR_BUILD_WORKER_IMAGE: Visual Studio 2015
-    #   PROJSOURCE: 6.2.1
+    #   PROJSOURCE: 6.3.1
     # - PYTHON: "C:\\Python36-x64"
     #   PYTHON_VERSION: "3.6"
     #   PYTHON_ARCH: "64"
@@ -66,7 +66,7 @@ build_script:
         appveyor DownloadFile https://sqlite.org/2018/sqlite-tools-win32-x86-3250100.zip
         7z x sqlite-tools-win32-x86-3250100.zip
   - copy "%APPVEYOR_BUILD_FOLDER%"\sqlite-tools-win32-x86-3250100\sqlite3.exe %SQLITE3_BIN%
-  - set PATH=%PATH%;%SQLITE3_BIN%
+  - set PATH=%SQLITE3_BIN%;%PATH%
   # setup PROJ.4
   - if "%PROJSOURCE%" == "git" git clone https://github.com/OSGeo/proj.4.git proj-git
   - if not "%PROJSOURCE%" == "git" if not exist %PROJ_DIR% set BUILD_PROJ_STABLE=1
@@ -85,7 +85,7 @@ build_script:
   - if defined BUILD_PROJ cd build
   - if defined BUILD_PROJ cmake -G "%VS_FULL%" .. -DCMAKE_BUILD_TYPE=Release -DBUILD_LIBPROJ_SHARED="%BUILD_LIBPROJ_SHARED%" -DCMAKE_C_FLAGS="/WX" -DCMAKE_CXX_FLAGS="/WX" -DCMAKE_TOOLCHAIN_FILE=c:/tools/vcpkg/scripts/buildsystems/vcpkg.cmake -DCMAKE_INSTALL_PREFIX="%PROJ_DIR%"
   - if defined BUILD_PROJ cmake --build . --config Release --target install
-  - set PATH=%PATH%;%PROJ_DIR%\bin
+  - set PATH=%PROJ_DIR%\bin;%PATH%
   - set PROJ_LIB=%PROJ_DIR%\share\proj
   - if defined BUILD_PROJ cd %PROJ_LIB%
   - if defined BUILD_PROJ curl -O http://download.osgeo.org/proj/proj-datumgrid-1.8.zip

debian/changelog

+python-pyproj (2.5.0+ds-1) unstable; urgency=medium
+
+  * New upstream release.
+  * Move from experimental to unstable.
+
+ -- Bas Couwenberg <sebastic@debian.org>  Sat, 22 Feb 2020 08:04:48 +0100
+
+python-pyproj (2.5~rc0+ds-1~exp1) experimental; urgency=medium
+
+  * New upstream release candidate.
+  * Drop patches, included upstream.
+
+ -- Bas Couwenberg <sebastic@debian.org>  Fri, 21 Feb 2020 06:55:18 +0100
+
 python-pyproj (2.4.2+ds-4) unstable; urgency=medium
 
   * Add upstream patches for PROJ 6.3.x support.

debian/patches/proj-6.3.0.patch

-Description: Add support for PROJ 6.3.0.
-Author: "Alan D. Snow" <alansnow21@gmail.com>
-Origin: https://github.com/pyproj4/pyproj/commit/82f01bdcd0cabd4a43b956db66c4b50d2c3351fb
-
---- a/pyproj/crs.py
-+++ b/pyproj/crs.py
-@@ -247,130 +247,6 @@ class CRS(_CRS):
-                 LENGTHUNIT["metre",1],
-                 ID["EPSG",8807]],
-             ID["EPSG",16015]]
--        >>> print(crs_utm.to_json(pretty=True))
--        {
--          "$schema": "https://proj.org/schemas/v0.1/projjson.schema.json",
--          "type": "ProjectedCRS",
--          "name": "NAD83 / UTM zone 15N",
--          "base_crs": {
--            "name": "NAD83",
--            "datum": {
--              "type": "GeodeticReferenceFrame",
--              "name": "North American Datum 1983",
--              "ellipsoid": {
--                "name": "GRS 1980",
--                "semi_major_axis": 6378137,
--                "inverse_flattening": 298.257222101
--              }
--            },
--            "coordinate_system": {
--              "subtype": "ellipsoidal",
--              "axis": [
--                {
--                  "name": "Geodetic latitude",
--                  "abbreviation": "Lat",
--                  "direction": "north",
--                  "unit": "degree"
--                },
--                {
--                  "name": "Geodetic longitude",
--                  "abbreviation": "Lon",
--                  "direction": "east",
--                  "unit": "degree"
--                }
--              ]
--            },
--            "id": {
--              "authority": "EPSG",
--              "code": 4269
--            }
--          },
--          "conversion": {
--            "name": "UTM zone 15N",
--            "method": {
--              "name": "Transverse Mercator",
--              "id": {
--                "authority": "EPSG",
--                "code": 9807
--              }
--            },
--            "parameters": [
--              {
--                "name": "Latitude of natural origin",
--                "value": 0,
--                "unit": "degree",
--                "id": {
--                  "authority": "EPSG",
--                  "code": 8801
--                }
--              },
--              {
--                "name": "Longitude of natural origin",
--                "value": -93,
--                "unit": "degree",
--                "id": {
--                  "authority": "EPSG",
--                  "code": 8802
--                }
--              },
--              {
--                "name": "Scale factor at natural origin",
--                "value": 0.9996,
--                "unit": "unity",
--                "id": {
--                  "authority": "EPSG",
--                  "code": 8805
--                }
--              },
--              {
--                "name": "False easting",
--                "value": 500000,
--                "unit": "metre",
--                "id": {
--                  "authority": "EPSG",
--                  "code": 8806
--                }
--              },
--              {
--                "name": "False northing",
--                "value": 0,
--                "unit": "metre",
--                "id": {
--                  "authority": "EPSG",
--                  "code": 8807
--                }
--              }
--            ]
--          },
--          "coordinate_system": {
--            "subtype": "Cartesian",
--            "axis": [
--              {
--                "name": "Easting",
--                "abbreviation": "E",
--                "direction": "east",
--                "unit": "metre"
--              },
--              {
--                "name": "Northing",
--                "abbreviation": "N",
--                "direction": "north",
--                "unit": "metre"
--              }
--            ]
--          },
--          "area": "North America - 96°W to 90°W and NAD83 by country",
--          "bbox": {
--            "south_latitude": 25.61,
--            "west_longitude": -96,
--            "north_latitude": 84,
--            "east_longitude": -90
--          },
--          "id": {
--            "authority": "EPSG",
--            "code": 26915
--          }
--        }
-         >>> crs = CRS(proj='utm', zone=10, ellps='WGS84')
-         >>> print(crs.to_wkt(pretty=True))
-         PROJCRS["unknown",
---- a/pyproj/transformer.py
-+++ b/pyproj/transformer.py
-@@ -736,23 +736,6 @@ def transform(
-     >>> xy = lons+lats
-     >>> '%8.3f %8.3f %8.3f %5.3f %5.3f %5.3f' % xy
-     ' -92.220  -94.720  -90.370 38.830 39.320 38.750'
--    >>> # test datum shifting, installation of extra datum grid files.
--    >>> p1 = Proj(proj='latlong',datum='WGS84')
--    >>> x1 = -111.5; y1 = 45.25919444444
--    >>> p2 = Proj(proj="utm",zone=10,datum='NAD27', preserve_units=False)
--    >>> x2, y2 = transform(p1, p2, x1, y1)
--    >>> "%s  %s" % (str(x2)[:9],str(y2)[:9])
--    '1402291.0  5076289.5'
--    >>> from pyproj import CRS
--    >>> c1 = CRS(proj='latlong',datum='WGS84')
--    >>> x1 = -111.5; y1 = 45.25919444444
--    >>> c2 = CRS(proj="utm",zone=10,datum='NAD27')
--    >>> x2, y2 = transform(c1, c2, x1, y1)
--    >>> "%s  %s" % (str(x2)[:9],str(y2)[:9])
--    '1402291.0  5076289.5'
--    >>> xeq, yeq = transform(4326, 4326, 30, 60, skip_equivalent=True)
--    >>> "%.0f %.0f" % (xeq, yeq)
--    '30 60'
-     """
-     return Transformer.from_proj(
-         p1, p2, skip_equivalent=skip_equivalent, always_xy=always_xy
---- a/test/test_crs.py
-+++ b/test/test_crs.py
-@@ -895,10 +895,6 @@ def test_to_dict_no_proj4():
-     )
-     if LooseVersion(proj_version_str) >= LooseVersion("6.3.0"):
-         with pytest.warns(UserWarning):
--            assert crs.to_proj4() == (
--                "+proj=ob_tran +o_proj=longlat +lon_0=-10 +o_lat_p=30 "
--                "+o_lon_p=0 +R=6371229 +no_defs +type=crs"
--            )
-             assert crs.to_dict() == {
-                 "R": 6371229,
-                 "lon_0": -10,
---- a/test/test_datum.py
-+++ b/test/test_datum.py
-@@ -1,12 +1,20 @@
-+from distutils.version import LooseVersion
-+
-+import pytest
- from numpy.testing import assert_almost_equal
- 
--from pyproj import Proj, transform
-+from pyproj import CRS, Proj, proj_version_str, transform
- 
- 
--def test_datum(aoi_data_directory):
--    p1 = Proj(proj="latlong", datum="WGS84")
-+@pytest.mark.parametrize("proj_class", [Proj, CRS])
-+def test_datum(proj_class, aoi_data_directory):
-+    p1 = proj_class(proj="latlong", datum="WGS84")
-     s_1 = -111.5
-     s_2 = 45.25919444444
--    p2 = Proj(proj="utm", zone=10, datum="NAD27")
-+    p2 = proj_class(proj="utm", zone=10, datum="NAD27")
-     x2, y2 = transform(p1, p2, s_1, s_2)
--    assert_almost_equal((x2, y2), (1402291.0833290431, 5076289.591846835))
-+    if LooseVersion(proj_version_str) < LooseVersion("6.3.0"):
-+        assert_almost_equal((x2, y2), (1402291.0833290431, 5076289.591846835))
-+    else:
-+        # https://github.com/OSGeo/PROJ/issues/1808
-+        assert_almost_equal((x2, y2), (1402285.9829252, 5076292.4212746))
---- a/test/test_transform.py
-+++ b/test/test_transform.py
-@@ -46,3 +46,8 @@ def test_transform():
-     assert_allclose(numpy.maximum.reduce(numpy.ravel(x3 - x1)), 0, atol=1e-4)
-     assert_allclose(numpy.minimum.reduce(numpy.ravel(y3 - y1)), 0, atol=1e-4)
-     assert_allclose(numpy.maximum.reduce(numpy.ravel(y3 - y1)), 0, atol=1e-4)
-+
-+
-+def test_skip_equivalent():
-+    xeq, yeq = transform(4326, 4326, 30, 60, skip_equivalent=True)
-+    assert (xeq, yeq) == (30, 60)
---- a/test/test_transformer.py
-+++ b/test/test_transformer.py
-@@ -261,11 +261,7 @@ def test_transform_no_exception():
- def test_transform__out_of_bounds():
-     with pytest.warns(FutureWarning):
-         transformer = Transformer.from_proj("+init=epsg:4326", "+init=epsg:27700")
--    if LooseVersion(proj_version_str) >= LooseVersion("6.3.0"):
--        with pytest.raises(ProjError):
--            transformer.transform(100000, 100000, errcheck=True)
--    else:
--        assert np.all(np.isinf(transformer.transform(100000, 100000, errcheck=True)))
-+    assert np.all(np.isinf(transformer.transform(100000, 100000, errcheck=True)))
- 
- 
- def test_transform_radians():
-@@ -496,13 +492,10 @@ def test_transformer__operations__scope_
-         "GNSS CORS.",
-         None,
-     ]
--    assert [op.remarks for op in transformer.operations] == [
--        None,
--        "Scale difference in ppb where 1/billion = 1E-9. Derivation excluded Cocos, "
--        "Christmas and Macquarie Islands but is applied there. See codes 8444-46 for "
--        "equivalents using NTv2 method. See code 8447 for alternative including "
--        "distortion model for Aus only.",
--        None,
-+    assert [str(op.remarks)[:5] for op in transformer.operations] == [
-+        "None",
-+        "Scale",
-+        "None",
-     ]
- 
- 

debian/patches/proj-6.3.1.patch

-Description: Add support for PROJ 6.3.1.
-Author: "Alan D. Snow" <alansnow21@gmail.com>
-Origin: https://github.com/pyproj4/pyproj/commit/5bd38314c58dd414da034fbe21b290cffa32e4a9
-
---- a/test/test_transformer.py
-+++ b/test/test_transformer.py
-@@ -551,7 +551,10 @@ def test_transform_group__area_of_intere
- 
- def test_transformer_group__get_transform_crs():
-     tg = TransformerGroup("epsg:4258", "epsg:7415")
--    assert len(tg.transformers) == 4
-+    if LooseVersion(pyproj.proj_version_str) >= LooseVersion("6.3.1"):
-+        assert len(tg.transformers) == 1
-+    else:
-+        assert len(tg.transformers) == 4
- 
- 
- def test_transformer__area_of_interest(aoi_data_directory):

debian/patches/proj-7.0.0.patch

-Description: Fix test failure with PROJ 7.0.0.
-Author: Bas Couwenberg <sebastic@debian.org>
-Bug: https://github.com/pyproj4/pyproj/issues/531
-Forwarded: https://github.com/pyproj4/pyproj/issues/531#issuecomment-588394294
-Applied-Upstream: https://github.com/pyproj4/pyproj/commit/cf2a1f92c7afa543ca152dab57d8944f0c85f689
-
---- a/test/test_crs.py
-+++ b/test/test_crs.py
-@@ -519,13 +519,22 @@ def test_coordinate_operation_grids__alt
-     cc = CoordinateOperation.from_epsg(1312, True)
-     assert len(cc.grids) == 1
-     grid = cc.grids[0]
--    assert grid.short_name == "ntv1_can.dat"
--    assert grid.full_name.endswith(grid.short_name)
--    assert grid.package_name == "proj-datumgrid"
--    assert grid.url.startswith("https://download.osgeo.org/proj/proj-datumgrid")
--    assert grid.direct_download is True
--    assert grid.open_license is True
--    assert grid.available is True
-+    if LooseVersion(proj_version_str) >= LooseVersion("7.0.0"):
-+        assert grid.short_name == "ca_nrc_ntv1_can.tif"
-+        assert grid.full_name.endswith("ntv1_can.dat")
-+        assert grid.package_name == ""
-+        assert grid.url == "https://cdn.proj.org/ca_nrc_ntv1_can.tif"
-+        assert grid.direct_download is True
-+        assert grid.open_license is True
-+        assert grid.available is True
-+    else:
-+        assert grid.short_name == "ntv1_can.dat"
-+        assert grid.full_name.endswith(grid.short_name)
-+        assert grid.package_name == "proj-datumgrid"
-+        assert grid.url.startswith("https://download.osgeo.org/proj/proj-datumgrid")
-+        assert grid.direct_download is True
-+        assert grid.open_license is True
-+        assert grid.available is True
- 
- 
- def test_coordinate_operation__missing():
---- a/test/test_transformer.py
-+++ b/test/test_transformer.py
-@@ -261,7 +261,11 @@ def test_transform_no_exception():
- def test_transform__out_of_bounds():
-     with pytest.warns(FutureWarning):
-         transformer = Transformer.from_proj("+init=epsg:4326", "+init=epsg:27700")
--    assert np.all(np.isinf(transformer.transform(100000, 100000, errcheck=True)))
-+    if LooseVersion(pyproj.proj_version_str) >= LooseVersion("7.0.0"):
-+        with pytest.raises(pyproj.exceptions.ProjError):
-+            assert np.all(np.isinf(transformer.transform(100000, 100000, errcheck=True)))
-+    else:
-+        assert np.all(np.isinf(transformer.transform(100000, 100000, errcheck=True)))
- 
- 
- def test_transform_radians():

debian/patches/series

-proj-6.3.0.patch
-proj-6.3.1.patch
-proj-7.0.0.patch

pyproj/__init__.py

@@ -47,7 +47,7 @@ CONSEQUENTIAL DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM
 LOSS OF USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT,
 NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN
 CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE. """
-__version__ = "2.4.2"
+__version__ = "2.5.0"
 __all__ = [
     "Proj",
     "Geod",

pyproj/_crs.pxd

@@ -14,7 +14,7 @@ cdef class Axis:
     cdef readonly object unit_code
 
     @staticmethod
-    cdef create(PJ_CONTEXT* context, PJ* projobj, int index)
+    cdef Axis create(PJ_CONTEXT* context, PJ* projobj, int index)
 
 cdef class AreaOfUse:
     cdef readonly double west
@@ -24,7 +24,7 @@ cdef class AreaOfUse:
     cdef readonly object name
 
     @staticmethod
-    cdef create(PJ_CONTEXT* context, PJ* projobj)
+    cdef AreaOfUse create(PJ_CONTEXT* context, PJ* projobj)
 
 
 cdef class Base:
@@ -35,38 +35,43 @@ cdef class Base:
     cdef readonly object _scope
 
 
-cdef class Ellipsoid(Base):
-    cdef double _semi_major_metre
-    cdef double _semi_minor_metre
+cdef class _CRSParts(Base):
+    pass
+
+
+cdef class Ellipsoid(_CRSParts):
+    cdef readonly double semi_major_metre
+    cdef readonly double semi_minor_metre
     cdef readonly object is_semi_minor_computed
-    cdef double _inv_flattening
-    cdef readonly object ellipsoid_loaded
+    cdef readonly double inverse_flattening
 
     @staticmethod
-    cdef create(PJ_CONTEXT* context, PJ* ellipsoid_pj)
+    cdef Ellipsoid create(PJ_CONTEXT* context, PJ* ellipsoid_pj)
 
-cdef class PrimeMeridian(Base):
+
+cdef class PrimeMeridian(_CRSParts):
     cdef readonly double longitude
     cdef readonly double unit_conversion_factor
     cdef readonly object unit_name
 
     @staticmethod
-    cdef create(PJ_CONTEXT* context, PJ* prime_meridian_pj)
+    cdef PrimeMeridian create(PJ_CONTEXT* context, PJ* prime_meridian_pj)
 
 
-cdef class Datum(Base):
+cdef class Datum(_CRSParts):
+    cdef readonly object type_name
     cdef readonly object _ellipsoid
     cdef readonly object _prime_meridian
 
     @staticmethod
-    cdef create(PJ_CONTEXT* context, PJ* datum_pj)
+    cdef Datum create(PJ_CONTEXT* context, PJ* datum_pj)
 
 
-cdef class CoordinateSystem(Base):
+cdef class CoordinateSystem(_CRSParts):
     cdef readonly object _axis_list
 
     @staticmethod
-    cdef create(PJ_CONTEXT* context, PJ* coordinate_system_pj)
+    cdef CoordinateSystem create(PJ_CONTEXT* context, PJ* coordinate_system_pj)
 
 
 cdef class Param:
@@ -81,7 +86,7 @@ cdef class Param:
     cdef readonly object unit_category
 
     @staticmethod
-    cdef create(PJ_CONTEXT* context, PJ* projobj, int param_idx)
+    cdef Param create(PJ_CONTEXT* context, PJ* projobj, int param_idx)
 
 
 cdef class Grid:
@@ -94,10 +99,10 @@ cdef class Grid:
     cdef readonly object available
 
     @staticmethod
-    cdef create(PJ_CONTEXT* context, PJ* projobj, int grid_idx)
+    cdef Grid create(PJ_CONTEXT* context, PJ* projobj, int grid_idx)
 
 
-cdef class CoordinateOperation(Base):
+cdef class CoordinateOperation(_CRSParts):
     cdef readonly object _params
     cdef readonly object _grids
     cdef readonly object _area_of_use
@@ -112,7 +117,7 @@ cdef class CoordinateOperation(Base):
     cdef readonly type_name
 
     @staticmethod
-    cdef create(PJ_CONTEXT* context, PJ* coordinate_operation_pj)
+    cdef CoordinateOperation create(PJ_CONTEXT* context, PJ* coordinate_operation_pj)
 
 
 cdef class _CRS(Base):

pyproj/_geod.pyx

@@ -39,61 +39,45 @@ cdef class Geod:
         forward azimuth and distance.
         if radians=True, lons/lats are radians instead of degrees.
         """
-        cdef Py_ssize_t buflenlons, buflenlats, buflenaz, buflend, ndim, iii
-        cdef double lat1,lon1,az1,s12,plon2,plat2,pazi2
-        cdef double *lonsdata
-        cdef double *latsdata
-        cdef double *azdata
-        cdef double *distdata
-        cdef void *londata
-        cdef void *latdata
-        cdef void *azdat
-        cdef void *distdat
-        # if buffer api is supported, get pointer to data buffers.
-        if PyObject_AsWriteBuffer(lons, &londata, &buflenlons) <> 0:
-            raise GeodError
-        if PyObject_AsWriteBuffer(lats, &latdata, &buflenlats) <> 0:
-            raise GeodError
-        if PyObject_AsWriteBuffer(az, &azdat, &buflenaz) <> 0:
-            raise GeodError
-        if PyObject_AsWriteBuffer(dist, &distdat, &buflend) <> 0:
-            raise GeodError
+        cdef PyBuffWriteManager lonbuff = PyBuffWriteManager(lons)
+        cdef PyBuffWriteManager latbuff = PyBuffWriteManager(lats)
+        cdef PyBuffWriteManager azbuff = PyBuffWriteManager(az)
+        cdef PyBuffWriteManager distbuff = PyBuffWriteManager(dist)
+  
         # process data in buffer
-        if not buflenlons == buflenlats == buflenaz == buflend:
-            raise GeodError("Buffer lengths not the same")
-        ndim = buflenlons//_DOUBLESIZE
-        lonsdata = <double *>londata
-        latsdata = <double *>latdata
-        azdata = <double *>azdat
-        distdata = <double *>distdat
+        if not lonbuff.len == latbuff.len == azbuff.len == distbuff.len:
+            raise GeodError("Array lengths are not the same.")
+
+        cdef double lat1, lon1, az1, s12, plon2, plat2, pazi2
+        cdef Py_ssize_t iii
         with nogil:
-            for iii in range(ndim):
+            for iii in range(lonbuff.len):
                 if not radians:
-                    lon1 = lonsdata[iii]
-                    lat1 = latsdata[iii]
-                    az1 = azdata[iii]
-                    s12 = distdata[iii]
+                    lon1 = lonbuff.data[iii]
+                    lat1 = latbuff.data[iii]
+                    az1 = azbuff.data[iii]
+                    s12 = distbuff.data[iii]
                 else:
-                    lon1 = _RAD2DG * lonsdata[iii]
-                    lat1 = _RAD2DG * latsdata[iii]
-                    az1 = _RAD2DG * azdata[iii]
-                    s12 = distdata[iii]
+                    lon1 = _RAD2DG * lonbuff.data[iii]
+                    lat1 = _RAD2DG * latbuff.data[iii]
+                    az1 = _RAD2DG * azbuff.data[iii]
+                    s12 = distbuff.data[iii]
                 geod_direct(&self._geod_geodesic, lat1, lon1, az1, s12,\
                     &plat2, &plon2, &pazi2)
                 # back azimuth needs to be flipped 180 degrees
-                # to match what proj4 geod utility produces.
+                # to match what PROJ geod utility produces.
                 if pazi2 > 0:
                     pazi2 = pazi2 - 180.
                 elif pazi2 <= 0:
                     pazi2 = pazi2 + 180.
                 if not radians:
-                    lonsdata[iii] = plon2
-                    latsdata[iii] = plat2
-                    azdata[iii] = pazi2
+                    lonbuff.data[iii] = plon2
+                    latbuff.data[iii] = plat2
+                    azbuff.data[iii] = pazi2
                 else:
-                    lonsdata[iii] = _DG2RAD * plon2
-                    latsdata[iii] = _DG2RAD * plat2
-                    azdata[iii] = _DG2RAD * pazi2
+                    lonbuff.data[iii] = _DG2RAD * plon2
+                    latbuff.data[iii] = _DG2RAD * plat2
+                    azbuff.data[iii] = _DG2RAD * pazi2
 
     @cython.boundscheck(False)
     @cython.wraparound(False)
@@ -103,45 +87,29 @@ cdef class Geod:
         between an initial and terminus lat/lon pair.
         if radians=True, lons/lats are radians instead of degree
         """
-        cdef double lat1,lon1,lat2,lon2,pazi1,pazi2,ps12
-        cdef Py_ssize_t buflenlons, buflenlats, buflenaz, buflend, ndim, iii
-        cdef double *lonsdata
-        cdef double *latsdata
-        cdef double *azdata
-        cdef double *distdata
-        cdef void *londata
-        cdef void *latdata
-        cdef void *azdat
-        cdef void *distdat
-        # if buffer api is supported, get pointer to data buffers.
-        if PyObject_AsWriteBuffer(lons1, &londata, &buflenlons) <> 0:
-            raise GeodError
-        if PyObject_AsWriteBuffer(lats1, &latdata, &buflenlats) <> 0:
-            raise GeodError
-        if PyObject_AsWriteBuffer(lons2, &azdat, &buflenaz) <> 0:
-            raise GeodError
-        if PyObject_AsWriteBuffer(lats2, &distdat, &buflend) <> 0:
-            raise GeodError
+        cdef PyBuffWriteManager lon1buff = PyBuffWriteManager(lons1)
+        cdef PyBuffWriteManager lat1buff = PyBuffWriteManager(lats1)
+        cdef PyBuffWriteManager lon2buff = PyBuffWriteManager(lons2)
+        cdef PyBuffWriteManager lat2buff = PyBuffWriteManager(lats2)
+
         # process data in buffer
-        if not buflenlons == buflenlats == buflenaz == buflend:
-            raise GeodError("Buffer lengths not the same")
-        ndim = buflenlons//_DOUBLESIZE
-        lonsdata = <double *>londata
-        latsdata = <double *>latdata
-        azdata = <double *>azdat
-        distdata = <double *>distdat
+        if not lon1buff.len == lat1buff.len == lon2buff.len == lat2buff.len:
+            raise GeodError("Array lengths are not the same.")
+
+        cdef double lat1, lon1, lat2, lon2, pazi1, pazi2, ps12
+        cdef Py_ssize_t iii
         with nogil:
-            for iii in range(ndim):
+            for iii in range(lon1buff.len):
                 if radians:
-                    lon1 = _RAD2DG * lonsdata[iii]
-                    lat1 = _RAD2DG * latsdata[iii]
-                    lon2 = _RAD2DG * azdata[iii]
-                    lat2 = _RAD2DG * distdata[iii]
+                    lon1 = _RAD2DG * lon1buff.data[iii]
+                    lat1 = _RAD2DG * lat1buff.data[iii]
+                    lon2 = _RAD2DG * lon2buff.data[iii]
+                    lat2 = _RAD2DG * lat2buff.data[iii]
                 else:
-                    lon1 = lonsdata[iii]
-                    lat1 = latsdata[iii]
-                    lon2 = azdata[iii]
-                    lat2 = distdata[iii]
+                    lon1 = lon1buff.data[iii]
+                    lat1 = lat1buff.data[iii]
+                    lon2 = lon2buff.data[iii]
+                    lat2 = lat2buff.data[iii]
                 geod_inverse(
                     &self._geod_geodesic,
                     lat1, lon1, lat2, lon2,
@@ -154,12 +122,13 @@ cdef class Geod:
                 elif pazi2 <= 0:
                     pazi2 = pazi2+180.
                 if radians:
-                    lonsdata[iii] = _DG2RAD * pazi1
-                    latsdata[iii] = _DG2RAD * pazi2
+                    lon1buff.data[iii] = _DG2RAD * pazi1
+                    lat1buff.data[iii] = _DG2RAD * pazi2
                 else:
-                    lonsdata[iii] = pazi1
-                    latsdata[iii] = pazi2
-                azdata[iii] = ps12
+                    lon1buff.data[iii] = pazi1
+                    lat1buff.data[iii] = pazi2
+                # write azimuth data into lon2 buffer
+                lon2buff.data[iii] = ps12
 
     @cython.boundscheck(False)
     @cython.wraparound(False)
@@ -223,46 +192,38 @@ cdef class Geod:
         float: The total distance.
     
         """
-        cdef double lat1,lon1,lat2,lon2,pazi1,pazi2,ps12
-        cdef double total_distance = 0.0
-        cdef Py_ssize_t buflenlons, buflenlats, ndim, iii
-        cdef double *lonsdata
-        cdef double *latsdata
-        cdef void *londata
-        cdef void *latdata
-        # if buffer api is supported, get pointer to data buffers.
-        if PyObject_AsWriteBuffer(lons, &londata, &buflenlons) <> 0:
-            raise GeodError
-        if PyObject_AsWriteBuffer(lats, &latdata, &buflenlats) <> 0:
-            raise GeodError
+        cdef PyBuffWriteManager lonbuff = PyBuffWriteManager(lons)
+        cdef PyBuffWriteManager latbuff = PyBuffWriteManager(lats)
+
         # process data in buffer
-        if buflenlons != buflenlats:
-            raise GeodError("Buffer lengths not the same")
-        ndim = buflenlons//_DOUBLESIZE
-        lonsdata = <double *>londata
-        latsdata = <double *>latdata
-
-        if ndim == 1:
-            lonsdata[0] = 0
+        if lonbuff.len != latbuff.len:
+            raise GeodError("Array lengths are not the same.")
+
+        if lonbuff.len == 1:
+            lonbuff.data[0] = 0
             return 0.0
+
+        cdef double lat1, lon1, lat2, lon2, pazi1, pazi2, ps12
+        cdef double total_distance = 0.0
+        cdef Py_ssize_t iii
         with nogil:
-            for iii in range(ndim - 1):
+            for iii in range(lonbuff.len - 1):
                 if radians:
-                    lon1 = _RAD2DG * lonsdata[iii]
-                    lat1 = _RAD2DG * latsdata[iii]
-                    lon2 = _RAD2DG * lonsdata[iii + 1]
-                    lat2 = _RAD2DG * latsdata[iii + 1]
+                    lon1 = _RAD2DG * lonbuff.data[iii]
+                    lat1 = _RAD2DG * latbuff.data[iii]
+                    lon2 = _RAD2DG * lonbuff.data[iii + 1]
+                    lat2 = _RAD2DG * latbuff.data[iii + 1]
                 else:
-                    lon1 = lonsdata[iii]
-                    lat1 = latsdata[iii]
-                    lon2 = lonsdata[iii + 1]
-                    lat2 = latsdata[iii + 1]
+                    lon1 = lonbuff.data[iii]
+                    lat1 = latbuff.data[iii]
+                    lon2 = lonbuff.data[iii + 1]
+                    lat2 = latbuff.data[iii + 1]
                 geod_inverse(
                     &self._geod_geodesic,
                     lat1, lon1, lat2, lon2,
                     &ps12, &pazi1, &pazi2,
                 )
-                lonsdata[iii] = ps12
+                lonbuff.data[iii] = ps12
                 total_distance += ps12
         return total_distance
 
@@ -293,35 +254,26 @@ cdef class Geod:
         (float, float): The area (meter^2) and permimeter (meters) of the polygon.
 
         """
-        cdef Py_ssize_t buflenlons, buflenlats, ndim, iii
-        cdef void *londata
-        cdef void *latdata
-        cdef double *lonsdata
-        cdef double *latsdata
-        # if buffer api is supported, get pointer to data buffers.
-        if PyObject_AsWriteBuffer(lons, &londata, &buflenlons) <> 0:
-            raise GeodError
-        if PyObject_AsWriteBuffer(lats, &latdata, &buflenlats) <> 0:
-            raise GeodError
+        cdef PyBuffWriteManager lonbuff = PyBuffWriteManager(lons)
+        cdef PyBuffWriteManager latbuff = PyBuffWriteManager(lats)
+
         # process data in buffer
-        if not buflenlons == buflenlats:
-            raise GeodError("Buffer lengths not the same")
+        if not lonbuff.len == latbuff.len:
+            raise GeodError("Array lengths are not the same.")
 
         cdef double polygon_area
         cdef double polygon_perimeter
-        ndim = buflenlons//_DOUBLESIZE
+        cdef Py_ssize_t iii
 
-        lonsdata = <double *>londata
-        latsdata = <double *>latdata
         with nogil:
             if radians:
-                for iii in range(ndim):
-                    lonsdata[iii] *= _RAD2DG
-                    latsdata[iii] *= _RAD2DG
+                for iii in range(lonbuff.len):
+                    lonbuff.data[iii] *= _RAD2DG
+                    latbuff.data[iii] *= _RAD2DG
 
             geod_polygonarea(
                 &self._geod_geodesic,
-                latsdata, lonsdata, ndim, 
+                latbuff.data, lonbuff.data, lonbuff.len, 
                 &polygon_area, &polygon_perimeter
             )
         return (polygon_area, polygon_perimeter)

pyproj/_list.pyx

@@ -165,7 +165,7 @@ def get_codes(auth_name, pj_type, allow_deprecated=False):
     ----------
     auth_name: str
         The name of the authority.
-    pj_type: ~pyproj.enums.PJType
+    pj_type: pyproj.enums.PJType
         The type of object to get the authorities.
     allow_deprecated: bool, optional
         Allow a deprecated code in the return.

pyproj/_proj.pyx

@@ -56,42 +56,35 @@ cdef class Proj:
         if errcheck=False and the forward transformation is invalid, no exception is
         raised and 'inf' is returned.
         """
+        cdef PyBuffWriteManager lonbuff = PyBuffWriteManager(lons)
+        cdef PyBuffWriteManager latbuff = PyBuffWriteManager(lats)
+
+        # process data in buffer
+        if lonbuff.len != latbuff.len:
+            raise ProjError("Buffer lengths not the same")
+
         cdef PJ_COORD projxyout
         cdef PJ_COORD projlonlatin = proj_coord(0, 0, 0, HUGE_VAL)
-        cdef Py_ssize_t buflenx, bufleny, ndim, iii
-        cdef double *lonsdata
-        cdef double *latsdata
-        cdef void *londata
-        cdef void *latdata
+        cdef Py_ssize_t iii
         cdef int errno
-        # if buffer api is supported, get pointer to data buffers.
-        if PyObject_AsWriteBuffer(lons, &londata, &buflenx) <> 0:
-            raise ProjError("object does not provide the python buffer writeable interface")
-        if PyObject_AsWriteBuffer(lats, &latdata, &bufleny) <> 0:
-            raise ProjError("object does not provide the python buffer writeable interface")
-        # process data in buffer
-        if buflenx != bufleny:
-            raise ProjError("Buffer lengths not the same")
-        ndim = buflenx//_DOUBLESIZE
-        lonsdata = <double *>londata
-        latsdata = <double *>latdata
+
         with nogil:
             proj_errno_reset(self.projobj)
-            for iii in range(ndim):
+            for iii in range(latbuff.len):
                 # if inputs are nan's, return big number.
-                if lonsdata[iii] != lonsdata[iii] or latsdata[iii] != latsdata[iii]:
-                    lonsdata[iii] = HUGE_VAL
-                    latsdata[iii] = HUGE_VAL
+                if lonbuff.data[iii] != lonbuff.data[iii] or latbuff.data[iii] != latbuff.data[iii]:
+                    lonbuff.data[iii] = HUGE_VAL
+                    latbuff.data[iii] = HUGE_VAL
                     if errcheck:
                         with gil:
                             raise ProjError("projection_undefined")
                     continue
                 if proj_angular_input(self.projobj, PJ_FWD):
-                    projlonlatin.uv.u = _DG2RAD * lonsdata[iii]
-                    projlonlatin.uv.v = _DG2RAD * latsdata[iii]
+                    projlonlatin.uv.u = _DG2RAD * lonbuff.data[iii]
+                    projlonlatin.uv.v = _DG2RAD * latbuff.data[iii]
                 else:
-                    projlonlatin.uv.u = lonsdata[iii]
-                    projlonlatin.uv.v = latsdata[iii]
+                    projlonlatin.uv.u = lonbuff.data[iii]
+                    projlonlatin.uv.v = latbuff.data[iii]
                 projxyout = proj_trans(self.projobj, PJ_FWD, projlonlatin)
                 errno = proj_errno(self.projobj)
                 if errcheck and errno:
@@ -111,20 +104,20 @@ cdef class Proj:
                     if errcheck:
                         with gil:
                             raise ProjError("Projection undefined.")
-                    lonsdata[iii] = HUGE_VAL
-                    latsdata[iii] = HUGE_VAL
+                    lonbuff.data[iii] = HUGE_VAL
+                    latbuff.data[iii] = HUGE_VAL
                 elif proj_angular_output(self.projobj, PJ_FWD):
-                    lonsdata[iii] = _RAD2DG * projxyout.xy.x
-                    latsdata[iii] = _RAD2DG * projxyout.xy.y
+                    lonbuff.data[iii] = _RAD2DG * projxyout.xy.x
+                    latbuff.data[iii] = _RAD2DG * projxyout.xy.y
                 else:
-                    lonsdata[iii] = projxyout.xy.x
-                    latsdata[iii] = projxyout.xy.y
+                    lonbuff.data[iii] = projxyout.xy.x
+                    latbuff.data[iii] = projxyout.xy.y
         ProjError.clear()
 
 
     @cython.boundscheck(False)
     @cython.wraparound(False)
-    def _inv(self, object x, object y, bint errcheck=False):
+    def _inv(self, object xx, object yy, bint errcheck=False):
         """
         inverse transformation - x,y to lons,lats (done in place).
         if errcheck=True, an exception is raised if the inverse transformation is invalid.
@@ -134,45 +127,36 @@ cdef class Proj:
         if not self.has_inverse:
             raise ProjError('inverse projection undefined')
 
+        cdef PyBuffWriteManager xbuff = PyBuffWriteManager(xx)
+        cdef PyBuffWriteManager ybuff = PyBuffWriteManager(yy)
+
+        # process data in buffer
+        if xbuff.len != ybuff.len:
+            raise ProjError("Array lengths not the same.")
+
         cdef PJ_COORD projxyin = proj_coord(0, 0, 0, HUGE_VAL)
         cdef PJ_COORD projlonlatout
-        cdef Py_ssize_t buflenx, bufleny, ndim, iii
-        cdef void *xdata
-        cdef void *ydata
-        cdef double *xdatab
-        cdef double *ydatab
+        cdef Py_ssize_t iii
         cdef int errno
 
-        # if buffer api is supported, get pointer to data buffers.
-        if PyObject_AsWriteBuffer(x, &xdata, &buflenx) <> 0:
-            raise ProjError("object does not provide the python buffer writeable interface")
-        if PyObject_AsWriteBuffer(y, &ydata, &bufleny) <> 0:
-            raise ProjError("object does not provide the python buffer writeable interface")
-        # process data in buffer
-        # (for numpy/regular python arrays).
-        if buflenx != bufleny:
-            raise ProjError("Buffer lengths not the same")
-        ndim = buflenx//_DOUBLESIZE
-        xdatab = <double *>xdata
-        ydatab = <double *>ydata
         with nogil:
             # reset errors potentially left over
             proj_errno_reset(self.projobj)
-            for iii in range(ndim):
+            for iii in range(xbuff.len):
                 # if inputs are nan's, return big number.
-                if xdatab[iii] != xdatab[iii] or ydatab[iii] != ydatab[iii]:
-                    xdatab[iii] = HUGE_VAL
-                    ydatab[iii] = HUGE_VAL
+                if xbuff.data[iii] != xbuff.data[iii] or ybuff.data[iii] != ybuff.data[iii]:
+                    xbuff.data[iii] = HUGE_VAL
+                    ybuff.data[iii] = HUGE_VAL
                     if errcheck:
                         with gil:
                             raise ProjError("projection_undefined")
                     continue
                 if proj_angular_input(self.projobj, PJ_INV):
-                    projxyin.xy.x = _DG2RAD * xdatab[iii]
-                    projxyin.xy.y = _DG2RAD * ydatab[iii]
+                    projxyin.xy.x = _DG2RAD * xbuff.data[iii]
+                    projxyin.xy.y = _DG2RAD * ybuff.data[iii]
                 else:
-                    projxyin.xy.x = xdatab[iii]
-                    projxyin.xy.y = ydatab[iii]
+                    projxyin.xy.x = xbuff.data[iii]
+                    projxyin.xy.y = ybuff.data[iii]
                 projlonlatout = proj_trans(self.projobj, PJ_INV, projxyin)
                 errno = proj_errno(self.projobj)
                 if errcheck and errno:
@@ -193,14 +177,14 @@ cdef class Proj:
                     if errcheck:
                         with gil:
                             raise ProjError("projection_undefined")
-                    xdatab[iii] = HUGE_VAL
-                    ydatab[iii] = HUGE_VAL
+                    xbuff.data[iii] = HUGE_VAL
+                    ybuff.data[iii] = HUGE_VAL
                 elif proj_angular_output(self.projobj, PJ_INV):
-                    xdatab[iii] = _RAD2DG * projlonlatout.uv.u
-                    ydatab[iii] = _RAD2DG * projlonlatout.uv.v
+                    xbuff.data[iii] = _RAD2DG * projlonlatout.uv.u
+                    ybuff.data[iii] = _RAD2DG * projlonlatout.uv.v
                 else:
-                    xdatab[iii] = projlonlatout.uv.u
-                    ydatab[iii] = projlonlatout.uv.v
+                    xbuff.data[iii] = projlonlatout.uv.u
+                    ybuff.data[iii] = projlonlatout.uv.v
         ProjError.clear()
 
     def __repr__(self):

pyproj/_transformer.pyx

@@ -404,68 +404,59 @@ cdef class _Transformer(Base):
             return
         tmp_pj_direction = _PJ_DIRECTION_MAP[TransformDirection.create(direction)]
         cdef PJ_DIRECTION pj_direction = <PJ_DIRECTION>tmp_pj_direction
-        # private function to call pj_transform
-        cdef void *xdata
-        cdef void *ydata
-        cdef void *zdata
-        cdef void *tdata
-        cdef double *xx
-        cdef double *yy
-        cdef double *zz
-        cdef double *tt
-        cdef Py_ssize_t buflenx, bufleny, buflenz, buflent, npts, iii
-        cdef int err
-        if PyObject_AsWriteBuffer(inx, &xdata, &buflenx) <> 0:
-            raise ProjError("object does not provide the python buffer writeable interface")
-        if PyObject_AsWriteBuffer(iny, &ydata, &bufleny) <> 0:
-            raise ProjError("object does not provide the python buffer writeable interface")
-        if inz is not None and PyObject_AsWriteBuffer(inz, &zdata, &buflenz) <> 0:
-            raise ProjError("object does not provide the python buffer writeable interface")
-        else:
-            buflenz = bufleny
-        if intime is not None and PyObject_AsWriteBuffer(intime, &tdata, &buflent) <> 0:
-            raise ProjError("object does not provide the python buffer writeable interface")
-        else:
-            buflent = bufleny
+        cdef PyBuffWriteManager xbuff = PyBuffWriteManager(inx)
+        cdef PyBuffWriteManager ybuff = PyBuffWriteManager(iny)
 
-        if not buflenx or not (buflenx == bufleny == buflenz == buflent):
-            raise ProjError('x,y,z, and time must be same size')
-        xx = <double *>xdata
-        yy = <double *>ydata
+        cdef PyBuffWriteManager zbuff
+        cdef Py_ssize_t buflenz
+        cdef double* zz
         if inz is not None:
-            zz = <double *>zdata
+            zbuff = PyBuffWriteManager(inz)
+            buflenz = zbuff.len
+            zz = zbuff.data
         else:
+            buflenz = xbuff.len
             zz = NULL
+
+        cdef PyBuffWriteManager tbuff
+        cdef Py_ssize_t buflent
+        cdef double* tt
         if intime is not None:
-            tt = <double *>tdata
+            tbuff = PyBuffWriteManager(intime)
+            buflent = tbuff.len
+            tt = tbuff.data
         else:
+            buflent = xbuff.len
             tt = NULL
-        npts = buflenx//8
 
+        if not xbuff.len or not (xbuff.len == ybuff.len == buflenz == buflent):
+            raise ProjError('x, y, z, and time must be same size')
+
+        cdef Py_ssize_t iii
         # degrees to radians
         if not self.is_pipeline and not radians\
                 and self._input_radians[pj_direction]:
             with nogil:
-                for iii in range(npts):
-                    xx[iii] = xx[iii]*_DG2RAD
-                    yy[iii] = yy[iii]*_DG2RAD
+                for iii in range(xbuff.len):
+                    xbuff.data[iii] = xbuff.data[iii]*_DG2RAD
+                    ybuff.data[iii] = ybuff.data[iii]*_DG2RAD
         # radians to degrees
         elif not self.is_pipeline and radians\
                 and not self._input_radians[pj_direction]\
                 and self.input_geographic:
             with nogil:
-                for iii in range(npts):
-                    xx[iii] = xx[iii]*_RAD2DG
-                    yy[iii] = yy[iii]*_RAD2DG
+                for iii in range(xbuff.len):
+                    xbuff.data[iii] = xbuff.data[iii]*_RAD2DG
+                    ybuff.data[iii] = ybuff.data[iii]*_RAD2DG
         with nogil:
             proj_errno_reset(self.projobj)
             proj_trans_generic(
                 self.projobj,
                 pj_direction,
-                xx, _DOUBLESIZE, npts,
-                yy, _DOUBLESIZE, npts,
-                zz, _DOUBLESIZE, npts,
-                tt, _DOUBLESIZE, npts,
+                xbuff.data, _DOUBLESIZE, xbuff.len,
+                ybuff.data, _DOUBLESIZE, ybuff.len,
+                zz, _DOUBLESIZE, xbuff.len,
+                tt, _DOUBLESIZE, xbuff.len,
             )
         cdef int errno = proj_errno(self.projobj)
         if errcheck and errno:
@@ -478,23 +469,23 @@ cdef class _Transformer(Base):
         if not self.is_pipeline and not radians\
                 and self._output_radians[pj_direction]:
             with nogil:
-                for iii in range(npts):
-                    xx[iii] = xx[iii]*_RAD2DG
-                    yy[iii] = yy[iii]*_RAD2DG
+                for iii in range(xbuff.len):
+                    xbuff.data[iii] = xbuff.data[iii]*_RAD2DG
+                    ybuff.data[iii] = ybuff.data[iii]*_RAD2DG
         # degrees to radians
         elif not self.is_pipeline and radians\
                 and not self._output_radians[pj_direction]\
                 and self.output_geographic:
             with nogil:
-                for iii in range(npts):
-                    xx[iii] = xx[iii]*_DG2RAD
-                    yy[iii] = yy[iii]*_DG2RAD
+                for iii in range(xbuff.len):
+                    xbuff.data[iii] = xbuff.data[iii]*_DG2RAD
+                    ybuff.data[iii] = ybuff.data[iii]*_DG2RAD
         ProjError.clear()
 
     @cython.boundscheck(False)
     @cython.wraparound(False)
     def _transform_sequence(
-        self, Py_ssize_t stride, inseq, bint switch,
+        self, Py_ssize_t stride, object inseq, bint switch,
         direction, time_3rd, radians, errcheck
     ):
         if self.projections_exact_same or (self.projections_equivalent and self.skip_equivalent):
@@ -502,31 +493,25 @@ cdef class _Transformer(Base):
         tmp_pj_direction = _PJ_DIRECTION_MAP[TransformDirection.create(direction)]
         cdef PJ_DIRECTION pj_direction = <PJ_DIRECTION>tmp_pj_direction
         # private function to itransform function
-        cdef:
-            void *buffer
-            double *coords
-            double *x
-            double *y
-            double *z
-            double *tt
-            Py_ssize_t buflen, npts, iii, jjj
+        cdef double *x
+        cdef double *y
+        cdef double *z
+        cdef double *tt
 
         if stride < 2:
             raise ProjError("coordinates must contain at least 2 values")
-        if PyObject_AsWriteBuffer(inseq, &buffer, &buflen) <> 0:
-            raise ProjError("object does not provide the python buffer writeable interface")
-
-        coords = <double*>buffer
-        npts = buflen // (stride * _DOUBLESIZE)
         
+        cdef PyBuffWriteManager coordbuff = PyBuffWriteManager(inseq)
+        cdef Py_ssize_t npts, iii, jjj
+        npts = coordbuff.len // stride
         # degrees to radians
         if not self.is_pipeline and not radians\
                 and self._input_radians[pj_direction]:
             with nogil:
                 for iii in range(npts):
                     jjj = stride * iii
-                    coords[jjj] *= _DG2RAD
-                    coords[jjj + 1] *= _DG2RAD
+                    coordbuff.data[jjj] *= _DG2RAD
+                    coordbuff.data[jjj + 1] *= _DG2RAD
         # radians to degrees
         elif not self.is_pipeline and radians\
                 and not self._input_radians[pj_direction]\
@@ -534,26 +519,26 @@ cdef class _Transformer(Base):
             with nogil:
                 for iii in range(npts):
                     jjj = stride * iii
-                    coords[jjj] *= _RAD2DG
-                    coords[jjj + 1] *= _RAD2DG
+                    coordbuff.data[jjj] *= _RAD2DG
+                    coordbuff.data[jjj + 1] *= _RAD2DG
 
         if not switch:
-            x = coords
-            y = coords + 1
+            x = coordbuff.data
+            y = coordbuff.data + 1
         else:
-            x = coords + 1
-            y = coords
+            x = coordbuff.data + 1
+            y = coordbuff.data
 
         # z coordinate
         if stride == 4 or (stride == 3 and not time_3rd):
-            z = coords + 2
+            z = coordbuff.data + 2
         else:
             z = NULL
         # time
         if stride == 3 and time_3rd:
-            tt = coords + 2
+            tt = coordbuff.data + 2
         elif stride == 4:
-            tt = coords + 3
+            tt = coordbuff.data + 3
         else:
             tt = NULL
 
@@ -581,8 +566,8 @@ cdef class _Transformer(Base):
             with nogil:
                 for iii in range(npts):
                     jjj = stride * iii
-                    coords[jjj] *= _RAD2DG
-                    coords[jjj + 1] *= _RAD2DG
+                    coordbuff.data[jjj] *= _RAD2DG
+                    coordbuff.data[jjj + 1] *= _RAD2DG
         # degrees to radians
         elif not self.is_pipeline and radians\
                 and not self._output_radians[pj_direction]\
@@ -590,7 +575,7 @@ cdef class _Transformer(Base):
             with nogil:
                 for iii in range(npts):
                     jjj = stride * iii
-                    coords[jjj] *= _DG2RAD
-                    coords[jjj + 1] *= _DG2RAD
+                    coordbuff.data[jjj] *= _DG2RAD
+                    coordbuff.data[jjj + 1] *= _DG2RAD
 
         ProjError.clear()

pyproj/base.pxi

@@ -9,6 +9,26 @@ cdef extern from "math.h":
         HUGE_VAL
 
 cdef extern from "Python.h":
-    int PyObject_AsWriteBuffer(object, void **rbuf, Py_ssize_t *len)
+    cdef struct PyObject
+    cdef int PyBUF_WRITABLE
+    int PyObject_GetBuffer(PyObject *exporter, Py_buffer *view, int flags)
+    void PyBuffer_Release(Py_buffer *view)
 
 
+cdef class PyBuffWriteManager:
+    cdef Py_buffer buffer
+    cdef double* data
+    cdef public Py_ssize_t len
+
+    def __cinit__(self):
+        self.data = NULL
+
+    def __init__(self, object data):
+        if PyObject_GetBuffer(<PyObject *>data, &self.buffer, PyBUF_WRITABLE) <> 0:
+            raise BufferError("pyproj had a problem getting the buffer from data.")
+        self.data = <double *>self.buffer.buf
+        self.len = self.buffer.len // self.buffer.itemsize
+
+    def __dealloc__(self):
+        PyBuffer_Release(&self.buffer)
+        self.data = NULL