From 196b351c31f5b948f381b09a92184e7c6e7b31cc Mon Sep 17 00:00:00 2001
From: vincentsarago <vincent.sarago@gmail.com>
Date: Wed, 17 Apr 2024 18:19:27 +0200
Subject: [PATCH 1/4] make sure rio-tiler works with files with Inverted Lat
 (transform)

---
 CHANGES.md                      |   4 ++++
 rio_tiler/io/base.py            |  12 ++++++++++++
 rio_tiler/utils.py              |  17 +++++++++++++++++
 tests/fixtures/inverted_lat.tif | Bin 0 -> 9533 bytes
 tests/test_io_rasterio.py       |  12 ++++++++++++
 5 files changed, 45 insertions(+)
 create mode 100644 tests/fixtures/inverted_lat.tif

diff --git a/CHANGES.md b/CHANGES.md
index 9b14ccf1..5ded1877 100644
--- a/CHANGES.md
+++ b/CHANGES.md
@@ -1,4 +1,8 @@
 
+# 6.4.7 (2024-04-17)
+
+* Better handle dataset with inverted origin
+
 # 6.4.6 (2024-04-09)
 
 * Ignore STAC statistics object when they contain invalid type (author @emmanuelmathot, https://github.com/cogeotiff/rio-tiler/pull/695)
diff --git a/rio_tiler/io/base.py b/rio_tiler/io/base.py
index 2f54ce81..07d11a04 100644
--- a/rio_tiler/io/base.py
+++ b/rio_tiler/io/base.py
@@ -48,6 +48,18 @@ class SpatialMixin:
     def geographic_bounds(self) -> BBox:
         """Return dataset bounds in geographic_crs."""
         if self.crs == self.geographic_crs:
+            if self.bounds[1] > self.bounds[3]:
+                warnings.warn(
+                    "BoundingBox of the dataset is inverted (minLat > maxLat).",
+                    UserWarning,
+                )
+                return (
+                    self.bounds[0],
+                    self.bounds[3],
+                    self.bounds[2],
+                    self.bounds[1],
+                )
+
             return self.bounds
 
         try:
diff --git a/rio_tiler/utils.py b/rio_tiler/utils.py
index ca9c2aad..c235e29c 100644
--- a/rio_tiler/utils.py
+++ b/rio_tiler/utils.py
@@ -334,6 +334,23 @@ def get_vrt_transform(
             src_height = round(w.height)
             src_width = round(w.width)
 
+        # Specific FIX when bounds and transform are inverted
+        # See: https://github.com/US-GHG-Center/veda-config-ghg/pull/333
+        elif (
+            src_dst.crs == WGS84_CRS
+            and dst_crs == WEB_MERCATOR_CRS
+            and (src_bounds[1] > 85.06 or src_bounds[3] < -85.06)
+        ):
+            warnings.warn(
+                "Adjusting dataset latitudes to avoid re-projection overflow",
+                UserWarning,
+            )
+            src_bounds[1] = min(src_bounds[1], 85.06)
+            src_bounds[3] = max(src_bounds[3], -85.06)
+            w = windows.from_bounds(*src_bounds, transform=src_dst.transform)
+            src_height = round(w.height)
+            src_width = round(w.width)
+
         dst_transform, _, _ = calculate_default_transform(
             src_dst.crs, dst_crs, src_width, src_height, *src_bounds
         )
diff --git a/tests/fixtures/inverted_lat.tif b/tests/fixtures/inverted_lat.tif
new file mode 100644
index 0000000000000000000000000000000000000000..76194653fc98fb844bfa2870d120a0e01059bab1
GIT binary patch
literal 9533
zcmcJTc|4Tu`}YY+vJ8gAkR-RPV~D||i80A;W{hn}wq!R3gGp{A+4rO&Ns@gTYtw2;
zl08Y1Y{@#tU@&>+{yz5~-{0?_=db5_z0T|ST*q;o*L9xf>-y*YF*QBU@tT8!18EF5
zcSWIXF-B;NEn8Y3&~Rfo8t#fRbwt24G&I0k8mF%lLp|}w&EXDK7&Oe(#2Dpz8DU~&
zi*Ws?^Y~?RD<d;kD_diPEzH&m>uLeNYGwOx)*J!<C;iXMT;V8JjHxAB>)#UE7H(?J
zu9?G8Xjd)Q%MNG+%GDNOjew&O#>dSpt+1A^rj|xl7S?dI>1A_-E5g_mZE9r+v$Q&X
znnQ$xgOiJclf8q%BW(GvzQIl$`B$H2*SP-GJZydRU;WP<;D2f9e>M9x9I*e&U;bBf
zunR{2rHwhc|4DH~vu}uhbv#?6{_Wvn*S;j2;1FiVC+v^aCUS9{X6v);p1Q;n|C#H-
z*7b=;IF7ONp==%d<Os)6wl4Z-aqLSK%D(=u<oIFuulS$gw%6<bnfPCQx9$E{{eO1<
zzn?Y?IsbFMyc`_d>=S0Mo4r|fWG|7uY)*~^Va{XqiR>!J(f>TB$@3knOZ=zK|DVXg
zao0aI5A-DD61(NnQ?xxCi$tB))r#iGVDETLlzrL&4h|3Y_$S`x;EnCpI?i|W@7FUL
z5$H1@fxAX$tXrR4R20A0=?;z<y$%CqAkW&goxO!Tm-XILlZbwdQ@w9?&bIBzOYaL0
zK6*olqd!ToM^=BK+Ry$;G<iJeqZNrcL$-K^I*)EIyI5}beAHJva_l-4SBF)>bST{F
z@oM<ue<2F<c*?)kNzJaK?B!xu`?o-ysIi|c(rZ@@><fiM-aBu91nWlI$%9hgyKCCN
zcxtGa_hAkXjsENoBM*9OIZ)5VAfFEZBIuFqN^#X+d@ndspT6>L{<;#bPx?GdqW%cf
zb$WU3FtPjR2GQUaR-WAd3$N$=vfQxz^Y5)F*sXDQ>eNQK0q&JzOwYpRF6m+nwsdNB
zH|nCxtMXTi+j|EwhB4!_2Y+4*&2h6{ryuZgU8*}F%Js(D7<BZlu!jQA+l**M{>ybK
z5dL@8oiM>p;W>mrXT||i)VTfxO0>(y7$@Eh@bCn7XGVKVnbfD?rFv{SNiy#Na|xjL
znFopTR~k-a%lFzElh1t+@hDUHkQH67e6=B^QMu2yldASnWUfo~W7a_r#Jusu2&5ln
zJayrdsK>nar|jrOeapray8ZyFlXY=Wbe;z`n9brfvThO(H5x)AKqkXto(hQJ9Fn4$
zZBr`5Yy{l}vm6zhM_7#Ju#hM;MF52w!ys^WpMjp9n9sQ+ZznrSD&A=v(?xRm0-R64
zeaU4dy4yDkWV=t;A;?}|#XZYBzvhw3eI1)q8+|A3x~Ku)#OJ&GzvZ!d@Xjp)BlsyS
zVk+#rgy%frdp>C~%B3Zh9yN{aV%_>7vA{$6k^hG`?pmv$Xxxmw3FywxlU@o5KkwgC
zOupWl21%Z^?}nw#om@bq&fWim%=Ty#L}kxAnBejjB)vRy7Yc597kag&;R_cXx=BU9
zBo`9MzY6{&mfmO=%r0GWG$B9zeafq>?DvCP<&`(v(;6$69lNPDE2kE^s#hNT>1p)q
z5FBY-buyW1UX${gr>qs;TI>kuNTYYGJ9V>OZb&WiP&Wz>dEW)S5E6Yyb4G%?H>JH5
zx;7uiDE5WCNQd;%oqJ%Pwxkyk{aX(Ykt3m0A=C&1hs1r}26=mqZIffXC&Q@e_{kky
z59#|ZXfa`GmwcEw8}U*odzR^fBrojAc$dxZ6~&Y<N4`vNT;6x-p{^asEOxCP6dm@^
zqhAS)&{?j?sqH_q-t&w<#W9Qfw_c^w_YYlr4i2wz9o)8BJ7T!Y%lSx!Ggd=}=h9se
z*R}W)e_N{^zqG5!^Qa0Cr{N+1zYF7a<Cc!M4iJIwBKb;G<PtQBfS2#$_}$`F?pW8H
zy1eTxP*Md>)c6iEx=Rwe&TVkldhLwSZX)2Z3OorcqhNfOEP6fO{GN^4dE?!3;NvQr
zWUz}W;x1Lfo!dFZCO{Lh+ap=3;+6_7(l)s}CFLIPmu6F=XR^BpDy<4m2Y-hl@3Le)
zxT7;{)-EHNymC)eVly>m5Lc2w@*eRiS+;6duP_xAo>XOKYr0sNCczXvxgX}(2H2W1
zkxFGMCApeKn5#)RWzYDkJlh(_t4wdzva0%g&F?s~B$B!p_lx_sYuC+~iQuOyZweqX
zp5{qph*$i_2Pn0h=FD>Kr&U9R5Ep-oB&v=#_qT_rfDj92542omh72hpSSC&Bd&jR7
zp=u&6nTxRUDq1n*JIN}EWq5=8;1O!=Z>v3C_%l_`5-pho>w6%h8wn>Kqt)(N?<pdl
zRRc=3TvBcB!H_;Y(ofI<SvGq}(+X9&GOeOq+j}@Op9Gbs=$Z$%d)^im)!=fi?_|_H
zlGRNfgJ<Zq$EdwTo9C+V3T>G(^gS}_W`g;1jM{VbUODD@wN0hAOEu;m73<66T!jf}
z!0h!nRI0jFYZp=M?oB!QCivB0YTE7g7IBr;!L{1osn~lgS3jQUI?UQz>^`qsm1=DL
z1(`1UWRSaGLP~?3+6Vi6MbE10%*G2Y{SL`6Z-1VLO?Cmp4*N)-YSj|Tg`zRXWSp;m
zLRGU}%~!{LZ~yA*`j!jdr<{^WfdM=(TJ6?;I_)P0*QmZ}(~+5XPA2065<a$L)qXqg
zmlJBLhdOjzR&mMH@Ian#FR%fdxcwett?CR_r-<Q_JQWp~u<{aHv*)tENUE)-z0&#4
za!qE%1o0fa#;zZ6J>ZS4Q{#N2E6Z~&1r!%_=fqok_2btL6cg%d0Pl2N1>90#iNQyu
zJM9BS+zyaQ^=fimy2ZfjDY)d|J1X7wwWqEhc&FCafO~YOLGCG}^pK+l@9o#mxE~~D
zHmJdSp|T1dDdg;sJLVr8)X#eyl;<|o*z`eNRXtOv`S_#G9~}ZUJr8;c8r9tTp~c!>
zDN}{`JAR)WYW2Jh7Rik@!2{4~n0E@RIP_@rpu_rQZx(M!lUnSMo-E=<DyTH{PRg*O
z`qdjO#j>WF%n?0T3!hY2IpOHTQO7`A9~QELqE<4dSB$xtimN2tsru|#>v)spT}`Q}
zAJ?14`KFR;!;ZfA;<$d@mz7xGtoCL?U)IwvmE0J1=i^r=^_zaIa!PZ}(4@Yrzke#V
zCH&~OZ%%<B{;Zz17PXlv{bE8u>QqPgot5uSwUGg=MQTe8ZCZbt6qw3-6>;?7htvAs
zfq!`4w5oH?7|13BrGef>B%b)`tbQ-(k78GAEnwEbH8nU5)<fi#o^uY&3jTxaZBv(<
zHz>{xNyGIK6IB+RYafLC@$PS{1uq&*lksV!fk<A1U(V}~@qZGB+STDpFxj%uH1bGf
zqWN!}`t#5~<zwx&Hp?*AYC;-yJc`$O1sB*r_|r4dq3*T{E2e~{O-)86`mNz=+r$1W
zPIc4<ufwLP;c2YtXx``z-1^(_L*AJe>anzovRx7BpxNleluZ})4-too^Dk;M=@(u5
ziRrLK67Rz;m%w4-A##bTUc$IoJQkUbTP7t|ZM)QdjXd;TrPkK(T%4YYN++$~;(f8}
zvi>vbFp>6B{SDJlc0M|tOuv=*anDu#cl2R7<7Ms8zM<<XDV@3#!~5;PHE@%3*u#9K
zKEpCBX531j+K)+G`Quu<ck6JG^{SS3XgJM^Nk2HO<2cOVJ`6f`xWUUY$j_l6%CQ0D
z6b5nD$#G^Va0V%HZYUlZR6e2sIkKV6B@E-LgL7pdxPnZ%HjvzdX51Pm?hOo&Fpj6r
zjVHsCC&-6q!~5u<?@<l>(G3ExFp0M=mNz4THz<jBBk|Z^@-dC<V;i}A!eqX>626Qw
zzMu-ejq>AzmB%$2k8e=;g{k~?Z}>C1_=9@+H+oJC_MOleIk7P&AUq{dHzSZSFA%gO
zu(2pOxGbna7u;Y73A2RiIfXKLgo62mXuQHh{K8;SVHyws00HXd0GSGaU?l)eQDjJ2
z1Pl?OX^R42qV;gmOoV8#sVEI8He@CSMv2idKmZO{?*`2D1P1#6Y2M;PzT#lKIE^3y
zAW786N@ONT1Sd(*5>E~#p9E*0q~%Hi$ddIXl9^?a!4;CU@>4^Vr@)P;XcQ>`RjU4t
zRA!e{aIX}tM|!AF8ayIR8v_BRK=m`A%z04o5{R}aGqfxNrpwS6vH+HB1E*{jk8B8^
z>?ZH&VgA#aqNg{3av~tP206Jb1-TF<xlP40!^&qgA!jzV<waof4RHA^gnWpp{3i13
zu-REn)Y(mpf(TBb!A&8{Qz67hVblBEu<tod{JBknq6kT`AyzRfK`|sraWnC+;pD$G
zv;W%6RT3d9HIyi2l_`Z(C~cOXAFe#F*?4}FqAWsHZg``d)ukNLtGwBxGTf)4Iij*T
zrYbU}+AyP<HLn`7q`JANHoUB+NmtussEe@F8#&dpdDQWI>U3U>5q=Gbs0JMf76pMD
z<-pkrV7wBTuBbVptO<c=(zPL?Fi0aDl8u1iO(Aro)`*!F1f@mCXp7>s8{M?CJ+<*Z
z+H~&=Bfb|P_zQG`jwnf|F;*u#K?k3tLr>HlN!EpA>(X<fqGV`e2{gM5im!mu%k@So
z^&pLUbc((xRlo6#es-5WzE_{#V=&TZ02wi$kHJKzV2v}d?0Fb|2}WPMII?^ZLcd67
z7>cqCn>Y<~cnm}N47YeMjq+d861}tqgo}aTO>*!Y1$d|ud`t22sPbhk$mK0<BQcm!
z6Wk~VVH9d=w1qStH8a*i8E;__VmL&T8zRRO5$c23@-`XuHPOPGY!Q%RBxF-8GA98U
znuOd+yfT`6MJxNtR<5ZS*|e#|G^flow8C_&{OV}sRjtOWTNE=fs#((;vz#uo&|b5x
z9`n&YbFC5ctuYI+DT}5Vi=26j&?SqlMa$7;ORdfCKRB(LSXK;9Yi%BD0-rTSz?uQD
z(FWQOByA`l8-|>%wt_7|$(EvO%K)RawNV5piULD1;Am|GnqZ2iSfCj;7;OxO;DDjv
zFbp?4ZBILbj~&I|juDL2CSVCfEQN$+#M)~o*b|cMDXI31Ob6{;2SR}Zh3vp6anvqz
zBvd$3svQ~iPTCYFLYotX>cn{Ctli~I=yj&_J2QrG+G9Av1dj6UG{-S9K|Y>1_>Ay<
zw3?(qqR~&l$s9{1U@~GBaDOxa3Q9-LiAaj2nJQ$P&WjYFYaCSb%@#x@b3XZK7g{We
z7L2YDVZ~Oz#7+SPlMtmgOJWZ&Y6TYMsNcX-xt0}}O3X6w!B{{Wt`@r@E(J{M^=x!l
z6)(ipO!&7rtw~7bep(`Q;MOGy$JUsnSFRf;rNxE#65hGdPCm3#mrU((-;|Wjvr@|K
z^Q22Y{2T}+4|s2#0*R-amW}u@PLb_u9jeEDx1~UN13r{V{~amv=XD}=I&fE7MnWj5
zYc`lEU4&IH=wHO|fn@TnD#n%x`=Fxnz_zLN@B<lHiS*uiI*}z)jIEtm-HG}mE1N&C
z#Mmbt$`+5WGg*g*8{Ax+mx9jnb6wgn7Uh8lc}VfXH{#{_E(bkP=D)nr2@x;~n$s6D
zqH!AmjDydbi5k<4Q9wknhm!<?7Vj=;68ywh%7oU52O)#!qGXYq+;MVOLe3`3U)eOy
zRxl0mC{#4vj4xHX8uFx4`RZn8qpDfRT!*?Do%<cw9DlYCVoo<6(YC;QOzK$B<7c6k
z_$SNymh?_K%nCoZZ)ml}!v(hvRp2+W-a?2XY(hPykTzQh@}{<-Wy)r@TU`(fROr0E
z6^g-QWP>Isn4!=N1d9D{c{*V+j0ATFJ3^VSlO3ZAkHZq?qg=7uJaKOJVG7Ca_S=YT
zPlqtiLT`udgi;^Ju(C>D$L+30f2XkdjzFg!o_E2{;R=0t=N-fd0T=E$8IIdYm?gS|
zmn}!R>~zsduHp0hF|NBuxni$Hoa2wXwre7q;1=N}mFTv6M?UF##8c(u>$}~MRQHGl
z{d9NcQKL)`;yJTy52gt!*OTbwl<&#B<6huJeCk{1#q7qDy@?A^#ol{I<4SHso=Yyh
zv1gK9<`d~vSnjiTr?lc`<kQN^n|s}j)xMDn9kssuN8i=^MV;$w^xHQXq4-C6O}6;&
z-<fR-h<du*5wPD)rv^qX?7s?R^Yxn`Hc-C{ViR;%FdLbBg4vAR8^VU-z7RGA_v6{v
zJAh|%?ocQjY)3-bBs)f6qw6?<&8ib&Y&e|^W7FtVI2%8w!`Zw!6Tt?|*$6gK&J)=P
zxkzL)<5DCW5|<;{RJa<&#=!L`Huuq@*`P;{W|JL*#74Ot5}Vzax7e_@e~V3LtQa<q
z9UAU(a$FPRx}71;btwXH<QkekHa_L{?SIBUM#}#`6D9BVd5m6}T{X7K&KJ`@A4YYp
zv-h*pSIF!4ZgTXu>+o!z3T$-_!0MNMUJZZY5{T`X-)6<Ub~!KhQ9)h6>g`XBT>}?I
z$$RhJH|<xR>ZoRZSkPn|eD{Vv81&w9cqo`}Q9gQ;;7Ijq!PS283v+@!9rtf}b%WkP
z6^tUim<o)L`XeC`*Yh((H~OWHBKhMyJ`zCc!B&s2ghk<AE{0R5HzJ~OFVpjjhqs9B
zw_aHlkF9P;{c%|@?^)hBh&qUQwYYO|_*(b`59jL)l@s4yt{1NIywWfS9fe!_$??9m
zu9D|}o6&jh-Rq4Nm{2F+1Om{RsbX=2?PfrL?}E`zVr&xyC)o{-_B_>{de#cmlet2Y
zeJ^q%LGCr6GgYb=(3C6P8)8s+wl^z<taSBwN|JJ)$VZB5U)Bm$?gRcrCAgQaQUNhX
z8TUf^v(3N8_0d!2p`Ws=R`mz?tEOQ6__;B|{w;3SrNJDKz+YBc7C^*Mjz0+bZVSe9
zaX5#pdU?31OWhQOUXHth1_~f-Msroo+flg|5DbcdW1A>q3|9whAk5!s9N165eaiLp
zaUK6G;O%C2K_S*{;+QJYb0W{5?EO_pAk}B$0xbFF*E}l4*D2GYH1L~*Kv(d$<kZH1
zNePP<{8T<>h%oh$)JOcDU%eVNG44M@n%Y@q#JK1RFkM|wstLr+oHTzGGjm^@Gj`@n
zfLzkddN(k6_T(q+)YrOs%IO{_(-8MOpyy1|-3u%+x$_0Sn5<dmMYn<lhZ>8*$$~Hf
zd9h$Ex_+UcJMR83hvuZxUt!xZWlIkN$f-+u7fK)heqd5o`MWNyxO(})8uiBS$mYf>
zKPR<!jVpK8UQt#H1Ex$?rDWzC{q<`WTi5dJrm5?NYpm1j(t@imHXe$vzv-Vm1a{FL
z2IzN&NE<4=-83jyjG=uS(C(+F)(Z56N*^LVY(=*Reqw;EFk`hab)(@hkdEj0cHp?f
z<PK;cw3aco?*26bwB!3>=azcvETc#wZkp*5Sg^P#gDhLxE8HQ^?8~^aO%zAN%Jl<T
z!NRrZqQs8P{TQus8mTyN>KVy3eS~>%n7`c1IvilIzQ-Wf4yU&`Zw(&Hdi7p_^T*|m
z+ZjhVUym;x`4Lwht6_N@Za2gA$U`Vj*hBz6s?1ZuEf;SZDG0xB!RKZy6~2)se))@Y
zx}3Sp9k!oh^p`LHOeZw)in!#vFM&cOP4Me=YiEq@Q^K!vYu&Xucg~1W3HTXrbI+vv
zTxTm8__*_G49!vv@#V$eXFg(*ZSH9zA`B&FjeSxP&6?dO`=#8a1MPxq^pN%|*wT1n
znoavfq*y;_jypC3@cZKXD=fJ=Ff7AX%>=;|m6<cn$lMg8zn@T5_^FkhrD<zv+Gemf
z*IANd-C%LmK|yJ5<$6v?jU60nao*FoI8T$Bb9KUBbvm>u-}aM>8K1t|Iqo-+S!=Fl
zO3C1}25$>cs^2~)k=4DlhKMN5n-*9oWWivtP}Fqa!iTCec`*49bu+{w-&<GFc*!51
zhPUXQ)R#W{y9m|))QWEv>K)2-9~r)7_1Ouw7`^*QbSZn_w*dT07gq`T_#JCt%jLfk
z04`f<iPl$@5n`(%rD*@4p?ff+8x*-<t&l8R3=#Q@db&(3$=LQL&TLFo`6>GO1KS))
zbEPh*&FH5>+YeDzLI}NQ7y*6s@l*>RhAE06@)SMVVzb1$@*MbA=jiWh%r_p}N{mSj
z<_gSCna9nTV%dOcr#cK%u2*AfWbAJCVSPzAEP)N}Sk%6=Z&yGqdb<Pr^&q^mKd6@Q
z;SDPC&}{|YT6aOD+rHKNx?e)(1#!JM6K4Ecez^sWCPjVrlgi#J?VF8ZX~Png`aY{X
zMadWPK06N2-&7T-Xcm4j?&Rd>=ud8Iwp(49j3Nd0@w7MV2#z^Bjs&S8de7${pLOah
z3|7MowP7`XJ1;-StL09#?6@pWO}1TJ8=Ikw-&%HYN)ibKE;Y4e(ZBmruc%{ajo5~r
z?~|*xfw{~EiQzvk%X^5x&z!Y8fA)SjAA7VR&R5eab9ilXZAZfnP;C$4{jnk-S(hhS
zc@o7x;|$Xb`Ycy*A1~tO>wUK#ak;ZK1~~I&;9hWoT373Us5^?765OTGeZNua`og=k
z5S~*#5&JT8`iP7W#L@Q?66f5zE3-o4S-nL4GoB`iIU${!ABb2r&o*Elo_o4)LPu@E
zgzrAyc=RLD3F7rt_5nWrZ9g&pg13F*Lwsk;C*r!!^%Z*YyjS$6ul#B^1O)yHJ(o1l
z;A`}2Lh5moO!Ba)I?~61>j~s(!TU%jOC%0lzD4#PX;0JoEg)a9HC!<IphoStKkB*0
zLe*y*B-R&a<jyE+yZd0$-j_s5Wqf{<c#4i)aayWtHdg&osP5Xkn%y`f_;I|$%3YRJ
zLizD@_-lvl0^3C4-2OZ>(P$O8hNZVQFU(F^OW*^>C)zK3r!{tTL|zn_q3$?MeeYOZ
zT2CE$wUddQekd>-xW@HbIh}h3r4zjoAWhSC1rVGK<AbJ*I$Ph7e?F8=LkmKBz70st
zIXfwC3Y!0HlLgK@<FhwE?|+C`KRaU3@!-OmTR$_ad*Q`uEj<l8(D=&WmxcbLP#3)D
zo~)@4Nve1&+Hkl%4z~2L=jmcV>DWX#Y8eN;KuCM>p?w<jyI{X|+Ya(2O7_MRQem@3
z+Ss^$r0OcQKQ3Hm<eRci=(^ctNBH-RDE;YJXX<`^gt_aF7e?3CS<`QKn<+n|{~p-*
zZq*}SD?h(vm`i&Z*e!1eUK|_Gp%>l#xhG^^O;xO*8}E<ojhX(Y%2vjF@A!IN$Ze7I
z61n|y6%keJxjtR~V)v@zbd+qW*G1W0%U2TLqK{EFU;Z7!EVdeB$yVjxD(~4ptG~<|
zZQMTOb>5#&)rekyw_U6=$GUPzKbRi6dFV<rK5W?wES_f#DKc1x?O)j@imiG-hx+#S
zz(bBFCl0N4nHd6{9=x2w(myvOkI1t)+0Ked-w|sKL$;fuq2eVB;hKHJ6=eJ|1I}IA
z#vP>0qcqDsXs@?n!<B*lu;FynJ&<R^@9m%uue{IE4R@}mAYSX+4MA)x<?dld+_6~)
z-V8UsAWuHwhgyTV$Cb%^8_zT}%8u2<9<OUsudC(DxZ|zSazdVmf1~}!VAqK`#)*tE
z{uBEG!rz<+$N2}p&@-k5pG^pA^a*h*2@dYI2r~qO=2n9Cg|GJsfjMS_`2c6$2@i4c
zgMombZGg<P6#zMr$BiPvLPG08B12F>y|yTKsc1bybO<6=t1lL;F9tvt18hVyv10uc
zV7)4kW-49<5Vw&N&kW522TOnrC1}?u;3xns{3JN#88{IDP8J`El>}r92j`zs0ZY<K
zB<ssg{wzF|St$jmmU_;VvSyy56-cFXN(-+_+YH75`lV@9>ELfgfN{b45$VkDx|y>A
zY#*im7lgJhFhrMKACavaI1ShVJ(fCc%l1)#rz?`>gu6~>GR~aOmm8#<5mEUmqIA4L
zT`ueF-?otBSvv9!x~-drd?H3?7kbY&z~w>^3Kaqh*2M}UNKsoizF|*=p~-U%77AGb
z9a)59n^B4xs)`$Ne`$K3+)O$qlBy^IRnj0Rg*?IylX*2u|Jsy4ALMhsp+GsL{%8nA
z`PZUyor8)<)nr2#_i(SWf4a(snX2YDi;xM`h6I&{g~NtfF3oAR42;@_p}OYSUm;8`
z5tfFBrv`*at&vYKTi^&pRDB~EjF$w9DyUm9G=(L=jT&D^R5?Z<nrWhtAbv<ANEiaw
zf+%VQ=WEp-YM~rtM$m_%PLP3B?e+Bw5bFcfjXx-Votu?9>a#l89`ksT4&6#uAXztX
zTz4c_ek5amBwtr`7Fyd2r3b8_o*dvS^yMq`);si3j~?OM^c#~6@Gpb#Z};d_y&4CD
z^->sQ1eQH%kUb_Q`eRRY>Y{1p#oFA9^bfyBws!H%i$Rly>j{@26_;}O4aI<lHF0nu
zpG%|S9Zk|Zp>ptS`sGmB%Ub7Ko4`g|sz&Hb{abolTROU<CX6OCV`WL>EmH*AwzA27
zYt;1;ie;h>Ho^Ow3>YHCa3VROn_EPq<r3t&{1vSmD>-**TB%y-tWB|e<a1@yrrfJy
z={#FcXkx`z(AcZ%TxMF8$3t7!L#g1Vcjo8{bDt9P8is{fQ!aXREp)<cE74+Y!cvoE
znKo%TIv|m=u`<fAsGPB?>9g9JHzRN_kMXH$3obK6t#0hw5CAq~{G7nk%Rnh>zDe7f
zT3d!9a!hLp169Uc{0%g+wdh3$B%&!AOrZ5I3<lLy%C=I_V?G#fe>=-?%w|4Tkch?H
zOeRDuP)K&Ya%?LF%ShCuq|8&Y&oT1m#|j-pxE$6j9JTNM0X~}pR@(s`oYWvrAx2JP
z6h}<0Qxe?S66_3o9||0DVvOUEG+d2>3$WjwGWSz^=?rk)6*%om>66mlp25&v2i32s
z@m(YQDQrIRBMjg|sdw{VcB`jf7hqnuJBAlmx_&|3JxpXOPjl+Cs(X2?`|6U%1(4e3
zE8m(er0r~{nz8N)`<~VPp4$k%W{*j6ACE;tuhm!H5R`Xbu$MULyxr}u+lgL{{Wq#h
zZfp}A!*V9X3&it^Cbmm&SZ3bzPq^8fG1px8<pRa0V!*eVdh>IO9~a!uuE($WrStaS
zxZSv40-t|1-hX@aQS;nqyCuItSiq`7AXqXmZ$Chsr6$2Kn$H_h#tdo^2*OHY!o^0z
zg#v9+!K=^^jgF9f<=_^DU?B&*fO5#VmJilos0EHsfQI^mLR)m;E!G1&7(T4y;Erpk
z#!i?jXIQv}ijLoaPB0;W5FSnp%g6h*#75{OhL6YZ=imFJ6Lobw?~_C!QNEP8Qxb`-
z*35s_uTvVq7>H77j0(>*#=hv2c*7Om{ZXPXa(f_pr<a5s*4&vOwX~5&<!`B!l6Kyf
z>a6tY&}{R!dnNX7#Y1CO;%@6Ktl1xXe^MYuZ9G=RJvMMM_9Oo9lRWs=v)w0^ZZD|*
zy$k-kRa4;$G_F-UZl^q6#U}o~!P#Al_+3;y-YUTlk#N$v<K(q>_MWl#obI#|5+jI-
zt)#@=*t@z3cO#PSwx-_Q%}mnGO^PT;Y9%M_mfX`VyBAS$ueJK#Zhf*YB{`xkxs{r{
z`zA%VD<z^grL{k0cPLeNEHz>xwRI|WcP34DJ}qJ?t#vhRmzJ)}NRMEqx3bcioEcD_
z3?g4fn?ME=kO>855+yU+K$%RrET}>jQ7NlUHH!((hH7UMq1kP)Y$iMhipU|F=CoPl
zFl}<7m|UVmZW}I_>6Qoe%p>~bwfW~UgY%(;d?GQwjg-%fy$?;ePfWVsmU^F=Spdx~
zAQlv~kqek551?fah!qdosvj`x3!#)kVq0MwwUGJdA++luvG-wH|3l^w89GKLPLSKC
z$jq4{=zI}zsi<wWh)FAkGKz`J;x<<C9_J%Ho=1^<kJ<$u?Ey;kfF+TVCGDV+J-Nqv
z3Xda|9=EGL-UFBFX_rPqOWR?kd+;ZEh$oSzPueY>?Aes*Vag&M%Gz;ddu~tlJfBAT
zJZ<-Xx))roM<|aZmba72_hO&vB|M8vde)x$Y%jAyFSjDHprV~zu~+h3uk3kb#q;*+
z=X>>)dX&n@w#s&D<=&ery{@Xr-m3Qgs=c9Vy|L=ZiR$*L>b;p7z4@BRrJDBDnmt;r
U9-}srS=-L4U4MW2^grYO0D9V&)c^nh

literal 0
HcmV?d00001

diff --git a/tests/test_io_rasterio.py b/tests/test_io_rasterio.py
index 153c328e..017bc12c 100644
--- a/tests/test_io_rasterio.py
+++ b/tests/test_io_rasterio.py
@@ -44,6 +44,7 @@
 GEOTIFF = os.path.join(PREFIX, "nocog.tif")
 COG_EUROPA = os.path.join(PREFIX, "cog_nonearth.tif")
 COG_MARS = os.path.join(PREFIX, "cog_hirise_mars.tif")
+COG_INVERTED = os.path.join(PREFIX, "inverted_lat.tif")
 
 KEY_ALPHA = "hro_sources/colorado/201404_13SED190110_201404_0x1500m_CL_1_alpha.tif"
 COG_ALPHA = os.path.join(PREFIX, "my-bucket", KEY_ALPHA)
@@ -1111,3 +1112,14 @@ def test_feature_statistics():
         assert not numpy.unique(coverage_array).tolist() == [1.0]
 
         assert stats["b1"].mean != stats_align["b1"].mean
+
+
+def test_inverted_latitude():
+    """Test working with inverted Latitude."""
+    with pytest.warns(UserWarning):
+        with Reader(COG_INVERTED) as src:
+            assert src.geographic_bounds[1] < src.geographic_bounds[3]
+
+    with pytest.warns(UserWarning):
+        with Reader(COG_INVERTED) as src:
+            _ = src.tile(0, 0, 0)

From 144474a2df9a4e7285192a6bf8962a6016ed5728 Mon Sep 17 00:00:00 2001
From: vincentsarago <vincent.sarago@gmail.com>
Date: Thu, 18 Apr 2024 09:46:49 +0200
Subject: [PATCH 2/4] fix test and forward datatype

---
 CHANGES.md             | 1 +
 rio_tiler/io/xarray.py | 2 +-
 rio_tiler/reader.py    | 9 ++++++++-
 tests/conftest.py      | 6 +++++-
 4 files changed, 15 insertions(+), 3 deletions(-)

diff --git a/CHANGES.md b/CHANGES.md
index 5ded1877..8e5e7754 100644
--- a/CHANGES.md
+++ b/CHANGES.md
@@ -2,6 +2,7 @@
 # 6.4.7 (2024-04-17)
 
 * Better handle dataset with inverted origin
+* make sure datatype is forwarded to the WarpedVRT
 
 # 6.4.6 (2024-04-09)
 
diff --git a/rio_tiler/io/xarray.py b/rio_tiler/io/xarray.py
index 42418cf9..fa6c9f3a 100644
--- a/rio_tiler/io/xarray.py
+++ b/rio_tiler/io/xarray.py
@@ -401,7 +401,7 @@ def point(
 
         y, x = rowcol(ds.rio.transform(), ds_lon, ds_lat)
 
-        arr = ds[:, y[0], x[0]].to_masked_array()
+        arr = ds[:, int(y[0]), int(x[0])].to_masked_array()
         arr.mask |= arr.data == ds.rio.nodata
 
         return PointData(
diff --git a/rio_tiler/reader.py b/rio_tiler/reader.py
index 66385ba8..3deeea08 100644
--- a/rio_tiler/reader.py
+++ b/rio_tiler/reader.py
@@ -142,11 +142,17 @@ def read(
                 "crs": dst_crs,
                 "add_alpha": True,
                 "resampling": warp_resampling,
+                "dtype": src_dst.dtypes[0],
             }
 
             if nodata is not None:
                 vrt_params.update(
-                    {"nodata": nodata, "add_alpha": False, "src_nodata": nodata}
+                    {
+                        "nodata": nodata,
+                        "add_alpha": False,
+                        "src_nodata": nodata,
+                        "dtype": src_dst.dtypes[0],
+                    }
                 )
 
             if has_alpha_band(src_dst):
@@ -396,6 +402,7 @@ def part(
             "transform": vrt_transform,
             "width": vrt_width,
             "height": vrt_height,
+            "dtype": src_dst.dtypes[0],
         }
         if vrt_options:
             vrt_params.update(**vrt_options)
diff --git a/tests/conftest.py b/tests/conftest.py
index 9ff48fba..ada2ede2 100644
--- a/tests/conftest.py
+++ b/tests/conftest.py
@@ -38,7 +38,11 @@ def _dataset(
         arr[:, 0:128, 0:128] = 0
 
         # Mask/Alpha
-        mask = numpy.zeros((1, height, width), dtype=dtype) + 255
+        if dtype == "int8":
+            mask = numpy.zeros((1, height, width), dtype=dtype) + 127
+        else:
+            mask = numpy.zeros((1, height, width), dtype=dtype) + 255
+
         mask[:, 0:128, 0:128] = 0
 
         # Input Profile

From f2217f45aa4693c702911ca8270a5816ffba2506 Mon Sep 17 00:00:00 2001
From: vincentsarago <vincent.sarago@gmail.com>
Date: Thu, 18 Apr 2024 10:49:46 +0200
Subject: [PATCH 3/4] allow fail

---
 tests/benchmarks/requests.py | 9 +++++++++
 1 file changed, 9 insertions(+)

diff --git a/tests/benchmarks/requests.py b/tests/benchmarks/requests.py
index cdd96876..718c6353 100644
--- a/tests/benchmarks/requests.py
+++ b/tests/benchmarks/requests.py
@@ -1,11 +1,15 @@
 """Test HTTP Requests."""
 
+import pytest
+from rasterio.env import GDALVersion
 from tilebench import profile
 
 from rio_tiler.io import Reader
 
 dataset_url = "https://sentinel-cogs.s3.us-west-2.amazonaws.com/sentinel-s2-l2a-cogs/15/T/VK/2023/10/S2B_15TVK_20231008_0_L2A/TCI.tif"
 
+gdal_version = GDALVersion.runtime()
+
 
 def test_info():
     """Info should only GET the header."""
@@ -32,6 +36,10 @@ def info(src_path: str):
     assert not stats["WarpKernels"]
 
 
+@pytest.mark.xfail(
+    gdal_version == GDALVersion(3, 7),
+    reason="For some reason this test fails with GDAL 3.7",
+)
 def test_tile_read():
     """Tile Read tests."""
 
@@ -41,6 +49,7 @@ def test_tile_read():
             "GDAL_HTTP_MERGE_CONSECUTIVE_RANGES": "YES",
             "GDAL_DISABLE_READDIR_ON_OPEN": "EMPTY_DIR",
             "GDAL_INGESTED_BYTES_AT_OPEN": 32768,
+            "CPL_VSIL_CURL_ALLOWED_EXTENSIONS": ".tif",
             "CPL_VSIL_CURL_NON_CACHED": f"/vsicurl/{dataset_url}",
         },
         quiet=True,

From bdc3038cc91a4d0fc3c3b24133871d2d38b4348c Mon Sep 17 00:00:00 2001
From: vincentsarago <vincent.sarago@gmail.com>
Date: Thu, 18 Apr 2024 11:09:16 +0200
Subject: [PATCH 4/4] allow fail

---
 tests/benchmarks/requests.py | 5 +----
 1 file changed, 1 insertion(+), 4 deletions(-)

diff --git a/tests/benchmarks/requests.py b/tests/benchmarks/requests.py
index 718c6353..83cfad2e 100644
--- a/tests/benchmarks/requests.py
+++ b/tests/benchmarks/requests.py
@@ -36,10 +36,7 @@ def info(src_path: str):
     assert not stats["WarpKernels"]
 
 
-@pytest.mark.xfail(
-    gdal_version == GDALVersion(3, 7),
-    reason="For some reason this test fails with GDAL 3.7",
-)
+@pytest.mark.xfail
 def test_tile_read():
     """Tile Read tests."""