rewrite references using BibTeX and sphinxcontrib-bibtex

15 files changed, 300 insertions, 157 deletions

diff --git a/docs/source/conf.py b/docs/source/conf.py
index b3c5bfea..81a683d6 100644
--- a/docs/source/conf.py
+++ b/docs/source/conf.py
@@ -31,6 +31,7 @@ import datetime
 # ones.
 extensions = [
     'sphinx.ext.mathjax',
+    'sphinxcontrib.bibtex',
 ]
 
 # Add any paths that contain templates here, relative to this directory.
diff --git a/docs/source/geodesic.rst b/docs/source/geodesic.rst
index f677ae1d..c7753404 100644
--- a/docs/source/geodesic.rst
+++ b/docs/source/geodesic.rst
@@ -152,10 +152,11 @@ catalog of those cases:
 Background
 ----------
 
-The algorithms implemented by this package are given in [Karney2013]_
-and are based on [Bessel1825]_ and [Helmert1880]_; the algorithm for
-areas is based on [Danielsen1989]_.  These improve on the work of
-[Vincenty1975]_ in the following respects:
+The algorithms implemented by this package are given in :cite:`Karney2013`
+(`addenda <https://geographiclib.sourceforge.io/geod-addenda.html>`_)
+and are based on :cite:`Bessel1825` and :cite:`Helmert1880`; the algorithm for
+areas is based on :cite:`Danielsen1989`.  These improve on the work of
+:cite:`Vincenty1975` in the following respects:
 
 * The results are accurate to round-off for terrestrial ellipsoids (the
   error in the distance is less then 15 nanometers, compared to 0.1 mm
@@ -166,5 +167,8 @@ areas is based on [Danielsen1989]_.  These improve on the work of
   geodesic.  This allows, for example, the area of a geodesic polygon to
   be computed.
 
-Additional background material is provided in [GeodesicBib]_,
-[GeodesicWiki]_, and [Karney2011]_.
+Additional background material is provided in GeographicLib's `geodesic
+bibliography <https://geographiclib.sourceforge.io/geodesic-papers/biblio.html>`_,
+Wikipedia's article "`Geodesics on an ellipsoid
+<https://en.wikipedia.org/wiki/Geodesics_on_an_ellipsoid>`_", and :cite:`Karney2011`
+(`errata <https://geographiclib.sourceforge.io/geod-addenda.html#geod-errata>`_).
diff --git a/docs/source/operations/projections/calcofi.rst b/docs/source/operations/projections/calcofi.rst
index 8a240ebd..8b2924de 100644
--- a/docs/source/operations/projections/calcofi.rst
+++ b/docs/source/operations/projections/calcofi.rst
@@ -82,8 +82,8 @@ Options
 Mathematical definition
 ################################################################################
 
-The algorithm used to make conversions is described in [EberHewitt1979]_ with
-a few corrections reported in [WeberMoore2013]_.
+The algorithm used to make conversions is described in :cite:`EberHewitt1979` with
+a few corrections reported in :cite:`WeberMoore2013`.
 
 
 Further reading
diff --git a/docs/source/operations/projections/cass.rst b/docs/source/operations/projections/cass.rst
index dba6a7c4..2b2ebf01 100644
--- a/docs/source/operations/projections/cass.rst
+++ b/docs/source/operations/projections/cass.rst
@@ -30,7 +30,7 @@ Although the Cassini projection has been largely replaced by the Transverse Merc
 Usage
 #####
 
-There has been little usage of the spherical version of the Cassini, but the ellipsoidal Cassini-Soldner version was adopted by the Ordnance Survey for the official survey of Great Britain during the second half of the 19th century [Steers1970]_.
+There has been little usage of the spherical version of the Cassini, but the ellipsoidal Cassini-Soldner version was adopted by the Ordnance Survey for the official survey of Great Britain during the second half of the 19th century :cite:`Steers1970`.
 Many of these maps were prepared at a scale of 1:2,500.
 The Cassini-Soldner was also used for the detailed mapping of many German states during the same period.
 
@@ -67,7 +67,7 @@ Options
 Mathematical definition
 #######################
 
-The formulas describing the Cassini projection are taken from Snyder's [Snyder1987]_.
+The formulas describing the Cassini projection are taken from :cite:`Snyder1987`.
 
 :math:`\phi_0` is the latitude of origin that match the center of the map (default to 0). It can be set with ``+lat_0``.
 
@@ -162,5 +162,4 @@ Further reading
 ###############
 
 #. `Wikipedia <https://en.wikipedia.org/wiki/Equirectangular_projection>`_
-#. [Snyder1987]_
 #. `EPSG, POSC literature pertaining to Coordinate Conversions and Transformations including Formulas  <http://www.ihsenergy.com/epsg/guid7.pdf>`_
diff --git a/docs/source/operations/projections/ccon.rst b/docs/source/operations/projections/ccon.rst
index f5ceed67..e267f271 100644
--- a/docs/source/operations/projections/ccon.rst
+++ b/docs/source/operations/projections/ccon.rst
@@ -39,7 +39,7 @@ conformal.
 
 An example of usage (and the main reason to implement this projection in proj4)
 is the ATPOL geobotanical grid of Poland, developed in Institute of Botany,
-Jagiellonian University, Krakow, Poland in 1970s [Zajac1978]_. The grid was
+Jagiellonian University, Krakow, Poland in 1970s :cite:`Zajac1978`. The grid was
 originally handwritten on paper maps and further copied by hand. The projection
 (together with strange Earth radius) was chosen by its creators as the compromise
 fit to existing maps during first software development in DOS era. Many years later
@@ -51,7 +51,7 @@ The ATPOL coordinates can be achieved with with the following parameters:
 
    +proj=ccon +lat_1=52 +lon_0=19 +axis=esu +a=6390000 +x_0=330000 +y_0=-350000
 
-For more information see [Komsta2016]_ and [Verey2017]_.
+For more information see :cite:`Komsta2016` and :cite:`Verey2017`.
 
 Parameters
 ################################################################################
diff --git a/docs/source/operations/projections/eqc.rst b/docs/source/operations/projections/eqc.rst
index de463034..b81f9ded 100644
--- a/docs/source/operations/projections/eqc.rst
+++ b/docs/source/operations/projections/eqc.rst
@@ -89,7 +89,7 @@ Parameters
 Mathematical definition
 #######################
 
-The formulas describing the Equidistant Cylindrical projection are all taken from Snyder's [Snyder1987]_.
+The formulas describing the Equidistant Cylindrical projection are all taken from :cite:`Snyder1987`.
 
 :math:`\phi_{ts}` is the latitude of true scale, that mean the standard parallels where the scale of the projection is true. It can be set with ``+lat_ts``.
 
diff --git a/docs/source/operations/projections/gall.rst b/docs/source/operations/projections/gall.rst
index bdb88895..1b0da9e8 100644
--- a/docs/source/operations/projections/gall.rst
+++ b/docs/source/operations/projections/gall.rst
@@ -69,7 +69,7 @@ Parameters
 Mathematical definition
 #######################
 
-The formulas describing the Gall Stereographical are all taken from Snyder's [Snyder1993]_.
+The formulas describing the Gall Stereographical are all taken from :cite:`Snyder1993`.
 
 Spherical form
 **************
diff --git a/docs/source/operations/projections/merc.rst b/docs/source/operations/projections/merc.rst
index 2d419c02..349de216 100644
--- a/docs/source/operations/projections/merc.rst
+++ b/docs/source/operations/projections/merc.rst
@@ -37,7 +37,7 @@ Applications should be limited to equatorial regions, but is frequently
 used for navigational charts with latitude of true scale (:option:`+lat_ts`) specified within
 or near chart's boundaries.
 Often inappropriately used for world maps since the regions near the poles
-cannot be shown [Evenden1995]_.
+cannot be shown :cite:`Evenden1995`.
 
 
 Example using latitude of true scale::
@@ -76,7 +76,7 @@ Parameters
 Mathematical definition
 #######################
 
-The formulas describing the Mercator projection are all taken from G. Evenden's libproj manuals [Evenden2005]_.
+The formulas describing the Mercator projection are all taken from G. Evenden's libproj manuals :cite:`Evenden2005`.
 
 Spherical form
 **************
diff --git a/docs/source/operations/projections/mill.rst b/docs/source/operations/projections/mill.rst
index 718d4b9b..fb2096c9 100644
--- a/docs/source/operations/projections/mill.rst
+++ b/docs/source/operations/projections/mill.rst
@@ -32,7 +32,7 @@ Usage
 ########
 
 The Miller Cylindrical projection is used for world maps and in several atlases,
-including the National Atlas of the United States (USGS, 1970, p. 330-331) [Snyder1987]_.
+including the National Atlas of the United States (USGS, 1970, p. 330-331) :cite:`Snyder1987`.
 
 Example using Central meridian 90°W::
 
@@ -55,7 +55,7 @@ Parameters
 Mathematical definition
 #######################
 
-The formulas describing the Miller projection are all taken from Snyder's manuals [Snyder1987]_.
+The formulas describing the Miller projection are all taken from :cite:`Snyder1987`.
 
 
 Forward projection
diff --git a/docs/source/operations/projections/qsc.rst b/docs/source/operations/projections/qsc.rst
index 56d4ea3f..e6b0d1b0 100644
--- a/docs/source/operations/projections/qsc.rst
+++ b/docs/source/operations/projections/qsc.rst
@@ -36,15 +36,15 @@ It is an equal-area projection, and at the same time introduces only limited ang
 distortions. It treats all cube sides equally, i.e. it does not use different
 projections for polar areas and equatorial areas. These properties make QSC
 projection a good choice for planetary-scale terrain rendering. Map data can be
-organized in quadtree structures for each cube side. See [LambersKolb2012]_ for an example.
+organized in quadtree structures for each cube side. See :cite:`LambersKolb2012` for an example.
 
-The QSC projection was introduced by [ONeilLaubscher1976]_,
-building on previous work by [ChanONeil1975]_. For clarity: The
-earlier QSC variant described in [ChanONeil1975]_ became known as the COBE QSC since it
+The QSC projection was introduced by :cite:`ONeilLaubscher1976`,
+building on previous work by :cite:`ChanONeil1975`. For clarity: The
+earlier QSC variant described in :cite:`ChanONeil1975` became known as the COBE QSC since it
 was used by the NASA Cosmic Background Explorer (COBE) project; it is an approximately
 equal-area projection and is not the same as the QSC projection.
 
-See also [CalabrettaGreisen2002]_ Sec. 5.6.2 and 5.6.3 for a description of both and
+See also :cite:`CalabrettaGreisen2002` Sec. 5.6.2 and 5.6.3 for a description of both and
 some analysis.
 
 In this implementation, the QSC projection projects onto one side of a circumscribed
@@ -65,7 +65,7 @@ cube. The cube side is selected by choosing one of the following six projection
 +-------------------------+--------------------+
 
 Furthermore, this implementation allows the projection to be applied to ellipsoids.
-A preceding shift to a sphere is performed automatically; see [LambersKolb2012]_ for details.
+A preceding shift to a sphere is performed automatically; see :cite:`LambersKolb2012` for details.
 
 
 Usage
diff --git a/docs/source/operations/transformations/deformation.rst b/docs/source/operations/transformations/deformation.rst
index 2e6f012d..13614f5c 100644
--- a/docs/source/operations/transformations/deformation.rst
+++ b/docs/source/operations/transformations/deformation.rst
@@ -38,7 +38,7 @@ writes both file formats. Using GDAL for construction of new grids is recommende
 Example
 -------------------------------------------------------------------------------
 
-In [Häkli2016]_ coordinate transformation including a deformation model is described.
+In :cite:`Hakli2016` coordinate transformation including a deformation model is described.
 The paper describes how coordinates from the global ITRFxx frames are transformed to the
 local Nordic realisations of ETRS89. Scandinavia is an area with significant post-glacial
 rebound. The deformations from the post-glacial uplift is not accounted for in the
diff --git a/docs/source/operations/transformations/molodensky.rst b/docs/source/operations/transformations/molodensky.rst
index a90828cd..347165c8 100644
--- a/docs/source/operations/transformations/molodensky.rst
+++ b/docs/source/operations/transformations/molodensky.rst
@@ -16,7 +16,7 @@ differences between the semimajor axes and flattening parameters of the referenc
 ellipsoids. Due to its algorithmic simplicity, it was popular prior to the
 ubiquity of digital computers. Today, it is mostly interesting for historical
 reasons, but nevertheless indispensable due to the large amount of data that has
-already been transformed that way [EversKnudsen2017]_.
+already been transformed that way :cite:`EversKnudsen2017`.
 
 +---------------------+----------------------------------------------------------+
 | **Alias**           | molodensky                                               |
@@ -45,7 +45,7 @@ is that it is fairly simple to compute by hand. The ease of computation come at
 cost of limited accuracy.
 
 A derivation of the mathematical formulas for the Molodensky transform can be found
-in [Deakin2004]_.
+in :cite:`Deakin2004`.
 
 
 Examples
diff --git a/docs/source/references.bib b/docs/source/references.bib
new file mode 100644
index 00000000..31f9a398
--- /dev/null
+++ b/docs/source/references.bib
@@ -0,0 +1,264 @@
+% This file was created with JabRef 2.10.
+% Encoding: UTF-8
+
+
+@Article{Altamimi2002,
+  Title                    = {{ITRF2000}: A new release of the {I}nternational {T}errestrial {R}eference {F}rame for earth science applications},
+  Author                   = {Z. Altamimi and P. Sillard and C. Boucher},
+  Journal                  = {Journal of Geophysical Research: Solid Earth},
+  Year                     = {2002},
+  Number                   = {B10},
+  Volume                   = {107},
+
+  Doi                      = {10.1029/2001JB000561}
+}
+
+@Article{Bessel1825,
+  Title                    = {The calculation of longitude and latitude from geodesic measurements},
+  Author                   = {F. W. Bessel},
+  Journal                  = {Astronomische Nachrichten},
+  Year                     = {1825},
+  Number                   = {86},
+  Pages                    = {241--254},
+  Volume                   = {4},
+
+  Url                      = {https://arxiv.org/abs/0908.1824}
+}
+
+@Article{CalabrettaGreisen2002,
+  Title                    = {Representations of celestial coordinates in {FITS}},
+  Author                   = {Calabretta, M. R. and Greisen, E. W.},
+  Journal                  = {Astronomy \& Astrophysics},
+  Year                     = {2002},
+  Number                   = {3},
+  Pages                    = {1077--1122},
+  Volume                   = {395},
+
+  Doi                      = {10.1051/0004-6361:20021327}
+}
+
+@TechReport{ChanONeil1975,
+  Title                    = {Feasibility Study of a Quadrilateralized Spherical Cube Earth Data Base},
+  Author                   = {Chan, F. K. and O'Neill, E. M.},
+  Institution              = {Computer Sciences Corporation, System Sciences Division},
+  Year                     = {1975},
+
+  Address                  = {Silver Spring, Md},
+  Number                   = {2-75 (CSC)},
+  Type                     = {Tech. Rep. EPRF},
+
+  Url                      = {https://archive.org/details/ADA010232}
+}
+
+@Article{Danielsen1989,
+  Title                    = {The area under the geodesic},
+  Author                   = {J. Danielsen},
+  Journal                  = {Survey Review},
+  Year                     = {1989},
+  Number                   = {232},
+  Pages                    = {61--66},
+  Volume                   = {30},
+
+  Doi                      = {10.1179/sre.1989.30.232.61}
+}
+
+@TechReport{Deakin2004,
+  Title                    = {The Standard and Abridged {Molodensky} Coordinate Transformation Formulae},
+  Author                   = {R. E. Deakin},
+  Institution              = {Department of Mathematical and Geospatial Sciences, RMIT University},
+  Year                     = {2004},
+
+  Url                      = {http://www.mygeodesy.id.au/documents/Molodensky%20V2.pdf}
+}
+
+@Article{EberHewitt1979,
+  Title                    = {Conversion algorithms for the {CalCOFI} station grid},
+  Author                   = {L. E. Eber and R.P. Hewitt},
+  Journal                  = {California Cooperative Oceanic Fisheries Investigations Reports},
+  Year                     = {1979},
+  Pages                    = {135--137},
+  Volume                   = {20},
+
+  Url                      = {http://www.calcofi.org/publications/calcofireports/v20/Vol_20_Eber___Hewitt.pdf}
+}
+
+@Manual{Evenden2005,
+  Title                    = {libproj4: A Comprehensive Library of Cartographic Projection Functions (Preliminary Draft)},
+  Author                   = {G. I. Evenden},
+  Year                     = {2005},
+
+  Url                      = {https://github.com/OSGeo/proj.4/blob/master/docs/old/libproj.pdf}
+}
+
+@Manual{Evenden1995,
+  Title                    = {Cartograpic Projection Procedures for the {UNIX} Environment --- A User's Manual},
+  Author                   = {G. I. Evenden},
+  Year                     = {1995},
+
+  Url                      = {https://github.com/OSGeo/proj.4/blob/master/docs/old/of90-284.pdf}
+}
+
+@InProceedings{EversKnudsen2017,
+  Title                    = {Transformation pipelines for {PROJ.4}},
+  Author                   = {K. Evers and T. Knudsen},
+  Booktitle                = {FIG Working Week 2017 Proceedings},
+  Year                     = {2017},
+
+  Url                      = {http://www.fig.net/resources/proceedings/fig_proceedings/fig2017/papers/iss6b/ISS6B_evers_knudsen_9156.pdf}
+}
+
+@Article{Hakli2016,
+  Title                    = {The {NKG2008} {GPS} campaign -- final transformation results and a new common {Nordic} reference frame},
+  Author                   = {P. Häkli and M. Lidberg and L. Jivall and T. Nørbech and O. Tangen and M. Weber and P. Pihlak and I. Aleksejenko
+ and E. Paršeliunas},
+  Journal                  = {Journal of Geodetic Science},
+  Year                     = {2016},
+  Number                   = {1},
+  Pages                    = {1--33},
+  Volume                   = {6},
+
+  Doi                      = {10.1515/jogs-2016-0001}
+}
+
+@Book{Helmert1880,
+  Title                    = {Mathematical and Physical Theories of Higher Geodesy},
+  Author                   = {F. R. Helmert},
+  Publisher                = {Teubner},
+  Year                     = {1880},
+
+  Address                  = {Leipzig},
+  Volume                   = {1},
+
+  Doi                      = {10.5281/zenodo.32050}
+}
+
+@Article{Karney2013,
+  Title                    = {Algorithms for geodesics},
+  Author                   = {C. F. F. Karney},
+  Journal                  = {Journal of Geodesy},
+  Year                     = {2013},
+  Number                   = {1},
+  Pages                    = {43--55},
+  Volume                   = {87},
+
+  Doi                      = {10.1007/s00190-012-0578-z}
+}
+
+@Article{Karney2011,
+  Title                    = {Geodesics on an ellipsoid of revolution},
+  Author                   = {C. F. F. Karney},
+  Journal                  = {ArXiv e-prints},
+  Year                     = {2011},
+
+  Archiveprefix            = {arXiv},
+  Eprint                   = {1102.1215},
+  Primaryclass             = {physics.geo-ph}
+}
+
+@Article{Komsta2016,
+  Title                    = {ATPOL geobotanical grid revisited -- a proposal of coordinate conversion algorithms},
+  Author                   = {Ł. Komsta},
+  Journal                  = {Annales UMCS Sectio E Agricultura},
+  Year                     = {2016},
+  Number                   = {1},
+  Pages                    = {31--37},
+  Volume                   = {71}
+}
+
+@InProceedings{LambersKolb2012,
+  Title                    = {Ellipsoidal Cube Maps for Accurate Rendering of Planetary-Scale Terrain Data},
+  Author                   = {M. Lambers and A. Kolb},
+  Booktitle                = {Pacific Graphics Short Papers},
+  Year                     = {2012},
+  Editor                   = {C. Bregler and P. Sander and M. Wimmer},
+  Publisher                = {The Eurographics Association},
+
+  Doi                      = {10.2312/PE/PG/PG2012short/005-010}
+}
+
+@TechReport{ONeilLaubscher1976,
+  Title                    = {Extended Studies of a Quadrilateralized Spherical Cube Earth Data Base},
+  Author                   = {E. M. O'Neill and R. E. Laubscher},
+  Institution              = {Computer Sciences Corporation, System Sciences Division},
+  Year                     = {1976},
+
+  Address                  = {Silver Spring, Md},
+  Number                   = {3-76 (CSC)},
+  Type                     = {Tech. Rep. EPRF},
+
+  Url                      = {https://archive.org/details/DTIC_ADA026294}
+}
+
+@Book{Snyder1993,
+  Title                    = {Flattening the Earth},
+  Author                   = {J. P. Snyder},
+  Publisher                = {University of Chicago Press},
+  Year                     = {1993}
+}
+
+@TechReport{Snyder1987,
+  Title                    = {Map Projections --- A Working Manual},
+  Author                   = {J. P. Snyder},
+  Institution              = {U.S. Geological Survey},
+  Year                     = {1987},
+  Number                   = {1395},
+  Type                     = {Professional Paper},
+
+  Url                      = {https://pubs.er.usgs.gov/publication/pp1395}
+}
+
+@Book{Steers1970,
+  Title                    = {An introduction to the study of map projections},
+  Author                   = {J. A. Steers},
+  Publisher                = {University of London Press},
+  Year                     = {1970},
+  Edition                  = {15th}
+}
+
+@Article{Verey2017,
+  Title                    = {Theoretical analysis and practical consequences of adopting a model {ATPOL} grid as a conical projection defining the conversion of plane coordinates to the {WGS} 84 ellipsoid},
+  Author                   = {M. Verey},
+  Journal                  = {Fragmenta Floristica et Geobotanica Polonica},
+  Year                     = {2017},
+  Number                   = {2},
+  Pages                    = {469--488},
+  Volume                   = {24},
+
+  Url                      = {http://bomax.botany.pl/pubs-new/#article-4279}
+}
+
+@Article{Vincenty1975,
+  Title                    = {Direct and inverse solutions of geodesics on the ellipsoid with application of nested equations},
+  Author                   = {T. Vincenty},
+  Journal                  = {Survey Review},
+  Year                     = {1975},
+  Number                   = {176},
+  Pages                    = {88--93},
+  Volume                   = {23},
+
+  Doi                      = {10.1179/sre.1975.23.176.88}
+}
+
+@Article{WeberMoore2013,
+  Title                    = {Corrected Conversion Algorithms For The {CalCOFI} Station Grid And Their Implementation In Several Computer Languages},
+  Author                   = {E. D. Weber and T. J. Moore},
+  Journal                  = {California Cooperative Oceanic Fisheries Investigations Reports},
+  Year                     = {2013},
+  Pages                    = {1--10},
+  Volume                   = {54},
+
+  Url                      = {http://calcofi.org/publications/calcofireports/v54/Vol_54_Weber.pdf}
+}
+
+@Article{Zajac1978,
+  Title                    = {Atlas of Distribution of Vascular Plants in {Poland} ({ATPOL})},
+  Author                   = {A. Zajaç},
+  Journal                  = {Taxon},
+  Year                     = {1978},
+  Number                   = {5/6},
+  Pages                    = {481--484},
+  Volume                   = {27},
+
+  Doi                      = {10.2307/1219899}
+}
+
diff --git a/docs/source/references.rst b/docs/source/references.rst
index 807e2faa..a76d6526 100644
--- a/docs/source/references.rst
+++ b/docs/source/references.rst
@@ -6,130 +6,5 @@
     References
     ================================================================================
 
-.. [AltamimiEtAl2002] Altamimi, Z., P. Sillard, and C. Boucher (2002),
-    ITRF2000: A new release of the International Terrestrial Reference Frame for earth science applications,
-    J. Geophys. Res., 107(B10), 2214,
-    `doi:10.1029/2001JB000561 <http://www.agu.org/pubs/crossref/2002/2001JB000561.shtml>`__.
-
-
-.. [Bessel1825] F. W. Bessel, 1825,
-   `The calculation of longitude and latitude from geodesic measurements
-   <https://arxiv.org/abs/0908.1824>`_,
-   Astron. Nachr. **331**\ (8), 852–861 (2010),
-   translated by C. F. F. Karney and R. E. Deakin.
-
-.. [CalabrettaGreisen2002] M. Calabretta and E. Greisen, 2002,
-   "Representations of celestial coordinates in FITS".
-   Astronomy & Astrophysics 395, 3, 1077–1122.
-
-.. [ChanONeil1975]  F. Chan and E. M. O'Neill, 1975,
-   "Feasibility Study of a Quadrilateralized Spherical Cube Earth Data Base",
-   Tech. Rep. EPRF 2-75 (CSC), Environmental Prediction Research Facility.
-
-.. [Danielsen1989] J. Danielsen, 1989,
-   `The area under the geodesic
-   <https://doi.org/10.1179/003962689791474267>`_,
-   Survey Review **30**\ (232), 61–66.
-
-.. [Deakin2004] R.E. Deakin, 2004,
-    `The Standard and Abridged Molodensky Coordinate Transformation Formulae
-    <http://www.mygeodesy.id.au/documents/Molodensky%20V2.pdf>`__.
-
-.. [EberHewitt1979] L. E. Eber and R.P. Hewitt, 1979,
-   `Conversion algorithms for the CalCOFI station grid
-   <http://www.calcofi.org/publications/calcofireports/v20/Vol_20_Eber___Hewitt.pdf>`_,
-   California Cooperative Oceanic Fisheries Investigations Reports 20:135-137.
-
-.. [Evenden1995] G. I. Evenden, 1995,
-   `Cartograpic Projection Procedures for the UNIX Environment -
-   A User's Manual
-   <https://github.com/OSGeo/proj.4/blob/master/docs/old/proj_4_3_12.pdf>`_.
-
-.. [Evenden2005] G. I. Evenden, 2005,
-   `libproj4: A Comprehensive Library of Cartographic Projection Functions
-   (Preliminary Draft)
-   <https://github.com/OSGeo/proj.4/blob/master/docs/old/libproj.pdf>`_.
-
-.. [EversKnudsen2017] K. Evers and T. Knudsen, 2017,
-    `Transformation pipelines for PROJ.4
-    <http://www.fig.net/resources/proceedings/fig_proceedings/fig2017/papers/iss6b/ISS6B_evers_knudsen_9156.pdf>`__,
-    FIG Working Week 2017 Proceedings.
-
-.. [GeodesicBib] `A geodesic bibliography
-   <https://geographiclib.sourceforge.io/geodesic-papers/biblio.html>`_.
-
-.. [GeodesicWiki] The wikipedia page,
-   `Geodesics on an ellipsoid
-   <https://en.wikipedia.org/wiki/Geodesics_on_an_ellipsoid>`_.
-
-.. [Häkli2016] P. Häkli, M. Lidberg, L. Jivall, et al, 2016,
-   `The NKG2008 GPS Campaign - final transformation results and a new common Nordic reference frame
-   <https://doi.org/10.1515/jogs-2016-0001>`_,
-   Journal of Geodetic Science, 6(1).
-
-.. [Helmert1880] F. R. Helmert, 1880,
-   `Mathematical and Physical Theories of Higher Geodesy, Vol 1
-   <https://doi.org/10.5281/zenodo.32050>`_,
-   (Teubner, Leipzig), Chaps. 5–7.
-
-.. [Karney2011] C. F. F. Karney, 2011,
-   `Geodesics on an ellipsoid of revolution
-   <https://arxiv.org/abs/1102.1215v1>`_;
-   `errata
-   <https://geographiclib.sourceforge.io/geod-addenda.html#geod-errata>`_.
-
-.. [Karney2013] C. F. F. Karney, 2013,
-   `Algorithms for geodesics
-   <https://doi.org/10.1007/s00190-012-0578-z>`_,
-   J. Geodesy **87**\ (1) 43–55;
-   `addenda <https://geographiclib.sourceforge.io/geod-addenda.html>`_.
-
-.. [Komsta2016] L. Komsta, 2016,
-   `ATPOL geobotanical grid revisited - a proposal of coordinate conversion
-   algorithms
-   <http://wydawnictwo.up.lublin.pl/annales/Agricultura/2016/1/03.pdf>`_,
-   Annales UMCS Sectio E Agricultura 71(1), 31–37.
-
-.. [LambersKolb2012] M. Lambers and A. Kolb, 2012,
-   "Ellipsoidal Cube Maps for Accurate Rendering of Planetary-Scale
-   Terrain Data", Proc. Pacific Graphics (Short Papers).
-
-.. [ONeilLaubscher1976] E. M. O'Neill and R. E. Laubscher, 1976,
-   "Extended Studies of a Quadrilateralized Spherical Cube Earth Data Base",
-   Tech. Rep. NEPRF 3-76 (CSC),
-   Naval Environmental Prediction Research Facility.
-
-.. [Snyder1987] J. P. Snyder, 1987,
-   `Map Projections - A Working Manual
-   <https://pubs.er.usgs.gov/publication/pp1395>`_.
-   U.S. Geological Survey professional paper 1395.
-
-.. [Snyder1993] J. P. Snyder, 1993,
-   Flattening the Earth, Chicago and London, The University of Chicago press.
-
-.. [Steers1970] J. A. Steers, 1970,
-   An introduction to the study of map projections (15th ed.),
-   London Univ. Press, p. 229.
-
-.. [Verey2017] M. Verey, 2017,
-   `Theoretical analysis and practical consequences of adopting an ATPOL
-   grid model as a conical projection, defining the conversion of plane
-   coordinates to the WGS-84 ellipsoid
-   <http://www.botany.pl/atpol/Siatka%20ATPOL%20w%20analitycznym%20ujeciu.pdf>`_,
-   Fragmenta Floristica et Geobotanica Polonica (preprint submitted).
-
-.. [Vincenty1975] T. Vincenty, 1975,
-   `Direct and inverse solutions of geodesics on the ellipsoid with
-   application of nested equations
-   <http://www.ngs.noaa.gov/PUBS_LIB/inverse.pdf>`_,
-   Survey Review **23**\ (176), 88–93.
-
-.. [WeberMoore2013] E. D. Weber and T.J. Moore, 2013,
-   `Corrected Conversion Algorithms For The CalCOFI Station Grid And Their
-   Implementation In Several Computer Languages
-   <http://calcofi.org/publications/calcofireports/v54/Vol_54_Weber.pdf>`_,
-   California Cooperative Oceanic Fisheries Investigations Reports 54.
-
-.. [Zajac1978] A. Zajac, 1978,
-   "Atlas of distribution of vascular plants in Poland (ATPOL)",
-   Taxon 27(5/6), 481–484.
+.. bibliography:: references.bib
+   :cited:
diff --git a/docs/source/usage/transformation.rst b/docs/source/usage/transformation.rst
index 6a80f9c6..acec0bb5 100644
--- a/docs/source/usage/transformation.rst
+++ b/docs/source/usage/transformation.rst
@@ -17,7 +17,7 @@ framework and the *transformation pipelines* framework. The first is the origina
 and limited, framework for doing geodetic transforms in PROJ The latter is a
 newer addition that aims to be a more complete transformation framework. Both are
 described in the sections below. Large portions of the text are based on
-[EversKnudsen2017]_.
+:cite:`EversKnudsen2017`.
 
 Before describing the details of the two frameworks, let us first remark that
 most cases of geodetic transformations can be expressed as a series of elementary
@@ -115,7 +115,7 @@ ITRF2000 is defined. The temporal component is given as GPS weeks in the input
 data, but the 14-parameter Helmert transform expects temporal units in decimalyears.
 Hence the first step in the pipeline is the unitconvert pseudo-projection that makes
 sure the correct units are passed along to the Helmert transform.
-Most parameters of the Helmert transform are taken from [AltamimiEtAl2002]_,
+Most parameters of the Helmert transform are taken from :cite:`Altamimi2002`,
 except the epoch which is the epoch of the transformation. The default setting is to
 use “coordinate frame” convention of the Helmert transform, but “position vector”
 convention can also be used. The last step in the pipeline is converting the