From 840b697f037f1a1d75865809816adb43def48415 Mon Sep 17 00:00:00 2001
From: Mary Isabelle Wisell <mwisell4707@sdsu.edu>
Date: Thu, 20 Nov 2025 02:52:54 -0800
Subject: [PATCH] WIP: Add medical image denoising demo for anisotropic
 diffusion

---
 .../data/download_cbis_ddsm.py                | 112 ++++++++++
 .../medical_imaging/data/sample_mammogram.jpg | Bin 0 -> 21912 bytes
 .../demo_anisotropic_diffusion_medical.py     | 205 ++++++++++++++++++
 python/demo/medical_imaging/requirements.txt  |  24 ++
 python/demo/medical_imaging/utils/__init__.py |   0
 .../medical_imaging/utils/image_processing.py | 141 ++++++++++++
 6 files changed, 482 insertions(+)
 create mode 100644 python/demo/medical_imaging/data/download_cbis_ddsm.py
 create mode 100644 python/demo/medical_imaging/data/sample_mammogram.jpg
 create mode 100644 python/demo/medical_imaging/demo/demo_anisotropic_diffusion_medical.py
 create mode 100644 python/demo/medical_imaging/requirements.txt
 create mode 100644 python/demo/medical_imaging/utils/__init__.py
 create mode 100644 python/demo/medical_imaging/utils/image_processing.py

diff --git a/python/demo/medical_imaging/data/download_cbis_ddsm.py b/python/demo/medical_imaging/data/download_cbis_ddsm.py
new file mode 100644
index 00000000000..bb4f108df28
--- /dev/null
+++ b/python/demo/medical_imaging/data/download_cbis_ddsm.py
@@ -0,0 +1,112 @@
+import os
+import numpy as np
+import zipfile
+from PIL import Image
+from pathlib import Path
+import argparse
+import shutil
+
+def download_dataset(output_dir="data/cbis_ddsm"):
+    output_path = Path(output_dir)
+    output_path.mkdir(parents=True, exist_ok=True)
+    
+    os.system(f"kaggle datasets download -d awsaf49/cbis-ddsm-breast-cancer-image-dataset -p {output_dir}")
+    
+    zip_file = output_path / "cbis-ddsm-breast-cancer-image-dataset.zip"
+    if zip_file.exists():
+        print("Extracting dataset...")
+        with zipfile.ZipFile(zip_file, 'r') as zip_ref:
+            zip_ref.extractall(output_path)
+        zip_file.unlink()
+        return True
+    else:
+        print("Error: Dataset download failed.")
+        return False
+
+def prepare_sample_images(dataset_dir="data/cbis_ddsm", sample_dir="data/samples", n_samples=5):
+    sample_path = Path(sample_dir)
+    sample_path.mkdir(parents=True, exist_ok=True)
+    
+    dataset_path = Path(dataset_dir)
+    image_files = list(dataset_path.rglob("*.png"))[:n_samples]
+    
+    if not image_files:
+        print("No images found. Please download the dataset first.")
+        return
+    
+    print(f"Copying {len(image_files)} sample images...")
+    
+    for i, img_file in enumerate(image_files):
+        img = Image.open(img_file).convert('L')
+        max_size = 1024
+        if max(img.size) > max_size:
+            img.thumbnail((max_size, max_size), Image.Resampling.LANCZOS)
+        output_file = sample_path / f"sample_{i:02d}.png"
+        img.save(output_file)
+        print(f"  Saved: {output_file}")
+    
+    print(f"Sample images saved to {sample_dir}/")
+
+def create_synthetic_test_image(output_file="data/synthetic_test.png", size=(512, 512)):
+    # Create clean image with geometric shapes
+    img = np.zeros(size, dtype=np.float64)
+    
+    y, x = np.ogrid[:size[0], :size[1]]
+    
+    # Circle 1 (top-left) - bright white
+    cx1, cy1, r1 = size[1]//4, size[0]//4, 60
+    mask1 = (x - cx1)**2 + (y - cy1)**2 <= r1**2
+    img[mask1] = 1.0
+    
+    # Circle 2 (bottom-right) - slightly dimmer
+    cx2, cy2, r2 = 3*size[1]//4, 3*size[0]//4, 80
+    mask2 = (x - cx2)**2 + (y - cy2)**2 <= r2**2
+    img[mask2] = 0.9
+    
+    # Rectangle (center)
+    img[size[0]//3:2*size[0]//3, size[1]//2-40:size[1]//2+40] = 0.8
+    
+    # Add light Gaussian noise
+    np.random.seed(42)  # Reproducible
+    noise = np.random.normal(0, 0.03, size)
+    noisy_img = img + noise
+    
+    # Clip to [0, 1]
+    noisy_img = np.clip(noisy_img, 0, 1)
+    
+    # Convert to 8-bit with high quality
+    img_uint8 = (noisy_img * 255).astype(np.uint8)
+    
+    # Save with maximum quality
+    Image.fromarray(img_uint8).save(output_file, quality=100, optimize=False)
+    
+    print(f"Improved synthetic test image saved to {output_file}")
+    print(f"    Image size: {size}")
+    print(f"    Normalized range: [0, 1]")
+    print(f"    Noise level: LIGHT (std=0.03)")
+    print(f"    Features: 2 circles + 1 rectangle")
+    print(f"    Quality: HIGH (no compression)")
+    
+    return output_file
+
+if __name__ == "__main__":   
+    parser = argparse.ArgumentParser(description="Download and prepare CBIS-DDSM dataset")
+    parser.add_argument("--download", action="store_true", help="Download full dataset")
+    parser.add_argument("--samples", action="store_true", help="Create sample subset")
+    parser.add_argument("--synthetic", action="store_true", help="Create synthetic test image")
+    parser.add_argument("--all", action="store_true", help="Do all of the above")
+    
+    args = parser.parse_args()
+    
+    if args.all or args.synthetic:
+        create_synthetic_test_image()
+    
+    if args.all or args.download:
+        success = download_dataset()
+        if success and (args.all or args.samples):
+            prepare_sample_images()
+    elif args.samples:
+        prepare_sample_images()
+    
+    if not any([args.download, args.samples, args.synthetic, args.all]):
+        parser.print_help()
\ No newline at end of file
diff --git a/python/demo/medical_imaging/data/sample_mammogram.jpg b/python/demo/medical_imaging/data/sample_mammogram.jpg
new file mode 100644
index 0000000000000000000000000000000000000000..b2fa341e976aafa94eefed71722a6f85e7a36f35
GIT binary patch
literal 21912
zcmV)-K!?Bo*#F=F5K2Z#MgRc;000310RRC1+W<oV0s#R50RaL50RaL60s;aA0|Ej9
z0tEyF0|W*I1_lNN1_lNP2?z!U2?quS2nz@a3JMAe3I+%Z3=0Yj2?`1d|G)qX2msRn
z;sF5>0RO}Q9{>OW1pxs80RaI300000000010s{mE1_uZU3Jd?l0JRVR0s#X90t5pE
z1q1{D00Dgg0s{a95d{(<F+m0sQDFxnaTGFvk)a49Lb1UWQqk}vGh%Xr2?`b$7#SKR
zCMPH<Dl|1VH#j*uLqtVJM@UIZQ&d$}S6EqEV`OD!XJ~0^b98lfcX)YvgM@{Khlq)a
zl9QB`mY0~BnxdnmrKYE-sj9NGw6(UkxVgH*!^FkL$H>Xb($mz{*4NnC+Tr5k<mKk)
z=;`Y5^7Hid_V@Vt`v2Mh2mt{A06zfy_V5k6Y4h4eZDT#p+O5o+wh$wb6SYV^a5MgK
zUh&``7+UIDYUtD1OJrpMzSNNx4%i`p810;A+dj1<nstrU!MCxQoXr6U&clB3jGmd~
z<EN?oJ1+x#Gt%@}ttE+D?b?HtiI<(A5}@<ck~5t5JX7`E6HU{v?k(rkCAEgx)J%{g
z120mneX`ic<5#t{oM=|?X;-jbp}0O{E0N`<a(8Y%n5!BWiZtmU7H0i5DQCSZn@a<3
z;I=Xc(*r)$=^qcR^=nqU)I35Nts5;H#K3uyY&;YF(Z_FE*YP)mCH~8~)!~}rRh9%2
ztViWQLB@H=!NB};Sej?SUkl$^w3BJ!q=yp1h)WBMfJ%TyIM05>@li$KTQ3Vu96ENM
z4asvHkTjcq%%z9`0x(VnbLt1JdDq4N0EW`*lg+EGnr>EPl3T~Y1c%%ddgtGvCcIwb
z!Co5Bb!(WcQ6os#?GhxB;9-@E@TWK#>(2x6u0zFsHSn#ZvMs84fJ_-4-ZoG+?#txo
zXgK=twd0;I@ivjKYD|XJWw<CzZ5Ji5<P0(605LfAHNbeE!IxUbp!S+(qb#w;!xEGL
zuyA)c+rc==#(HC=b6P)xE%f(MHSCtRLQ@1ybcKWD=WYk6E7!5Djcda#VRRiPh2e=$
zn-k+O1B?;Zf;#88_RUS<yDe+OYL|L7g}D1Re5H~xnV6|Oo`aE$`+j<!61nk5iom>|
zSX(QYe|IO@NI-<I%v=(3J9Xs!GhT~jqUej{Tx(Ynwa=DHMz6G%P(T<t#xv0Rk80@S
z@a^TMn-007nN{PJp>~o&6$jnG{0TS*sU2}zn*RWe^;?e)M`Fmcsbp7@>H;0qlBDA#
zo|*bCdsf|ori<Y>J{8cSS?077yWIkOvH=f<>7II;;=DzD<IPK4hgP`O($#e3R-WDn
zlr_O8cgk_{jP*SZ?|Spk9%_0%y{WD4kgSFc{{U&1AP2}DdFV3Uxy5q#&|A#+R(d+L
zFiXtWu*tGdlw@rOoMk}Ap%|_bXq6_LT(GxxOKHH3<v?3$0YdiNaJ_qeM058ks!MSs
zlIeEn%;d5VpDE8u{T%VsIyQ@P`lg2!j1b#`l1U44=V9S<-;!&k_;agU8(1!Gu8f7`
z%VrFzLF`8a<Ye*5uKxhT`c|ppyQ#mktc*TNkF$Ad11gQoc{sx3AC+|e3e#Upfi0ym
zGr=&Ld0tX{k_q`wzdp4;jC3pgYQXAN8cDjEaAKb4Vfho2fHBWc(>20tc0MHV?4k>u
zHfX%#<y=8?jjTvJ2OWCnJ-GC(Z-t%*lz4$Pm;x*@+^D&RM-940ED0cF@y-vmd$yP1
zd#G(SC@pX81&c$4ww3Vlf~V#j_aODG9Uor2*DW8y4XNB2<zJoJ0=s@z$`CSo5tClW
z@UP+SwWsMZwY{7U5e&hVx*eeH+tdS)J&h^w-M5G)i&mDw)=#!G+_Nbd7~=qQ*CX=x
zuQd3r4uaaf)|+)2!7SckOGS)?EOU;&-_E{J_~YRnCLLPuS<>#tX0}s0$fP=rxj#S-
z2kGCB9`QegHOqTJ_Ow{-?k(8I9KeHyQbMpC07)m`JmS2gQNFX&lxtV|Olfy*B#=Ui
z2wkIoG1EPB&wpMk2UYPQeQwCjq^jGZ31#x+0I?*T_4gx;^~XV3^J<n_?ye<G5UI6w
zXc{wuuB7hUyMvHG&U<=SQFn1+;tf7{b)6<#7m`59a*dESIc88t9FQ~kih|d}dV92~
z9)&dag4@juv2s}zj?y!e#yj*U-nD!;;Qb!WAh2y0G0d{dD>{KA$Z!Avbl{RdKWg97
zygi~{=~7!?Guuv!Adk$Ai$x(Vl5%*@Jx{+CzXqW#){ApzqAo10tx(-MA>4c9WOO+n
zmTJAmw=S`3CA9iu3#71i*_@QxhdY4fasdS344<u1@b8W@;oJS!5<+C&E%6~?8mPhf
zPDV4y<DsroUh%#4?ZuV3yE3p<KWM-K#u#j1!=MYF!`N3pd8pWGme&?q%wW7t<-*Yd
z9zn?^za<&^)|QLl-4|Ta-X?+<X^<HsVi8FSdCKy{5&>Sgz^rXY;q3O-M%PKUc$tV>
zozQ%WoBnJbyn1oa)=VD+b$ggn;dJPnNTGwaqQJIBbA!|AUc=#U63wq^b}>t>BF8Mo
zNhG;#+lXRGIp{fK_*biF{{R?t%WEO0MKosPh*1{Y1z-Rzlh|O6Ph3}{d_VZv;tPA2
zqOyocu*O0bKngg@g>LyE=hu;5r!}Sap<^$JuI#0^mPU}WLc?r?w*&$)k;flO((q2L
zuIX1F*ta?r{HSf6p}1tv&6B$r?VOH#RV`meu-E6gj6o~Q6ux1J70D+ZTRpMQ;A>05
z{sYx~4`b)(Fto02RI!j8WCkh;_2lzjoAAf`DjS3SqXd6yy@Pbb$-FO8aDKTs=dae2
z!9D`fH18U{%v#jf07(0mICy3(MtC_rI(Gak>2CmddL0K#xSI0G+Tu$o5#x=UC<|wH
zdJZ}RUVnddW7TyH8%vQIcvQ*eG=+B_SYy+l=U+R1(fa+p{{V@0_>jqKB1h%f?NEHi
z+AwqMPJa{lc%P3nRZRi)D6W=CS{<u8J`^9j$~QpV^X*;*@h;m`FC%!uMTGw4k>Ya7
zz=p<Gr$g7?ys!H<9b-_uiu+{Y8P!DYLaStuGw4^h*13-tSx2ULhD|TXhs;Zz$Q~27
zXO>5nbC5te`r!4f?K;awv$2F5$)~uA<o^Jpq=8fuyPs}OKM!10-Ct3Wb)W1w^*PvH
zMv+;<6op201QzE91-j#$a&cWYou%H}*$5=n<@;O>mWU~1AY>>`I~?(!L&ZU^d_%mm
zxNCT{sTJdhl%0rNJ}?FdIOKmnt#kH1Ki72|T|(c)a74DU+$xDu<Zo~Y0Y_1QLv_zU
zE3WVdie=X?g2i)j96m-?6N4LV`@4?=?ym<t<E3{=qT6fM=foOy*=gJeZe1nXQ3~OM
z4xn?4;~$-I-w?FN@W8efj}(@tYHeZ?SvUlP>xKXfde<Q|Hn&%^&pdJ2OAKtr7X_G{
zlaL1^JC8oS1`Ts7;lB`Ry1ds~<WCC-*<uLfE9KaH-1EnOdi5^{cz)wnv$wdiwhwm~
z63Hh4Natw>zB&A>(EK~9=}_L>-D$Vb+y=II+S)`=t08QbIp?+k=j&W`_1vp_Ft*VI
zk%?Q&ks@8Jx!{g^V!XFb)Frpk(iPSLnHc%XGK{AqArDe;Kp<oAtG*radg+tgHT~O(
zw9B|y6EkjL%Yqn=hXnqEjw{@L34C7DB-G4K^G0PjiZHna@HUM1B#!?8&o%Y`0K#92
z+JB0?HkKDRc9E2gn(8?dBsm}x?AY!-YrnGipP*@y+GzHgr<J}_9FT_I(qLzt{v7eo
zYLeUHHjSxxo$b6PO{Aee(gUtHee9=lPX{02U7e<bdE#vj=GHCn$qY=3_J<4mz`rYh
z^}Aqatyb`6sigcYvzF&vvlGi?OT&HTw1E7qNIYZPkF9s#34B3iuWNF6*H5=g80=Lf
z7Re9`D0Rs>?sK1@uSM|&xvt(>+<#?V&1o6ferQUeEH@4S<YW`c^sW=cTEzMum|sV?
z`&9E>6}v?TG4(vRBN)N!o(+8O@sbz5*X8?Ol-CUyR{N|@2+1Ix2L$nv>0V3Y*Oo61
z-0AS5iRYgs=Y--h<`Ksj9CrG9*M#_^!n&2tx?x+%j_a4q!ctZtoB_vC&vV<_xi21k
zI@Yu+T}nrISIw3}*H+^{XDr1*Jvktboch-RscC}339a=jBYx^##l}-;7$+wqf--s&
z$*OuSw}<VKX1&p6c_o<0lSZV4V#R<wbBrDjJet38f8o6gTeiE^BU_u#v&SMx`*6e1
zA8(hh_iIXP+us)J65VN<iLs4i4fcdc$YMH@-GMv-T*r%VB-0{mYs&|_lqM%TLv9!<
zbI%y-$G==xFL$9^+|32r=@LU>J1Tt6=o2^v@;ZVKPgCi|X&(&hF?e?4`!8CG+7<aN
z%Vk;i<Q(I-%1FjDgI=HT{{Y0c5L$?!R*<!@jEo=N79_Le_ZU5WIOr+A5PV`rpJwcl
z-&<UlXr{FTa9xCKBhUkr{`GPacmr3t)g!sFf*W{d%tX9^yO&UL*PtMBefZ6EkbEzD
zj}*ysd8IwX59Na#jko~2xLl6Klw@*neQR?<zSFdsrn;KO;?>eQ78Q8|7hr`+BOvDk
zkb9o>&}ljjhIG`{?R?`a!tBvW3XDV;;lb_lo(Ek0Y1ZBx@YJ^P%CO%kWCSkxc;lxA
zzHHHaA>tXdo1JnCMuFgGc?<2`wSizUo`jLe^r`#>;V%(uaa_l)*(1+A&~>|zkhmLo
zAaR^yBR#zhdPl>%+wBKllSkJyrxHSrErKA#e0BijVQ_QTp|4o*-;QC?tYSJP!Yr~Y
zB-2RSHZf*mpl1giJ^ELDYvLu;*H65A+3m`^%F{D~2~&_6hebbLYtTF`uH46cYpdB^
zNUc4}xOpYR5;7kcJbc{!1$xcD#;r$4)NPXH`4rvYEKGNW#zyAn87DuNy>-d)OI*|q
zzMplVTfUof12nU^U89CL$7~XM`qx|V_WuA!vV(k=GAcrgxuUfIHd7b_xA{C*NAV})
zkAytf{5=iD#AnC`<1(%z+q(lE_)+h+yl-0g!+m2Uk!U)@q8ow)k~hnaa-fc!FwS|%
zu5aQujwjb_i>|n$;LM?zD04B(XP$PiL5lM)5ons<iL|XsT|ps^+7)GKWBbn`-SZKg
zk&fR=^9@JBntq!UH*(!O<(=da9i#)?05O~s$87V`o8wl}ZG0};289aF%^v2CPrjtT
z4144j<Nez4y<_3@_d489VRsy5r4AySa|qP+Cmz7@?mrt>egwOk*56yQfEzn?kljTa
zozv&gb|-M>kJ7QcP2gV&Nv6*vQo3p|qO(J^F^x_MQH<o{Cy!d=Z*>;YwE1+4Pd+Pa
zRNo{Zjlhk>4lpyvO!PhZtnFvT9v|^DBWi+MQu2A8XgdQrEx|b8{J{P`)RWuID2fw&
zve}{tip4MlE;2F(bJY6-is>$VIpJwMGH&j*R+c7=IzX9F#d5ef>(d_mW16?{;$1&b
zj(9G7N0ez<V2Io8a(-+Ac*24^_0D_qKZkrXVl1z&l2#Mh#_As9{ZvSDrw80|kUxj%
zSpF!`yi=jDnJ1L(Wtucuz{;Eq?!X*lJfBhCy%)is6}26HJ0Ckwxt2+86%yVRP;$%5
z?mVafbJwY^&rI+JjSiiwj}~enWr2Ru*6{TppP5yMV}N^gH3qAt`Mw`T()Ae3mn0S?
z0P+V-qt^^aO2oF)z9(B-$848!ff@w>-l#YO6~*{p!Wyo<WpE(4hWhRXb8!@`7s_m8
zbR7UaD}O}M_01<up4#(5Yl)jG(Z<}Wka3RPcI5pl3QaFk(mYt+B#s!_+hmeObCfJd
z<fs7j3y#&(L#J9xtNp)L(dClvX5AgM9TXgr#DUk2ab0epy3T=PZ0}(k?oyFLSOj25
z+N1;Sdh=WU7V(dWtTf=ack{<B#4<;31a8Gq^4Me7p~v7o>$>nIw~KCXV6v9$X4!^s
z1Lf{yBRvlRK|e#!TJ)*)UlC}!F|eGQ$>eV?NSQ{}Vn`)hpyMQuL0+k$d_dGR&kmS0
z!yICGWQ|$eIR}tQ1oX%pit`T^-hW_P-bbq5M{{#`5i-mMJh$E1fgLl`t#Aj#n!b~9
zHHGA|-IkG~f=g)!?p6L=Kli<IJqM+9S{H_378dbqcBwRz7ebrJE%N8nsmD0ur9<N@
zt2=!z?#?NowzV*oXTc7{so)-`x$Dq=d^d*tJEdIcP+Wh)OAFfvBHzgy(cGMVXTQ_h
zvaKz2oma-QTpuLc+aZb4HD)f#%yMzYf6scY@dnSrnkI&g;p<JV(h>HJ#D4Dwg&p}L
zlkZ+lajD&D+SK<EYO+Z++)S2Hg1a_3Jh3My2apF$)+O${rr+zrJAFP$0Y{2S4a%X6
z5&&fGP5~L`k4$rm;)r}NrRqaNW|~nnc#=6{Q=BQoWqah2n&iG9cvcuJ66Odk<k>%&
z6iDvHR5<~h0mgXij<v{Xwh-vreXYfn+_A*-6)SG15(A&Rk<?&w?OR?ix=kh_5t3_*
zC{aj=Ib*+|03W-zsp;R^sA@Ks+AfJViEbrWM0W2TydQ4M0L72p1fDQd9=ulH!as;w
zzOA4}?+I8u=aNX6e6CQI17YOybNCwcuLEd<#FquNopx!NCOe@jzDpo%l*TYuYV+2*
z?LNsgeQG&$-w(|tq`x7RiZC3W%h-IQfIlkh^^Hfv+SSL|BzF?|ParV945X4uZti#k
zZ(e&>Ipb!6Yx_$m*xSWoe88mqz~GIC2cY98-@Rp}hOHt87q-W6&R@*Q`tjDbei-;(
zNyL_RmL5!sttMu;Azbs@zDXQ<k9zIAOK+%YcaZ7n1*B3NeTjULfw+CvYz%hk$Gvlx
z-T`e|<~!N6<9m~XaIAPaY~+r>0mg88*P-}h;pLx(E#Z5hUFLP%(OZTwC@qqB<EiK0
z@T;--IxSztlIZ&4ODu9JFuZ(>kly~l`g>LFBjGd+crS^VELP#8h0wO*7=Sn(PavOD
zir=yDgxC7qSJoDe8Y$W@B0dlvFaQIR039*k*1Os7ukI}8(w!RJZR9dD%MLvCC`Vo~
z$8W7O#P_xu7l+eLd%LK2d2MM7BY9#$CoDO~Ptv>x<M+g6wAO8|G&x>Gh04toPq@ds
zIV2u3Mkw${jPLJ(5G0pRE2xUz<~38$;AD*T2lMS+2BGm+QjT@e?zGI!6uYyGAz%(z
z5AR`f&U2c=)_x{O;y)0zlm=9qRU+cnCc;F5J-fSZIXyd8eekvLZkuH#z3e|{nGC1R
zK(8E}i~u?l!Q-H;j}v&CP4M*aYW^6rV{nAee;R<T<`r?cwtj49rBc-H?|f6FC!c9*
z%OQ$Kk(+kqBi9@bqaNAsShpSsSuSn=0JE&5x4lTl)@H$H{{UDzaC;A4eQT5Whu|L)
zTIrKZe{T$p46o$Mk&sEi$nAgu9QzT}*Pi&VN;WpOl1n|vmBCqMg@WbXh+;wKqXp}X
zlT<t{1+Jxc4w>R9t|e!ci}_EzQZjH0H-14qax>{$`h-3W)O8C*o@Z!mvL}`}kD0J!
z;fU$M&rX>4rPRC?;eQlp^XeJ^4I{Sb;&dU11=O)8xg>Bgz~-}WJ{N1cM}%BIhNhnC
zNtAi8$GNvQ?W}pmPEYybtNb+hR}Q79NFGbuJKV~!+M^Jp91;P>19klC+B^Z_Ee}bx
zw*JnwnoEcub;Bayt1dCP6VE*7jP9<vJTZR<ho-X@TPBWUJj)#9@r-fL%ir=9fu(p;
zPK!{BQNO-qhA8BlYTfQAN?@zwIOV`N^gZjKwbA@PdnLD>EJk=4$bw?VIa7AxasW91
zp1rDP`*?CApGOOylzpUh{6%W`Q^%SPgJ-5S-iD2}aV&{0*(G8Q+zj_8_58X^lus8e
z#-nAY<WF{Ed=OYE!3WR|8`~N0TUuU~rTAX{O)ag~&E5uSX0~G9Vq+aSRX70S+O#ws
zQrpB$9ks2%f)#y|I4Yo><#v!cJ-O&{?_D07s7qtwCh*<NZE%8VOC*beQG2QARP?|(
zsr*f+&pO9xY(<o+&e6uV1c0FMtVzHl1l7+ATj<w7SS94MY7&rUb(?94#(VMnJ;zM@
z){eEH_>WGC(hF?M=0rZ$?wMeoxFq0=AHtd9TYn9BJHheXUdwlK(#|2BEs1#{Q=UL4
zj19Ttr{P~Sd|>bw_Qs}{8cMM8uH`EvK%Y55$Ojw(F^pBe3~LrTK9BZz*2yl}v%Oyd
zP!`~xxxne~?^1ta>bD+Wh&4S?Qqn7SWpKbg8zebWd;Gol1lKvPYPLEqq?g(x&v%q0
zB(d+pFd>+n6UKAzf%U5J_=iZoj%^0{BnCMfdoUPP8<-qmk-;mT`R$t3y7+ypUtJ?<
zkqLg&9>^LHcRHS+px}%Q<ocTJyfN^yShUMcF5QtOn{k#)gc(te2vO*8GxZqaslUVP
zjSE(aFA?5Osc9LI-7UE)#h<egoB}cFRD64_SX*h5+&-Nu%)5+JHgP0*QMdv*$mi&z
zrFj;;;9VBVM$uW~xVdRTxdvs)3fqPWB;*`$J!_8Dd>nisJ;t}?$8V>|NM&3!#Jd2!
zM+1?boC;;f!`rQUQZ~p!T|7BOw`I$Zn3ZwU7~>f|M_S*}JSPQ)p=)8T*iF4H!--)9
z3kAsl3}e4N@I`kKL7@!@2<#2y<jEda0b@k~0iL7yx_VaU!~XyXcu!K)ZnS+$%HBpx
zW(rXmQQRm$d$A=)Bm1V8#z`)b;xVP^wlHZ?#KU{X9)8il8D9Jq$@cu}j;-+<QPlMr
ztP@I(;@V^eNuw;PTP#X|co@LQ>yM>v+4w_6zJVZ@cbjsCOG{KdFbr5JAo4!FO=fCV
z_dXSo;JnkG;yErERw?A#BiadY2+j{o52bUb#u^L<$o|NP?D=xH+#GsTJ{|DykLJ{i
z@hU@xfGXulBcT}NV;-N*yZ-=(I?k$6TTNy?pXqKShCCLGpzFAS$2@%vZB2LLJx1+s
zA6K@#j^Hu~rHF)yK;w70$Q9W555sr*DYd@4v;rISf9&^@0V60*<}JoS#(nFtZ5zZf
zUPG!iyb@b)%u~hCfC(UV#~D9SSW|f4#UJpN#AUGE9OyGVhZtgU^8?AoJ9Z=1w`Y<u
zs?Ve8_aQu`Qs|L}k%Mw@z;HMfv1{=L+ru`p8=)}3^4;4nuwav%1B~bO&wBGO5_qdu
z)MvfbG~3A_x4ECni_BfSP`4-3oB_|Za$20RS?U*8o*vXAMU7b?bp_0_Hgb9C_lG})
zc@D3j_@hpq<lkMv3|CIyYC(`#<L|ECbDzT^nc?3YYgdZC9lO;hOLmQ#4>C=kBC%eW
z^yjbViu6lA4cf`z`@8KQP_~j+wO}_maV7`N@5#?jzl~>jOT(A?-SSxtGA}MmkawvI
zh{zj92fhgUV!eaF`ag!XoA_*O?w&ZUqGHoPyS&%=l#jp+XWqVx{g{3S_;=z@#LYLt
zn%CKF^w*U`&fqI4$@y2;p+8#pAKUl#2GTwod^2AWcwN)%s{(|yN8Zy&bA$zY4tv+;
z?~Ohl_|wDwADdUzV3sK$4H5|ma|B1`<a!(sJJ*jtfUGrt7!MWrf5S-xv#R~F_YawK
zo!C>r<wqkuy-jEQQM84u8&^#_N1o0=B!H+Y0{p6a=ef^Nb6K7M(ONlfV7rlGcRy)c
zHJlexL1H#%)tA$Z*4(;^okL8s*7Xpv!Di9ii$E4tQg8tt*aNSBt!;QN*Hc;b3)@6e
zwz50N=1M+cysD0Rlat5P*0t^Th`cuem87<RE+uBQfle9v1^R)%oYhTdSJf_JE8-hA
zTbL(zg6#6%SkD;t1aset=e$9qc)w1H-$sx3usn}2xm6te$0WD6=5Pgf{tL9c(wF;g
z?^SjXTcI$o-x=M4I|G&<_t2_dHikR78WoOd1kk*!vdp8&$Sg1c<nlfKQJ%Hp_cxv;
zisnea;TX0%T#W7mztX;umr&MM!tpkxqNkk{CE3~|8C6l&8RQJ`IO$Vp+E0hIJv`q&
zsP?ay4-(_&a7NzzXLmf~o@+Yb-rrsafTS=+RkufRi9;zZ&#rs^KU(w+dsp!V*M~Ii
zKF(-mja_4sQo+_WP{gk8xjbW_;=4Z&YqEHQ!&drck#uzZHwGDZVZ6A?l^aR?GJXF5
z3RB>?uP&q0ber3Uw6WWZ$uRFO*5I%EJoNUY({CY95M4`crJ2kzt}a>1l2CG^9Wrn-
zLGS5ZhlwowA#w2KwI`Ao%(AreTgJopM2?JelgJ#`hj{Z%u)Mi_CqvXDcb9NfRLV(>
z-Eu|%AYlC~$u;dDTKLLKyKgQ;nBi@rm*fw-0O#?^$9l=swe410lV`5m#RQ66z}{E{
zkC-y^*Ez`g9916+_<H-qqDT@uSkaxn&vPV%jij6qr;szy4!EsJ^o#!h3F=EVoQoxm
z$}}svbDg<kl5w}UIO=P%@Q;SBZ7w3ey3}WmJDHic#V`!Zla=IyfTxUl;=L=t-W}4T
z*KN;-B%K{2EpR1OQ6Vk6VPZJIB#wE{y?qPtS6IFHdE!CivF6`BoXAD6Uz#N$n4FH7
z&-uk_e$*ef&%}R;v+1^f+w9h*PA#K+6;4%41Hi`|{#Elw#9t9#>&tbd>yFNsWnJb*
z;K#U(5<1}R2mb)CUU{m^;hh9No21DzJ)ERLFid1)8B>4)fEZ(+*1U7a+MVZzH8y}u
zp&5Mi2-~EUoE%2@7|+duGsXz&YmT{w{A*I`t*Sz;AlM>_f}P!y1Cla2k^DKUo*&ak
zhvxqPgqWG-n)wTPa=B%|EJK{&^dswv=&d!qH^eu$EvrLw2t%^X#H<!VxxoX!@b=?(
zsIIS0{?v<6kQlAT)XK+oA;#s}K|M3y7|(j>yhR1(kEBa=uh~a+aeHi@TiHq&4l|Mg
z;PMVSe=&h^FN|({8>VU2@+6IEZHeyg@Jnt71xVu{;QlxQtoU<U@m8YISYKJL%Se2=
zWESdhRBq1PdJJc`6{#+tp?I%EX4N!SeL^`0lF^cjBx8*22a-YK+M|<0l*q>0NPtI^
z<%4SMcJ}8y@!qqvzZ2^oBJixA+7L)Qxa3dX1BTjLcOdbC2e|-M?+;yT_TF5YwY|AX
z=3Rk?K^qBT+?<T|?~(MaogYwbM$+U<Z5tvt%Ef+dm?L)YL5y?#E2>Q|MY*!INVNkH
zSP{$~YZC4hvnl5a5B~tJwrw>XUORgWODUj^*6mcpg_k=QJ$dR#&pwsw8ZG>mx;fIt
z>?d`o^4{74%z%PG{7)Qv`c-WcN}IzEB-R#iO$Evv#GiSW<|U4Ew2nW?u0!Hyi5d+m
z;%gn&v8<%qgpd`JD%*j_J3+y(G}8Y7w<NuJ?n=CpOcqO?peyH+++=j=T(^xh&kXo>
zE6p0q#E?m-3!RD`i<t0t72VD=lh@bUvhI964zd=<Rf-F}i{-ImcNYMX%k7YUZ$n$Y
zAMl;eh-QM$$HQ}8v@)!Aw>G3XPM8NBHh_OIT{Y*!U0TOZhS*w38z2E%NP2GQ2{{F_
ze!qo!55vEN_c0V*LrRq2TS^Qvd1b~r<O9@v$3IhE$*+7Y@HUsLS;?Sk2{r7t(am%(
zR416~q-{_SKAGblmC*b=@OGE1*iCMWZ}zo`UQk<!Bnz};cNqP7t`FlPXucn@8eXTP
zu$EaCIWAxyakrd~dJu8|?fBP@>%Ix`h2@RMh&(;zPovv0GHwB6JnjI0calf>74x6O
zFNxFYwpY<<dhM#lR6Mq5qkjP9at{Y8IX?B|4dZ_jYC7vCy5DB<r$&k8!z-xEY#o5)
zXY0*&8XtpntGylWZMBayOL6AOEOP=OW-M`nGH}QAt&MxadgZT#KERhoL$yJlXjBQE
z+uI|80ul!RZmtb{Mw$oMG#x};-AdO=1kZrQKp|TPpa2iAA6o0@_?--TgSF<OBD!Kk
zD<lM#OmIgS!RhUvdf5K}gx1ecvRe&TR+8TeUgktxM*Q)~1A&F(>+4-7hkR9Gsm~6J
zX?GEb5Yk+{dFYI>o(~{4GCjR(*t|913+p{1ILvW0GTSMFOLScMji0=5an$pkxcn>Y
zzYbuE#qA}9?f%IaiYZJ4ypl0v)qaMuBG7ECB3ly`x$;y;nG5bX<DPoe%{$?W-)Ky>
zQju)(05(JwS91f?xWGTHSMf)LzS-fX{>s)1iJxl%s)`?H-Wx`E&+EyrKW##N9s}X~
zEB^qrJVApstAf$F$tpqWFn>{7mYyQ<2CHUt>#a%x(5b{DEWi+O3FADkZopST;j3#M
zRV^pEx@#7OFg()1a!Bj6oQ2MK_Z92k4E`c%{vYvGpNH;rnPsz;N=I~KBuY@=md^kj
z9`);2*FGQAv_#eHtf7YD=?sChf>0+Q2I>npPpRjM@b8Nr9q{(CC)yuWmdsl-3qYST
zgNEya(<FNil~2Mx8@+@#hjjf?3%tH=BoU2~Hq3*A+<qL_9r0)2t?$FFYHfSOFheie
z;K(pXg`LSOxON<8o;d@pdDfNU3mAMvlE-3_UCD9fItK~lA>c0K*pd!;$v)Nf5A1a>
zgmf#(*ZvXp@{yFQIbs;F+E}pb&*m%W4~}26YiNHCCh`7=7V#w>bXJQXWRoC}6k{L^
z{dpDTz8v_1pAYouVAeG|i)+_`nO1aPES3i-dx8(k_=@B-;i>q$Si12<UL{j_^E`iN
zk+@JYNzQYU#Mf!?OX8-BqWF?aC@<rl-rg@P>JzR?I6UXj921l4RW%PB=@9EWrO%9f
zM-HTka`HyY7s%OWIpF#?UfHhz_|N-DTix5=XgUhRYv-(%_a0BmsUFs0xDJi}*ROi`
z%i=GHV6@iLPSN5|AV}(p92X%$&R8A+-Hh?yk4oAi>V7M;*{n#h#cdYIk{DJ+Z~+@~
zdVBMl?|cd1FSFlSh_xB*)IlecWxx2jB;y=$^LNF1jpg@;v>%7p);=eO^j&$w?EBB(
zBLR8FHu3oL)aJaa$GR?&tms2his6(+DkM|wEWvi1j-33=d9H6z)2wgB&9vH*%Q{93
z5<*$_LJ0v_jErY*(1BeBjp4g}KS?finXjW8l6za6EMp<M5y{C1@D=I54{W>_r(Qv)
z_<~jh-X!W_QMq64;GAO{f&DqJb<(_0np}&e>H55h8zRn&D*|CayM{mn1K062qvL-U
zPZGm<YipRyWu6}_vbvIfZ)1;N@vc65uQGR%T_s)F$z~`4Cp?^P;=NYiN7Hm?Uk`Yh
znmc1`R@rT#$qom~#kw8|{A-Bur-i&TZvts6dpFu)-10`z074Z5g&j{?@(KJYyIsAy
zmgre8nC%&aXU;4MJBjJQ993w%Re7fPowWToN{V#2n`H6`fhVa4cR&VDPIG}*qtdM;
zlG|C-uh<ur{L5(3m7JE~V4i?u>0Yhj%ZIVFOPi=9jUbJ|RwaWm$iZx3eKYQB)xHL5
zz8?6VmwFz#H`(s3<-`)Sn}7hY+kw{@7_K|xC&L|k!oD+;O*+(mT!_)a(LO|+a4~~{
z$A0;)E5d#})O<guUyVtUZ6?{m%^VV9Oty2m{a5)Ke~7*$+x$Ygi%W*#<$@6YO{BPX
z0a)N}VbiBt<uAM=V_~nQ%$gfY@{$QjRG71pI2h@Hf$3du!@n7L!pg?$!`>g31c^ap
zl?h#hZY_hH{{VRB+P#xT_=)j?Ux!xS9-CA1F4O)zq?+l91Y&m~9WoD3Z^F6j$wBab
zqI`R;+NInm{{UsTWN$Y)CvZIH8%MQx)|KKv6lwatqpE3oMcvK3(yT0Fm4OUf<|~7e
z+ecoVYc6Y#ANXbqtzW}-Q(H$M65c};;Q5iBNx=l<6=T!AY50fZ=ZCy0q(y6aqm_9P
zP5Zgr#kX)kz&vmX#}(#YI`MymykRt#dS#88+h0b;=@$UXcORD{spRw1@~;^2u7RQ0
zS<4Na?Ht#!nInOwRdfV&$84Mr&bjXjOFWixX>iSOK64yEOc|ml<lJ~a#7;5m>t3Jm
z^4{hfiS)g9SluKSnKw%}-#9-m?Bg8d01kH!m2={DrKsxdC5`s1k~A@yV3ysOEJXut
z2LqhoW9oftgVj7usX-iaYFb3o#$aEzTchneL4Ax*rW?1Rtxt!L>6(sp5i4n;R{{u`
z)z(xXbsc!Y<JcPPUg9lVO>GLpPd75CXrz%O9RV3^`e1%1KBQMq;afX<9byDXA(%xH
zKiZ0wQskB(AL2M5^WU27uf8RCWp${pF7){%g53mc!Lz*jXV8#w{&VPB0ACd7Z!A|V
z(8VJB#Nkj706>3;1<prfgY>HMX_l^SzqPIx%R6@ZV~lmLZt&-UE;S7WJV73|ru)gm
z!67V0OAr7flk3Oxu0LMzrH-{7-KDM8`<cm(Sy*iHVcNY3JdEe_uQt^D6RcTYy}jE(
zJTP1|lTN9-%vBlg2VwyHDy;q&*C)0MVWL581Id&Mo89EWJhvSQ&tAVue0~+yHMk?S
zwy;BQZmtk3gpd;0X3ug4c=q<KL*i?11!>UOURcQ{>3p*?D3DC7+?GC=<2dIz=Dp|O
z-@(ZIW8r(R5a}0ETPzE{+Dk||nCFr5p1rg6uQ>SeW2{{3!%vphWW9za3Ott1FmOv}
zJ$Gl<>0U2+W2uX~%R3giR+UxTD55Y*=KvfMc+WnSr{de+4tRdf$5qg7Wfw9xF!^9*
zj?i~|WCMd)7k)8&JvLkUqm06G@`P}#z_9>>&TtM7upYJ4{>{E6@atS%UTQ2Qp4!oZ
zp<ZMeV5E!yKnIhL(x%om4-uaf%+}hK<Y*b(TH3L~DJ1YdqkuaPYPYUh-0JW}cc*KY
zQd<pzXS$1Xs^frgIs@Es#tl*Mhs1vlcs->o@r0V?qmtR9A=nhbVsZZH?b5iP8hlc*
z(I!~*i;YGga_w|u${pPle|HK!e}vbVc*n)|y3<W@cc#f2UPYLtnMp^?fp7@ueM#s~
zy?Keg*R+2YPBgfRp|tX&j^%g8CBbeso_PGp_r+MY)qG9jD~aL1@ZGJ%mXYlj2s0pH
zNLI(N0At_2Dt`xf3&a=N#qH0EZltqvtFVw7GW&TrADQ&W(z-7b>7v5qr-W@^P_&kI
zbQ>6D_!u79$8XY`<12p<vVj!n)~?&3Xn+PWg1d8%nSH;)yx+w0_^?}xWYQ!S0RusI
zkQlbtQaJklp0!F16I`~owZ7A?g4{?PM8U??!t~GJ064Du!oDTFg2hdy(RPzFTwF+Z
z0@*w>9u8O#209w=J{Ne?#2z4fT^PcL&Xl`I5V23a*RVauZhCgFNt?i$Z;3TKSv2^a
zai1|DlWNN0P8XgsH+IL~uD8S5wzuK^Is7%JU0Pe5g(_2a=?LL*gPspPsU-1@)vS$d
zZ6k9TE4?jF8ysYF(!Pf9Mzi7#8sAmZr||5NT}GR1<np%yes>2al0XC>!nym+e$Lax
z5NS6WBh7Vnb-HF%18Nb}^NfN9YXinN_PSNo&6bx6UN)dvZ6YAB=t(4c4Y>Ss+OqFk
zREtbD#^x_RL||rA+pw_Zw{wgx4{Uz5kD!Qxfp2uO-rU`>jun8A+gY~^<174*`po#B
z;yaHIPR|aGsS#z0+2L^V>^)m)`?xsZV!aRGzwL*w{5|lUy6IN4$ql=<2JM^q&KPrp
zoDAl(ejxm2@rRD?f8jB0BHXhysT5EiWpWAmi5zst<~!G&YM1vm9x*z8p&B8&we!5F
zMiMszk+@_Y0pp&P=kVR#c*n(Z>C)=*uz6S!A_8JY7#PPGJ8In8KBuUQuLmQmym0wD
zB`dkg=PIY`gWEaBYMaBFezW0Qn7laBDZ2AnS*{;$!L^l_BZkg62l`j6d_D0cj;pDk
z2&$!`TIsS2$>oeJ?)W9P78n4ZmpysyT9IGFVD?&V=AAP`aAaA^x7|*v8>UCy&(zn9
z{9N#EuO_=A*-aEkuo&WrxKdagZ3~{62iq0SFZQ+GnKy@Tu3)(-8no9jw%j<#`E!i@
z;C~W*E0xsBT|;H1i`fb+Y!MbLl_pGN`T?|&*CW)5zo+~lx6v-HHH|XsK6ROn5XcKJ
z&GN8PIXUgrW35l2cuMvk_?TR1amtsdt0ZwDQiPlq&rEbZPgB>eYg+hoQqdj!GpR*u
zboWq{gn8xFK_?t_`=xV}{OPye2h^=&miETQZAv&*h&D+$KA&_FIQKOm@Rp;gX^Ep)
zT10{82Xyfg7-t#BAa(7G6YMJ8x4~^y?PAgVO>Ym`nWQ$c!xl#_*axNuR{sEi{VT?P
zBYYdX_;EAb+3KrykxV66R!|z;klPi;I2;~+hPmA;PZMbO7h_4c63zClq9%|OMj(;^
z03Vc+c^Ulc(0&zI>w@Vaw!de%g768VA1kQ~$WF(%>&<%Sg{0JV4-zfS-msYqv_>RU
z+W1^!83Tjgy>?oyzN<9Ty`|g|+<&iWpOO^s<K-YP$}x{mTEMyS*NPV6L2aYh`C-qO
zwTT;Z$36J3sC-A`{WjZ1X*@4)D@_rXadC+oc0tL{1CfE;S0N1A)~NGpT2;hk%EvL8
zH*Ug$2RwpDA5mOe*)D^u$9a3G!fp(WA&I3292aa1k_K=y-lx0N;WpAtRxHf&EPEGl
zBnBn04<CoHsOR`=;w!yA)5Ufggi-IeW|Z48$`QMOI}?GQ=Ui{amb$XgUI=vIH0)WU
zQS-E!3&#V304LM0df=_E{8ge_!nUmKZ9My=m+e7JxI4DD9R@M)_}650>-&8f?{t~>
z++2;XYV(57KbCM%5rPJLdsbwN;+<Kgw$vUd%cBLIpDe@6BsO^CkJMHlguE5vy+>G{
z>gsrwF$faRZ70lrVpp7Y0D5&b4~Zhx{55?(u9}-(c^A$0Zzz0`LR>gIx|72$J@cCG
z{tf&;(!5_|sCcKuQsVAxH<OZD-vI3(oa58ilj~lyp!_EApND)r-Xdv>`7t?_(2{|$
z-zHlP)RWWNxqG;-d^%^;JV`XSHx~}gWp3oE6M_KXoDM+e(-p#ATiNQ`)#b(Q<o5Bc
zx9-wWxJJi0j1YLq+xh)zhrk+b#;TL(+Ek4JE9SXk(J|zW*z7<AeFaz7eiq!@U&B9z
zbik4$gE(?RDJXK=Ph1?4>zdQE@D25>l3Z(%hlO${oyjG43=*JpB)3mcD&C{;y(RFY
zYkn1zXW8<T1|yaJ>gSB*dFp!BJ(ihst<GV*(L?!=4<bNznAZif&N0RZ*A=-p!+#Q4
zvT3V*_LwyQ1-uA(RygI0an%0+o-3lzeidG6+F*S|T#J{ERu*iyBLui4ejR?Em7A#O
zEoG-&+|0LAGRTpoyYQ$!UAQ^O13h`_YlQge;W#aPKX<O%=n{Rt;Z<5G;TV+0K1K)m
z4o9_o*P;Ia!fU2@vdc`>;Q>VNB&@IVG^8D@#Ev`u0=-t=!^GYl(r$I6hkMVNqDbLT
z#c))*_W*uJKEkwq1^Ba6ZEc$S%{%b%%Q459AW}fj;nTh`UcKXsyDMAkJGm|H(i?>f
zbrXyNHsk<1=Pl2_=a+e`nC|6{)55AESr|kxoGy6%E7@!{Ju6w&rqfdD7_Me{3drMY
zOh)XG2d5a$b6RknBTVpvT#Ku?=7L7^t(Mtx`5+PxP%+MN`Bxv~8-EJvclQ_XY4(R!
z5xn8cc>&isVU7UCIQ7k6@Q#tCYqx^aP1L?(?u!)2KX#<HT|v(o&*fduhBP}}Qqk;g
z^%<Z<gUGjpCf_wic;l(!xPKn}GEEt-CGi}HWZmUPBa3TiVFRJ*$j{ciN5(dCS>1nN
zrIXzuU7|!n4+P}%o<Tk9pKlpIgY@}s8qmnf$NoPQmQX?Ij{W!?WO`H{9r%}NeQ7+J
ztmY!}nBk4$ZR+40F!n#zyRQss16$HG=Z9GGhS5FJMSx12{H@cTPyYZ~tJ>;%K9eQ>
zhks~uC=ch0a@bySp<A9u9sd9d^e={bOcxPbLE+6&Y*uMyX(JJB-k^d+nH}&CrcHXM
z#D5d`UfV>qv+%}}*YQk}B&#vYD~#p#FVJ9-{8i-sA<(qT?J9Y-8|#sA1<uGK+q90k
zZP~#2vNP$@xos;&wwG6ZB29iqP~n*+T(e`Gmd-K3$MWW~wQY7S2gUGenhV|Be8xEK
z&&sWgsm?g#XydLr*Gu4!ANYKDj$imn9!I&npAm(TO2)+J7~`Vwee1S=5oxi0_8AiH
z%gSbDWg8CFW0F+xSEn44pGu1I8y^+9xzRNPZ*FYdM{OtFnexQ9A6`52T7DMzec_9X
zQzoRVkcpHFhK<4G^Vo$P3h#7u@Xn2Yb9Hfj0$Rr!Eu?L33-2J{=bixP`R23!A$&=;
zmd8z#O<Q=;SR~S{cx)+bFgWKq9C{3MUP<F`9_jkByp|e_(7mZMO67K}jm~%;nIGg+
zH0wP-#(F)(SMHB<Z4#@>-~hN$wYlnWNXNLYAL1v$PYqh$KB1#ss?8*ZS>?bBw2}!K
z&rIWvD<@R&q$>c9`t@4SGRGRB+DQSik&K+@IplHq>G&(bI)$C>t&!7Gcu`fLNgHb@
z$3uWKfLjCH`qx#d{5jUVKv8bBIT|FKVHA>y7#q0&^aJwUS+_nIw1VPDZM;o<(W8}(
zmfT>dJRZK4?3$0m+igAM@h^uOS*-2eU$gvxfxA5cIV1Rev0B#<M=VyWXJQ#5c6-U?
zeB*^GgkvA)y>Om5wbbqBxS1iD7IrTfi3Zm=8T~lMD>Fg7Zx3sc-)b7Q#k}&vwi5sh
ziZa<;_vbx51$sw@d{=+r>y&6z+z}+3<(Tf{92{o=d-_uN?&9yoJ{N6y?Vs&%ypnv6
z@`(8u<0Jd0lfbVg{?3jKYSs;A*v8T!W)CQ3#v3Cir*Ez+h4{L@4xUTh5w4-NH#lXD
z2nIp4xC4$o@z)2^xi1Mt;VVfjE+^66RgjV9$q{m)bI#M#o;f@f{41>R=D!4!-{|Qx
z87=Llk~0#G#O^Jgzc(C=XQn&XK`)H_D|e-()bKsfTbS12jzqE?0;&!V9OoT4HNW7`
z4{4Xvrm1VJLlvduj4kcwm<^A+*C+Dtj`gvjYI1ld;uYL+*!{W{jabVp#F-a?gzN`Q
zXVSe9r|=((e0J=bR1n+UM7yo7$@wxc$8a8}j8~F;anf(REAc9Lyc?k)g4RI5#6x8=
zQ>o{J$od-MVAQoOQd?_HKT7hvw+C!C*j`tYkFM{0kF7gHx6qr!PdxSxuZaXl&Du%H
zP<HhlKc7nLEPOd>rrNO5yx5}q63uaCJC?@$kCbPO^ZD1L4~HHm(|ieQb*1RfZ6y0{
zZQgYOnifLJRe##Z>&6d$@_1+CPK)79ZFKJx_)2S|_GML412JYN=g7#xCqBI@&Y9v(
z6UEwX!pmtCcf#XsyS&&j1M;ZHs8BdL?Ocr7Eyb;}`E$uRvuRc-cP~{XSxG030qNfz
zb6kIplHGK8wCi}xhFKZsU8}i~zlD1C+mr9_Qs_3?Ri3pRUv9UK(raDP{{H~xkN0B%
zyAhBvgWkH0dq>kXTS@LT#+G={1`L4WMFF=1jxoVKdK%)S({!t?L{DWE%UZ14ZEcv0
z#;OO*dVL4hyGZm=rb7CxdfTk=#0(J29!FA1=LejgJ*yAKUM;x1w)-{AD{UqkISeul
zn;7y*&pk=@&$VYo;k_utzHPjY`$p#sK4aYVucx%12Aj)=@a$_Z*}liS%?*JZokw1P
z3}+pUQPsX5O=qY>rrYY10UUd-SbVDGiN`^f>s4?56ho>>b8S3rDzrtCWo)kSyc`qB
z$mz!>yx+w>9Pu58h^+iL)=u_G13k+$c-(SI<I}k{sp3sj!`>LsQ&EmVBRY9U%9Id;
z01Sl#AdZCB2c!H~y3_4$G-k1$7;cf|iB<OrHxdG#dLBJLoh7cTplZ-T4epuc$m_O8
znoaSZIb5Dd2RN=b;_rqu``usdmX~U_i)_HlAsZuAQ^KCuBv+hx4K?^KCA*d>MXStV
zS>{r|C?|1N$EhHmGsSu|vuUwN`r^q2s;8Odn{))mfZ+0TjNo<dbI7i3HFdjvTF%o?
zdy96w&9nL2i#JZZW97leu0?NnN5qfe%SbNnbvQ&ORg&M$OtPJ%s*T4yvvKK)wc%fg
z(rPys7T5aH8^(ytBAJ4TKw><wdBzVU`gX4W0K@+P7p=6ZQ%b*UX{L@}IcHV`I)d3=
zPH+hTdg7({sryR&Fq6dg{uaB?#P=IRERFLJ$C20n0M+@{&XM391IHS0n@GC0k~41{
zlgAvP3}H#&`kc0L-m|rThuRkRGHH+J+0HXPx!^`0YK1-iW1MHPtY^}$b**hxXDug|
zVIfhM$N`VeLiNA{uTE>wE&e0;>%%@Dn?~1fq)U0A5D6}!+!iJIS8hf|6F9-gAXhWu
z-G9Rt+CAj@CY2`JyG4!erXjqxT;Zc0*dL$0N1|)Cmhf21d!&e()<t%jQa4DwNf_kg
z1L@DLV)(DddS%tr$zvf9{k56AnLz~|^4wz?B%h^ur;5BwaSx1blStL1@|8%22VBW1
z`GLc7FfpIQj<v0(>l#Is`b!L}G?xmgG;&4>iFppnx#y?p#%r<AJXbyblDY<!9720(
zrH0$dk7t(|1%gONC$C<oimxuOC7!nn$7LOy5GHOMNJA+tobpa^a4Lqi@mlg!{=&P~
zE$*&ljlRzksUBJ0Knzr8k`Dyu`PV6bsOZ=FY!)qevb<MJ_Y7qv56p9b4{i@5pdOXm
z+5Bm?vb2iKcztmcYJ>MmdJJ_d&3hk)zi3TT?kjscD|8mCJ9*a%>R6C6cqgBfAJV-m
zO|+ZG*A}-{w(A^MQWXx_K1@olob&2<=eX-#ocd3Oth5#TDjS!x8`9vBy#^TMsK<N^
zAJ)9{#r8TkiDN8aVjeXRNbQN0$ZTf^jP%dFdC!dg8eUuI_cpf@#E6lfF3nq=q~NGL
z@w@WRy?Cajpz8Jl>dNlQ*3qLt5rXm)JP=PkKQd3ID?9!Y-9y9LD#p5@yfGs@L&*RD
zcX9w7j~O`NdsS^`#20=m(qqzNv{okA%JW+7Wgw$3Ap5sKDyEg8>qkU`!S+5(mp2YR
z&^)pPU;@m*@tp8G{*`9a!yYQtZ2th^5tc~4&XZ>L;J5|RWFggil1M!d@~U1d@Li{d
zn%>I7c~M&AOBAL;QI$yq?8i9B{6;I4y0o?NWY>wQ>bDk7&Skb}<K86%o&t=4$Rz&&
zy^hstKZd$@h81F-z)LGyD}QLLj2suqz&ZXL0(0wGUNP}Kw}!krb-OF+?k2Y~DlXBM
z+@~sd&!#(j@mO9O@rBpKi_3o$>TyRM%u*O$S)X$h=m+R~de^CVKjM6`_@djxdQPIr
za1@r6w*+<zxmB~*q2zkwwR=B={wZ8*mSWdi5r-;*Nz?!zlb$&Owk2&c!%uj0TYHR$
zK6LSNu_%6X)Un4o{{U<8tgS~@@m8`VxxU)gQQDZDP#2TS1mq5KIw&U{>OU4}I%kY+
zB3l^Fnj@H~v<ktU$pK3G4yW4{k*jF;`aE&zisURMNsynM6O!2=@(Jtt3ZF*!b9JrT
z+Um~~cJofb%E1{bUnDRmKj#_kSiU9rS#@=(>QY^4SZzLh+RUUCF_P*=G0^ANx%$-|
zOTv=*@acL|i8T<5Ew$<Ze5?$x;NTo&bL;I{{v_~+iM$naHN}*(3FqjW$O_~iPBI7b
zK9!$7w|{n!>GE5!xwwpvI2@doCDf7+QbuzDAQj*0Rpanv7kca3i<E*!C8cScWr}XW
z<F_NHO{T0`cs=zic2TH9018pjU>K2z0g8}0;~jrm(?1J*T^aK2E$vz*;x_Vu%Jx5%
zdN0FI2}5(IO{)a+QfVB?Y||8E`A--hQSE|j+kP0us%SbtopAHR4a|?2ww6`r3_%(F
zbJDg?#O)>Y>HIaIK#k?9nOWR}upYjnuYX@^#Iy0r%N)hwF!_-(wXD(nxXB^6ckl06
z`qXz`Bh$?Cu!{Cr-Z*@@7y_dV2^@fV1Jv<e1@T8oys_~|m3^qZws!Cx<`Un$a9H3T
zG0Opt_2&218eQh7Z0+Wq-D2`%5vf%S2>>31_v^<4(xcHlA>us`O_4OlV<q%tLRaQs
zp=BNT0Fi;|fr{+D9$rnQc#g|g@g>6Bc`%TWn0|SQQ}Z(p=L~937QMZetY1faQ*kkA
zLjVfJ7~XKW&QC%`aK8~O{8E;edeznRX$9rXX|^cNW9KRZ1C=KkJ-c+SSHgNer=@t3
zKlM`E!i^(aggA}W`WDZoPCH_+c%NJGABZod{>RjsD1)-Z#yz+LAe<h(0MF97jcdYs
z?bJ7L>oN(#w#bRN;aQ(6xF5rT>T54q(Hp~;$eIPDBIPA#HsG)&xycFu9=II$tQ`kf
z@y)iZw{2x4?{9Jq#$zCb`|4Lc`mfX0y*uHb#4i-+c6SMK(8NB{;H-#7;MpgyPQ87-
z>)QMW;{9V)Um^72=1PX)7BRRmX6ec2pMER1()?|Cr)$SWlTSbciKU!Nz+!Wbc-l@g
z^cC0LXc6k#s`!Ln%JV@IHPS}B2F7x%3F(hs=9#Yek4cMI`z5}TvMjctw!C1HOn8$8
z`sC-4-}9;M{A=O6-8xMxPqk2DoPC}alVDz@h&=WjbKl>YFT5kC+IVkJ)vb)t+^b~^
zxOn9SaunyUB<Gy+Ycu{8UdwpyEbX%I0<?l=^7k>oZLQZQc0Vy(FT@=N#^UO8dL{EE
z4ICar<)I*gRQAd0Seo)b8*0*cV(m0NB3G5y8)w1F$-MA-j>99`mqqdYlHE&ftZ9+n
z&m<vIW_<6E<nl)g)RH<JSAFpR07vnEi&}EFFxjF;%?y_0V29u~2i@Rx{{RZ=qwp7s
zVL=z#E)LO*p~&OYo|Je}=ECbp7kXX%$#roEl*U&Im|_7Rr{*(Wo#AbF#d@}oi+g@y
z5XyGJ<b%Bre@<|J15oNZ=96!%E}wa89gWrfo3OZ*f+T?FYVGfXjP=eytbY>CquMQm
z8m6(bNIqLgN|q9{2gvP_j^3HBYUkn&_KSN4r*2vlf@9_@V8%lz`=ckW20pc!@j`7$
z{Jluw?eAe}Hh{S;$@vKD#c>`Fwo7aGmHasc#A;m4=7YU9<S7}(c|P4gF1Nyd7PQyz
zBzvnjis2(IX&WKNasUbk@Zh)S?^}B3!};|c9P3vae0Fk$cnOFAyFxZwq31uPc{abG
zS=(PpJH;iutC+3sA>Ffp#J9f~&sxI0llwzkywG&%FC{i^8tU_B?q4cSRCEWN{uR<`
zu<7<5AGp(W>-cUaxl(3VLAjZ-SPnCuhXc~N+b;zV5lR<PDytl(MHb9ehyjiP>^|Tl
zJ&kMXUJtyNP>W1ytkNYdwd9D5_6p2<j(-qplz$B#*I2vOXVL!CBwN19V@>MH6m!AG
zLFe!vg|YC$Eq2(Yt%bxBJS^r!VH$-QZg%|%_at#q#qfVuHd4uRVQ;ky%<lZ{Bjy<-
zAK_ESzrA%n3)XIQDQ|Qr?Y7%m$nsk4!#Ox5S8y3o*CVGjzhfVWHHl=u(`Slhc-5}M
zVM^{|!2n$DUT{u2kZaT~^u2P^L$tWIw}#3qi*1Ui;N+4?!0Fg^KbLs)dtEZZf9&0A
zC|#}L2WEbGU9*BS)MuVMij!K<G)-J1L#ZXSh18-fM7SMJPCi};9Qu3LPNz?>(qmb#
z4EB-Slz|!yvmL*9`w)6(`By9P3sr+v(k<@nk&6t*uWu;=5DOMo1HN!`+nVHjVP!g8
zz2&Tt+<l@=t1CG9NK!pI5#Jp1UU%akglnT~>X&M}omH;ojv<2CBml#q&Ob9ooBKOi
zUu%P0(HV?(*pgW8RPAPLk<&Tnl6|Y!{t|db`@<M{rj41{sc4B2lu^p$^#`8axW#()
zwu5CgymG#!V;N)w#+bsAIu1Hl&AL8?KZmqgZnZ0^ArVLAqemd;XvQ34lh(Wc00elx
z`e(V7Ef*2nf1^scY@M8p@thKRXB6Lv8upnswP$gx#Lkdf;ZKyxsCJOxXO6l2zb<1^
z@jiv(D@L<bGFq8@?LJl;+7?WL8;<$Ka~~6QoE<Yz(wSwpk<bYD?<o7Z4fqUVy4elh
ztEDBjt$hf4D<lvjL-zn~Tc=!-IOO82LE%3N+FXei`Ux(fbq^BCKw?HoU_S8YI6Mvo
zW@wEri{neniwnh$G>zhSi-EkWHV}Jt{5Y=0T{gkICl#)>sKU20rs!mk=HRCUbJwXS
z+Pufd{u{EA_QqQ|Mq+`TVs{E<1+YtfdHR~*?qu;rzM*-f!)Vr$HQQzh3Cpi5xL^^F
zq3fS|==2{6UTOBxTWWSl(Y#`3675Ia!G~VpkLB%J(|B^~%UFPUVjUnTFD#>RX>q(K
zxIM=ee@l=yopn9_m@VFB^5T*e1<^N?>5g*WkJ7tu1NcK%lS-O9#X}|B%FE`k+N2OU
z7(DjJsjVw-1g4j+HMXZB2`0A_$1HJ((2SB-AI123RozeFriHBN5Z<n%=gzk!UIkts
zlGw*g1pws#0Q$JyXW&{#4WIU1;|H4TyowwK+QIgN&!?}a=UsdYu6TdMQfXE;MlOrz
zGv#*-oB_vd4tnCdjUHbUd2m`ao10@J8<^F2OksK-a@_kG=zLFmYjt@dU&CiCF{S{G
z25=An2^`?(k=p~MW7uj7;u0j1bV(yBBUtgBxgZ9|9kZIz8aAS~R<Eh+5vc{)L}<h(
zI01p`Fc0TTZ3U&^k4UzSytW0dBe(=Jl|0}9j2@p#>3l)pyKN>Igf^DQZ&az8IQF!o
z89Q(fely1yu47yHHXUP{O*YwKHz3aC%aWMqo}szNPo;CV-WZ4^x?PNo98rl5F}aHp
zaHDa@UNQCKJXX{;b6eO$6_%ePeD4j!X#hZ^jN}|+uY7l{cx?V72$V;Hp+?g34CCLP
zwei)@jx`(j;MO!<T3d^GHkONV#IqbTGXu~A#xg6@z7P1ZZahb6qO@V1H^!nD;U*vK
zkA89VuHWKsj1HcD%jEg&E%)EYG&l;o4>%ogb6$Djty14qZ|u!J{=mU;xo7(<qcNva
z2=_mhD_7xcc3SU>E_8N`$!BjVnPp?NGK>w_Jx5-Z+iQLUm-}V)n4^$DFeKWCYJrtF
z8Rr~#sPxZ=IyZ{-7<Ek@0Bn!#6q_nn9B)zFcRAxcR9-LmeQn`b*5|_7TC&FZEawdD
zL+&S@8yFm)dbgrWVSA%2*E(}6`BEU-o6Lj<;A5QONzGF67ld_<a%(H*(s)~F_N3AP
zV{b6-PBFnb&wr(6c<<q5{f*;aTiMM$%AnqoAQ@fV54XO1Rv*M44{LfZo2JXIWJDx)
zc;P}BkmRV~=jG}8S66FsrRq8=>CLCfa#WM%Ns;9e2wl17aNIxdp4CHFA7j;RbnRi$
zZBN{HlN;q2DxiP>IVAoy)A(;%ztSxf+D{BqNjn!w3lfo*Aok?oW~}(WZw=q+@#{Li
zuoNH>h#B%w7UbY_{(b77hrC^&Hmm)iX|LTmUC9J751K)^e9R975Oa#&i&pT?i+V3b
z`^fiEiDZsd%0M&nfx-O_cr~>R-k&|hQ)ybF&viP?k_J^DEU=UyZd1`)Dn`}^7(7z=
z-pbFzaYJEga}l`Nx$bTs0l?k4M|_e$9<_y{>9>>WK6Q-jVscm`g<{S_jjBhdq3i`&
zn@ZNLBaS^U!y*}^o8(}?7aO}{uW{|Vt!thj)O3hplTqH0;o3sOgdus)bI(4N&;J0z
z-JLFPHD<6|pDGoRn~4U@C?MmJ#twa}BVYK};unc@n_WKcdwaj1O)S6avQ%StUJFya
zKjL2$TI!QFsP__CE1x}i-mKp*-OpwO=hC|!KSR-V>o{gJMI2Jh@-Q2LS6p&WrsG|9
zkv5BO3~hU)%{81Zk?FHLGVPEq3CQcuudYo?FN|Zik)Rh4VHl|Sm~`Wy99QNCh5TC%
zjTARCY6vZ4<yGYm=5Rp-gV!e*K9$(~IQW5n(qCR%!6b7@aP4a(Xg4q<9OEDilfm{C
z)ZF;@#hT^fXxIA4Np(AEB9dE@Nh(3wFunem>+M$m0J5OaE%dv68%=`RH!T^6FYhL7
zu?GjSJmmGjuG8U9gM3Wd+iE(E%d9GaWE;s{gz=Dh1K3x$cyCtH<BdE$b*INQyizhc
znR3WE>CY;Cc{#1Ud&BlxWzY7`mhtJ)1Xvlv4XUJ_ybhrHds5x_3N2H^V_3brwvr&l
zouX_qWdX+0dF1=?Ts4M;r@pdv4-KOU*r*7SY%!h8zbMZf;=KdI9vAUs*YGUbJZ%bx
zXO+NV8C;GBKyb(MseE|?TUtqJa`LU*(lWvFk`x6v`G-#6{{X#SS>jzcQ@yi|HZWwq
zz>{wX0+F<C&JU)0)f>$lQL>cX>sOEq2r=at*Y~qA1bo=w0qg$&)}&7i=#66zyRAj%
z%jOkZd)bIW3BYhVoK!w5u<>QoQtG;1on~XhE5!iEmH_F`as553$2Gqf_^R;8k!i8{
zntjUK#$qfMNXn86Wd8t<O6q<Zc%EG;X0jq#f}j#PfU_i^j0Pt-1RVVZdMAKdKM-nf
zV|cLnh}+(GAP9iz!R|p~D`Qjmg=eDplv_zCMYM$OmL|xT<-i1b@t>j3(>F!<k!!AN
zHxpgkO>%6aRgM>j<xG+MM~(poy{F-mcsgx7c!N`fd2JiLoN+<r{M-^)DfGzCPAfw1
zTk!q*N2o^jQOOj+O|s>2@{yJYsmTL@+PwGTGHJH9g{4z6v}B~QGf2@4HsILE=NKPT
z>0GXvs%l!6fqMD}*ugftce0SvFV0(@hZ{fu;<9uf8<oD7OT6&C^k&I6#d491Gm(LT
z&<v^V(y88fpIFp&;<~=;Tc?S)M2bk0JrH`}XY{L{5AoKoWp!~4w1#MH)QL>Wb~yl?
zdJkXon(wsF9Qd17(Zsh>LaT0L5k~G4&SRD=IO&g<@uZhon&L*dy77FG%Oe(&HGUC~
zBn~<QUzm0t1)IY5msU~g6H67+Y*<1Tk~Sfnao4X09@wtW!P<L466y2mj|JzNe7Gfi
z?F>gb2Ll+u2RJp`cuT|9RyQeiscO(%T*Gqgj%3f6F~Iy!TyzH&*m%B69}&xEcOA@1
zEp2?ONMjN!V3J7587<rLuVnB}muvljXLRvT3Po-sQ90aP{{VRUfI#d0D|1cMbw3O|
zR=y&!lIaeiD70|~!m7+lf<V9kdiKXP-C0WxwIq6Or*6;*M)G8QMvkq3IO=hej-Q=7
z#-0w;ZcA%AeyS}EqP&9N?RZ5{2qZ5c_0K)47sGm7UMst_)bv}4t|wVMse&LIF}UX#
z=Q-`|UZ<dHQ9<FUCZ6dn<Vjm>Vlv_~3w-|o?4E=1HOSm*T9=6YT_&9*WspYYT#^<?
z<c2$j0Re~k6~=r^(|kc?6Y2UrjIA_)d9xVzDGQRMkUqS186SpDx#2P5bdOrrB4aEO
zi2R7;03I+%9nS3Ij+F(^gYTfQmKTH1Hh61)EeUL#fV}aLOAZI{)$LnH@g9kzthRbZ
z+ug_l&g+Fx+gK{M9-{@ijo$U*I);I$c!t|i(!7>{O)4}pOs9XCxKq#$2sz`>)gKmW
zQG6)yloI{1uBVk<5)~tKviEk$&pZM3HAlkV6-!w^i1m5q)9w;6FzK~8`A1yf=kVsO
z%j3tBPlo1rrj~1&8couF5)Zmb>N@@d@vlJe-^9H$!Zr!3M`mT2iIL((;4T%E9R1)j
zIsFZHJ`eq&b;#j@(D-)FHMxWR8FR3|-XnJ;bRW~LXLzULMchsk><e!&ourVc$(4sa
zM<<-|*ZEZXmyT?#v}cJt%gE<aCxainUO2};-o}2Wx$SSp*M1hVZ6XVE7`u^@UN97f
z!k(OUUUS%-aBB7TtF3A<wceR-TI~@Wu$^%}Fam+b)5d<STz#L4JXd?-xU~&Z<~gI5
zK$1ZU?<fmDLNGI&ay>xpT?P5O)$LgLYet^gTo{a4Dp(QAFx!H7;QRX5ckuH<oNDiN
zsM%c3V2MO;9Jo~l*x-VAz$c#g=DMSCVG%J(VGO)5YjQ|F=e>S#_~t0%*6k&Y0<&DI
zmR=4wCx68Lb+d7&*=gFh`W3Rr`@*;bw$O3;*Rbf?CXz2;L8rt5MG~M3_d&@ZdSr_9
z4NpSRq0}_EqSIkeWKeLNo)3PscuMl?%$6pYvPG=<QLFsGBIgT(+PzBG8it0@-T@+!
z86c1Eb*rBYNBSd5%vrag%dt2nOq^$@L0u-P8ec8iULcf5X}JT8;O_(Z6UA5XR)eJ3
z>(N7LYvn7vOv@aMpp%^RKGo@(b%mT>E%|NKi6Z&c{$P5ZhmNAGc#gu(+Wn=swn*4Y
zm@%$8XVd9iAB_OF)5}5vo$_=YdY*kxy?H3UySiz0ON)kS#nC|7zzRUf$4}C&ctYZ7
zEG)do+-KBCl(5NQH{fK~b>c~-@@IfZ!4fi6LB<J)<o$WC9`Qz-V)jCFY@yT0ky&^l
z$>j0e*Cp|L!JZdB432lv<7Zo-iq6bPUoplxK7dz`_|r|XvezCvYh_gy&f!SN-Gk3H
z%3Nv}*BZ67b6v>8<lwpfD6)T0E9g%J=(<txo{cfEk8zqjoxbL}t2K&G2Ty+X@fnh2
zjU+KHRY@EUIQz%*r19>ZV|jIM_cn;cwZA2Oa;NhERwL;ed|o3GrM?>~+1|&Hzz+38
zTS2aPyHlT5W}VeZR%SmgCLm)SI2D88_m@ewj!OyTjf66=K4CoPIpmtuDz=_2u#VJu
zOC^+YNF(`oKlhG1@_#z_?}b`Uoj->YR=BiDE@QTLNo4tmVID!}+@H?6Z-sUixAyOG
jbrE;Cw04;8BY6wKI5pGV_(Q`2>QFRGfAx+09`*m(XDnZ3

literal 0
HcmV?d00001

diff --git a/python/demo/medical_imaging/demo/demo_anisotropic_diffusion_medical.py b/python/demo/medical_imaging/demo/demo_anisotropic_diffusion_medical.py
new file mode 100644
index 00000000000..ba7fb0a68e7
--- /dev/null
+++ b/python/demo/medical_imaging/demo/demo_anisotropic_diffusion_medical.py
@@ -0,0 +1,205 @@
+import numpy as np
+from mpi4py import MPI
+from dolfinx import mesh, fem, io
+from dolfinx.fem.petsc import assemble_matrix, assemble_vector, LinearProblem
+import ufl
+from ufl import grad, div, dx, dot, sqrt, inner
+from petsc4py import PETSc
+import sys
+import os
+from basix.ufl import element
+
+sys.path.insert(0, os.path.join(os.path.dirname(__file__), '..'))
+from utils.image_processing import (
+    load_image, save_image, create_image_mesh, 
+    image_to_function, function_to_image,
+    compute_psnr, compute_ssim, edge_preservation_index
+)
+from data.download_cbis_ddsm import create_synthetic_test_image
+
+class AnisotropicDiffusion:
+    def __init__(self, image_shape, kappa=30.0, dt=0.1, comm=MPI.COMM_WORLD):
+        self.image_shape = image_shape
+        self.kappa = kappa
+        self.dt = dt
+        self.comm = comm
+
+        self.domain = create_image_mesh(image_shape, comm)
+        
+        P1 = element("Lagrange", self.domain.basix_cell(), 1)
+        self.V = fem.functionspace(self.domain, P1)
+        
+        self.u = ufl.TrialFunction(self.V)
+        self.v = ufl.TestFunction(self.V)
+        
+        self.u_n = fem.Function(self.V)
+        
+    def diffusion_coefficient(self, grad_u):
+        # Compute gradient magnitude with small epsilon to avoid division by zero
+        grad_magnitude = sqrt(dot(grad_u, grad_u) + 1e-10)
+        
+        # Perona-Malik diffusion coefficient
+        c = 1.0 / (1.0 + (grad_magnitude / self.kappa)**2)
+        
+        return c
+    
+    def setup_variational_problem(self):
+        c = self.diffusion_coefficient(grad(self.u_n))
+        
+        F = (self.u - self.u_n) * self.v * dx + \
+            self.dt * c * dot(grad(self.u), grad(self.v)) * dx
+        
+        a = ufl.lhs(F)
+        L = ufl.rhs(F)
+        
+        return a, L
+    
+    def solve_timestep(self, a, L):
+        A = assemble_matrix(fem.form(a))
+        A.assemble()
+        b = assemble_vector(fem.form(L))
+        b.ghostUpdate(addv=PETSc.InsertMode.ADD, mode=PETSc.ScatterMode.REVERSE)
+        
+        solver = PETSc.KSP().create(self.comm)
+        solver.setOperators(A)
+        solver.setType(PETSc.KSP.Type.CG)
+        
+        pc = solver.getPC()
+        pc.setType(PETSc.PC.Type.JACOBI)
+        
+        u_new = fem.Function(self.V)
+        solver.solve(b, u_new.x.petsc_vec)
+        u_new.x.scatter_forward()
+        
+        return u_new
+    
+    def denoise(self, image, n_iterations=20, verbose=True):
+        self.u_n = image_to_function(image, self.V)
+        
+        a, L = self.setup_variational_problem()
+        
+        for iteration in range(n_iterations):
+            if verbose and self.comm.rank == 0:
+                print(f"Iteration {iteration + 1}/{n_iterations}")
+            
+            u_new = self.solve_timestep(a, L)
+            
+            self.u_n.x.array[:] = u_new.x.array[:]
+            
+            a, L = self.setup_variational_problem()
+        
+        denoised = function_to_image(self.u_n, self.image_shape)
+        
+        return denoised
+
+def demo_synthetic_image():
+    print(f"Demo: Anisotropic Diffusion on Synthetic Image\n")
+    
+    # Load or create synthetic test image
+    test_image_path = "data/synthetic_test.png"
+    
+    if not os.path.exists(test_image_path):
+        print("Creating synthetic test image...")
+        os.makedirs("data", exist_ok=True)
+        create_synthetic_test_image(test_image_path)
+    
+    # Load image
+    print(f"Loading image: {test_image_path}")
+    image = load_image(test_image_path, normalize=True)
+    print(f"Image shape: {image.shape}")
+    
+    solver = AnisotropicDiffusion(
+        image_shape=image.shape,
+        kappa=20.0,
+        dt=0.05
+    )
+    
+    print("\nApplying anisotropic diffusion...")
+    denoised = solver.denoise(image, n_iterations=40, verbose=True)
+    
+    # Compute quality metrics
+    print(f"Quality Metrics:\n")
+    psnr = compute_psnr(image, denoised)
+    ssim = compute_ssim(image, denoised)
+    epi = edge_preservation_index(image, denoised)
+    
+    print(f"PSNR: {psnr:.2f} dB")
+    print(f"SSIM: {ssim:.4f}")
+    print(f"Edge Preservation Index: {epi:.4f}")
+    
+    # Save results
+    os.makedirs("output", exist_ok=True)
+    save_image(denoised, "output/synthetic_denoised.png")
+    print(f"\nDenoised image saved to: output/synthetic_denoised.png")
+    
+    return image, denoised
+
+
+def demo_medical_image(image_path):
+    print(f"Demo: Anisotropic Diffusion on Medical Image\n")
+    
+    # Load image
+    print(f"Loading image: {image_path}")
+    
+    target_size = (512, 512)
+    image = load_image(image_path, normalize=True, target_size=target_size)
+    print(f"Image shape: {image.shape}")
+    
+    print("\nInitializing anisotropic diffusion solver...")
+    solver = AnisotropicDiffusion(
+        image_shape=image.shape,
+        kappa=20.0,
+        dt=0.05
+    )
+    
+    print("\nApplying anisotropic diffusion...")
+    denoised = solver.denoise(image, n_iterations=30, verbose=True)
+    
+    print(f"Quality Metrics:\n")
+    psnr = compute_psnr(image, denoised)
+    ssim = compute_ssim(image, denoised)
+    epi = edge_preservation_index(image, denoised)
+    
+    print(f"PSNR: {psnr:.2f} dB")
+    print(f"SSIM: {ssim:.4f}")
+    print(f"Edge Preservation Index: {epi:.4f}")
+    
+    # Save results
+    os.makedirs("output", exist_ok=True)
+    output_path = "output/medical_denoised.png"
+    save_image(denoised, output_path)
+    print(f"\nDenoised image saved to: {output_path}")
+    
+    return image, denoised
+
+
+if __name__ == "__main__":
+    import argparse
+    
+    parser = argparse.ArgumentParser(
+        description="Anisotropic diffusion for medical image denoising"
+    )
+    parser.add_argument(
+        "--image",
+        type=str,
+        help="Path to input image (if not provided, uses synthetic test image)"
+    )
+    parser.add_argument(
+        "--kappa",
+        type=float,
+        default=30.0,
+        help="Gradient threshold (default: 30.0)"
+    )
+    parser.add_argument(
+        "--iterations",
+        type=int,
+        default=20,
+        help="Number of iterations (default: 20)"
+    )
+    
+    args = parser.parse_args()
+    
+    if args.image:
+        demo_medical_image(args.image)
+    else:
+        demo_synthetic_image()
\ No newline at end of file
diff --git a/python/demo/medical_imaging/requirements.txt b/python/demo/medical_imaging/requirements.txt
new file mode 100644
index 00000000000..677b01f645c
--- /dev/null
+++ b/python/demo/medical_imaging/requirements.txt
@@ -0,0 +1,24 @@
+# Core Dependencies 
+# (Based on Last Stable Release: Dolfinx v0.9.0)
+fenics-dolfinx==0.9.0
+numpy>=1.24.0
+matplotlib>=3.7.0
+scipy>=1.10.0
+
+# Image Processing
+pillow>=10.0.0
+scikit-image>=0.21.0
+
+# Visualization
+pyvista>=0.43.0
+
+# Data handling
+pandas>=2.0.0
+kaggle>=1.6.0 # Unique to the CBIS DDSM Dataset
+
+# Testing
+pytest>=7.4.0
+pytest-mpi>=0.6
+
+# For DICOM
+pydicom>=2.4.0
\ No newline at end of file
diff --git a/python/demo/medical_imaging/utils/__init__.py b/python/demo/medical_imaging/utils/__init__.py
new file mode 100644
index 00000000000..e69de29bb2d
diff --git a/python/demo/medical_imaging/utils/image_processing.py b/python/demo/medical_imaging/utils/image_processing.py
new file mode 100644
index 00000000000..8bb6245d0d4
--- /dev/null
+++ b/python/demo/medical_imaging/utils/image_processing.py
@@ -0,0 +1,141 @@
+import numpy as np
+from PIL import Image
+from dolfinx import mesh, fem
+from mpi4py import MPI
+import ufl
+from scipy import ndimage
+from skimage.metrics import structural_similarity
+from scipy.interpolate import griddata
+from basix.ufl import element
+
+
+def load_image(filepath, normalize=True, target_size=None):
+    img = Image.open(filepath).convert('L')  # Convert to grayscale
+    
+    if target_size is not None:
+        img = img.resize((target_size[1], target_size[0]), Image.Resampling.LANCZOS)
+    
+    image = np.array(img, dtype=np.float64)
+    
+    if normalize:
+        image = image / 255.0
+    
+    return image
+
+
+def save_image(image, filepath, denormalize=True):
+    if denormalize:
+        image = np.clip(image * 255.0, 0, 255)
+    
+    img = Image.fromarray(image.astype(np.uint8))
+    img.save(filepath)
+
+
+def create_image_mesh(image_shape, comm=MPI.COMM_WORLD):
+    height, width = image_shape
+    
+    # Create unit square mesh and scale to image dimensions
+    domain = mesh.create_rectangle(
+        comm,
+        [[0.0, 0.0], [float(width), float(height)]],
+        [width-1, height-1],
+        cell_type=mesh.CellType.triangle
+    )
+    
+    return domain
+
+
+def image_to_function(image, V):
+    u = fem.Function(V)
+    
+    # Get mesh coordinates
+    height, width = image.shape
+    
+    # Interpolate image values onto function space
+    def image_values(x):
+        values = np.zeros(x.shape[1])
+        
+        for i in range(x.shape[1]):
+            x_coord = x[0, i]
+            y_coord = x[1, i]
+            
+            # Convert to pixel coordinates (with bounds checking)
+            px = int(np.clip(np.round(x_coord), 0, width - 1))
+            py = int(np.clip(np.round(height - 1 - y_coord), 0, height - 1))
+            
+            values[i] = image[py, px]
+        
+        return values
+    
+    u.interpolate(image_values)
+    
+    return u
+
+
+def function_to_image(u, image_shape): 
+    height, width = image_shape
+    
+    # Get mesh coordinates and function values
+    mesh = u.function_space.mesh
+    coords = mesh.geometry.x[:, :2]  # Get x, y coordinates (drop z)
+    values = u.x.array
+    
+    x_pixels = np.linspace(0, width - 1, width)
+    y_pixels = np.linspace(0, height - 1, height)
+    X, Y = np.meshgrid(x_pixels, y_pixels)
+    
+    Y_flipped = height - 1 - Y
+    
+    # Interpolate function values onto regular pixel grid
+    # Use 'linear' interpolation for smooth results
+    image = griddata(
+        coords,
+        values,
+        (X, Y_flipped),
+        method='linear',
+        fill_value=0.0
+    )
+    
+    image = np.nan_to_num(image, nan=0.0)
+    
+    return image
+
+
+def compute_psnr(original, denoised, max_value=1.0):
+    mse = np.mean((original - denoised) ** 2)
+    if mse == 0:
+        return float('inf')
+    
+    psnr = 20 * np.log10(max_value / np.sqrt(mse))
+    return psnr
+
+
+def compute_ssim(original, denoised):   
+    return structural_similarity(original, denoised, data_range=denoised.max() - denoised.min())
+
+
+def edge_preservation_index(original, denoised):
+    sobel_orig = ndimage.sobel(original)
+    sobel_denoised = ndimage.sobel(denoised)
+    
+    # Normalize
+    sobel_orig = sobel_orig / (np.max(sobel_orig) + 1e-10)
+    sobel_denoised = sobel_denoised / (np.max(sobel_denoised) + 1e-10)
+    
+    # Compute correlation
+    numerator = np.sum(sobel_orig * sobel_denoised)
+    denominator = np.sqrt(np.sum(sobel_orig**2) * np.sum(sobel_denoised**2))
+    
+    epi = numerator / (denominator + 1e-10)
+    
+    return epi
+
+
+def add_gaussian_noise(image, var=0.01, seed=None):
+    if seed is not None:
+        np.random.seed(seed)
+    
+    noise = np.random.normal(0, np.sqrt(var), image.shape)
+    noisy = image + noise
+    
+    return np.clip(noisy, 0, 1)
\ No newline at end of file