[Insight-users] mutual information again

Tue, 7 Jan 2003 17:49:08 +0100 (CET)

  This message is in MIME format.  The first part should be readable text,
  while the remaining parts are likely unreadable without MIME-aware tools.
  Send mail to mime@docserver.cac.washington.edu for more info.

--202006272-1637312257-1041958074=:7438
Content-Type: TEXT/PLAIN; CHARSET=US-ASCII
Content-ID: <Pine.LNX.4.44.0301071748031.7438@g1-095.mcnet.azn.nl>

Dear itk-users,

First of all, thanks for your useful answers to my questions so far.
Not all problems are solved, yet, so here a question on the use of 
mutual information using itk, again!

To get "good" results, I am still working on a solution for a particular
problem. Whenever the two image volumes that I want to register with
mutual information, do not have "enough" overlap, an exception is thrown
when I start the registration (Exception message:
> itk::ERROR: MutualInformationImageToImageMetric(003C6D20): All the
> sampled point mapped to outside of the moving image ).
Somehow I understand this problem, but what to do when I have to images
that are already aligned up to some level, but where the volumes are
not overlapping sufficiently (eg when I have two cubes with the 
center-of-mass points on top of each other but where one is rotated 
wrt to other, along some axis, by eg TENS of degrees; A FEW degrees does
not give a problem!). A more practical example is shown in the attached
jpg-file where you see a MRI (grey) and CT (colour) slice of medical image
volumes that I want to register. The MRI is rotated (with about 14
degrees) wrt the CT, so that I do not have a good overlap between the two
images. Moreover, the CT and MRI are already aligned reasonably. Note: The
reason for this rotation is that the MRI has a rotation in the
index2physicaltransform and the CT has not!

Personally, I see a solution by making a mask where I set all grey values
of pixels that do not overlap to zero and then, when using mutual
information, NOT using those pixels with zero grey value in the itk mutual
information program. But how to do this? Is it somehow possible to exclude
pixels in the registration calculation using
itkMutualInformationImageToImageMetric. Hopefully someone has an answer to
this question.

Note: To my opinion it would be very useful to be able to exclude some
pixels in the mutual information registration: then, not only areas where
no overlap is present can be excluded, but also parts where there is a lot
of noise, or parts that you do not want to use due to deformations (of
organs or tissue) that are present.

Thanks in advance,
Regards,
Jorn.

--202006272-1637312257-1041958074=:7438
Content-Type: IMAGE/JPEG; NAME="ct_mri.jpg"
Content-Transfer-Encoding: BASE64
Content-ID: <Pine.LNX.4.44.0301071747540.7438@g1-095.mcnet.azn.nl>
Content-Description: 
Content-Disposition: ATTACHMENT; FILENAME="ct_mri.jpg"

/9j/4AAQSkZJRgABAQAAAQABAAD//gBYQ1JFQVRPUjogWFYgVmVyc2lvbiAz
LjEwYSAgUmV2OiAxMi8yOS85NCAoUE5HIHBhdGNoIDEuMikgIFF1YWxpdHkg
PSA1MCwgU21vb3RoaW5nID0gMAr/2wBDABALDA4MChAODQ4SERATGCgaGBYW
GDEjJR0oOjM9PDkzODdASFxOQERXRTc4UG1RV19iZ2hnPk1xeXBkeFxlZ2P/
2wBDARESEhgVGC8aGi9jQjhCY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2NjY2Nj
Y2NjY2NjY2NjY2NjY2NjY2NjY2NjY2P/wAARCAFRAVEDASIAAhEBAxEB/8QA
HwAAAQUBAQEBAQEAAAAAAAAAAAECAwQFBgcICQoL/8QAtRAAAgEDAwIEAwUF
BAQAAAF9AQIDAAQRBRIhMUEGE1FhByJxFDKBkaEII0KxwRVS0fAkM2JyggkK
FhcYGRolJicoKSo0NTY3ODk6Q0RFRkdISUpTVFVWV1hZWmNkZWZnaGlqc3R1
dnd4eXqDhIWGh4iJipKTlJWWl5iZmqKjpKWmp6ipqrKztLW2t7i5usLDxMXG
x8jJytLT1NXW19jZ2uHi4+Tl5ufo6erx8vP09fb3+Pn6/8QAHwEAAwEBAQEB
AQEBAQAAAAAAAAECAwQFBgcICQoL/8QAtREAAgECBAQDBAcFBAQAAQJ3AAEC
AxEEBSExBhJBUQdhcRMiMoEIFEKRobHBCSMzUvAVYnLRChYkNOEl8RcYGRom
JygpKjU2Nzg5OkNERUZHSElKU1RVVldYWVpjZGVmZ2hpanN0dXZ3eHl6goOE
hYaHiImKkpOUlZaXmJmaoqOkpaanqKmqsrO0tba3uLm6wsPExcbHyMnK0tPU
1dbX2Nna4uPk5ebn6Onq8vP09fb3+Pn6/9oADAMBAAIRAxEAPwDmaKKKACii
igAooooAKKKKACiiigAooooAKKKKACiir+lyLH5u6CGXOP8AWpux16UAUKK2
nl80bWt7YDr8kCqf5VH9nEhJWKP6BRmgTaW5k0Vsiwz1iVfqBUw0+Hui++FF
PlZlKvTj1MyL/kDXP/XxF/6DJVMAkgAZJ7V1MVnbrYyr5anMqE5A64f/ABpq
wRIMLGij0AxRGGrMpYpLZHOrazPnEbcevH86nGnSbWJZcgcAd63di/3F/Kk2
J/dA/CtORGDxU3scvRWrLbxpIyrGhAPdap3kaJs2DGc54rI9FO6uVqKKKBhR
RRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFF
FABRQASQAMk1MlrM/SMgZ78UEuSjuyGirqaa5zvkVfTAzVhdPgUHIZvqen5V
ag2YyxNOJlVIlvLIRtjY56HGBWqkKR42oBxjOKcpwcCtFR8zCWL/AJUZyafK
2NxVfXnkVoWNiIw/zFiSOgxUgHGamtv4uvbpQoJGEsTUfUd5MY6KD9RmpecU
q5z2pduTzVJdjFyb3Yw4PWnHPpS7R2OPrTzGcdPyFAhUBNnJ/wBdE/k1RbSK
nVD9klGP+WifyaoCjDoP0qIrVlS2QFR2oxS7WA5BpQDx61dibmfexjzgxzyK
zL8KPLwc9a3bqFpFUjqDWLqSBPLwCOvX8KwmrM9bDy5qaKNFFFSbhRRRQAUU
UUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFACoxR
ww6g5rbBBGR0PQ1h1r2jBraMnjjGPpxVw10OLFrRMmpQaSgVojgHcUgUZyKA
KDxV3YgIyeKsWwJ3YHp0FQDFWrEkb8DPSnuIkUENyDTwm4jgj8Ke8iSDAIz1
4oWXaNuwn3qWhism3qaYkrE/dpzS4PKr+OKerc8Jt96LgPQFrWTqMuvb2amB
celTqcWz8H76/wAjULSYPXFTFasqT0RG7sgz057VXLMWzmpf3k3y4PHNK21V
2lDkdxWkUSV5cmNvnGfpWLqgYeVuz3xx9K22dexP86yNaUDycY/i6D6VnWja
zPQwj0cTMooorA7QooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACi
iigAooooAKKKKACiiigAq5Y5ZWUdjnr/AJ9Kp1YsDi5Ax94Ef1/pWlKXLJGN
dXgzV6gcc0UdBiitHueSAPtSnntQKXFNCG1oadj95n2/rWeQas2rqgfLhema
cU3sFrlySIAcbs59KFQ8DihJo5SQpzioWnSOXoSQavkbGky15IPp+dL0Pei2
uI5Y8uDuzjpTS2enH0NZSTW4yfCm1k4Od6/1qo6DdxkcelXEU/ZHYnjcv9aq
TnDjHp6Uo7sJbIetyqnhMHH96oXcsWzwCfSnAbuCSajnjJQgZ46AVoo3ZK1I
ABkgc4rP13P7jK4+9/StFdyAA/Kao+IFJWBsYwSOmPT/AAp1leJ1YeVqhjUU
UVxnphRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAU
UUUAFFFFABT4X8uZH5wpycUyigTV1Y3ic0lRWzeZbo3PTnPtUuK33PFkrNoP
wpSQoJJAFGPSobzm1kAOCBn8ua0hG7sJK7sVL65ZhJEgG3j5s8+vFV7GN5S+
FzjGSxqsNz5xW1o6qsrGRRtIG0e9eo0qMND0XajBNIhila33ssbYHBPbNXbe
RZoBIcAnOfarF3FvilSOMFcZ+uOf51lRxkZMuUUEYGcDNQpKor9R04rEvdJm
rbuu1sMQB39auxyDGQFPGelY1s/73DPvz0xVxUG0c4wawqQszkq0+V2NPzD9
kkOGGHUYzx0PSqzHcc4p0T7bGUAE4kT+TUzduAJGK5LWbMZbIfF9/oelExYM
AMgEU9BlunamT4wMMAelXF6krchZfnxntVPWFMlvgnJUEjPbpVx88BT264qt
doCoUnIOQRVSV0y4vlkmc7RQylWKkcg4NFcB7K1CiiigYUUUUAFFFFABRRRQ
AUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQBpac+YSueVPT
0H+c1c61mac+2cqTww6Y6kf5NaY6VvDVHk4iPLUYuahmjd2BX071Jn2p2a2X
u6mCdncxm8hpzC0eNpIyPatfTox5fyqfl4561QuLJzdNLGcE8jHY4pNMleQy
pJIezAH9a7Z/vI3TO6tN11zGyJEzyM1DfATQ7WbGGyAMf571JJkwg5yP9kVR
mKCTJU59CKxprW5z01Z3I2KxKFAGcdSOakjla3Rtzs+ccelQzkT7QnG3PamK
hWbbJyB613K0lqerGca0eSa1Ogs3U6U8gBG6ROv/AAKkB3c4/SqKSiHS5jnp
NHj/AL5ery4KhlPDDI4rzakOWT9TzsVRdKVmWohlz06fWmygfMMD64pyna2D
jPtTJmKkYGaiK1OVEDYzjn8Kr3iAbM479qsPksMY6VWuh9zPvVsZh3kfl3LY
GAeR/n61BVzUQCUYA+hP+fxqnXFNWkz16MuaCYUUUVBqFFFFABRRRQAUUUUA
FFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUASW7bLhGzj5hk
+3etoZrBrbik8yJG4yRk4rakcGMjsx/SjHNGPWgjNbnALVU2PlwTSQqNy8gD
0z/hVnpzU9rht/HpTjNwd0XGTi9Cnb3BjtAtwdoJIHuKlES4yNpBGck8Yp99
a/alTadrIe/TFVJIpYpdgEsgQAA9FxiuiLjU1TszphBVNnZksyM0bCIgMf4g
artayl9zOGPrjFSrdJFGTKhUknAA61MojaNZJWwpGRxjFWpSpjftKLs1YitC
xsLkOgz50YA/CT/ChJbgXCg4EYwMe1SYRbS5aI4UzxEc/wCzJVcyugOV4PSk
5RipVJ7f8BHu4aEMZSVOa1N5GDsSAeKc8QZeSPXAHNVrNi21nJ+Ze351e2gx
HG38q533R83Vg6cnHsU3QBgACeOuaqXeMoF96vSYQ+gqjdNwhBJ69qRmZ97G
DbE85XnAFZlbLASKwOQGGKxq5qy1uehhJe60FFFFYnYFFFFABRRRQAUUUUAF
FFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFamnMzW2OMKxA
/wA/jWXV3THxI646jOfp/wDrq6btI58TG9Nmjj1ox9aOTR0611HkhU9oCS4B
54qCrNluBfbx0oAlwO5x+FB69cj6VYWItxtpGtwAc7s49qh6DRlPpsRChHZT
nk5zxVZrKZZCq72QY+YsB+laxgXO4kihYATy1dMa0l5nRDEzi9dfUzI4h/ZU
5LOoM0eGYDniTp+dQuZJBliNmevSuijiX7I+VVwJVIGO+GoW2MmWWNBng8cm
sajVa6ltc9TD5s8OrxgtTNsDJGFHZD90gjPOa3GRQu4flmqA2gg4HHtUzO3l
8J29abtZJdDya9Z1pubWrGzY3/d7dqpXgxs4x1qwd7nIBFJqcezys98/0o2M
TMbnGaybldtw4znnPT1rYK8nH8qz9RTa6MT1GOnp/wDrrOrrE68LK07FOiii
uU9IKKKKACiiigAooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKACii
igAooooAKmtG2XUZx3x+fFQ0ULRkyV00bxB7GjFIjb0DjgMAaXj0rsPEejsF
X9Mi8zzeSMY/rVDGa1NFz+/x/s/1poRdABOADn6cUN8pxj+VSkBUG44pMjaO
1FtRmZeux2MvHYnNMjDOqk4OfQVav4la3dlTJUbsk9Mdf0zVTTZw0bD7zofw
A/zmuhRvC6NkrwuaEAMdqwK4IdeAPZqXzAP7wPpUpBkgcbdvzL/I1D5fzgEj
j071zx3ZEtkQ4lJwApPpUkkUnlfMFHTpQxlXkMBn2oLtt6DNNoggWLI759qj
1dceV1PXr+FXEJIOcdao6x/yx4PfoPpRYLGccD0qrqC7rbIxhSDVnAzzTZk8
yNkGORgZpSV1Yum+WSZjUUUVxHtBRRRQAUUUUAFFFFABRRRQAUUUUAFFFFAB
RRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQBrWD77ZRk5U45qyaz9LbmRCT
2IH8/wClX8iuqDujx68eWo0GK1tDOBPyB93+tZNa2hf8t/8AgP8AWrMTSZMj
5uarzwn59m4HHABOKubccjqaQqMYOefemMxLmdraHZdE7ZMjA7jv/OqOk6lF
aPKZkARlH3Bkkg+59zWxqGk/bI0WNgJFbOWJ6d/6VRXw5J9r2GTMAx8/AJHf
j161yVnWvaD0Po8vll6w7jWfvPf/AIBfTVoJ7ZmjR8GRE5wCCd3+H61J5pPT
ipI9ItYbdo41O0spbJzuIBGf1oe2ESFgSMdzW9CUXHXc8jE+y5/3WxGqMpBJ
J/KoyhLkH19KmViVGOgPU8/r0ppbLEHJ+nNas5GKsJA7Vn6uCPKyf739K0Gy
enH1qhrDn9yCP73T8KnURmYHejGPT8KccGmcZ60wMi4TZO64xycAelR1b1FN
swYDhh19SP8AIqpXHJWZ7NKXNBMKKKKk0CiiigAooooAKKKKACiiigAooooA
KKKKACirfmad/wA+l1/4Er/8bo8zTv8An0uv/Alf/jdTd9gKlFW/M07/AJ9L
r/wJX/43R5mnf8+l1/4Er/8AG6LvsBUoq35mnf8APpdf+BK//G6PM07/AJ9L
r/wJX/43Rd9gKlFW/M07/n0uv/Alf/jdHmad/wA+l1/4Er/8bou+wFSirfma
d/z6XX/gSv8A8bo8zTv+fS6/8CV/+N0XfYCpRVvzNO/59Lr/AMCV/wDjdHma
d/z6XX/gSv8A8bou+wEdi+y6XnAb5Tx/nvitfFZyzaerBltboEHI/wBJX/43
WuJbNlDLBPgjI/fD/wCJrelJ2tY8/Fw1UiGtnw//AMvHA/h/rWdvtP8AnhP/
AN/h/wDE1qaI0J8/yoZF+7ndID6/7IrW77HJyruabe4/KkA2807cR/yz/wDH
/wD61Jnd1j/8f/8ArUcz7ByruJuAOc4pryHt0qTAzny+fZv/AK1NkKAcq3Ho
3/1qal5DUV3Fm802+6ADeSvGcCqj6hFbRqlywMhBPyHcOtXEdDbsCpxnH3uv
6VU8m0u5VxASUHUHt78VNPd8y0Nko2XMPmgEgHTPsKriJ1lzjjPtV/AHy4fj
jO4f4U0xx7gWV/pv/wDrVak+xjyruVG3f/rrM1bB8rBP8X9K32jh4JU+3zf/
AFqyNWW3j8rzI5WBzjEgGOn+zRzPsHKu5jYBFAO3gDNWN9mf+WE4H/XYf/EU
nmWg/wCWE/4TD/4ijmfYOVdzN1Eb4g23lT19v84rOrfnNm8LKbe5PGcCYZP/
AI5WX5mnf8+l1/4Er/8AG65qt072PRwr9y1ypRVvzNO/59Lr/wACV/8AjdHm
ad/z6XX/AIEr/wDG6yu+x0lSirfmad/z6XX/AIEr/wDG6PM07/n0uv8AwJX/
AON0XfYCpRVvzNO/59Lr/wACV/8AjdHmad/z6XX/AIEr/wDG6LvsBUoq35mn
f8+l1/4Er/8AG6PM07/n0uv/AAJX/wCN0XfYCpRVvzNO/wCfS6/8CV/+N0eZ
p3/Ppdf+BK//ABui77AVKKt+Zp3/AD6XX/gSv/xujzNO/wCfS6/8CV/+N0Xf
YCpRVzzNO/59Lr/wJX/43RRd9gKdFFFUAUUUUAFFFFABRRRQAUUUUAFFFFAB
WvZvvtUPGQNvHtWRWhpj/K6HHByPX/PStKb1OXFRvC/YvVbsbmW38zyzjdjO
BVPp1qWA53YX0rpPLLpvbhuDKce2Kb9snXpM5H1qLH0H1ppxjqKYywbyfr5s
n/fVNa8uCCDNIMj1qIBcc8mggjuPwFAFuC7mjtJHZ2lxIgwcDjD1at9UXALx
sPpg1nov+hTZ4/epz+DVCQQRtcgVK1bLk7pHQSTxO52yoMjuCTT9wZcjcR64
GKw+gyCPzppmdehwR7U0Rc3wc8g1j65keRu5Pzf0qFLmXH38/gKr3JLBcnpn
pTAr8nrik288ClxQce5pgJyPesaVPLlZcHg8Zra49KzL9cXJP94A/Tt/Ssay
urnZhJe80VqKKK5j0QooooAKKKKACiiigAooooAKKKKACiiigAooooAKKKKA
CiiigAooooAKKKKACiiigAqzp7bbkD+8CPp3/pVanRtskV8Z2kHFOLs7kVI8
0WjcAqaHjd0xxUG7J6YqSDJJ7/Suy3c8WxMQWoAVfc+wpScDGKZ9eKYDxnvj
8DSUi4Bz3NOqW7ATIP8AQpc/89E/k1VshztAOQewqymDZS5/56J/JqiHyk44
71NN6sp7IB90A5HvTSF9yaGOevFKTxx1+laEiAYpk3Rc+9OPXkUycH5aTBkO
M+1BXjilIOKYc7hgnimgQYb1/IVU1JP3aNnocY/z9KuqTnkVHdpvtZBnHGen
pzUT1VjSlLlmmY1FFFch7IUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAF
FFFABRRRQAUUUUAFFFFABRRRQAUUUUAbFoQ9rGcY4x+XFW4B978KzdOcmFky
flOfzrRt25aumM9DxqseWbRPj1JpNvOc0vJ7UHA607t7GYf56UfjmgU0kk43
DFHqBYU4spfXzE/k1VyCeT+lTICLGXH/AD0T+TVCA+SOQPanBWbKeyFGFOQv
FIWOeBxUkUTzEpGCSBnFaEOmRFFySCQCeKskzArH7tLewPAI94IzmtKSwKOB
GSVx6UzxCAPs+P8Aa/pUtiMbFNI9CfwpxpM+tK4CYOOp4pOCCCNwPXin5phX
A9apNMaMR1KOynqpxSVPeptuCcAbgCAKgrjkrOx7UJc0UwooopFhRRRQAUUU
UAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFAFnT32
z7TnDA1qR3CoXXjOARWJGxVwRxg1pQ28pJO3JOB6Yo5nsjmqUoOfNJloXLsT
1pfNc+tRpG6kkgY9qmVsADaKxnJ3tciqqaXuoVED884B6VOoC9APyphYgE/y
FCknnBrppcq0W5xTbepaX/jzlxx+8T+TVXMcjnKqxFXbWEyWrbshWkXn8GqV
VEI2Lg5Oeea6Ibsl7IntLRbSQvu6jHpU6uhkwDk5qu5YKNx3c9hSx7yw+QgH
uOKbZJbdwp6fnWPrrBvIwem7+lXpQy8kOfxzWDrU26eLa/y88D8KcY8w1G5F
ik4psYIBFOPSokrCClHSkApcUkDM/U4z8kmP9kn+X9aoVr3yb7VuMlfmH+fp
msisqi1PTwsr07dgooorM6gooooAKKKKACiiigAooooAKKKKACiiigAooooA
KKKKACiiigAooooAKKKKACiiigArbsXWZF3HkgZ+tYldN4TfEU/K5Vl4PbrW
1F2Zy4qN4XLUdlLJlUXPHqBUiaRKeroD3BOSK0EnlcndlRTDuEhIDM3XnitX
Si3dnnJva5AulqoIkfJz/COKlQQRtmBQTjBBH/1qkEZfl+v50/APCEk/Si0Y
iuMIJikyQPmXgD2NNVVA+7n3NPdY1gkMjYO5c/kagBU/6qT5fr3qE22ypbIf
5ZXksaa8yRoSznj0Xmqkty0Ckyk+wA5JrPk1GYyKVQjvt65rop0ZS1HGm5F2
aVppcKT5eQNxyCB34qHUUhHlbBzz1H0pLK/+0M4YgEcjI7VLqhUpblRxhu30
qqkeVWsE48uhQxSYpaKwsZCdaWiikAxgGUqw4IwawyCCQRgjqK3j0rHu02XL
9cE5596iqtLnbg5atENFFFYHohRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRR
RQAUUUUAFFFFABRRRQAUUUUAFFFFABXR+EYw63hOeNnT/gVc5XU+CgCl/kZ4
j7f71OLs7mdWPNBo6FVVTkAflSbYzICfvZHajzFP3SM+9QyzxqWBb5wOg9a3
TbPHZPKURsNjOPSqdxds6AIu0565qGW+gQbriVkPTFZ39qx7iIYi59+Aa0jR
nLZFxhKWyNTzVFo5m2DLqMn6NWNdaoqzGKJACDjO3gipT9pvrKYS4WPzoyFU
dBtfP64qm8MNrMoxvbhgMf59K6KFGN2pbnRCknZPVlmeNpZxNI7H5gSvTpV0
24hAYR5yeCDzWa11Iwz5YPoB/WtiwiupzG86lY8cq1a1Nk76Dr0507KehBaa
Yxl8wrhGJ/CnaxD5KwKDkfN/StjAQBVGB6CsrXc/uMgj73X8K5JzctzjlJyM
mijmioJCikLAUxpPSmotjFc8EVm34G5GBySCKu5OO9V7xQYTgY2nP9Kc4+4z
ooPlmihRRRXEeqFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUU
UAFFFFABRRRQAUUUUAFbPh67Nr9oxEH3beueMZ+lY1bPh0J/pG9gv3cZ/GgD
Y/tn59nmfODgjb0qlfamWZwPvEcNitO5tLRmaRotzOeSy47Y4/KqU8llAS6r
GCf4VHWvSpOm7WTbPNcacZNcrbMf7O08paUtnHfjNTwWgW6K7gBz1PQVdCyX
UJ8iIgHkdsVK2jz+cHLJ0zub19K6lU6N2LU+jdiZHEGmyRRqWYMoyvPJD/4V
AthM8hFw2Bjgoa17G38qCQAcMy5P0zT5EKtlF7deK41V5XJRMXVcLOBFBaRw
MWRcDH8RzVuNwRgEYFO+Vjxim+WqNuJ/DNQ5N7mEpOTvJ3EkIJ7/AIVma+Wz
AGA/i/pWjJKEPHAPasfXCSYN3P3sfpQoisZ24CmlwO4+lMoCjrjFUkh2EJy2
cGj6gilxig9M4qhjcnsp+ppGQyIyk4yMdKXJ9KXAPbNDKTtqY9FSXClbhwfX
P51HXntWdj14u6TCiiikUFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABR
RRQAUUUUAFFFFABRRRQAVs+HrSO6+0eY23btx+tY1dP4Ls4bv7b5w+7sxzjr
u/woAmOn3M7bLd1Tg7s55q3aaHbwIplJdwMkH7oP5VvXilbcmCJQ3rgDiqKJ
I43SEZ+tddKpJQsnY8/EylGVl1IpFk+7GikY65PFWEZTkLyakVOOP0pQiqcq
oBobucZGN3zem4fyNQzyMjgc9Owq2fuHI7jmomZgPkGRUQWrHLZCAFOQQaqs
0rylcbVJ+8KshlfgYzTXHHArVaEojCbehJ9yKzteAxBgf3v6VpZYHlwvsayd
Xff5XzA43f0p3C5lnik+opc01iKpDA/Q036U4AnoKTbj75xQ5KO7KQhOPeky
xPQ49QKRpAv3B9CarStvOZJcYGMLXPLFRWkdTpp4eUtXoR3mwyBkYMCOcetV
6kk8vaAgOfWo65+Zy1Z3xjyqwUUUUFBRRRQAUUUUAFFFFABRRRQAUUUUAFFF
FABRRRQAUUUUAFFFFABRRRQAUUUUAFdj8P13f2hg4/1fb/erjq67wEqt9v3F
hjy+mf8Aa9KAOwntYp0CyLkA5GDjms4WsaHoQVPrWszBRkg/gM1RuCqux6Dr
zxWlPexx4uPupjBjHFGaha6jTgnP0qpLqC7f9Yo57HNbWPONHPyE+4qGSdUO
CKzZLwG2lb5n+dB6A8N/hVN7lgeMKMdKlaN3KadkahulTlVwf9qqU13nd0PP
QVRaUDG5j+JqFrpF6YodSK6lxo1JbItvckn7uOKrXLs20u3r7VA92x6cVTv5
3+TBABz0qXW7I3jhJdWWHeNc5bOOwqF7uNfugfWs4sx6k0lZupOXU3jhYLfU
uPfMegqFrhz3qGiotfc3jCMdkOLsc5JptFFFiwooopgFFFFABRRRQAUUUUAF
FFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFbXh3XF0b7Ruj
L+bt6DpjP+NYtFAHZx+MoN37w3OMfwxJ/jUEviKC8nk2I4GBhnABPHoOK5Og
EjpTTad0RUgqkeVnQS6jDnLkkjioX1WHnAf8hWLk0UOcn1MY4WmvM3l1NTpk
7hWwJox09Vf/AAqn/aSnr5n4AVDF/wAga5/6+Iv/AEGSqdZpau50KMVsi99r
i7q/5UhvEPGHx9BVKirGXRdxD+F6guZVl27QRjOcioaKACiiigAooooAKKKK
ACiiigAooooAKKKKACiiigAooooAKKKKAEooooAKKKKACiiigAooooAKKKKA
CiiigAooooAKKKKACiiigAooooAuRf8AIGuf+viH/wBBkqnRRUrdgFFFFUAU
UUUAFFFFABRRRQAUUUUAFFFFABRRRQAUUUUAFFFFABRRRQAlFFFAH//Z
--202006272-1637312257-1041958074=:7438--