From db45aaa9eb272ecc29bf31421dfad287216992bc Mon Sep 17 00:00:00 2001
From: Greg Lucas <greg.m.lucas@gmail.com>
Date: Wed, 23 Oct 2024 11:21:12 -0600
Subject: [PATCH] ENH: Add HEALPix and rHEALPix projections

---
 docs/make_projection.py                       |   3 +-
 lib/cartopy/crs.py                            | 109 ++++++++++++++++++
 lib/cartopy/tests/crs/test_healpix.py         |  22 ++++
 lib/cartopy/tests/crs/test_rhealpix.py        |  36 ++++++
 .../test_global_map_HEALPix.png               | Bin 0 -> 9844 bytes
 .../test_global_map_RHEALPix.png              | Bin 0 -> 8657 bytes
 lib/cartopy/tests/mpl/test_mpl_integration.py |   2 +
 7 files changed, 171 insertions(+), 1 deletion(-)
 create mode 100644 lib/cartopy/tests/crs/test_healpix.py
 create mode 100644 lib/cartopy/tests/crs/test_rhealpix.py
 create mode 100644 lib/cartopy/tests/mpl/baseline_images/mpl/test_mpl_integration/test_global_map_HEALPix.png
 create mode 100644 lib/cartopy/tests/mpl/baseline_images/mpl/test_mpl_integration/test_global_map_RHEALPix.png

diff --git a/docs/make_projection.py b/docs/make_projection.py
index 8f50fd559..b0aca511c 100644
--- a/docs/make_projection.py
+++ b/docs/make_projection.py
@@ -95,7 +95,8 @@ def utm_plot():
                   'OSGB': 4, 'LambertZoneII': 4.1, 'EuroPP': 5,
                   'Geostationary': 6, 'NearsidePerspective': 7,
                   'EckertI': 8.1, 'EckertII': 8.2, 'EckertIII': 8.3,
-                  'EckertIV': 8.4, 'EckertV': 8.5, 'EckertVI': 8.6}
+                  'EckertIV': 8.4, 'EckertV': 8.5, 'EckertVI': 8.6,
+                  'HEALPix': 9.1, 'RHEALPix': 9.2}
 
 
 def find_projections():
diff --git a/lib/cartopy/crs.py b/lib/cartopy/crs.py
index 600f8ac24..d3e4fbd7b 100644
--- a/lib/cartopy/crs.py
+++ b/lib/cartopy/crs.py
@@ -1832,6 +1832,115 @@ def __init__(self):
         self.bounds = [-5242.32, 1212512.16, 1589155.51, 2706796.21]
 
 
+class HEALPix(Projection):
+    """
+    Hierarchical Equal Area isoLatitude Pixelisation (HEALPix) of a 2-sphere is an
+    algorithm for pixelisation of the 2-sphere based on subdivision of a distorted
+    rhombic dodecahedron. The HEALPix projection is area preserving and can be used
+    with a spherical and ellipsoidal model. It was initially developed for mapping
+    cosmic background microwave radiation.
+
+    """
+    _wrappable = True
+
+    def __init__(self, central_longitude=0, globe=None):
+        """
+        Parameters
+        ----------
+        central_longitude: optional
+            The true longitude of the central meridian in degrees.
+            Defaults to 0.
+
+        globe: optional
+            An instance of :class:`cartopy.crs.Globe`. If omitted, a default
+            globe with a spherical radius of 6370997 meters is created.
+        """
+        proj4_params = [('proj', 'healpix'),
+                        ('lon_0', central_longitude)]
+        super().__init__(proj4_params, globe=globe)
+        # Defaults to the PROJ sphere with semimajor axis 6370997m
+        w = (self.globe.semimajor_axis or 6370997) * np.pi
+        h = w/2
+        self.bounds = [-w, w, -h, h]
+        self._threshold = w / 1e4
+
+        self._boundary = sgeom.LinearRing(
+            [(-w, h/2), (-3/4*w, h), (-w/2, h/2), (-w/4, h), (0, h/2),
+             (w/4, h), (w/2, h/2), (3/4*w, h), (w, h/2),
+             (w, -h/2), (3/4*w, -h), (w/2, -h/2), (w/4, -h), (0, -h/2),
+             (-w/4, -h), (-w/2, -h/2), (-3/4*w, -h), (-w, -h/2), (-w, h/2)])
+
+    @property
+    def boundary(self):
+        return self._boundary
+
+
+class RHEALPix(Projection):
+    """
+    Also known as rHEALPix in PROJ, this projection is an extension of the
+    HEALPix projection to present rectangles, rather than triangles, at the
+    north and south poles.
+
+    Parameters
+    ----------
+    central_longitude: optional
+        The true longitude of the central meridian in degrees.
+        Defaults to 0.
+    north_square : int
+        The position for the north pole square. Must be one of 0, 1, 2 or 3.
+        0 would have the north pole square aligned with the left-most square,
+        and 3 would be aligned with the right-most.
+    south_square : int
+        The position for the south pole square. Must be one of 0, 1, 2 or 3.
+    """
+    _wrappable = True
+
+    def __init__(self, central_longitude=0, north_square=0, south_square=0, globe=None):
+        valid_square_pos = [0, 1, 2, 3]
+        if north_square not in valid_square_pos:
+            raise ValueError(f'north_square must be one of {valid_square_pos}')
+        if south_square not in valid_square_pos:
+            raise ValueError(f'south_square must be one of {valid_square_pos}')
+
+        proj4_params = [('proj', 'rhealpix'),
+                        ('north_square', north_square),
+                        ('south_square', south_square),
+                        ('lon_0', central_longitude)]
+        super().__init__(proj4_params)
+
+        # Defaults to the PROJ sphere with semimajor axis 6370997m
+        w = (self.globe.semimajor_axis or 6370997) * np.pi
+        h = 3/4 * w
+        box_h = w / 4
+        self.bounds = [-w, w, -h, h]
+        self._threshold = w / 1e4
+
+        self._boundary = sgeom.LinearRing([
+            # Left edge
+            (-w, box_h),
+            # Cut-out for the north square
+            (-w + north_square * w/2, box_h),
+            (-w + north_square * w/2, h),
+            (-w + (north_square + 1) * w/2, h),
+            (-w + (north_square + 1) * w/2, box_h),
+            # Right edge
+            (w, box_h),
+            (w, -box_h),
+            # Cut-out for the south square
+            (-w + (south_square + 1) * w/2, -box_h),
+            (-w + (south_square + 1) * w/2, -h),
+            (-w + south_square * w/2, -h),
+            (-w + south_square * w/2, -box_h),
+            # Left edge
+            (-w, -box_h),
+            (-w, box_h)
+            ])
+
+    @property
+    def boundary(self):
+        return self._boundary
+
+
 class LambertAzimuthalEqualArea(Projection):
     """
     A Lambert Azimuthal Equal-Area projection.
diff --git a/lib/cartopy/tests/crs/test_healpix.py b/lib/cartopy/tests/crs/test_healpix.py
new file mode 100644
index 000000000..01c1ebb0d
--- /dev/null
+++ b/lib/cartopy/tests/crs/test_healpix.py
@@ -0,0 +1,22 @@
+# Copyright Crown and Cartopy Contributors
+#
+# This file is part of Cartopy and is released under the BSD 3-clause license.
+# See LICENSE in the root of the repository for full licensing details.
+"""
+Tests for the HEALPix projection.
+
+"""
+import cartopy.crs as ccrs
+from .helpers import check_proj_params
+
+
+def test_defaults():
+    crs = ccrs.HEALPix()
+    expected = {'ellps=WGS84', 'lon_0=0'}
+    check_proj_params('healpix', crs, expected)
+
+
+def test_central_longitude():
+    crs = ccrs.HEALPix(central_longitude=124.8)
+    expected = {'ellps=WGS84', 'lon_0=124.8'}
+    check_proj_params('healpix', crs, expected)
diff --git a/lib/cartopy/tests/crs/test_rhealpix.py b/lib/cartopy/tests/crs/test_rhealpix.py
new file mode 100644
index 000000000..552488ddc
--- /dev/null
+++ b/lib/cartopy/tests/crs/test_rhealpix.py
@@ -0,0 +1,36 @@
+# Copyright Crown and Cartopy Contributors
+#
+# This file is part of Cartopy and is released under the BSD 3-clause license.
+# See LICENSE in the root of the repository for full licensing details.
+"""
+Tests for the RHEALPix projection.
+
+"""
+import pytest
+
+import cartopy.crs as ccrs
+from .helpers import check_proj_params
+
+
+def test_defaults():
+    crs = ccrs.RHEALPix()
+    expected = {'ellps=WGS84', 'lon_0=0', 'north_square=0', 'south_square=0'}
+    check_proj_params('rhealpix', crs, expected)
+
+
+def test_square_positions():
+    crs = ccrs.RHEALPix(north_square=1, south_square=2)
+    expected = {'ellps=WGS84', 'lon_0=0', 'north_square=1', 'south_square=2'}
+    check_proj_params('rhealpix', crs, expected)
+
+def test_invalid_square_positions():
+    with pytest.raises(ValueError, match='north_square must be'):
+        ccrs.RHEALPix(north_square=4)
+    with pytest.raises(ValueError, match='south_square must be'):
+        ccrs.RHEALPix(south_square=-1)
+
+def test_central_longitude():
+    crs = ccrs.RHEALPix(north_square=2, central_longitude=-124.8)
+    expected = {'ellps=WGS84', 'lon_0=-124.8', 'north_square=2',
+                'south_square=0'}
+    check_proj_params('rhealpix', crs, expected)
diff --git a/lib/cartopy/tests/mpl/baseline_images/mpl/test_mpl_integration/test_global_map_HEALPix.png b/lib/cartopy/tests/mpl/baseline_images/mpl/test_mpl_integration/test_global_map_HEALPix.png
new file mode 100644
index 0000000000000000000000000000000000000000..273f4d6ab9be62a858a94b44e6f6a4149e02499d
GIT binary patch
literal 9844
zcmd^Fhd<T-_rH`quMr`8&&bLQ*SJ>p-VzbnD<mT_%U1bVSxHFt%*r02WM*Y#@66xv
z{WpF+9&#TK-1~l?*E!F1-Y>N@m5K0a@DT_E5n2VM1K+!@zHqVO{}Q85|G_srR~2Iq
z1cJo+>I<V#w%`#0!Oo6G$?5sN+sN=U(mS~LCL2n44fC36$e)nM+^**CSmYvd9=&E5
zrB5gNRl7IrH%)7-UJ{It^|2k~KQbXRd2L2SX~rowC7+1u`+5s${P55DY3RVq=X1AP
zczrVm|N8y$vl=+O+w0VHdensorhYLTot$^gD(wv+LF;X@GDbocjQ{U18ni+R#~gdQ
zy2?DM{ryh%x$NxiO(h5$AK2TEj1=7GyY<3CNJvPf#K2S{#@)~Qx_0rHu;VD_?)G-p
z!QY>z5_GpdS5$B@G6pv^Hm)vqC#Dc9M<gXt>>eKGx3sj-aBQC({OxXQ3ss7zH}_!s
zN^nnG8`;{{Hn+7EUsEH_nXKC0*=ZudiwIWD3(3r6ctOGt=k@o8T#+6#8(YVBYYap}
zVq&%)ce%^7#Lh{P9=BQq79uh>He9C&WoAZqc784xaO_G!Nog$cQE1dxPPK<=%51Ll
zTYI~auxPR>){`esYHDk(*T>7ZcX#D8B)lBMSX)KNf>id+mDQNKn9Nv`RrPJziiSf}
zrbly?<69>tCTtQ4^|*cg{oyqYO7G9Y%8G*^8ZKjQWd#@B58mrq=Cg~tk+^<-em7EC
z+1atHs;b1@7b6FbMD!DAcp{VN!o$O<CB1ngafQ+eoqm22aeB<1j2}xY_5*3N`aMrY
zRW&3xmn9N6L(->8MU}aWseLR@jY|DLshg-U`@#RR%gTmq6589_TY7r%QDL(|rzL&6
zZ#dDh6z_A^Tl&)Pc+qR7tG|FcrL$e#q>qnJ*8`!Q<6Dy6qsGd|LH%z;9wwu@3=Y<R
z>Bh{>cSH|42zT}OD;gP5N8;9dZwS@13JMCAm6szICTcx65WnlUniNVc8n=g)8JJ{c
zn>*1u3)cl%OQl?Yw+I>2q^m#NT~A?*aP5nX#NFA7AYq95tEi+D5;@cqcq;VISy%gW
z-_G%l?Nwpo>7~h1pUgIy$9T29Ao4xksAfv-om*J2bw}IS+Rpv{jW9Gc%y@bOiF^)!
z91>||XBR#wdr5b%Ulr|`#HgsKNG<x9nnuheM$%`~YVbelmEl~RkdP49{T1!S#YILj
zF{)S^k?(J}#PV37r1bRk5_WMshM%l@Q}`H}QOe5Ee_f_(w9RlhO*utG*_JrDx#gix
zp!;I;)KW1KySs>Im*tH3ilr8_AuH0PvDdW?M?abWe9whm^|(A<h(DN-WDaR=mStvU
ze)aBMjA^yY*WO-&xw*L;$?rRHb@cSY%F1pxobC3eN%^^_OEqOP=u%ASb5p(%ex&oq
zh>3~GsMb9%j_{tb9AX5mqZ8NOjygU$x$s)YW_*KbzwA-=nHv+?^N&KVh%h2PeuwSl
zORmdzV~!%&#B{hkJUlTpBKUCLd~nNJa3VrScZ7vUZ4yRt<)WgZZa?TE&5+_oTAv<n
zAcCP$=3wBAzBnNN3}(i6L{r(e1Y<bNG}Kz4A3Hn8{?N?PDf+#%r0?Q9cJ^;i4h3VQ
zXkkI0laq4~^{RIXmxLrjoxKgtJ{k&0ZB8Z^H&$$DJMrZaW5V+lFJYH}8wLgwI?RRl
z_U}X<K6em^48k0K!l)Bki0ker-p~-kfmT=NU36!YlG1TKubZ1(T3oc+`_rE$WQU)m
z%63D*BIZy3TT|Y8J2$zoq_!VFs`}be)6<RT4vmeCRnX|KFmA5rCtmOK^XJ=NlDm3(
ze)Lvl)>~4|GkPob1RFUA$YHbeb$hN&+sLS}!Xgx#2th5thFqk+>(1KHcrGVQj7n8S
zyV0Q%<rZUu&QEf)vhGbam6eeyL=bnpAewxgbe-2&Crl|eM%OhSAfvI7%xmp^Fm9nG
z!;tPHeWts2?{@X}DwvqiK)D1xSF|eBY_#Q6Z*<9?@9SpZ^b8J?6dRPZ{w-~eYkYb%
zKae4rJW-vPNWSRqm#j;c9d4GbmBgd0sYwh72oSY1;Z*uOawu1s#Eqf(;)I9p)~&gv
zrSKBd>Zn;mLozvHj)K&C`!W<!>ts4oN?JoV%^Hj%Z{B%mwf*~R&9loUxl<7(Tbpi=
zdJi0jU!TSK@8ADDBi(lBwllBjxixbO6~?adhLN9Nv+rh_WWAlBoa$`M)RdX9^Mu~+
z&pN|Wi*Xobd;$WVsmpwgxA#zG6&1yvsyd%fyQb@XI(D_D>%2OQ*frmAermRE!|qyg
zCq>RLEbvS<8x$*EUYx^dYCS?<q}Vb>ZYGhk5THin^ro;Al9Tm|&hhn>$pYtR8wq@U
zDZWwD5f+ZwRoDA$(YLm?{@&WM%qP^@tqAz%4t@Scq}~oGw?7@$+A6R9PQ1KN@mGx-
zi;k}DP6w4`OIzDDRG4+!3%T536?%md1KZA+n}9;kV`F2#USaMRhhfa1prEd<u9msB
z2sza8&*qEp(NR5dNy+H)<i5{~fBwwxuZ|2kR#MV0bSH9jEKwMejF0~iHA+}qd>}5)
zIk!TsD<4Bhs&RdImMcTcyO>-5HHuqbUG3WXM)X34*nOX0l}<m08vHy)iw#K<IqoYp
z2b~X<S+@xdeueqLNNdI)Uk^$0*-@(!AVtE}>E=sfA_(cFs8M0VxyrQrxAO}MSh%@+
zvWVm4dxo<Wa59W0s+@}xHD11aiA=m)9?YZ-+GTS++P>H3Q=OyPC0V9DuFr{fcNbb;
zPhyPXV(Oe^Oo(>V%#fhGpD9zKBb8Sa!}{XE&E~VuMx^P4$6L>$n@1cgZPxku`B%Q@
zMb}&$?<_!Xbu6(77O*~~M6Km){i?v(TOR!R$aO?HmEXKQ`Ug-)Xa!lG<W_S2gQ^cO
z1X4;44xGX?M5B%(WaQ-KPiAET0Ep*zfzN5Zck^(mZ{u$=bhe$h`TF`dW4$$5QcHOE
zj*ynel(gfwT~9mNn@I8pp_!EU@AkJ_R^x8HKHscwc_Y$7Wr=WeYbm$wkexO*?jT~)
zZy`l<>;6Y0+hml0O%%MoKEHLiF&Xjux7CinMWa6@!u9VDzMXo55|iX#skN+vZy7t+
zG(H#<UtLfJ_ndLX!<%26{SuOr6l%sR><2qP8dto^%)B(D{BRLaUQl44JYc&wx0k~x
zZ8?d3vm{)DZ&HZF<f@d5-mW<+?L()&no+wn!tZ7A6?YO>mvI_SxQ9v^wg;bKZZ+cc
zrC&oSV$h0VS_WdBZC`WUWpzE^juN_qlzE&sb87Ib=f#v+%g<`@cTIs9OxIi&I*E4t
zt6Zi{%w%O`=+93M+7`aWA8ZcNC8Q~m>Uvf-dHucdM#L#!3@;_;3+u@uXZL$_A}bFu
zx5wiQ@lUf<SdWA%ju*(zL@B(VrT$MW6zbXh_^5KWbDG(5ekxw$FT!5+D75uS6PG5R
zK`Hj1{!s6El-biI%Cmc6x~St=+@r{>%FiF(XtBvCqHqm4e{U_Z%`}gvQFA&!e)T;E
zi<AtL@D{#$8lo=^BVZp%N`7s993$Y!+|AK@4s-h(Lyapl4!)urt%j$tmX=m9kn7s$
z4FhalH5+$1^}ED!#bd2AGxUKc`$ePPA*^bt{K-iyckYlA(Fk+0gr%lZiMmWhwzS-1
z<=`;n@~%{h4z`OJ-WNcJ-OqUClT}t63u3QrU<RY~VC8xBj6Z>D9$@Jx(spGi``YL;
zeY{7NCJlQ7lKxBQXBC6N-{Ob6`-z+*w`q{r<QQ+ot`!?#>@8nIB3o95Ti!mwoUX$t
zF}<eo7IR>rBOF_ftd{O#X2yu-CSIOdToO0_5mzGt^F2+hwDA=U0ICdkOFcckyia*t
zY=X6nuV24*7|sdrJ~{c~Jb9zh|KP*Ew|up0Al1RI3VR6>z?%h_Xa!#LI*MwSX(5A^
z<BR!>=Ij47l@scuqn-BFX(@khCQ;-r+8=MvZ66(#-VW~_80dV(q~O_PTJbQ}TFz0p
zC6S}0#FSIQhj=dV^hUhbzrkn03*Vj}`&un)GK$`ajN#RkV^`<Y)RvbIKHShAe~<Sc
zCk5vTkH`p_HR;AR<WNTtW!x>veDb$L5g&zcE92p!lV=R@6A6U8*7VdfB*q?`GROu=
zQQ#{kuboFnMfHq^nb&)3EI$m7h@j=sYYcd{Lp;GR@t!4e_Tzck!D<*F1`{725pW=;
zwp;Z{)Pv*2#QSAFTNZ$_NKjgBAM}f}#vI))F9I_TunDQP@_Z{zob!V6pavg2%uJbU
zkY5wum3!>iIc!eJ1q3wIoV|V0RHS1T{OG?*DvMXZxbpm|EUc{US=jAS;?FJ`d@Bfx
zlti&VtUuKfc6m9`WS;FXMaEi@kkm3UVN$@ZacpF4{PpWs3{==m<I}%%3#%gqzW<g|
zJG^{0ChSr?B>ndbXFMFHY6L-VY=2i1SPih$K24#>_x)j!W~|YZ`)DNp3tYcM!13+c
z+S=iG*$bnkK=Y=6bfs7tMO4Qf7m&2J;aup(RI`;>&3>*eRWiy47OWvD{)g-6x1!ED
z(oYYy<mE$ZHs9U#SOJ(EEx6C{M#AgE@%dD(2Ny4IjJt=&%g7$}y@lVu@2#$`?)Qt$
z4t)*9;Q~d)yZocXRBU%|kA*o0r~@A#pP4zQuTTACyNzV^;{0^J4zQAgr<mqYg+^aY
zI!yu#siLlq0y)rbSpevpTL0;G>n5Pxc)gFn{QP`-NPR>EPSDwoaz!4!{5mQOoyZPK
z+H$L0?8ywTx0joj*I1EMS!JanoQJ1`0liFsk^R$HrBIM6`gQl;3aY57UDI|`_yZoH
z%ymv)P_x!7PBHq%-0o8ESks?(9Up?DsRY(ofByUT&+f-vs98fpLr;Gbos{noZ#n)0
zJG*S<c=}h}53uc!gBg+$VzU9d6x>WqFm;2<nKFT-q@;}hDX$97!X02Dv!V8sqMdIr
znL}b@iGt@{6j_vFS~jQaWBmzlv?ElBALA1ej{2>@RpL51IbFYg{rAXy*^WZ(+_I-f
zE@x+ZgElkg$%Ynk-{NldCxkqB_%H`pU5k@q9GEq)YRvHKHrRokDFp=uuaxtn=1a=d
zu}|i-sIaUoMyY^fQY<Vi-|gSeHD`rvcNgMiLn31=+Vp7PlY>q-#4;$XM+>!WzTQ_^
zYrecV6`I7G-Y6+40UOcJ+%=S~0Bq>iD8v%QZS^$-Yrd9)m34*J!`(fB_e{)rg5u`Q
zn@Z7Gh)))cQP4+U!9T%4imf%N?a7*|$OfH>=jG>rhoTfW3NSD*)Va=u9ACnut|Z=<
zVPejay*LhyRN)o=pvNuYxk6fE+5F$k(=BOc?*HaTM@ABbX*4QrU(b&g>Fyu=v>X9r
z)!p0M>u{QJ$3E`(2wD&?`1tZk^SKfVH8V4#S85@>bINC0HPn20{sbI$*2Dfc!~gb{
zi?Xt^`kL+t2{{a9$+0GOJ9I=*`p0Mf_n+oum2)oWMI9LPuED{0XiJsiu^$bXHB8jd
zh151);NaFbHiqlGoZc4{&@R=0?CW@b4d?g5!Zi*K4#YFVIhWa{Tfct&vZ?y;;K05+
zK02Bph2Q)`HLtF@%*hH1DRT3sL{uRe?h+AmYxQhUu}{YAvrEvHF4J`{XKzXQYE8K;
z^(0IDwSV|<Eou5=Nr~IuAB_yJgoFf7r58O*?jOx+$*2S@P+K#Nu@hY%;xAWue<X3~
z4!uW{=ibMxw&}npMPd`tgn^^T0)@lKoY5p}ZEZb&a>8z}#hsjX+d4#hLBf3zx6XSb
z%42!prdM;Yi)}%1u~Yg$5De)&h1Aba<|D@5^X+(~tHijyU0uP|)naDm=C9MzB4a??
zAzu*R`Y<y6tS?QlP`~&_WE4&+9<I;;k2H0A*V%o7oM8vo-NkN45sFoX;xQY5#zM<x
znMhpci7yz)Y|tX}4F*77JO&@{Sz9wb9C-Iq)OjNE7E?w>#?D%Csm)w#xKb<+Ob#jx
zGzpDYGv-2PEIC`cqz^xWQC5~gH(wpue7V~zaNt5qPahW?j6ig>ir>Bab!?1M%w>x7
zr{B@Ej~uG3y87Fd?)I`1y6JzgmUoBn-U^=88e#9{7sh;c0Ra-zs>gC#TCc|K+Gb`f
zIPw5TN`KAY1ZNZ3#aVKB0xb4Ili^b`k6{$c!HoZ3zV7t_Fd&M4r$<{Pl9Kg8lz5iU
z&MU|9Wgv_(t98d}ndS2CJ_qd8k742B>SkRgI3uC@-v2x!;}&2_930>Z^es@PsJoBP
zPh(}ynm^0SBOqy1*b`r?v*-4v3M4A!G9{}jlF&=tR+YNAxIm%ymouNy*gebJ+20@a
zGwg=sfs2PHIv_BROG4rem~^wvs_g9S_gb7QzrTjQHz+OwU17&f_JZ&xbJLRtfDVg)
z|Ed)g6^&Nd-)Cm7gL|sk1UtolF<t3o=zqMOtvxWzd)3;bNfamslM(|C9-*g)k7MYi
zX&@R&`A95+4L#xN;SpM2Uq3*R-T4rp!ET<a(??A(=*(9zoQRm1-$&XyQ$8H$JeiO3
zv$@aK3>I>=&+ZKT2KC07YO66Xxt6<c$+nxjJDGhv2?2qEr>8KxdKx2t=3S3`?YNAT
z6l_UJNw~FdoTnqoK|wOm4`B>J#~4T?vQW1mPCO13^*w{IH&NB7$$P_OXDh*uGP}=k
zSGnNxf)8OhmLL*Bi!|UJ-)oXyYu8H*q9Gtz6Q-7Vc40=B2HAq7`dd+=XBYqME93ZL
zjLg8U>v2;O-;xLeg?Z(G!0fSs-G2qvc%vs8EZzUS3A=`nLm<2vC67U*{_i~6TOu?>
zqQE5s!V0MID(APFl?UL~t)Dz0QC3zaq!sO3DRO;eL`tr!r^l)(&00ASf<*wRkFdA5
z7k)fmW^(7HgD?%R$(KTP?mz46F*^(2+6D%QQP;)9s@=amq-_OXm953u^Tk~6Dp$Ko
zKN1sj_ZnY$asRt1^dPe4tTmiKVnc1K`LgMlIWsfpq^M_^35&qShF;1S!F*ofu)@N|
zhJ^})i!wYN$HT*meECw&-(Mn*?rsPU2|WT4#~>RgHsh=C`1eO%lad?5!^4cs0cShk
z#^NeiPud3t2dR_=L>-3itsDIVqW%rJk5Qw)NYllMQ)rJ|GSpeHwVu1jiGfq2+b@Hd
zUdrL~h3WD@#x>*|kV5i|_vTcsO>Ou1&Ld%(hKrNG5`4?@-`YXZO8IUhF3*-_4-7o7
z<KWnW4F;~gV<KhN?<D`lgzuJwXXM&wkwRMpiI@Xzy1Eh^2`~Y~39plrDK}=DGe4Qv
zKO7zZ{{4GoVGW7C-IdMeGp|cZ;WJgz*49>HKs>awQ-oKnsl6{F3kwS)fk>569b>hJ
zhlfdNX@wR2lTVs8P#vKKP}LQN%VQj`jfrFcNtALN9rCn=^2Uwt;MjRDgM4z=@*}9&
zuz*#H-7x^g5{*ER4hSr+|EhRx*nu~x>H+$$VgI`tE*_p>AvQiIe-iX=uI;xt_RPTT
z#-kbX!Fr9iq7p=*d9q9T+122uQ5_Zfo3qUf1@|*=q$cxUCgBt?vM++HK%syvjaQwa
zX;BOT+ZdqWPk@Xax)a#GfSGrkX}AklW*{j_MMX8y_*CqcOhDS<=Csq0C|LjeZ}AKx
zAb+|iCh8^xRgjdFia>O%IrB3O00_d~>qcv{&9eVaPS!q3y8garQuX*|>@CR@lWG^f
zAO;#5)dI~7W=l&;8Y$n~UhCu9X8C&Dkc+CZOGoTl$z8NgPQIO*o_+*p1;Sr1FA$`l
z<M;P?{h6|vRnC)2>FMb?zGH0S;%O6=PRc-<cO1?^{veUB<+E@1?Buz+x@N&}>E<fM
zu7D{Je)OFcTsm#Xw{N;uR*axFU%!1zi-m)eD4)$pLP9c8=fzWER+|9KisnxpNnlf3
zfr9TZc9UlW-cJ__kBq$f7(SCPS-_Rki_5v8+9NE|4V=cj2%FCrV}_RcPb~?k_`}rM
z@2mw_u_(udMMOl%{l)jtSf!<V$PLwxx9d$|;^b86qfZWUqSoGdhXGjkk)GkU)xBC{
zp-PijIKh|3WeL-JYHa{mBM>aC{V1U=v1l~0&(gNL>%C6RXVH)bG2i`9pFSx?3*5f_
z%4MeEk=C;^Acd>@bbts4g?c~i0NzMZBJh1oEHHonf6{LO_;C@G6ck!8i~VVWq4M(b
zfNGi)^z@pL9^@cdlwLs;U0YM51>tymYaj?j2gp+AW%l;=uHbzC{AgAe^vDWOf#^Rg
z(OWSe-qKsVr??(0{?d5WSdKe1KAr>^frwv=R#im>7uEno{GLZE!wPh7`8~1l2nmH9
z)swTcvJk+)!;+6WL9VnnhchAaN;3_LJ0JtxhC(;$HZ_%TMdslf>1P9kFRrelp^op5
zkB_<YSRpPHc3)Hhhm$fq-_zZl4L=$JydvWElnp#d+MH=rgS=fh@Wk`vcPL3@BqwYR
z90c`W0F>%*Qh07o-4&$l9ULrxekIr5DRX)b@)wCrt_>1Ok9n7ou{wKsA?CF?l@1Yq
zhaCuP&-HP^+;7e|<i41c+$FZ9a2)%fh(@DFKj`!LFT{#DfbLA5kdc&pySuyF<C#~=
z@3T3@6sx2FlUL$8eaC*_75HiM=AaD!qfHIKIbo|71T%Am#tXwYetUnANNsR0D`1!Q
zx91e#0>5}HYXSn12Lc%CRlAr+#l*aa{erQ*eC_ub7@9@*T20#Smgm>FL`24ppVNv+
ze{YLU=^%p`j9Sp<I&>{hWyZEk`&Dxm?od!rw2X|9hvJZ=m<jOk5X8mBIZRgZfxf#+
z(nyh_kH@XP;NvQAadWo<>I^kL^)WIvU4SkglWektOb)VT@xT*de)Bq3@1x<L$o+!M
zAbGG5uBV628C$8m#v$Gt6K;m&9-rM7lt^!iT0WMO!|>jkA*waA$out&$Qr_vvd^D^
zE~t9UmmvXWWOj?XH7C6ZT^D|p^u@Z2DX#ofLOIMfNndGFsi%j1S<c~eUHgW72I!`f
zH6AzP)feMbRfM}akUuIy@=)qp_vZN6*Ky(t7Z~Kk(9mlzuisgB)YEQ12Mo^;cgvab
zf`LGzI{e&o=0)+DI>BnV1Jk1)=OKF}K-?6-xj^u@R|r>+?K}yBq+Cr^PcLcU^U#o%
zNv%7Z#v5Vu`n-**5h>o)2EqEW;N3eq8dyeDKK_*ubbjzZy`o7_E?8b(&IH&c@$~TH
z%o0d#SR0YRF;XCrawydK=@ADIGN^sh<qIJkSYmwqaY!rX^1ev7;5&SLPm=P**5yUK
zjg5`^j6!fE*6T&c<|%JP-OrFvgx<UAV*jONLv?yxh$dHLFaO1h<C&r3=%rOrC=eF;
z;31IOJFI?yf$w#S93VSpkH2pPG6JP1b=OiKEa25>BUa#DI|)m4${3fjh8<Sc#)=b2
z&Nn8jS2KKOvcd4^fXePC`<bH{-PO^dxIBh~>W~hC#e?(ObiEkX_3QuhR5=n}YiVpU
z&!J@6uC8U8`n<?tNWxc&OPg{az0+1#k6C<A2}>Gq(Ww^x-S5j)txts4u=0QeM;$BE
zBB;7xT45}-p{Z5@ND4wiFi$``cmMeDgV*iX<m3qG#z(NO`Vw^R5AqNjLXcQMtS%q`
z_)xo)j}c5?9LJ&zN-Z>Si0QJxszcd%&-X{}qfj6PS?U5$tI{ZAU}9=LmZK&eD@{zw
ztT1Zp>Y{Wz)>v3rLf6;VVe{25%wavs+)QsVYVvVV=9CZ>1|9J?fPspN<3|qcAQSsy
zjF9#N?~t|*Kr=u+J6~}a(n?<RXvv&yD+p==1p*XVfj8*t>$|Qc^)N&AOFqG!Y(&d{
z(oYF#SZ`2Ma~jxQaYeCTzas2FAo$PQn}+7^c)K}2pW~s>-w%Q|Muxn;eAuxP_J6Gq
zaQEKo_l1NMR7te0{Mye_fq^aHvPt~-tFSyW<gU1Q8Dl4ER#G&$IpRCMKwZH8%4_ig
z0puL$qmTQGILRG)clZrbdz#uuw8+TF+RNAR!!be66uOFUUg2sFPif>)$IQ${2YLHw
z>GMzAI;GEpF3xHP97L^m?U;oo*7tN^`x+8SN=xgsYpahvmvIDe0jn9oZox2Ka%A$U
zeCC%g!VO@<C4N=ZdHuaQb%7bDtnlbH$o*|uku%CQkAHo3x4b7{b__;4`k%G6yuB_N
zEiFNw&dd9<=M<`WpF34_c4hcYs|Z!|zOyF4=1gE^vRdG9Yqr$qA80Gj_A7%1IdyMi
z{8-0g2tq^+b|rIPb&K@2_x5-^+`{>Z)KYkn1aK!a?6DfJ6nwW;zd3S4zFqZXHV&4v
zjlBoc&EL!4M8mEPxc)l$u5T~NIkI$$)Y);Y*yV!XFrocO-pQZ=DlDFEh>Vt)8U}FG
z!nPTa_3ZM(VZ1C`o1lA1+V#&X)QVW19(ONqV%?LJ1XRnYJK`b5^|Ftwl+q3-tlC9-
z)Z(sLnh7#~yARx`i^@Ov^7s$+F`+FiEYjniznBtGL)#GFMcg64G*Ti{U?|&V<>JCU
z7dh!SqCN>u=c#AYN<V$zg6?C(vBHo;?5XPf!YOVL>ISSXyUOUnI4V6QFtB6b^+4++
zpW5=H`?zqdnBX^6bjEN#NTFlU`V5@9YC0LGHjd~%J{&7KSTd^WOyd+WqoJN8%&lM>
zxf{3EWK{AK+NqCnKQ9|tm?}TwWx}zthF?{&`G`{}g$8}fk)jM@!m=ok=C5Iu&~jR8
zYQ*@s(eb~3*|1cGq>Cb=;=`EZ8KhxpcRkoawmB!AGQP9*^!#%Cc@Q*17VJe=ewH%)
z=}Qk4wC#-CBE-YRZ2^39Jy_FkJgrEgyGKNN{q`r;&I_+G4s>@lS28^<UwC;1BirBs
zWN_j>n{>o<;=$lyDXX8%D#5Driiu!V|B@u9qN0NOse?n|F(w_PU;ubGH<Y<WMGcnq
zCy8i9QPDT}ey^_L!2%rSi-nO9feHhC7X3FoEbPiP%m?vu;b<vb$N0d}{hCGjw#OfJ
zR#w(D*1t?DWHdBUAY7xj8c6Bub|vPE#vDI_#4`IyVJ;lTe<yEt+;sO2kTq;1F%huT
z^Gbng5JAK&gccQpH9zDA9GdTj{#nen6M-+4m{tx52!ON!9o>U5$|#3;gocd^5N8<t
z+1JTQT1`z&S1{Xtm#B@2N(}>FeQ4>Oo#$r!0(*$gZ{JdlWg*f9K{qeq?d{FZ(}14<
zss<-B@bD&V=E?(458W+EINCk`^xaG}Owtn_J6lS%Y(IwNkCz)qtL-{Qk7~}CV@Jy4
zqImi2djQ60TNoOLA0N|?O@!wvoPP4VWmpFNQw+c;tgIBtS5LP|SsUzmz?UieEUuwJ
zszSe#mK&|AitqUIQ*37}Z5#ZXcSne5k=`>*fSz3QRPXKY4_C^4Zc?J;?Jb&<n+s1s
zI8<Q^Gu5*#W`@h@kK}8R&#p{l$zAJvBSHX>D&}<UO?HA&PI(nG8tP|lZ;x~J@BsH=
zj}{9S#zoW1rlux|<{xadt9{dTyo@UIi;M5YMATGN%8aODqCUN3t3OpSPlW!AW?^-T
zf*!hY{rL1@U44Bt$t)lm!bp=CLyk;VR`%9i@KR-9$D_UHdA24-!-#WFSp5#y)6;f;
zu-*-(;m6eq&;{GcOWa0CKu}RbhZ(~9q+tg;1>_6Rdi0SiC1+r)V+b8)Ug<Ch3~MK+
zXoxV!IPV5i-6PF9{zFVFN(ei+tCZa_j1H4+63v_8o}~2Z0i2|W5?n<BJiQWaxTLxx
zkBfPBPGCeHk%OO*=FaJ;kNgnqlt72sIyC==8wb%9?AV~%-3QOuTyd3LDdM&m<YIYO
zo5F)wola(uk-$-Xq5@CsVDFg)q2*XwZG-BbJYMVVKQ9XC?WpmLNMF1T#<Ko@`{MT{
Z#yhKTLkX&gb9jslfmYN+eUi5d`9DPmu)P2P

literal 0
HcmV?d00001

diff --git a/lib/cartopy/tests/mpl/baseline_images/mpl/test_mpl_integration/test_global_map_RHEALPix.png b/lib/cartopy/tests/mpl/baseline_images/mpl/test_mpl_integration/test_global_map_RHEALPix.png
new file mode 100644
index 0000000000000000000000000000000000000000..17c836c760652d674bcd68e112328754a9c5eb04
GIT binary patch
literal 8657
zcmd6N^;cA1xc8xv2I=nZlpeYp5s(^CN(t!@0qGJDknU~}kdRPX7?4&P6eOfUI^NB9
zz3cvR|AEU|!z@@cXU^Ww^L%Pgl&-c4E*2#g1VOlJs*3vHdGP*$i4J~O8`tE3rw6X8
zChib~55IpPmC2SmKoCQSnxdScPu6a>ubJWW)qreI-QXTpOP&?IdU;xF?Sq<f4hK`)
zj5+5ITTg`CFH*1_99t4IK0J5GXf0?dYDq-qh{ooP#A1)3?~xEfAdatWGPru&C|qr+
z|M_-Voc+rXKU;9mK{*q$3w#$EiHXyb95ek}GziX3L5kTE`2`tHPH0Sz!-R|-3RkoW
z=8whXCI}}CM(gHhtb>s&qyArf(F^xyjH$!<hMXgbk5lTKnv~)wMQEETV#BQ6f(O^#
ze8ik%8xfiYc8__aq@)s(lY1HALw#$=dYmb$!w5w|D78Y$-VURsjj^GHTn(R=9wSJT
zQLc-SY2{}fL#NqCRPmrCk$_csoPZ6QFH5vgLS`oRa?aN>Nmov0A}9g#+e3#uJ<?U@
zpOU(qDcl4lX3UC;B?x*kF|wMPngk@*tz|m<N%Vwov=-vxa~fyjVrSz*P0bkFQ%KV_
zae>Fx3Gm$@3=uO|{c#$-_q+)m85q#mAa(+^1hQZ{4lnQH_k72(Rv%;{VXplA@V7qH
zO_YPf{_GhT_>H-AQezggv-q;9ogW*TzLSw4ad0FQ73J{+-E79|B9w6EP@o)_CX9<E
zqw`u*cp8>KEzEJIO6r|qMVy`Lf_{~8C$(k3C5IzQ8Sa-)dsnLZBQQVhhfoj-*4W4m
z{hs7tqGc00-Epk7FLF7$gM&kvWmo5cG+BTvwtMhZHUts95XrJ^M@u~wOI)$X{xZYt
z#jiuB(+*RPB@eEC+AKehM2EVf0nyxkej3?Ev;Vbb7uQ&_d`M~-m%Bs|A*pB29BeO4
zSyWU6L6pMwSXNR&KG0{#Bgn7<szbfG3<%A?!G~=d<1ofqP;%q-6L2fj_w~Wx`KJx)
z#quv-atl-5q?%S6*-VybC8!r4oKOcnvFr#Id0~rmPWa^*dX}Inzp^$(l=~(l)@k18
z*WD6L^pn3ug*72>o^P$1qG+fziaH*qmEs@c9do7x#c7GZ3i8?w>Q$tthw2(6W6H|-
zJvMdmlUU>AH4zRAojg77Xf}`<QN}P}4MZFZ=ylXY-_SL4Z;@<AvmUUqv2{pIQObMh
ziwYjCxw_W`pDs1&CFGVi#YOn0WAW$0sHs=<b#?!An&RTn;qCm?n@%RF`y>-InqAAl
zxy>c)Xh=sT=B#+Ow^->+g&kTB9gEbPnGd6K5PvDtGkx`Q4*6*_a;6x{*RNtoI~cBL
zFbU@OR}c{r^;fpsvO2Hb0LU|6zau(ROhLeVkK3&CwF!aBt5*VlMzdo!H|-lp;(gT_
z9x^3;!-T$FU!Qav7(_{Uvm!^nd>Q%4CD3O-hOd(zLJg{_u69`M4nqxk>Esml$gwG7
zZ>dAWv8hcn=jwEq37V>dEA+ff`bJxwfG~);a7-Z5N>hkv<kBuEXe1f#446~%qM?zn
zk8o&qnh$s17=`!8Qb0i|Tu-oG$;-=+eSEHg$P}y9G^%`3iw+f)lw=ezBZNXtGL8HC
zzQj_%mOEj}=wlTm1XS3lXj>{y;{|At6}X*8-h0&Y3JS0OAShIe^)+*4XiBxSyAkO^
zi76>zJrVdQ&@({_Mn=Y!zkjha@6eI2Z*E46r2cU9O+8Mop>WUN_n)dS{G^6P&~YcD
zCnq<1Zen60f2{7^*WXV>_+Z&niyiVf{H272g~fS{$XAF^EgUC*MkQ9C7R^7e>Gy-_
zE@HvjDQSCvE*Rc36z_eyW74tT+1Xk5al*($TwFY1c-__8`+m0_f7Ha0S@`5HcxY!!
zgyl^<e%ND3V!)+)#)UgrX24cZR22Q&QJ9ia@EImj-hSs*%>rJ_yXC+;^D>^%9BHm8
zQqGj<{x~WXXtR0w4O5uBUb&N52tNs2H2C#TFEWmlxjJkz-m#_^KZbha;^WCIF3-<>
z#yv3bs8u2Rt%>62d(Ulc-{<D$@*jDEz0#XErq@(Tt4_)>dw~{YT{i95=s1lvk|FYN
z#$bzmw%mF4Ba8D=y9wNHcdiKqQdCyPf^()|!QpUrZth-CH3DjCoMQc`w{EQL6#WE<
z8H`!)HOw}z-Dp~AKTJ3ZI?k7<hiogC%ghBSDJ0a?)Vk-(1>E&yPpgekAg}$UxUw>K
zZsR&=b#=A-#%p(Ob!WE0eyQCz#r6k93A1gb9B48^8k)Gi7%~({=3?!kS)0d$z}t&b
z0+nL@{8zvG%bqs#f{!zWV^BjAlapcJ&D%WWe0(06Hos{(ADfvWWnf^q?~Nk;a)LLP
zCk~qFf?*q~<URbJRC{nX^E-8Qbw)6lrgc`q1Y!~r8PM2Lhre3+-#SB>y!=}Wznq*L
zMMdI=4;4@_@QTis11)Ql?rzSO4Sx5>g`GRie_EHpG~lp(Vf}A+o(#N1F;^yFtja`)
zn3$MeR+i3tZ-MaFreo5dkxcBQfOS1Hv)F_Lf{wbiwP!7@tqcvJ{-Blc7^I}AyKYXO
zrc<0fZFZ)}H$)i3my?uTrF8nA)mWSZoD8)*29r-fsJAhiU9J*On^Gd~G}q|%f~CIF
zVX_32jC1uWg6_?r$o<Xy4S6IXfeZ->3f>ookB@J4sYEL?PnVs5ge1&n;N9bl8PFaT
zS%GF|)adxM+^_Mmu@&_7^;0I{IXsBALVh+j3>g`jpY=8o78uv(hrjdgJ{5zZs(RXt
z=h#F!|H(<g%8I_td;cA{7y}<4J{lTYWL#Xee&>GsF&;GyO=50tEg3B#;o)YOb&eyV
z8LrM6($plfu(&wfUUqnR=(8F|HulNsDY)jz_B3|a-L2&Nq$E1=Kph<&Y6%a_H*enP
zvJ<cS>8{}#faaX8j0y=s`Pk4<Z}k;fQALGGK!ET`twrDA?*Rd?9b-6+k)|dNw32H1
z=IiEop`NjEBxrd>MMYS(QP12Q<=yQSUyVgaGys(7Kk12y!Jrs6h%|xot<v1FuE4;Z
z<t{pCWoIY0y<KXOG{g6!fG3+u6;r$2ciQoHo^Eb##apE%B@?oD*9!QwQhA^bx`0^9
zrybv`7o!FlRT;|}lxQ+Q_eId>Op(6?+TxeC1k;cW2`1rqrc0+9w7HQRIP)-~t4o&0
z<YVaN>8=LYLoOjKBMQ&IGc4J+(~eCE&c?>%94YFT;H9Dbq?iVT@#cWZ`V5VYRWB|s
z*4EcqObbu8CL=_h<_x!(841VRJpW=sS{cGoc(hXNMTfs)_$Xge^+<UIZvv5V@)d(F
znXoTI<ORE^=m$~f1=XdcCDfpBd^+AX;vCh_I|18M<<DGPbf7Y$YExoTQf08A)gPS@
z4KaO~MQTZ`e%<6kwm3A>evO*#D=SvWFqf_yuVP8`sK^@2F6|dHRfMcel{QqzJJuG<
z?Tp+O?c|`YX;Qv|j?)#7;q9@pIM5tnO{SR3uXx#8SO_{UXYBs(-zao+bcOT7QcP*4
z!Tq->M&*W$0asoho0>=*Gx<8tyBS~F<j5c-^z0GeTEK{)qM<z))~+zXIjtHr4DS=J
zuZS)!Wg{XY0@XrL6_6(#_F|zz%tQ~audhQF!=_n&kZ&*wk_Oku|9ew_h75rb`3gGU
z^=!|4eALhP=iFR(PY(*1Xx5MK1jPDfX0D=kLpEDH|B`}@3I~8?wcHs{Kvg}Yz?-U2
zZiuzaXVD(XYu*Np)yRBETh05vv-8r&-=B)wk~8Jk56i$lP~CDe-bB@vm6f2kZ;`?M
z!T|?%cU)h;ZgiR_1$qby`uF#Tjj__`x7fjNNo<EhT(!v$T!3s`^52&Z4G(9OlzfUA
zc(4+@f>B+u_UBLc<m6<_enz<p5<)<~(S1W3F38WwSTMH*a4*WV$&saMUXYA3E<7BU
zszeNb>a(|tX|pqp(_-6W7nUzz0@O-%a#_RV<&xNRNW1Muv(#ML55o^SG*Wq=A91j=
z=g~-d>H7Hi6hEx)?_FNbg>(%Ki9r)qcc&16F7_>Yt$A1@O@hq8czLzo6_z?bIhiz?
zD{BciYckaWeVF|Iy*%jM!o|<t2H>5oEiDPn&0;i<T)*V0soKop1+i*psk*hcFSP%t
zwbTJ9!wP{#DT6&Av#6pJany~cl`aE_CI54^C*ohulgjBvm*q~*S_dI?GCm7tNZ9x2
znYyNCk+}OGOs~10KgSkY-M5a`hvTos?1ByD)dt)lRo>LpDEIRNxZhsu2In;iHa@;s
zYc`9TmR7MvhhL&1qkw&!2|%`v%k9dyd0;qH8*E1kZf~v(2Hw$xxg~ocnxB0mLWWsb
zSR`j>J5~HL#5_1S$12uWD^!jv;CcFqTm+vq5-It_Y@*zNts|9JQZf@ro^pV7jeh?e
zDW&~{on}9Vy`zy-*4EZW4Vr6oP}rSs7Sq6epOvMFOC`!IDJf~%^(G)ChliJ!*TzO$
z`<?D7=oa(+r9>u`P=7o;Jk%gCy~?w*vqie>>|9(LQ(AuuU01#!A;XUUj0mY+N*L3l
zyY9}3bWp-iHpk`tpFWx1U+(HcxC0Dhh51{>$Hm2&pbEbH$?x{+6^D1Wv1hup|0mtk
z5Dfh7MR!k65_6v2^>w9Y=S9H`9!1kLO@K*^TwIkZdp3ss23!db*|c>5W|EMSQnGvp
za$)o4_dWkGGGgK5vq2tLr%sHF>KYn_2M4aEt!~8@uVrOrIeB?Sqj^YP0S-lm4S(Pf
za%^*JHDiw`EiKI%16}(7x;Q`A?eO~->O&>^GOw<VKc9k_jO>Z?u`maRnk#Ro>uR^#
zahKOF<;}^YP9wYr{pjF;o{Ni%q4w@(KTw(Hg+q>Ma4fmNTfiXYpeH?UE<Fp^G7}PD
zpS^Y&EiEnGJUv+$3ywFwKZGvVQ!Uw<m=x@dMC`}%7A&VEJh#WLFOG);OlVVFIxoT@
zU0d6d@yVHgV$j%Vw&ZY}`|%$w<h%`AwKEjtAW|OFA~IQVag*A8TPG(iwPLe>bh2-<
z9H+|ie;`>l+7tE-3>b91x!DA7?dPqa)v1O>M#^062stja2!w@)Q;OUH8b^i!?k?tf
zQloo!d-DJS<$Q2ANXVp&d8ly;_ehP2!PZLS|61k+P&65^rS_{=uavb?BR99FE8p%q
zJ3A*QB^7$#UH<!fKO)e{=>P#?wuEOiG=@R<nm%-~%`@ip>(_%LBXPI4S1w;mXY1fS
zWo2cmkz%$8a(JaX5M*Jl+2mc9+gydIP0Tt?uqj88Y?*+$&ELlWA>N#?V%+aAL`*;c
zgH$<FCIGIZB87&ACg_(Jff>|;!61<rGqI2^kbC6-A%yKlsSD?dGv|R`%F>&hJ{*<Y
z-yBIKaLBaG4LD6FxCM}qdbKl{$aFspibW)`m6cqrt*v_2)&&F`J@~TB(AwsvDnLDX
zvnWORfXkD7U@I8*#_|=yq_n;n@DLan06M$x4(G)-J%E&9fzIDDtO20dHf&TkFEK*s
zxVuIXy}gPRk1Ib|B$dN^!cl{O*kzkrm~C`OqF)L|!Fb=%5s)4n3h3D2=b5bEiG7>f
zIu@aHRlP>Za)-b7*LS}M;xmid6}brtkB+~~VhgpmaR~~(kBp3tj<%ou_+0+HSR)lE
z4?R;;%Ehk0W%i4|?>^1olBd0Wuf?-CU7!>rs8ut;o}vx}N1^ZWA3e|qW#0Qs>;04V
z!hQ||y}e^VuxLOhULp6$$Y)_>^g$qo;|sEYAZ46#Er*~@KaStkZgV;CHU&NT{VVqc
z-rLg#RGp83PfN}+azrle=RtABCD{=ZX%F<3JT4VExiMVIJfjT>-@L)_#}1f1CWtvv
zHS`Lo4Pny@^=@PVB&@Ik1oqnYwmR5+5{Uy)u;K>~SuxK_wK8p!+W}Rno}QkbfWe?-
z={N#@qtSU$B04fM53spjwP_Pr^5aXhXU~e{3MLluadAmxQg<phFJH3+dR(6EyEGo9
zgrZq9M#sebBfa(9919kv#C>jO_jiPen1o~;R9+|O^mXEQ*}JU0#m{G|zB*Zt5}|u|
z=8%@oKHZ&HVbjSeDmN(A2hdOsRBUoq*7aZr8Q=`VS_|onN3ROx4n-(&!$Lzb6;_Jb
zXr>)Ssx7-L4unzc12R*nKWVF~(yOSbxOsV%0BWG>Amzoy?(GEtf?X5Q`wNjHohjk9
zqYC=92vDCs0Pm(_!^&x4Ag~3Sj6nlCfQ@`R&*<%arhbKIeXsI&c6Npk=7#;?ejPyX
zCue3>lGVJBNNKn5-xXQC4Z|io`I9Phdv0@iI6(JcEC=p_WlYb?%Ec7}=~fz5<3!gU
zkJ`u19sp3Ch$i9M1a<T{+0@Ht@;_Qr17QRK(%gg5ih^<bds~GXM9ODzJ?QnJND4l@
z?&i1K^aB%`uC=r@Hg2^Uz>AKJEqXruVLUhR)<IXctfCI4i3mbMQQW{#H*h0nlag9^
z_P6=huV0uDppx35@1t4b2Ag6O^D=;T1FlX<ZqY?WMZ3Fufo0SK{kCHI>4b!wTy;2=
zuQY~?k8*^TH1fgt#eVhUN~4q<Dk0n0CD}J%zLFCXtPJNpj)~T@dHDKkyKdE?lhI&d
zQGi6xY;2g4E$#jSO{vD4n)Lbe=cVP1Mek)6NPd!VtcH^o_)i*fH$`$7km8NNnF(r1
z3V{77s;jG;qZynm)tLltsP>EzDCr*i7(T#HIsmkb^YZ9Jl=C7RvyXtL(g*jT)p>AK
z7%xLKL_<K3CnjKgSz(Xa0p<!dh*Hd1Gw|kLs%st>6ru!LoT0wH9#v{C2Q5xnq1k1*
z1l({2wCiYUvs?e@^*KLM-bkhx6|v<H22n33BU)cy-zkt5AU9S%G^#OE15;59;6Xoz
zl;<S0&VHN`SURLQ8p#hpAn}BL2Mq)^weZI?E?<fXVw>Q!gsd!T$BCi`0Kh1oUqtLZ
z!mz1na-9Ff1DxECycV#ntk7fHrV3ucbj?RD#cLwY3y<8qz1h6Iy|Hj`ipIut<<@|!
zzHi$6s8fJ`o4e7d^i?zs4Gm0yE<40O9KEy(P9rdHx*ebklyWn5|ASxvu(-9grCg^u
z{M4mG321WkKIvD#P?6SM@PtUt_Loc1ut_(?ZP~XNjl8!`%QwRZ0A!fKlQfI<$AOOy
zWlj6ZG006IG*+mBzx}iRH&fg>wd@-+dr|U-ivly@RTu`9(^F-7OnzsJ-chQ5Cw;)0
zl@~jVDX_xcIxlx>8I)?@S7>`IA0>`j{JopII~qPWkZC0XHvRXE$=}-vh*^#vjICnP
zVpXxbp`M;>ESfd?>E5CS%%qn7@z)<eQcWA};jR&Lx^LgU9UK@KSkd;bwde@#*lH`L
z2LvWqq4iiUUso3fgSBc4fzDxhJtLx_ad!g{McU^e2;?wIN=hkex!|&foFpPnbIBle
zwFecNu22IBf&K0Eg>sVyYf_>iXEBfwA-^{M{CN+|6E1~dsLb89dq?WwtAy0lC(6Wb
ztyva+TP(<p^RGNt`MRz@+MBkq`~$%cQ8z&Q%}E{E=fP<q$;rgqlco4bkC)L(=g&GE
zCRq#TJl2OuLavC%$H(vC2{@2qx;N*+5a#adC`D_-&p1kFN9d@ix)v5WD}L+@48al}
zoAkoXz&WZI85tEkO`M#V80_y?{;p>AZJduZl48tWSlDrj-Oa<JYQnv(fzV~C-NDeG
z7z^!wd2d5PlztbNl)P_f5b8Mp72~lnN&`PPWh1u1fsw0&Ob9sqcYdckYBMu4Q-Gb^
z&dw?gIsel*AdvtDg!x2i=qaY(<`~2Fa+hpzQPE}?7@Qy7czUFf>BbGVxJdl=EY<S)
zc|g7<#KqzM+nO9;KCLosims}v+UyR)Ch1oLImb#q)pwHyTf)45j`lBJJS|S&Qju$}
z!G$_bCNxe)#FlmC<bvhCfADBBZKFE^O^bP^uG5?G%e*otb;T}pYRV8e2g7Ym-~FYI
zgQKHvAOu#|)<TtHNYi_B02RRs$5*Dx4gUKq<UG$ITnGr{w{bwogn*$?d&xBdNE-z*
zZgFL>wzfW)eq#9*{N&mBpY6qYF9(o82fNM7+l&5SW%XcwHM1Vo54{_z3>z;pLeg==
z6$JacJ=auePo1u*EaEtoClhcvHrrsQ<!k_EezW6jiR&z<=FY~>j?IUw)nkiDTv9R#
zY=+Zv=i=Jw7P+_^GZM+w_X-?<F2EkKf<o85V6KrvATEwKR)N8$+yHSQAkH+w=U711
z+qh*~_^y*&pU%5<kGr!(#m1sS)H1Kh0okBH=T}!u;^Guwm@)D2L|X$SO0=`FAXOcm
z+L>dbC{!kucc_G`K#=weS38nxIg%U{&CrA3?mj#H_MQc7k+wFgi;Iina%V?QAiDtw
z<aFYw#7dA*&^N=tNF@RB2u#!WwU~#vEd2ZlbImTbWFHa}<(__enZGqzIte_HQpd>O
z>jzZZlCc|d2vt2jLZm!EkDFn`&jc^`KJT*w(UTK+>+3j{7i?lz8l`HS+vv8go&h{M
zE;)ZtR#p}}J^hy=Jo<`gXR2c$bREE&?+s;sh`HaQ(UwvHh<5LjwXM09q9W$6a;Vt7
zu~p`uKR{#a!>Rwefh6L&HZV6c)BHtIF9-^XN_7wD_cBLE_GaoOlVP}Ez@^96Q6S9|
zYHRsM(uE>G$`hq3iJ`C0+(woiX%BoSs$r*E6+oqK{!gQyrLlW)&Mz*;-M>qhl6`t}
z<F}e?QL0VT5m?sR3X-sK>3S+EP`rV7x>AtnK51cG1-DHMcreIICo&=;VrBWPkIb@4
zI38>z7V@NcztEr0r@eH)F3OT+3JIV{Oh79AR4(no3kzmzYilb&9+Hh-Q}J6tAbGZO
zak=}cJ?QD-5n5PS*q1vt-{$!Oy0?j~XXfVShrq$Df{wpWjq3c*obSVg<z+fx$VhoC
zEi7n3_8hj)nKIPU|C{4c@39`}ql42^QoK*ka$d&<4DAY;hX74kR9ws;E<XFKci0eL
zhHVkGrKP0~mgXifueb_C<%^ms@_zd?Eqb8sM4lMo*a#z76FU6Q7)3>A%w`2*_4xmq
zR)N+^6Oz;eLd~R3vlt||ds|ykInw?aT%WPj)YRmFs{i=$BW*JW;$hG%s@vhOz5yA1
zkDZ?as2CVbl9E^Tb6=&zrG$nTffbJi0V)7=GH{|0qF($~lhuF0R)BO=h$#VB=~aMs
z2L}gVfEHib-;W3FDC^unj=fQC>_|jG5p{EQW&@)8G5{!HqK2U?ydiLIl$c66VM?vm
zRS*-V@>%u+*#z?T7@!2YhljC&zwRF9>xu<{kq=U*me&_<>2?o^h(f>npYD7G2oUk}
zr!i1^^Uv0;a)@F-@wOm)3@_q}%j|$i)(|-A%e5q((DiX(ZoKvuO!D_<CZ72EGlEDS
z6ACyU5%W_akS*6Zq&NmyIEax;Ct9#cxd|*zx^6FHB57&e(Srcbu*=Hk0-scF$T>7P
zm=_s|d2b)4E1&2@ZUDy&%q8FN5H;_vMkjDi#|X9`{O1W}Enu;Gp4gskU;~CT=mVc+
z&4%TJ_4Bd`amvZ;E%PIwTW4Y<!C4Tcgoi!Ff1$cD;$1hmK60bV1CSryNA4>PO-0T7
z@BM2G^!10FY%i{^uygxPg&JUKLv`fAjW1_Gf{=ykeFlI6AWw=Q-(4gjQ>Y!tXf=WS
zySHb+?laig*=5KE%77Js?E1gJqEGbCY@hXM$U<iT6{tyYVIdO)<4myzl=}0tHy<Mt
zQ(|HwQo3e1hd~K6rhUD;)IkO5g7rMJvr_;eDM-nomG$+<<vLv9<v1|}gCBUFf-_96
z_5}$LuA$>nhRTJYr}&lt*Ves3A!KU>w6^l1T0ehK7$T(c#Sc8FVlBN6=-bVSVr=Nk
z&=A(MxeM(-6n%`SM+@%XkzgrgguEsnQG>i#UxRF4xj6>TEl*DwhYq5U_0hL6<M|30
z1cd$Zbn$=)jUHit=br{|EikHtySU&>N`7{fmIXQQLdz=@5L^g<e2!>%r-aEYBSQ<$
zXbuh!%V&<wKYR7fyL80?xsWC4Jvule@LiCUApvyR*S<b909Dmql(M(yOpw+2!D_!C
z_<RL=5C<s$7oMJ2R;UthJyT`U?Vx78sZ9m)C*i0#L_?GkB4m^G#Cf3w33l~9B0Su>
zFPa1jx*wVOc`7I&Jv|C2Cq@uP3kV2+%$1#!Q{K(($fiEqcSzu$VMqS3xP{Lu8Zr#@
zJ)^(Bf5$ot2;<PrTU>*|Qu%yt(P6>C!NNWVR&K7Y^i)(-l(KKyzYcgl{;!;jjs(!#
zh<`0Dk6s`4k)=I)#R!TL4*q<()J_3%Rb>$A9E_{bt^SUeRo2%h0`lQOx#49Rk6n%}
zx9zZ#tu1RLCRpC{y}sVwm5E{va^P4%=xgIOpJqoELpb>R%X=pqFtOx)*5n#t|8smF
zL#PVU?Qk55?qiBrCAa{x$)}eTp#8&M{HT$3HXv92Vy#>cbWM>yM^*KRm`G2dvF7u%
zU8lz0-_bndto<MPTVdeL2At<0fi(aBAE^)RgxQhPu)0_y%^eSVYN&DJ+R#xUq|DQl
v=y8sVtOW{<acrXg&whL0e}ALlJaD71&A@r%DKYqt2BfB>tym*(_4fY(ZIU%)

literal 0
HcmV?d00001

diff --git a/lib/cartopy/tests/mpl/test_mpl_integration.py b/lib/cartopy/tests/mpl/test_mpl_integration.py
index 6b528b982..6b0f0a9b9 100644
--- a/lib/cartopy/tests/mpl/test_mpl_integration.py
+++ b/lib/cartopy/tests/mpl/test_mpl_integration.py
@@ -169,6 +169,7 @@ def test_simple_global():
     ccrs.EqualEarth,
     ccrs.Gnomonic,
     ccrs.Hammer,
+    ccrs.HEALPix,
     pytest.param((ccrs.InterruptedGoodeHomolosine, dict(emphasis='land')),
                  id='InterruptedGoodeHomolosine'),
     ccrs.LambertCylindrical,
@@ -185,6 +186,7 @@ def test_simple_global():
     pytest.param((ccrs.RotatedPole,
                   dict(pole_latitude=45, pole_longitude=180)),
                  id='RotatedPole'),
+    ccrs.RHEALPix,
     ccrs.Stereographic,
     ccrs.SouthPolarStereo,
     pytest.param((ccrs.TransverseMercator, dict(approx=True)),