Hi, I'm running into a cryptic error. My CRS in the geopandas dataframe seems to be properly set, however I cannot convert it to any other CRS.

>>> print(gdf_PLU.crs)
epsg:2193
>>>gdf_PLU.to_crs(epsg=4326)
---------------------------------------------------------------------------
RuntimeError                              Traceback (most recent call last)
<ipython-input-92-46b32e6c8012> in <module>()
      1 print(gdf_PLU.crs)
----> 2 gdf_PLU.to_crs(epsg=4326)

~/anaconda3/envs/hthf/lib/python3.6/site-packages/geopandas/geodataframe.py in to_crs(self, crs, epsg, inplace)
    384         else:
    385             df = self.copy()
--> 386         geom = df.geometry.to_crs(crs=crs, epsg=epsg)
    387         df.geometry = geom
    388         df.crs = geom.crs

~/anaconda3/envs/hthf/lib/python3.6/site-packages/geopandas/geoseries.py in to_crs(self, crs, epsg)
    283             except TypeError:
    284                 raise TypeError('Must set either crs or epsg for output.')
--> 285         proj_in = pyproj.Proj(self.crs, preserve_units=True)
    286         proj_out = pyproj.Proj(crs, preserve_units=True)
    287         project = partial(pyproj.transform, proj_in, proj_out)

~/anaconda3/envs/hthf/lib/python3.6/site-packages/pyproj/__init__.py in __new__(self, projparams, preserve_units, **kwargs)
    356         # on case-insensitive filesystems).
    357         projstring = projstring.replace('EPSG','epsg')
--> 358         return _proj.Proj.__new__(self, projstring)
    359 
    360     def __call__(self, *args, **kw):

_proj.pyx in _proj.Proj.__cinit__ (_proj.c:1170)()

RuntimeError: b'no arguments in initialization list'

I can't find this error message in any other issue. From similar issues it seems some data files could be missing (if so, how can I install them; I'm using the binary that is dragged in by conda-forge geopandas). Or is it something else?

Connected to issue #152

This depends on: https://github.com/OSGeo/proj.4/pull/1208 Also, one of the tests fails, but should be fixed with: https://github.com/OSGeo/proj.4/pull/1205

There are a ton of changes - mostly to work with the new API. Feel free to make comments, ask questions, and be thorough in the review.

Dear sir or madam, My name is Ger Smit from the Netherlands. I have problem with pyproj / proj. I use Raspberry Pi, Python, for making screne from data with satpy. No problems til now. I install satpy, try to run python script and get error. I I was in contact with @djhoese about this and he referred me to you ERROR: Minimum supported proj version is 6.2.0, installed version is 4.9.3. For more information see: https://pyproj4.github.io/pyproj/stable/installation.htm

Then when i run the script i get error :

Traceback (most recent call last): File "/home/pi/a-satpy-testen/cursus_seviri.py", line 83, in scn.load(scn.all_dataset_names()[1:]) #HEEL belangrijk File "/home/pi/.local/lib/python3.7/site-packages/satpy/scene.py", line 1163, in load self._read_datasets_from_storage(**kwargs) File "/home/pi/.local/lib/python3.7/site-packages/satpy/scene.py", line 1183, in _read_datasets_from_storage return self._read_dataset_nodes_from_storage(nodes, **kwargs) File "/home/pi/.local/lib/python3.7/site-packages/satpy/scene.py", line 1189, in _read_dataset_nodes_from_storage loaded_datasets = self._load_datasets_by_readers(reader_datasets, **kwargs) File "/home/pi/.local/lib/python3.7/site-packages/satpy/scene.py", line 1214, in _load_datasets_by_readers new_datasets = reader_instance.load(ds_ids, **kwargs) File "/home/pi/.local/lib/python3.7/site-packages/satpy/readers/yaml_reader.py", line 945, in load ds = self._load_dataset_with_area(dsid, coords, **kwargs) File "/home/pi/.local/lib/python3.7/site-packages/satpy/readers/yaml_reader.py", line 1109, in _load_dataset_with_area ds = super(GEOFlippableFileYAMLReader, self)._load_dataset_with_area(dsid, coords, **kwargs) File "/home/pi/.local/lib/python3.7/site-packages/satpy/readers/yaml_reader.py", line 841, in _load_dataset_with_area ds = self._load_dataset_data(file_handlers, dsid, **kwargs) File "/home/pi/.local/lib/python3.7/site-packages/satpy/readers/yaml_reader.py", line 713, in _load_dataset_data proj = self._load_dataset(dsid, ds_info, file_handlers, **kwargs) File "/home/pi/.local/lib/python3.7/site-packages/satpy/readers/yaml_reader.py", line 689, in _load_dataset projectable = fh.get_dataset(dsid, ds_info) File "/home/pi/.local/lib/python3.7/site-packages/satpy/readers/seviri_l1b_native.py", line 458, in get_dataset self._update_attrs(dataset, dataset_info) File "/home/pi/.local/lib/python3.7/site-packages/satpy/readers/seviri_l1b_native.py", line 582, in _update_attrs actual_lon, actual_lat, actual_alt = self.satpos File "/home/pi/.local/lib/python3.7/site-packages/satpy/readers/seviri_l1b_native.py", line 610, in satpos semi_minor_axis=self.mda['projection_parameters']['b'] File "/home/pi/.local/lib/python3.7/site-packages/satpy/readers/seviri_base.py", line 703, in get_satpos proj='geocent', a=semi_major_axis, b=semi_minor_axis, units='m' File "/home/pi/.local/lib/python3.7/site-packages/pyproj/crs.py", line 436, in init super(CRS, self).init(projstring) File "pyproj/_crs.pyx", line 1738, in pyproj._crs._CRS.init pyproj.exceptions.CRSError: Invalid projection: +proj=geocent +a=6378169.0 +b=6356583.800000001 +units=m +type=crs: (Internal Proj Error: proj_create: Invalid b value)

Please can you explaine to me how to install or solved pyproj / proj in combination with satpy. I try a lot withoud result.

Thanks in advange @gersmit

• [x] Closes #411
• [x] Closes #408
• [x] Closes #403
• [x] Closes #426
• [x] Fully documented, including history.rst for all changes and api/*.rst for new API

Code Sample, a copy-pastable example if possible

A "Minimal, Complete and Verifiable Example" will make it much easier for maintainers to help you: http://matthewrocklin.com/blog/work/2018/02/28/minimal-bug-reports

Now that's tricky because I cannot do it. The seg fault very much depends on the input data, and it's part of a complex app. But with some input files 2.3.0 seg fault.

Problem description

Under 2.3.0 I can get the following using python -X faulthandler

Fatal Python error: Segmentation fault

Current thread 0x00007fa0f79c4700 (most recent call first):
  File "/opt/conda/lib/python3.7/site-packages/pyproj/crs.py", line 303 in __init__
  File "/opt/conda/lib/python3.7/site-packages/pyproj/crs.py", line 434 in from_user_input
  File "/opt/conda/lib/python3.7/site-packages/pyproj/proj.py", line 145 in __init__
  File "/opt/conda/lib/python3.7/site-packages/geopandas/geoseries.py", line 304 in to_crs
  File "/opt/conda/lib/python3.7/site-packages/geopandas/geodataframe.py", line 459 in to_crs
  ...
  File "<stdin>", line 1 in <module>
Segmentation fault (core dumped)

I've downgraded to 2.2.2 and I get a Python exception of the same lines, so it might be the same thing

Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  ...
  File "/opt/conda/lib/python3.7/site-packages/geopandas/geodataframe.py", line 459, in to_crs
    geom = df.geometry.to_crs(crs=crs, epsg=epsg)
  File "/opt/conda/lib/python3.7/site-packages/geopandas/geoseries.py", line 304, in to_crs
    proj_in = pyproj.Proj(self.crs, preserve_units=True)
  File "/opt/conda/lib/python3.7/site-packages/pyproj/proj.py", line 147, in __init__
    self.crs = CRS.from_user_input(projparams if projparams is not None else kwargs)
  File "/opt/conda/lib/python3.7/site-packages/pyproj/crs.py", line 435, in from_user_input
    return cls(value)
  File "/opt/conda/lib/python3.7/site-packages/pyproj/crs.py", line 304, in __init__
    super(CRS, self).__init__(projstring)
  File "pyproj/_crs.pyx", line 1308, in pyproj._crs._CRS.__init__
  File "pyproj/_datadir.pyx", line 18, in pyproj._datadir.get_pyproj_context
  File "/opt/conda/lib/python3.7/site-packages/pyproj/datadir.py", line 99, in get_data_dir
    "Valid PROJ data directory not found. "
pyproj.exceptions.DataDirError: Valid PROJ data directory not found. Either set the path using the environmental variable PROJ_LIB or with `pyproj.datadir.set_data_dir`.

This is the same environment as I've detailed here: https://github.com/conda-forge/pyproj-feedstock/issues/47

Environment Information

after downgrade to 2.2.2

System:
    python: 3.7.3 | packaged by conda-forge | (default, Jul  1 2019, 21:52:21)  [GCC 7.3.0]
executable: /usr/bin/condapy
   machine: Linux-4.15.0-1037-gcp-x86_64-with-debian-stretch-sid

PROJ:
      PROJ: 6.1.1
  data dir: None

Python deps:
    pyproj: 2.2.2
       pip: 19.2.3
setuptools: 41.2.0
    Cython: None
     aenum: None

at 2.3.0

System:
    python: 3.7.3 | packaged by conda-forge | (default, Jul  1 2019, 21:52:21)  [GCC 7.3.0]
executable: /usr/bin/condapy
   machine: Linux-4.15.0-1037-gcp-x86_64-with-debian-stretch-sid

PROJ:
      PROJ: 6.1.1
  data dir: /opt/conda/share/proj

Python deps:
    pyproj: 2.3.0
       pip: 19.2.3
setuptools: 41.2.0
    Cython: None

Installation method

Conda in a Ubuntu 16.04 Docker image.

wget --quiet https://repo.continuum.io/miniconda/Miniconda3-latest-Linux-x86_64.sh -O ~/miniconda.sh && \
    /bin/bash ~/miniconda.sh -b -p /opt/conda && \
    conda update conda -y && \
    conda config --add channels conda-forge && \
    conda config --set channel_priority strict && \
    conda install \
        geopandas numexpr bottleneck

Conda environment information (if you installed with conda):

proj4                     6.1.1                hc80f0dc_1    conda-forge
pyproj                    2.3.0            py37h2fd02e8_0    conda-forge

     active environment : None
       user config file : /root/.condarc
 populated config files : /root/.condarc
          conda version : 4.7.11
    conda-build version : not installed
         python version : 3.7.3.final.0
       virtual packages :
       base environment : /opt/conda  (writable)
           channel URLs : https://conda.anaconda.org/conda-forge/linux-64
                          https://conda.anaconda.org/conda-forge/noarch
                          https://repo.anaconda.com/pkgs/main/linux-64
                          https://repo.anaconda.com/pkgs/main/noarch
                          https://repo.anaconda.com/pkgs/r/linux-64
                          https://repo.anaconda.com/pkgs/r/noarch
          package cache : /opt/conda/pkgs
                          /root/.conda/pkgs
       envs directories : /opt/conda/envs
                          /root/.conda/envs
               platform : linux-64
             user-agent : conda/4.7.11 requests/2.22.0 CPython/3.7.3 Linux/4.15.0-1037-gcp ubuntu/16.04.6 glibc/2.23
                UID:GID : 0:0
             netrc file : None
           offline mode : False

pyproj/geod.py - add a parameter return_points to Geod.npts() to determinate the return type; improve return type hint pyproj/geod.py - add a parameters initial_idx, terminus_idx to Geod.npts() pyproj/[geod.py, _geod.pyi, _geod.pyx] - add inv_intermediate_by_npts functions that uses given buffers out_lons, out_lats (and out_azis) buffers instead of creating new buffers pyproj/_geod.pxd - add interface for geod_lineinit pyproj/[geod.py, _geod.pyi, _geod.pyx] - add inv_intermediate_by_del, fwd_intermediate_by_npts, fwd_intermediate_by_del functions, similar to inv_npts but using forward calculation test/test_geod.py - add tests for new functions

Related PR

https://github.com/pyproj4/pyproj/pull/825/

I am attempting to install the latest version of pyproj in order resolve the same issue as described in this closed issue 1113

Installation method/steps

• installed pyproj through conda in terminal, see below.

(base) Richard-Pelgrims-MacBook-Pro:site-packages data_science$ conda install pyproj
Collecting package metadata (current_repodata.json): done
Solving environment: failed with initial frozen solve. Retrying with flexible solve.
Solving environment: failed with repodata from current_repodata.json, will retry with next repodata source.
Collecting package metadata (repodata.json): done
Solving environment: done

## Package Plan ##

  environment location: /opt/anaconda3

  added / updated specs:
    - pyproj


The following NEW packages will be INSTALLED:

  pyproj             pkgs/main/osx-64::pyproj-1.9.6-py37h9c430a6_0

Environment Information

• pyproj version you are attempting to install: 2.6 to resolve issues mentioned in this issue
• PROJ version : Rel. 5.2.0, September 15th, 2018
• Python version: 3.7.4
• Operation System Information: Darwin-15.6.0-x86_64-i386-64bit

Addresses #226

This will need a good review.

I am unsure about what the CF unit mapped to. It is either units or vunits. I am guessing vunits at the moment, but it is only a guess.

Others I think are just another maping to datum, but not sure. Maybe there is a vdatum in PROJ I am missing.

• geoid_name (OGC WKT VERT_DATUM - ex GEOID12B)
• geopotential_datum_name (OGC WKT VERT_DATUM - ex NAVD88)

cc @justb4

Via https://github.com/geopython/GeoHealthCheck/issues/240, looks like with the latest pyproj release on pypi there are install issues:

Downloading/unpacking OWSLib==0.17.0
  Downloading OWSLib-0.17.0.tar.gz (149kB): 149kB downloaded
  Running setup.py (path:/tmp/pip_build_kralidist/OWSLib/setup.py) egg_info for package OWSLib

    warning: manifest_maker: MANIFEST.in, line 3: 'recursive-include' expects <dir> <pattern1> <pattern2> ...

Downloading/unpacking python-dateutil>=1.5 (from OWSLib==0.17.0)
  Downloading python_dateutil-2.8.0-py2.py3-none-any.whl (226kB): 226kB downloaded
Downloading/unpacking pytz (from OWSLib==0.17.0)
  Downloading pytz-2018.9-py2.py3-none-any.whl (510kB): 510kB downloaded
Requirement already satisfied (use --upgrade to upgrade): requests>=1.0 in /usr/lib/python2.7/dist-packages (from OWSLib==0.17.0)
Downloading/unpacking pyproj (from OWSLib==0.17.0)
  Downloading pyproj-2.0.0.tar.gz (408kB): 408kB downloaded
  Running setup.py (path:/tmp/pip_build_kralidist/pyproj/setup.py) egg_info for package pyproj
    ERROR: Cython.Build.cythonize not found. Cython is required to build from a repo.
    Complete output from command python setup.py egg_info:
    ERROR: Cython.Build.cythonize not found. Cython is required to build from a repo.

----------------------------------------
Cleaning up...
Command python setup.py egg_info failed with error code 1 in /.....

Currently don't have CI runners to build/test them. There are options to potentially cross compile, but the wheels wouldn't be able to be tested. Not sure what the best way to handle this is.

Related:

• https://github.com/rasterio/rasterio-wheels/issues/75
• https://github.com/geopandas/pyogrio/pull/64

Hi,

Could you please update the docstring documentation for the Transformer class? One minor typo is proj_from in the from_proj method, but the most important ones are the itransform and transform methods.

They both state that the input should be in the order of x, y when they actually are in the order of y, x. And the output is in the order of y, x. That got me quite confused for a while.

Thanks

Related to #1011

With PROJ 9.2 coming in March 01, 2023 ref that means that 3 minor versions is 9.2, 9.1, 9.0. Currently thinking that pyproj 3.4 will be the last minor version to support PROJ 8.2.

I copied this line from setup.cfg when updating the conda-forge recipe: https://github.com/pyproj4/pyproj/blob/main/setup.cfg#L12 But the actual filename is LICENSE_proj4. These should be made consistent.

https://github.com/pyproj4/pyproj/actions/runs/3596938993/jobs/6058180954

=================================== FAILURES ===================================
___________________ test_transformer_group__network_disabled ___________________

    @pytest.mark.grid
    @patch.dict("os.environ", {"PROJ_NETWORK": "ON"}, clear=True)
    def test_transformer_group__network_disabled():
        with proj_network_env():
            pyproj.network.set_network_enabled(active=False)
>           trans_group = TransformerGroup(4326, 2964)

test/test_transformer.py:869: 
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ 
pyproj/transformer.py:197: in __init__
    super().__init__(
pyproj/_transformer.pyx:215: in pyproj._transformer._TransformerGroup.__init__
    _Transformer._from_pj(
pyproj/_transformer.pyx:563: in pyproj._transformer._Transformer._from_pj
    transformer._init_from_crs(always_xy)
pyproj/_transformer.pyx:614: in pyproj._transformer._Transformer._init_from_crs
    self._initialize_from_projobj()
_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ 

>   raise ProjError("Input is not a transformation.")
E   pyproj.exceptions.ProjError: Input is not a transformation.: (Internal Proj Error: pipeline: Pipeline: Bad step definition: inv (File not found or invalid))

pyproj/_transformer.pyx:324: ProjError
------------------------------ Captured log call -------------------------------
DEBUG    pyproj:transformer.py:197 PROJ_ERROR: hgridshift: could not find required grid(s).
DEBUG    pyproj:transformer.py:197 PROJ_ERROR: pipeline: Pipeline: Bad step definition: inv (File not found or invalid)
____________________ test_transformer_group__download_grids ____________________

get_user_data_dir_mock = <MagicMock name='get_user_data_dir' id='140700317128784'>
tmp_path = PosixPath('/tmp/pytest-of-runner/pytest-2/test_transformer_group__downlo0')
capsys = <_pytest.capture.CaptureFixture object at 0x7ff7586c11f0>

    @pytest.mark.grid
    @pytest.mark.network
    @patch("pyproj.transformer.get_user_data_dir")
    def test_transformer_group__download_grids(get_user_data_dir_mock, tmp_path, capsys):
        get_user_data_dir_mock.return_value = str(tmp_path)
        with proj_network_env():
            pyproj.network.set_network_enabled(active=False)
>           trans_group = TransformerGroup(4326, 2964)
E   pyproj.exceptions.ProjError: Input is not a transformation.: (Internal Proj Error: pipeline: Pipeline: Bad step definition: inv (File not found or invalid))

pyproj/_transformer.pyx:324: ProjError
------------------------------ Captured log call -------------------------------
DEBUG    pyproj:transformer.py:197 PROJ_ERROR: hgridshift: could not find required grid(s).
DEBUG    pyproj:transformer.py:197 PROJ_ERROR: pipeline: Pipeline: Bad step definition: inv (File not found or invalid)
____________ test_coordinate_operation_grids__alternative_grid_name ____________

    @pytest.mark.grid
    def test_coordinate_operation_grids__alternative_grid_name():
        cc = CoordinateOperation.from_epsg(1312, True)
        assert len(cc.grids) == 1
        grid = cc.grids[0]
        assert grid.direct_download is True
        assert grid.open_license is True
        assert grid.short_name == "ca_nrc_ntv1_can.tif"
        assert grid.package_name == ""
        assert grid.url == "https://cdn.proj.org/ca_nrc_ntv1_can.tif"
        if (PROJ_GTE_91 and grids_available(grid.short_name, check_network=False)) or (
            not PROJ_GTE_91 and grids_available(grid.short_name)
        ):
            assert grid.available is True
            assert grid.full_name.endswith(grid.short_name)
        elif PROJ_GTE_911 and pyproj.network.is_network_enabled():
            assert grid.available is True
>           assert grid.full_name == grid.url
E           AssertionError: assert '' == 'https://cdn...._ntv1_can.tif/'
E             - https://cdn.proj.org/ca_nrc_ntv1_can.tif

test/crs/test_crs.py:622: AssertionError
=============================== warnings summary ===============================
test/test_datadir.py::test_get_data_dir__from_user[str]
test/test_datadir.py::test_get_data_dir__from_user[Path]
test/test_datadir.py::test_append_data_dir__internal[str]
test/test_datadir.py::test_append_data_dir__internal[Path]
  /home/runner/work/pyproj/pyproj/pyproj/datadir.py:38: UserWarning: pyproj unable to set database path.
    _global_context_set_data_dir()

-- Docs: https://docs.pytest.org/en/stable/how-to/capture-warnings.html
=========================== short test summary info ============================
FAILED test/test_transformer.py::test_transformer_group__network_disabled - pyproj.exceptions.ProjError: Input is not a transformation.: (Internal Proj Error: pipeline: Pipeline: Bad step definition: inv (File not found or invalid))
FAILED test/test_transformer.py::test_transformer_group__download_grids - pyproj.exceptions.ProjError: Input is not a transformation.: (Internal Proj Error: pipeline: Pipeline: Bad step definition: inv (File not found or invalid))
FAILED test/test_transformer.py::test_transformer_group__download_grids__directory - pyproj.exceptions.ProjError: Input is not a transformation.: (Internal Proj Error: pipeline: Pipeline: Bad step definition: inv (File not found or invalid))
FAILED test/crs/test_crs.py::test_coordinate_operation_grids__alternative_grid_name - AssertionError: assert '' == 'https://cdn...._ntv1_can.tif/'
  - https://cdn.proj.org/ca_nrc_ntv1_can.tif
============ 4 failed, 832 passed, 5 skipped, 4 warnings in 22.53s =============

Going to need to do a git bisect with PROJ.

It would be handy to create a Geod object from an Ellipsoid object. E.g., this is how I would expect:

from pyproj import CRS, Geod

crs = CRS.from_epsg(4326)
geod = Geod.from_ellipsoid(crs.ellipsoid)

or is that too convoluted? Is there a simpler way? A few current working methods:

Geod(ellps='WGS84')
# Geod(ellps='WGS84')
Geod(a=crs.ellipsoid.semi_major_metre, b=crs.ellipsoid.semi_minor_metre)
# Geod('+a=6378137 +f=0.0033528106647475126')
geod = Geod(a=crs.ellipsoid.semi_major_metre, f=1/crs.ellipsoid.inverse_flattening)
# Geod(ellps='WGS84')

and a bonus curiosity, probably a floating-point precision error

assert Geod(ellps='WGS84') == Geod(a=crs.ellipsoid.semi_major_metre, f=1/crs.ellipsoid.inverse_flattening)
assert Geod(ellps='WGS84') != Geod(a=crs.ellipsoid.semi_major_metre, b=crs.ellipsoid.semi_minor_metre)

from_proj was added based on pyproj 1 patterns of using the Proj objects with the transformer. However, with PROJ 6+ and pyproj 2+, the CRS object should be the object used when creating transformers.

from_proj does not have all of the kwargs that from_crs has and that was mainly because from_crs should be the method used and there isn't much point in supporting both. Additionally, the Proj object has some additional overhead as it creates a Transformer and a CRS when initialized.

I am thinking that removing from_proj will remove confusion.

Compatibility: from_proj was added in the initial version of the transformer in pyproj 2.1.0. from_crs was added at the same time, so if users modify their code to use from_crs instead it will also be backwards compatible.

Does anyone have good reasons why from_proj should stay?

Problem description

The +R_A option is not being honored in Pyproj, whereas it seems to work fine using the geod CLI tool.

This problem happens if I use pyproj.Geod() or pyproj.CRS().get_geod().

Code sample

My Python code:

import logging
import os

os.environ['PROJ_DEBUG'] = '3'
logging.basicConfig(level=logging.DEBUG)

import pyproj
# DEBUG:pyproj:PROJ_DEBUG: pj_open_lib(proj.ini): call fopen(/.../lib/python3.10/site-packages/pyproj/proj_dir/share/proj/proj.ini) - succeeded

geod1 = pyproj.Geod('+proj=lonlat +ellps=WGS84')
print(geod1.f)
# 0.0033528106647474805
geod1.fwd(-70, 40, 60, 10_000)
# (-69.89851793919364, 40.04498640681145, -119.93473805953444)

geod2 = pyproj.Geod('+proj=lonlat +ellps=WGS84 +R_A')
print(geod2.f)
# 0.0033528106647474805
geod2.fwd(-70, 40, 60, 10_000)
# (-69.89851793919364, 40.04498640681145, -119.93473805953444)

Equivalent CLI invocation using geod:

export PROJ_DEBUG=3

geod -f '%0.4f' +proj=lonlat +ellps=WGS84 <<< '-70E 40N 60d 10000m'
# pj_ellipsoid - final: a=6378137.000 f=1/298.257, errno=0
# pj_ellipsoid - final:    ellps=WGS84
# 70.0447	40.2273	-119.7864m

geod -f '%0.4f' +proj=lonlat +ellps=WGS84 +R_A <<< '-70E 40N 60d 10000m'
# pj_ellipsoid - final: a=6371007.181 f=1/  0.000, errno=0
# pj_ellipsoid - final:   R_A ellps=WGS84
# 70.0448	40.2282	-119.7855m

Expected Output

I expected the output from geod1 to be different from geod2, and that geod2.f should be 0.

I was also surprised to see that the output from the geod command line tool was different from the output from pyproj.Geod.fwd. Presumably the difference can be ascribed to different options when compiling the CLI tool (installed via MacPorts) and the Python package (installed via the wheels on PyPI).

I further expected to see "trace" logging output comparable to what the geod CLI tool emitted, but I suppose that's a separate issue.

Environment Information

MacOS 12.4 on ARM64.

Versions

pyproj info:
    pyproj: 3.4.0
      PROJ: 9.1.0
  data dir: /.../lib/python3.10/site-packages/pyproj/proj_dir/share/proj
user_data_dir: /.../.local/share/proj
PROJ DATA (recommended version): 1.11
PROJ Database: 1.2
EPSG Database: v10.074 [2022-08-01]
ESRI Database: ArcGIS Pro 3.0 [2022-07-09]
IGNF Database: 3.1.0 [2019-05-24]

System:
    python: 3.10.6 (main, Sep 19 2022, 14:40:49) [Clang 13.1.6 (clang-1316.0.21.2.5)]
executable: /.../bin/python
   machine: macOS-12.4-arm64-arm-64bit

Python deps:
   certifi: 2022.9.14
    Cython: None
setuptools: 65.4.1
       pip: 22.2.2

Installation method

Installed using Pip in a Venv:

python -m venv venv
./venv/bin/pip install -U pyproj