Vascular Modeling Toolkit / List vmtk-users Archives
Thanks for helping keep SourceForge clean.
Screenshot instructions:
Click URL instructions:
Right-click on ad, choose "Copy Link", then paste here &
(This may not be possible with some types of ads)
You seem to have CSS turned off.
Please don't fill out this field.
You seem to have CSS turned off.
Please don't fill out this field.
Briefly describe the problem (required):
Upload screenshot of ad (required):
, or drag & drop file here.
Please provide the ad click URL, if possible:
Mailing Lists
Brought to you by:
vmtk-users
You can subscribe to this list .
results of 74
(Page 1 of 3)
Hello Luca,
Thank you for the answer, it is helpful. I actually need to combine
images from the same geometry and from different geometries at different
points in my work. I have some higher quality MR scans that are on a
smaller domain and would like to understand how to combine it with a
lower quality image of a bigger domain. So getting some more details
would really be nice.
Best regards Adrian
13:53, Luca Antiga wrote:
& Hello Adrian,
there is a way of combining level set files. After all, they are all
& image files, so here's how you do it:
& in case they are all produced from the same input image (i.e. they
& have the same geometry), you can go
& vmtkimagecompose -ifile ls1.vti -i2file ls2.vti -ofile ls_combined.vti
& If you have more than 2, you'll have to call another vmtkimagecompose
& in case they are produced from different input images (hence they
& have different spatial extents) you'll first have to resample them
& a common grid us let me know if this is your case
& and I'll provide you with more details.
& On Mar 25, 2013, at 3:08 PM, adrianrh wrote:
&& I am working on segmenting water in the brain for simula research
&& laboratory.
&& I am using the command
&& "vmtklevelsetsegmentation -ifile image_volume_voi.vti -levelsetsfile
&& level_sets.vti -ofile level_setsOut.vti"
&& to continue working on a project.
&& what I would really like to be able to do is to make multiple
&& levelsets
&& files and then combine them. L
&& vmtklevelsetsegmentation -ifile image_volume_voi.vti -levelsetsfile
&& level_sets.vti -levelsetsfile level_sets2.vti -ofile
&& level_setsOut.vti
&& is it possible to do this in the vmtk software? or possible to add
&& Adrian Hope
&& ------------------------------------------------------------------------------
&& Everyone hates slow websites. So do we.
&& Make your web apps faster with AppDynamics
&& Download AppDynamics Lite for free today:
&& _______________________________________________
&& vmtk-users mailing list
&& vmtk-users@...
Hi YiQiang,
the angles you get from vmtk are explained in ,
"Computing bifurcation geometry" section.
The ICA angle is the difference between the two downstream inplane angles, while the planarity angle
is the difference between the upstream and downstream out of plane angles.
Hope this gives you a hint
On Mar 28, 2013, at 9:10 AM, qiang zeng wrote:
& Dear VMTK users,
& In my works, I need to calculate the ICA angle and planarity angle like the Stroke 2005 paper. but I didn't get the result from the Paraview.
& Does anyone know a way to do this?
& Thanks in advance!
& Best wishes,
& YiQiang Zeng
& ------------------------------------------------------------------------------
& Own the Future-Intel& Level Up Game Demo Contest 2013
& Rise to greatness in Intel's independent game demo contest.
& Compete for recognition, cash, and the chance to get your game
& on Steam. $5K grand prize plus 10 genre and skill prizes.
& Submit your demo by 6/6/13.
& vmtk-users mailing list
& vmtk-users@...
this is what I'd expect from -seedselector idlist, not -seedselector profileidlist.
Are you sure it was the latter and not the former?
On Mar 29, 2013, at 2:07 PM, Anouk Post wrote:
& Hi Luca,
& I already tried that and it didn't work. It gave me only a partial reconstruction of the centerline (a very short one).
& On 29 March , Luca Antiga &luca.antiga@...& wrote:
& Hi Anouk,
& -seedselector profileidlist -sourceids 0 -targetids 1
& (hoping you have a recent enough version of vmtk)
& On Mar 25, 2013, at 7:36 PM, Anouk Post wrote:
&& Dear VMTK users,
&& I am using VMTK centerline to calculate the centerline of an stl file as follows:
&& vmtksurfacereader -ifile /Users/.../Phantom_Best_Catheter_now.stl -f "stl"
--pipe vmtkcenterlines -endpoints 1 -seedselector openprofiles --pipe
vmtkcenterlinegeometry -ofile /Users/.../Phantom_Best_Catheter_now_output.dat
&& For all the different input files I use, the seedselector option always leads to me pressing 'q' '0' '1' and then it works. So I close the viewer, and I think I set the source id as 0 and the target id as 1.
&& Is there a way to make this an automatic process so I don't have to type q, 0, 1 each time? I've tried using seedselector idlist but I'm not sure how to use it (it does not give me the right answer when I use -sourceids 0 -targetids 1.
&& Does anyone know a way to do this?
&& Thanks in advance!
&& Best wishes,
&& ------------------------------------------------------------------------------
&& Everyone hates slow websites. So do we.
&& Make your web apps faster with AppDynamics
&& Download AppDynamics Lite for free today:
&& vmtk-users mailing list
&& vmtk-users@...
& Best Wishes,
& ------------------------------------------------------------------------------
& Own the Future-Intel(R) Level Up Game Demo Contest 2013
& Rise to greatness in Intel's independent game demo contest. Compete
& for recognition, cash, and the chance to get your game on Steam.
& $5K grand prize plus 10 genre and skill prizes. Submit your demo
& by 6/6/13.
& vmtk-users mailing list
& vmtk-users@...
I already tried that and it didn't work. It gave me only a partial
reconstruction of the centerline (a very short one).
On 29 March , Luca Antiga &luca.antiga@...& wrote:
& Hi Anouk,
& -seedselector profileidlist -sourceids 0 -targetids 1
& (hoping you have a recent enough version of vmtk)
On Mar 25, 2013, at 7:36 PM, Anouk Post wrote:
Dear VMTK users,
& I am using VMTK centerline to calculate the centerline of an stl file as
& follows:
& vmtksurfacereader -ifile /Users/.../Phantom_Best_Catheter_now.stl -f
& "stl"
--pipe vmtkcenterlines -endpoints 1 -seedselector openprofiles
vmtkcenterlinegeometry -ofile
& /Users/.../Phantom_Best_Catheter_now_output.dat
& For all the different input files I use, the seedselector option always
& leads to me pressing 'q' '0' '1' and then it works. So I close the viewer,
& and I think I set the source id as 0 and the target id as 1.
& Is there a way to make this an automatic process so I don't have to type
& q, 0, 1 each time? I've tried using seedselector idlist but I'm not sure
& how to use it (it does not give me the right answer when I use -sourceids 0
& -targetids 1.
& Does anyone know a way to do this?
& Thanks in advance!
& Best wishes,
& ------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& vmtk-users mailing list
& vmtk-users@...
Best Wishes,
-seedselector profileidlist -sourceids 0 -targetids 1
(hoping you have a recent enough version of vmtk)
On Mar 25, 2013, at 7:36 PM, Anouk Post wrote:
& Dear VMTK users,
& I am using VMTK centerline to calculate the centerline of an stl file as follows:
& vmtksurfacereader -ifile /Users/.../Phantom_Best_Catheter_now.stl -f "stl"
--pipe vmtkcenterlines -endpoints 1 -seedselector openprofiles --pipe
vmtkcenterlinegeometry -ofile /Users/.../Phantom_Best_Catheter_now_output.dat
& For all the different input files I use, the seedselector option always leads to me pressing 'q' '0' '1' and then it works. So I close the viewer, and I think I set the source id as 0 and the target id as 1.
& Is there a way to make this an automatic process so I don't have to type q, 0, 1 each time? I've tried using seedselector idlist but I'm not sure how to use it (it does not give me the right answer when I use -sourceids 0 -targetids 1.
& Does anyone know a way to do this?
& Thanks in advance!
& Best wishes,
& ------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& vmtk-users mailing list
& vmtk-users@...
I do understand how I can open VMTK from MATLAB, but not how to actually
run a VMTK 'script' from MATLAB. When I use the system function it just
opens pypepad. Do you know how I can actually let it run a VMTK script?
Best wishse,
On 29 March , Luca Antiga &luca.antiga@...& wrote:
you'll probably have to use the system() function in Matlab to call vmtk
& the way you call it from the command line.
& It should work.
& PS: a couple of random references on the system function
On Mar 25, 2013, at 6:06 PM, Anouk Post wrote:
Dear VMTK Users,
& I am using VMTK to calculate centerlines of stl files. This is part of a 3D reconstruction process for which I use MATLAB. I really want to incorporate the VMTK part of the process into my MATLAB code.
& Does anyone know how I can call VMTK from MATLAB and execute a VMTK script?
& Best wishes,
& P.S. My VMTK code is:
& vmtksurfacereader -ifile /Users/.../Phantom_Best_Catheter_now.stl -f "stl"
--pipe vmtkcenterlines -endpoints 1 -seedselector openprofiles --pipe
vmtkcenterlinegeometry -ofile /Users/.../Phantom_Best_Catheter_now_output.dat
& ------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& vmtk-users mailing list
& vmtk-users@...
Best Wishes,
you'll probably have to use the system() function in Matlab to call vmtk the way you call it from the command line.
It should work.
PS: a couple of random references on the system function
On Mar 25, 2013, at 6:06 PM, Anouk Post wrote:
& Dear VMTK Users,
& I am using VMTK to calculate centerlines of stl files. This is part of a 3D reconstruction process for which I use MATLAB. I really want to incorporate the VMTK part of the process into my MATLAB code.
& Does anyone know how I can call VMTK from MATLAB and execute a VMTK script?
& Best wishes,
& P.S. My VMTK code is:
& vmtksurfacereader -ifile /Users/.../Phantom_Best_Catheter_now.stl -f "stl"
--pipe vmtkcenterlines -endpoints 1 -seedselector openprofiles --pipe
vmtkcenterlinegeometry -ofile /Users/.../Phantom_Best_Catheter_now_output.dat
& ------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& vmtk-users mailing list
& vmtk-users@...
Hello Adrian,
there is a way of combining level set files. After all, they are all
image files, so here's how you do it:
in case they are all produced from the same input image (i.e. they
have the same geometry), you can go
vmtkimagecompose -ifile ls1.vti -i2file ls2.vti -ofile ls_combined.vti
If you have more than 2, you'll have to call another vmtkimagecompose
in case they are produced from different input images (hence they
have different spatial extents) you'll first have to resample them onto
a common grid us let me know if this is your case
and I'll provide you with more details.
On Mar 25, 2013, at 3:08 PM, adrianrh wrote:
& I am working on segmenting water in the brain for simula research
& laboratory.
& I am using the command
& "vmtklevelsetsegmentation -ifile image_volume_voi.vti -levelsetsfile
& level_sets.vti -ofile level_setsOut.vti"
& to continue working on a project.
& what I would really like to be able to do is to make multiple levelsets
& files and then combine them. L
& vmtklevelsetsegmentation -ifile image_volume_voi.vti -levelsetsfile
& level_sets.vti
-levelsetsfile level_sets2.vti -ofile level_setsOut.vti
& is it possible to do this in the vmtk software? or possible to add it?
& Adrian Hope
& ------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& _______________________________________________
& vmtk-users mailing list
& vmtk-users@...
Hello Greg,
first of all thanks for your patience.
Replies are inline:
On Mar 21, 2013, at 1:44 AM, Greg Fokum wrote:
& Hello everyone,
& I just came across vmtk and was wondering if I could you it for my work. My question is this:
& I have a tubular 3d surface ( length of colon) similar to the vessels you show in the tutorials except that there is not branching only bending.
& I would like to map this 3D surface to a plane (or regular cylinder) using VMTK so that I can then create try to subdivide my tubular 3d surface ( volume) in to sections ( alond the central axis direction) of equal volume based on mapping uniform strips on the mapped plane back to the 3D surface.
Note that we don't have full conformal mapping in vmtk: we can create a longitudinal harmonic map, but the missing part is the harmonic mapping in the circumferential direction. (I'm just mentioning this in case you are looking for a fully conformal map as in the literature).
& I see a two step procedure:
& 1) first create a mapping from the 3D surface to a plane ( or in fact a regular cylinder such that volume is conserved)
& 2) cut sections of the cylinder and relate them to sections on the mesh.
& Is that possible with tools in VMTK ? If not please help me with some suggestions. I can send a sample of my mesh in wrl/vrml format.
It's quite easy to create a mapping to the plane of a cylindrical structure in vmtk. Just follow this tutorial
and even if you have no branchings, just try running the scripts the same way.
For a mapping to the cylinder, then you'll have to code a bit. Actually a lot of the basic tools are there in vmtk already (differential geometry classes for the mappings etc), but it would require you to code, you can't just do it from the command line.
Feel free to seek more specific inputs on the list.
& thank you
& ------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& vmtk-users mailing list
& vmtk-users@...
I'm not sure about the problem you're experiencing. I'm developing vmtk on a OSX 10.7 machine and I've never
bumped into this problem.
Did you try holding the Alt key while pressing the left mouse button to move the planes?
On Mar 24, 2013, at 7:52 PM, Gawlitza, Matthias wrote:
& Dear Vacular Modelers,
& I am experiencing strange problems with the VMTK Image Viewer and the other render windows - if I load a 3D data set (e.g. CT-angiography) into it can zoom or rotate the images and I can adjust window levels, but no matter what I try I can not move the image planes.
& Maybe anybody of you has an idea how to fix this problem? I am using a Mac (Intel) with OsX 10.7 (middle mouse button is assigned via BetterTouchTool to my Magic Mouse). If I try this in a virtual machine in Ubuntu 12.04 [precise] the thing works (okay there are some other problems but that is another story...)
& Maybe anyone has an idea or has experienced the same problem and found a way to fix it?
& Greetings, Mat
& ________________________________________
& Von: vmtk-users-request@... [vmtk-users-request@...]
& Gesendet: Donnerstag, 21. M?rz
& An: vmtk-users@...
& Betreff: vmtk-users Digest, Vol 65, Issue 19
& Send vmtk-users mailing list submissions to
vmtk-users@...
& To subscribe or unsubscribe via the World Wide Web, visit
& or, via email, send a message with subject or body 'help' to
vmtk-users-request@...
& You can reach the person managing the list at
vmtk-users-owner@...
& When replying, please edit your Subject line so it is more specific
& than "Re: Contents of vmtk-users digest..."
& Today's Topics:
1. conformal/parametrized mapping ( 3D surface mesh to 2D plane
) with vmtk ? (Greg Fokum)
2. Re: vmtklevelsetsegmentation (Luca Antiga)
3. Re: Computing centerline geometry (Luca Antiga)
& ----------------------------------------------------------------------
& Message: 1
& Date: Wed, 20 Mar :17 -0400
& From: Greg Fokum &gregthom992@...&
& Subject: [vmtk-users] conformal/parametrized mapping ( 3D surface mesh
to 2D plane ) with vmtk ?
& To: vmtk-users@...
& Message-ID:
&CAFh_-ON1xGCGS9-0-tVuZifC8KHPoOyQNgiCobWdC-uXUG-pVA@...&
& Content-Type: text/ charset="iso-;
& Hello everyone,
& I just came across vmtk and was wondering if I could you it for my work. My
& question is this:
& I have a tubular 3d surface ( length of colon) similar to the vessels you
& show in the tutorials except that there is not branching only bending.
& I would like to map this 3D surface to a plane (or regular cylinder) using
& VMTK so that I can then create try to subdivide my tubular 3d surface (
& volume) in to sections ( alond the central axis direction) of equal volume
& based on mapping uniform strips on the mapped plane back to the 3D surface.
& I see a two step procedure:
& 1) first create a mapping from the 3D surface to a plane ( or in fact a
& regular cylinder such that volume is conserved)
& 2) cut sections of the cylinder and relate them to sections on the mesh.
& Is that possible with tools in VMTK ? If not please help me with some
& suggestions. I can send a sample of my mesh in wrl/vrml format.
& thank you
& -------------- next part --------------
& An HTML attachment was scrubbed...
& ------------------------------
& Message: 2
& Date: Thu, 21 Mar :48 +0100
& From: Luca Antiga &luca.antiga@...&
& Subject: Re: [vmtk-users] vmtklevelsetsegmentation
& To: RichSchmidt@...
& Cc: vmtk-users@...
& Message-ID: &B07DFE7E--BC3C-9F@...&
& Content-Type: text/ charset="us-ascii"
& Hi Regine,
& did you try inverting with vmtkimageshiftscale?
& vmtkimageshiftscale -ifile input.mha -scale -1.0
& If you need to re-obtain positive intensity levels for some reason, you need to pipe another shift scale for the shift part:
& vmtkimageshiftscale -ifile input.mha -scale -1.0 --pipe vmtkimageshiftscale -shift 4000
& Hope this helps.
& On Mar 18, 2013, at 10:39 PM, RichSchmidt@... wrote:
&& Dear all,
&& a friend of mine came up with a dicom data set of the lung. She wants to extract the trachea which is represented with a lower intensity as compared to their environment (since it is filled with air). We tried to use vmtklevelsetsegmentation, but it didn't work due to the inverted intensity distribution, I guess. I tried to invert the dicom data set and read it into vmtk afterwards, but I wasn't successful.
&& Did anybody have the same problem and could tell me how he/she solved it?
&& Thank you for any help!
&& ------------------------------------------------------------------------------
&& Everyone hates slow websites. So do we.
&& Make your web apps faster with AppDynamics
&& Download AppDynamics Lite for free today:
&& vmtk-users mailing list
&& vmtk-users@...
& -------------- next part --------------
& An HTML attachment was scrubbed...
& ------------------------------
& Message: 3
& Date: Thu, 21 Mar :50 +0100
& From: Luca Antiga &luca.antiga@...&
& Subject: Re: [vmtk-users] Computing centerline geometry
& To: qiang zeng &zengqiang2006@...&
& Cc: vmtk-users@...
& Message-ID: &F397CC31-B034--5BD4BB5789EA@...&
& Content-Type: text/ charset="iso-;
& Hi Qiang,
& the script works, I'm attaching the result and a Paraview screenshot where the Glyph filter has been applied to the Frenet normals.
& Note that I specified a higher number of iterations and factor than the default, to diminish the effect of noise on the finite difference approximation of derivatives.
& vmtkcenterlinegeometry -ifile wf--L-3-1-3.18_centerlines.vtp -smoothing 1 -ofile foobar.vtp -iterations 1000 -factor 0.5
& On Mar 20, 2013, at 11:47 AM, qiang zeng wrote:
&& Dear all,
&& When I am following the script "vmtkcenterlinegeometry -ifile foo_cl.vtp -smoothing 1 -ofile foo_clgm.vtp", I cannot get any centerline geometry. But
I am dealing with split centerlines, the script "vmtkcenterlines -ifile foo.vtp -seedselector openprofiles --pipe vmtkbranchextractor --pipe vmtkbranchgeometry -ofile foo_clcg.vtp" is OK for the same vessel (as an attachment).
&& Did anybody have the same problem and could tell me how he/she solved it?
&& Thanks.
&& &wf--L-3-1-3.18_centerlines.vtp&&wf--L-3-1-3.18_model_cl.vtp&------------------------------------------------------------------------------
&& Everyone hates slow websites. So do we.
&& Make your web apps faster with AppDynamics
&& Download AppDynamics Lite for free today:
&& vmtk-users mailing list
&& vmtk-users@...
& -------------- next part --------------
& An HTML attachment was scrubbed...
& -------------- next part --------------
& A non-text attachment was scrubbed...
& Name: foobar.vtp
& Type: application/octet-stream
& Size: 370380 bytes
& Desc: not available
& -------------- next part --------------
& An HTML attachment was scrubbed...
& -------------- next part --------------
& A non-text attachment was scrubbed...
& Name: Screen Shot
at 6.13.26 PM.png
& Type: image/png
& Size: 97833 bytes
& Desc: not available
& -------------- next part --------------
& An HTML attachment was scrubbed...
& ------------------------------
& ------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& ------------------------------
& _______________________________________________
& vmtk-users mailing list
& vmtk-users@...
& End of vmtk-users Digest, Vol 65, Issue 19
& ******************************************
& ------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& _______________________________________________
& vmtk-users mailing list
& vmtk-users@...
On Mar 22, 2013, at 9:48 PM, Evan Kao wrote:
& I am trying to segment a particular set of 3D-cine images with vmtklevelsetsegmentation, but the perhaps due to the shape of the vessel and/or quality of the images, the suggested settings of 300 0 0 1 isn't working well.
I've tried playing around with the level set deformation parameters, but the resultant shapes feel just as arbitrary as if I had simply used thresholds directly.
Can you recommend any guides or references that can explain in more depth (i.e. quantitatively) how the parameters affect the deformation, and/or also provide an explanation of the additional options available in vmtklevelsetsegmentation that weren't covered in the tutorials?
The level set parameters are those of itkGeodesicActiveContourLevelSetImageFilter and itkGeodesicActiveContourLevelSetFunction.
The main reference for level set methods is
J. A. Sethian, "Level Set Methods and Fast Marching Methods", Cambridge University Press, 1996.
but you can find a lot of material around.
You may want to take a look at my PhDThesis (good old times) at lantiga.github.com/about/
If you have noise issues, I suggest you to use a non-zero featurederivativesigma , which expresses the width of the
gaussian kernel used for performing partial derivatives for the feature image. It's expressed in physical spatial units
(e.g. mm) and it's roughly related to the minimum scale that you'll end up detecting in the image - i.e. the noise below
that spatial scale is going to be filtered out, thus not affecting the evolution of level sets (again, roughly speaking).
Feel free to post a couple of screenshots of the problem you're having, so that replies can be more to the point.
& Thank you,
& Evan Kao
& ------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& vmtk-users mailing list
& vmtk-users@...
Dear VMTK users,
In my works, I need to calculate the ICA angle and planarity angle like the Stroke 2005 paper. but I didn't get the result from the Paraview.
Does anyone know a way to do this?
Thanks in advance!
Best wishes,
YiQiang Zeng
Dear VMTK users,
I am using VMTK centerline to calculate the centerline of an stl file as
vmtksurfacereader -ifile /Users/.../Phantom_Best_Catheter_now.stl -f "stl"
--pipe vmtkcenterlines -endpoints 1 -seedselector openprofiles --pipe
vmtkcenterlinegeometry -ofile
/Users/.../Phantom_Best_Catheter_now_output.dat
For all the different input files I use, the seedselector option always
leads to me pressing 'q' '0' '1' and then it works. So I close the viewer,
and I think I set the source id as 0 and the target id as 1.
Is there a way to make this an automatic process so I don't have to type q,
0, 1 each time? I've tried using seedselector idlist but I'm not sure how
to use it (it does not give me the right answer when I use -sourceids 0
-targetids 1.
Does anyone know a way to do this?
Thanks in advance!
Best wishes,
Dear VMTK Users,
I am using VMTK to calculate centerlines of stl files. This is part of
a 3D reconstruction process for which I use MATLAB. I really want to
incorporate the VMTK part of the process into my MATLAB code.
Does anyone know how I can call VMTK from MATLAB and execute a VMTK script?
Best wishes,
P.S. My VMTK code is:
vmtksurfacereader -ifile /Users/.../Phantom_Best_Catheter_now.stl -f
"stl"
--pipe vmtkcenterlines -endpoints 1 -seedselector openprofiles
vmtkcenterlinegeometry -ofile
/Users/.../Phantom_Best_Catheter_now_output.dat
I am working on segmenting water in the brain for simula research
laboratory.
I am using the command
"vmtklevelsetsegmentation -ifile image_volume_voi.vti -levelsetsfile
level_sets.vti -ofile level_setsOut.vti"
to continue working on a project.
what I would really like to be able to do is to make multiple levelsets
files and then combine them. L
vmtklevelsetsegmentation -ifile image_volume_voi.vti -levelsetsfile
level_sets.vti
-levelsetsfile level_sets2.vti -ofile level_setsOut.vti
is it possible to do this in the vmtk software? or possible to add it?
Adrian Hope
Dear Vacular Modelers,
I am experiencing strange problems with the VMTK Image Viewer and the other render windows - if I load a 3D data set (e.g. CT-angiography) into it can zoom or rotate the images and I can adjust window levels, but no matter what I try I can not move the image planes.
Maybe anybody of you has an idea how to fix this problem? I am using a Mac (Intel) with OsX 10.7 (middle mouse button is assigned via BetterTouchTool to my Magic Mouse). If I try this in a virtual machine in Ubuntu 12.04 [precise] the thing works (okay there are some other problems but that is another story...)
Maybe anyone has an idea or has experienced the same problem and found a way to fix it?
Greetings, Mat
________________________________________
Von: vmtk-users-request@... [vmtk-users-request@...]
Gesendet: Donnerstag, 21. M?rz
An: vmtk-users@...
Betreff: vmtk-users Digest, Vol 65, Issue 19
Send vmtk-users mailing list submissions to
vmtk-users@...
To subscribe or unsubscribe via the World Wide Web, visit
or, via email, send a message with subject or body 'help' to
vmtk-users-request@...
You can reach the person managing the list at
vmtk-users-owner@...
When replying, please edit your Subject line so it is more specific
than "Re: Contents of vmtk-users digest..."
Today's Topics:
1. conformal/parametrized mapping ( 3D surface mesh to 2D plane
) with vmtk ? (Greg Fokum)
2. Re: vmtklevelsetsegmentation (Luca Antiga)
3. Re: Computing centerline geometry (Luca Antiga)
----------------------------------------------------------------------
Message: 1
Date: Wed, 20 Mar :17 -0400
From: Greg Fokum &gregthom992@...&
Subject: [vmtk-users] conformal/parametrized mapping ( 3D surface mesh
to 2D plane ) with vmtk ?
To: vmtk-users@...
Message-ID:
&CAFh_-ON1xGCGS9-0-tVuZifC8KHPoOyQNgiCobWdC-uXUG-pVA@...&
Content-Type: text/ charset="iso-;
Hello everyone,
I just came across vmtk and was wondering if I could you it for my work. My
question is this:
I have a tubular 3d surface ( length of colon) similar to the vessels you
show in the tutorials except that there is not branching only bending.
I would like to map this 3D surface to a plane (or regular cylinder) using
VMTK so that I can then create try to subdivide my tubular 3d surface (
volume) in to sections ( alond the central axis direction) of equal volume
based on mapping uniform strips on the mapped plane back to the 3D surface.
I see a two step procedure:
1) first create a mapping from the 3D surface to a plane ( or in fact a
regular cylinder such that volume is conserved)
2) cut sections of the cylinder and relate them to sections on the mesh.
Is that possible with tools in VMTK ? If not please help me with some
suggestions. I can send a sample of my mesh in wrl/vrml format.
-------------- next part --------------
An HTML attachment was scrubbed...
------------------------------
Message: 2
Date: Thu, 21 Mar :48 +0100
From: Luca Antiga &luca.antiga@...&
Subject: Re: [vmtk-users] vmtklevelsetsegmentation
To: RichSchmidt@...
Cc: vmtk-users@...
Message-ID: &B07DFE7E--BC3C-9F@...&
Content-Type: text/ charset="us-ascii"
Hi Regine,
did you try inverting with vmtkimageshiftscale?
vmtkimageshiftscale -ifile input.mha -scale -1.0
If you need to re-obtain positive intensity levels for some reason, you need to pipe another shift scale for the shift part:
vmtkimageshiftscale -ifile input.mha -scale -1.0 --pipe vmtkimageshiftscale -shift 4000
Hope this helps.
On Mar 18, 2013, at 10:39 PM, RichSchmidt@... wrote:
& Dear all,
& a friend of mine came up with a dicom data set of the lung. She wants to extract the trachea which is represented with a lower intensity as compared to their environment (since it is filled with air). We tried to use vmtklevelsetsegmentation, but it didn't work due to the inverted intensity distribution, I guess. I tried to invert the dicom data set and read it into vmtk afterwards, but I wasn't successful.
& Did anybody have the same problem and could tell me how he/she solved it?
& Thank you for any help!
& ------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& vmtk-users mailing list
& vmtk-users@...
-------------- next part --------------
An HTML attachment was scrubbed...
------------------------------
Message: 3
Date: Thu, 21 Mar :50 +0100
From: Luca Antiga &luca.antiga@...&
Subject: Re: [vmtk-users] Computing centerline geometry
To: qiang zeng &zengqiang2006@...&
Cc: vmtk-users@...
Message-ID: &F397CC31-B034--5BD4BB5789EA@...&
Content-Type: text/ charset="iso-;
the script works, I'm attaching the result and a Paraview screenshot where the Glyph filter has been applied to the Frenet normals.
Note that I specified a higher number of iterations and factor than the default, to diminish the effect of noise on the finite difference approximation of derivatives.
vmtkcenterlinegeometry -ifile wf--L-3-1-3.18_centerlines.vtp -smoothing 1 -ofile foobar.vtp -iterations 1000 -factor 0.5
On Mar 20, 2013, at 11:47 AM, qiang zeng wrote:
& Dear all,
& When I am following the script "vmtkcenterlinegeometry -ifile foo_cl.vtp -smoothing 1 -ofile foo_clgm.vtp", I cannot get any centerline geometry. But
I am dealing with split centerlines, the script "vmtkcenterlines -ifile foo.vtp -seedselector openprofiles --pipe vmtkbranchextractor --pipe vmtkbranchgeometry -ofile foo_clcg.vtp" is OK for the same vessel (as an attachment).
& Did anybody have the same problem and could tell me how he/she solved it?
& &wf--L-3-1-3.18_centerlines.vtp&&wf--L-3-1-3.18_model_cl.vtp&------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& vmtk-users mailing list
& vmtk-users@...
-------------- next part --------------
An HTML attachment was scrubbed...
-------------- next part --------------
A non-text attachment was scrubbed...
Name: foobar.vtp
Type: application/octet-stream
Size: 370380 bytes
Desc: not available
-------------- next part --------------
An HTML attachment was scrubbed...
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screen Shot
at 6.13.26 PM.png
Type: image/png
Size: 97833 bytes
Desc: not available
-------------- next part --------------
An HTML attachment was scrubbed...
------------------------------
------------------------------------------------------------------------------
Everyone hates slow websites. So do we.
Make your web apps faster with AppDynamics
Download AppDynamics Lite for free today:
------------------------------
_______________________________________________
vmtk-users mailing list
vmtk-users@...
End of vmtk-users Digest, Vol 65, Issue 19
******************************************
I am trying to segment a particular set of 3D-cine images with
vmtklevelsetsegmentation, but the perhaps due to the shape of the vessel
and/or quality of the images, the suggested settings of 300 0 0 1 isn't
working well.
I've tried playing around with the level set deformation
parameters, but the resultant shapes feel just as arbitrary as if I had
simply used thresholds directly.
Can you recommend any guides or
references that can explain in more depth (i.e. quantitatively) how the
parameters affect the deformation, and/or also provide an explanation of
the additional options available in vmtklevelsetsegmentation that weren't
covered in the tutorials?
Thank you,
the script works, I'm attaching the result and a Paraview screenshot where the Glyph filter has been applied to the Frenet normals.
Note that I specified a higher number of iterations and factor than the default, to diminish the effect of noise on the finite difference approximation of derivatives.
vmtkcenterlinegeometry -ifile wf--L-3-1-3.18_centerlines.vtp -smoothing 1 -ofile foobar.vtp -iterations 1000 -factor 0.5
On Mar 20, 2013, at 11:47 AM, qiang zeng wrote:
& Dear all,
& When I am following the script "vmtkcenterlinegeometry -ifile foo_cl.vtp -smoothing 1 -ofile foo_clgm.vtp", I cannot get any centerline geometry. But
I am dealing with split centerlines, the script "vmtkcenterlines -ifile foo.vtp -seedselector openprofiles --pipe vmtkbranchextractor --pipe vmtkbranchgeometry -ofile foo_clcg.vtp" is OK for the same vessel (as an attachment).
& Did anybody have the same problem and could tell me how he/she solved it?
& &wf--L-3-1-3.18_centerlines.vtp&&wf--L-3-1-3.18_model_cl.vtp&------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& vmtk-users mailing list
& vmtk-users@...
Hi Regine,
did you try inverting with vmtkimageshiftscale?
vmtkimageshiftscale -ifile input.mha -scale -1.0
If you need to re-obtain positive intensity levels for some reason, you need to pipe another shift scale for the shift part:
vmtkimageshiftscale -ifile input.mha -scale -1.0 --pipe vmtkimageshiftscale -shift 4000
Hope this helps.
On Mar 18, 2013, at 10:39 PM, RichSchmidt@... wrote:
& Dear all,
& a friend of mine came up with a dicom data set of the lung. She wants to extract the trachea which is represented with a lower intensity as compared to their environment (since it is filled with air). We tried to use vmtklevelsetsegmentation, but it didn't work due to the inverted intensity distribution, I guess. I tried to invert the dicom data set and read it into vmtk afterwards, but I wasn't successful.
& Did anybody have the same problem and could tell me how he/she solved it?
& Thank you for any help!
& ------------------------------------------------------------------------------
& Everyone hates slow websites. So do we.
& Make your web apps faster with AppDynamics
& Download AppDynamics Lite for free today:
& vmtk-users mailing list
& vmtk-users@...
Hello everyone,
I just came across vmtk and was wondering if I could you it for my work. My
question is this:
I have a tubular 3d surface ( length of colon) similar to the vessels you
show in the tutorials except that there is not branching only bending.
I would like to map this 3D surface to a plane (or regular cylinder) using
VMTK so that I can then create try to subdivide my tubular 3d surface (
volume) in to sections ( alond the central axis direction) of equal volume
based on mapping uniform strips on the mapped plane back to the 3D surface.
I see a two step procedure:
1) first create a mapping from the 3D surface to a plane ( or in fact a
regular cylinder such that volume is conserved)
2) cut sections of the cylinder and relate them to sections on the mesh.
Is that possible with tools in VMTK ? If not please help me with some
suggestions. I can send a sample of my mesh in wrl/vrml format.
&?xml version="1.0"?&
&VTKFile type="PolyData" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor"&
&PolyData&
&Piece NumberOfPoints="2774"
NumberOfVerts="0"
NumberOfLines="2"
NumberOfStrips="0"
NumberOfPolys="0"
&PointData Scalars="MaximumInscribedSphereRadius"&
&DataArray type="Float64" Name="MaximumInscribedSphereRadius" format="appended" RangeMin="1."
RangeMax="3."
offset="0"
&DataArray type="Int32" Name="EdgeArray" NumberOfComponents="2" format="appended" RangeMin=""
RangeMax="24"
offset="21068"
&DataArray type="Float64" Name="EdgePCoordArray" format="appended" RangeMin="0"
RangeMax="0.996"
offset="31360"
&/PointData&
&CellData&
&/CellData&
&DataArray type="Float32" Name="Points" NumberOfComponents="3" format="appended" RangeMin="23."
RangeMax="58."
offset="36396"
&DataArray type="Int32" Name="connectivity" format="appended" RangeMin=""
RangeMax=""
offset="66416"
&DataArray type="Int32" Name="offsets" format="appended" RangeMin=""
RangeMax=""
offset="66432"
&DataArray type="Int32" Name="connectivity" format="appended" RangeMin=""
RangeMax=""
offset="66448"
&DataArray type="Int32" Name="offsets" format="appended" RangeMin=""
RangeMax=""
offset="71676"
&DataArray type="Int32" Name="connectivity" format="appended" RangeMin=""
RangeMax=""
offset="71724"
&DataArray type="Int32" Name="offsets" format="appended" RangeMin=""
RangeMax=""
offset="71740"
&DataArray type="Int32" Name="connectivity" format="appended" RangeMin=""
RangeMax=""
offset="71756"
&DataArray type="Int32" Name="offsets" format="appended" RangeMin=""
RangeMax=""
offset="71772"
&/PolyData&
&AppendedData encoding="base64"&
_AQAAAACAAACwVgAApz0AAA==eJzt2vk/Vc//AHCue213IURlSUKShEoqNS97RcpSpKKSIpJki8qSUHYVUoosiWStpKQU2YoWS4lCCZX1Xjt9Xn3/je/7p+fj3HvOzJyZ17xm5vE4XssYLn5i3HDz0nP3zHXccESJsqrAlxtKNSS+52dxg3+08ssd97hB6fnLfS9fcMPRCDvZgxXccEO70VaGww2WYf4S9uPcoF9jf2HJGgrsFNlcWBBKgW63oImZSArscD69a1cFBWzyjDkPZiigMXw6pWc3D5zybKyfseOBn3FjCV+P8gBkroxvceYBW++KldQiHih9/OpX7nMeuOvv92tDFQ+82HY8d+FbHgjwW3LZ6RsPLPm5d7cdhQqb7R9EfNekwg2jLZ2pOlQ45/1GbZs+Fa4I/nx3CGWfXPB9/XYq1Mv1sVkWVNhmtkCdjTbqNZo8RW3egZj4Tiq4Rd0Z10Dfnix2yDalApfeuROfTahgK2puUbiVCiuWKwezUY7RV7dL2/D/uq028sZUgAseicvRWf/RkQbUp9p69Rg+3/7gUqQkelAlMOo6/n5BatsufXzO6tAihwNoWyKf7g8jKpSL+6R2oRYXuzwtDajQzdgjJkeooETeWj/bSIVw36vC/4zfbPvOG+0JoCedWUcFXRfNX37oQ2tO+zt8f82AZc9q1lJhv6LukadrqLCGuyD+5ioqLLV2HLyO5osNWN9Bb9roVs5HtbreyUmjzw9VpHGrUiEtzFvIHa9DCkLaetWoEPVlqr5LnQovzbZXUDSoEGp64jYP2nu77ZUO+rOvcY8emiVmm1CJCrhyy7xFQy4adJ/A+vV4rJ5uw/b48QjsNUbVZHTWJ6KlH2dm/3n6QKy1NbZ7Vo7asweV35vktAjfR9IuqNAOzVSuO3sQnZG/tk52PRWo8crzw9DVuif076NZUdtY9Wi6YIb0P4v7nzHfoW/vzk+W0KaCkPzpnA//lCn8fnozFW7RH8aGY3zYWapvj0CPSkhm6epSgbbhU38K9rvw8O/gOvTsYQOxelTi7vZAK0Mcj5f75K3R5I/MZ2PoRp+h3locv6idL7ekYHzIdwUkK2F8vR/viVuBchvOU16M465166Gzxg7s74vPoqIwrhTvM5tOWWI/3Lp+ffEuKjieN3+yFTX1kDJfuZsKKys2k1b0wKGbXCdtqBC53denAR2xnKas20+F3xa/+qttsd1bixheB6hQ8jjDWvYgFfZpjjVaoS/enCyab0+F3N3d1jdRw8yTZyePUmFDidPQB0eMm/O9Be1onuvFupNO+L6nzhv9RFcu3dW32wXvU6nZ1umK80E902OeOxUa9neXrkcrX+z84e2B7y17x/2sFxVUhBUX+6OTW27fM/Wmwvqs3RkjvlTgFR24NIauOMlft/ksFc47vVvYew7f97mtVD+6W6g1oSMA423DyjGvYIxLqZf+xahUtflTxQtUaE04EHUZzblk+fBMCBWcRgwyn4dRISZrEZRcxPibcoh7dIkKazfWXEwKp0KSQeF4N7rP9UXCoUgqRGweveCD7rVamKkQRYWwXQtO/zOGc0FlEn9/WtxZwI/XwnfCf3Dw+ntyQeY+9OXxY1pPo6kwWKUjd+Ay/r5wu817dGJNWnvKVSp0RfGcDoqnQrb5vXc/k6jwJ9Q2O/Y6FT7dyZj4cYMKO9doMkbRM+6a6TduY/wJZqwoSKOCb1zqxg3peJ2e//Qj6nX8WXYv+oPFJV9whwqvfDkiLdlU0HinGqaRg/N9vunPHWiAgtNR/1yc38qCH6ULqLCVX1rHAR2SPWls8gCfy9yoGPwQ58dZt1fkEc5/i+inT9BMVs/XLNR551LJuhKcN07RLW/Rpkcim4seU8F+1cTqf3b7LD9wFyW3Ow8U4P93nm7R7EFXtfcMnUWp98VDz6ANZq36dCzvC7P/QSDWd3KjnsJxrH/X/b3knxu3bAw+hb7+0Xnk3/VHk0uUOmynYq6NZS26Zz9/ryx6VSthx1F09/IlNXP52K6wmMSz+J56IQF7jNCRvMG3P7Efzm09wPmdhbLPHjuIxkx6zHfE/hz+bN4giSZ4fXBhoA3TDV1/kqnQsVPtyrlE7PeIl3ZTCVT4lWlr/zYO5407v/MgjnO0ydNKExxfiovoQoJxNGO9YHY4FPuLWJaPY5w5izeqVWH8Xdr3UdPzPBX8mW+sfmOcLtHNdmtHY11zEtgYx7umXcoCTmM8FO36e9oH56HQWrdjqOPe2GgR9JOvcvhenB+renrr9U9inlnxrcHiBMbJoUWHdHE+wasvA6bO2P6AWxH/rCvg/ztwDOc5f2ttAM7H1VIemXpHcN14oPpZC6VUP/f+jvPXc3+E1OdDVHjW6c87bofzKXe65ibOf8GH3nJmqNoaO6fifVS47Ssj0bkX89MOp6Ewm39xt0y9ZQ8Vyozyt+ZYUSF4vleREDqfrz9lHkrbPumjjJYuKBUfwnwTzl2WEYd56MHM2OM0FET+TEri7xVHZOz/5aMFKwrNZ/H3PP+Ynf9Mr2rcNI0+nbr54Q+WE/r+5qU2axxftfJFM+janLxAF5t/6y9/7Gls32puak4geuBG3c9/1+EvXAJL0aT4hkERzGv6AnyyLfh+qg2HFC1Qw6UJ0Xao+M4OSizmtTPK3Caf0Z4tggcsUA3VEyaHsH8uVD5WcjxMhbtu8sEx2G+5ibPNcdifHt5yftVoPWVT3hrs55MznYcMMb/tyXLQ3Yc60AU0bx+ngrH805HVOE7vThYvWY9mfE0bCv7nptSMWnT6adRrS8x/ctKTI36nsF8kTcra0S/Nfm8f4Xi/audUHsI8aBS95l05ut3XQ6IKLTS9G9SAcRGXrL+Pg3ETSqENvPCjwqZ9KXwVqNLN0sHPKJu0H6BinMk97r/YjnHWkSzoLexPhWNOC3kvBmKe3Kjuq4pxWR2xq8EE3XbuxfJraG7vh+VKmB+XsSfO9WNcD7baSLtifnx/t+pmHObDiQ2jRtfQutOz207gPHhjk+FCx/xmZqVZ5oM2C5U5KcdQwUR76tL6WCo8GffekYge0+rXccX543VVNDsZ819tmWt3M6poYqp6/ArmC+1RbRnMh/F1UxIPcb49U06r97uGz+3jHefHvPh20FLjHMo4lvjsI5p+oiVxLc7T8uOv5n2+iXnX6X6sVyoVFsvEfOxCrY7sv1GG89lD9JOzIObFmXMGa0oycV3I74tyxvn/getW8b88EVzA/WY95sPcrVwfrtzHfZ4ra4kJ5pPNTUqhoeiyvum8r5h/gtUvRlVjvgqoXhz1pRTj7Ajnz4ly3L/xSXqbVOB4msVudqvE9dPqq4NE9b/9Yt/qPY2YXxYVrZF/R4XDEw63O5uoYC1zqKC0GfNvcUhI9Scsz+DrPMl2zJPEJ47Whfuo8wsPdv/A95NPyDPrQaV/dfT2YXmrwivuD+F4mEv9uTZOBZctOWXX0VtGDy/+QLd5iD4lk1juAtnzqdMYf8av/H/N4LwWd3ueNUuFgdcxwW/QHW3qyjJ/MV9/oAmuQ/f+mKcEXDQIKHF+9pKHBstjMm+3oCLf952T56WB85/D1hHoAC1iYjUfDTaEsgM0BWjQKFv6+YMgDUIt63bLMGgwPXDxYy2TBm+ySzP2CNHA7f5j3aWiNMj6lNWxDQ1ILomMFaNBrmaseSKauPyGSxS6fYSlmIFyVmQencT7on8NfX0vQoPU9izJs8L4vLTynAjKazrdmYHlrj8zNZDIosGIsOJ7faw/cIcSJQbb4/9m/wwL25mzgt1KRRMKN8sswvdgevkbB3LTQGV/gmb5HBXmUWXVnLF/dqsaiq+bwH3FVhUeb+y/bukXtlVjuO/y6A41ZWP8CfSquWN/l/zOK/3wiwoiYslBUujhfSHvInsxf+Yy1N9+p8IRdduUhG7c93kqNWng+HWe7b+29huu+/6WVk87cL9jaxN7GvX7RZNkt1GhKP+GWzOOv1vSUR4T9LSeaNtFjAv5yJ0SORgn3wO/plbV4DxLWrolD+NKSK4x8z5qEnrr6Ty0N+64W/ATKlianMhm4DrL28RX74rr52HRs8mteRinoRc7m+/iujM9Ks+LcT55O+L4twzM27mGKz1xPriaOpOyFJy/3jyWd3BfcbMh+U8Dzif5NX77k3C9S7c/U+SK83G5mbtiP87by0V7sgJwnm/gMtqrhnLX5wQn4vwvPGgZIY15Yd8fA39LXPd2xSYcPYZ5RSGiQsoC3b1xZloa92PbVlHd/61vpGjm2jTu51aKnBC8jXnM3fSu8wDmyRF5vocJmDelrXIO5mPelRn84eWDefje0/Vmr3DdCasUb0vB9ebbk+aRJnMqfOZ5Gh6Npp+p4VLA/XDNd0aoJKqVZSuthfvlgXk3B+u2YP3KEUX2uJ/28H3wsAbPaXN+Vka5uC/33sm1+TJg+UeGE/TxHHTQSV5pfBOe1+Kg8QTa4CFveAr39QuGHhXxb8D8k+30vBH3/RIrx48sQPvivZ15Uf7p7VmmeH44+62T82T1v31dSMs3PM8o+VstNMDzzVPuHNNCPPcIF1ipXEUvXA3x3q6C+y3ZG5GLUZZHcKswysf9x6xzBa6f15cOnMJri+oN576jS1rSAobRi840fmt8Pogu0puBfuqRGfvno80pKUXoyzmpjZlY3yGXX8UD6Ptdr/ZexXboxlVHyeJ5iStV6W0bmjifuTcS22ny5EyJDJ6fOP6PV0ziOSmMS/H4PjwfFUjdk76MFr/jL2Lje20RvyMYie9ZaH7OeSH2x+YpqzczaHNkpJELnnccRrsCX2P/vbbfHduM/To3rWDQpIfrXV2J+BD2t0nBbzU69n+lqFB8LNp7NcCyEc+3tJ1udVU4Tt2qQ0k0PLf80cppnMJzsjavgTQP7gcUL389x4vrPnnx4rAgrt82TWsHOvA8UmjT4aaMcRGzarveF1yXZag/Fy7H/U6H2kzYMoyf78+N0wo9sT/1J/0foAbRm9Sm0K4rShd4MN4aTf5mhuH54XBSqORbXAfDPv4ISMD1zpOb6euL+7Q4mb2d/Lj/590+1tmA61K2+ud3l3Ae8Ms8TLDD/bmYZL3oLHqnJuAODdenn48dl87H9amutTHz9i3MG+nZ0iK4L3+2USBSEo3aKvhtyz18X4VzlZXobPu770ycp0z7WtHFqJPNuS8E15Hkg3fsGC+p0CKU9eDyGyqYp59LmdeA+VtpYuGvFnyfK9LeqbhOBO9fR/X+jPHB2Xm76Cs+55be8wHzzLXLWise/KHC1G6ZgDTMXz+zK1+aTOH+qO3ni1OoSqRMwzvMg7An5ogDlQaPXCikAPPp2Npv1/dj3k0qNnY0X0iDm9+jvZqVMZ8/8tf5pkqDioZ7MhPo7VKjEK51NEgZE1yfoUWDg9fcNxegz1wzvpzfTIPRHalOJltoQN9J2ahtQoPrKZqWh9FN00aHeHfT4LHF17SV9jRQrTVcEHGcBvOuDCeno4l/cvYX+9JArGDxh8EAGmx+ta4/PpgG+4Ie1G26QAM+7R12BqjRisr2oTAazKhYus+iuQFFnQ7hNOihd0V6oiJey/iKk2hg3mXlFZCMvyfqzz+YR4MpuVuJ19F7fS7upQ9o4OEgbGv3iAbprxterXpNgw5TywcnUXn1wmfCH2mQ3UnRkkWVcjzcxFtoUFpo+2DLN1xf1u0t5ftJgw/3R5PnoU/KyGzeEA20q4wZjSjXyqbIxdP4Xr9P/lD456ufUmbcvDAbm/4wApUVDFDoEOKFxw/mM5TFecGW0npeD7WS41y4JscLA1dv7WhV54VM48rd39A0+qlteet5Ie7QLYHvaC2l6/Pnbbzwp5LlsdqCFzoE47KPoGtPLBNZ68ILvQflLn9Fl7ff9o724AXNo+Y/j3jygqeg7sI0lCf8rtLfMF7gbgheR73IC47DNdmDkVi+l0A7PQrLYRsPvrvCC7ArRWLBVV4obfcuf3gN28Xifv0U3d1/UjLpFi+YFM2kUtJ4IWGsWE8EpR0biDYs5AWV0I+XMp/ywrZmg1ff0K9SP7vGK3gh4EdFI7zEchOFOJUfsbyD8Lse3R5RbHO5jRdOaDh8eYj+zebk9v3GcnzOWr/4wwuTFPv5uqO8wHv+qE44+vflc3/NSWy3lsoOe9T0utZrCSE+ULP3l1ggygduL1pCbOX4QCdOe76DOh/wn4t8GIbamffpzmnzQZhZqFzqdj747XXGU9aaD34emgQxJ7zOTForfZ4PHn0T+5t3gQ8Gja321abwwfVziZaNaMVrSv3NTD4oDj/+8G0pH0wULlAIqeODhcUf78m+wXKvn68sbOKDIccJa54RPvghYTm3Cm1WuOaTweED+XuGO54y+cF1+uJfTxY/hPFY+R8Q5YctQXnpccv4ITHeqzVRhx8mf3tr3ED9PB+7gC4/xAsp8Obu5YfNYglL7qN7dr1zWbePH5KG3X71ufJDZvha+a+n+cG/PsXqgy8/9F/JmjeNhm8I1Q0J5oe2io66wZtY39LlOmceYz05m+KFP/HD46Inkq6f+UGDOpHWMMkPXTa5uicWCIDXjmDPalkBEBT7+ChJXgCOahQuuUQEwNCoWdBEVwBUlTQuiBwVADL8fGWQuwA0rN+jtcNPAC6cLvL/ckkAfpf81szNEoDxtWBgUygAvhueUJzLBcDW98cjTqMAaHdkrPjaIgBileffxgwLwLHq16wpfkHIlRGI42UIQmDK9wVXWIKgESya3iYnCPJ5W5+c0xQEh9TLZQbaglC56sSXcVNBCC0/9bfqsCDU18gFzHcShHOdc64XfAQh5Kb5YUaMINi4RUgMxQnChHVz5JV8QZCxNOI/XC4II2Z/b25rEgR9WDvE/CIID+X1DV79EgS9X/0hnHFBkC3f7pZCo0NIZlf8d1TV+J7lPDodDqfs5m6aR4fKfZfWt4jQIXp/1EH5BXRYu4zrFLcCHRi0mI4gtOL2ziOr1ehgbbH98AZURDaHO1adDj98rZ/3adJBtMhHtUqLDgO1Oskn19Mhs3Q8RVyHDq1zXY0RBnSYi/MyrzGhA3hr3siwogMzLy+2EB2NN73UYEMHrQuTx+8cpQOX++Y5Jxc6FJz2VjY4Q4eXNvdn+sLpUNJhvCQhDtvbk0FNuIWeO7ngJiq8155jV4jtYoyM5pfQocvq5pL5L7HdU4ueHUetevUtD3ykwzPFh+JCTXRoOVl46ncPHT6PpNty/6JD7eHrr5x/02Fi+SUnmQk6vP0k2bVrkg4XfkBKHo0BNhJbl0eLMOBR/KDNQXEGzJxNjo9CXzAsE2lLGfDszncXfgUGvPIOMh5YxgBmpNErWQ0G2E6+cpbaxACHz97v7hIG3F8umK5pwoCT0waCBB3LNYuYs2bAZ68cETd7Bmwv3+2w3IEB0RTeLT5HGPDwZruolzMDSk5xdui4M6A+V2DQ8xwDWifc5xeGMyAi38t2DpVfGXvnTxQDom5dTSRXGLDhfOfvoXgGZIgrPOxLZMDqyO3rlVIYcPst7iGzGbBks0xqeQ4DpmILeFYVMMCQnnV/ZzEDXGqUe588YMBRn3WNDejk16OxuY8YwBM5m/8AtemQDTYvZYBS/9hvF3Rjg05jeBn2g4aHaDRqw+V94slzBvDf9X9pW8cA5azsCXf0sPa5x9pNDPCuVJw73sKAoJIMd75WBlifXSA+84kBAsWZGz51MqBi9Ohw5ncGDK8dWznQw4Dsb04lFf0McPzmoNOBDl5JCl7zhwHldtHS29DP19im2mwG/JR4q8A9zoAad4a/+hwD6OVjmd//Yrv87Bo+0phwOSKmw4GPCd+ik1+cRmue3z97W5AJX/b49f2zSvvMSzqdCc/Fj9xaK8SEpSBkpibKhEQvt1dnJJiQ0u99oU0a77d0O2Yqy4Qla87JXpVjAg+9qywD7Q4SddqqwISGw0YS+5SZEK49vtABXRuzlkSrMMGmNdMmX50JHyQ8dDevYcLOQO+sn2inhpsFYx0TvDnf6AqbmQCs5lMdhAkmD10rAnSYkBBRwd5hwIQLQ3tuOKOKZWvqZ9A2sSDrOVTj8NMGfyOsd4fj+VA0yS/FgLaDCYIfyzJ40aJs9uoC9N67PzXMnUx4y2ynPjRnwjM9Wc/VFkywjH/l32TNhNOHL9BK92B9FrS2F2iA4qx+5H4mTBxrir+Klrdvhat2TLDwPpGfgGqkF022oly3x322HmDCj6NTGbQj2K9/mLULXZkgZzKVrnqSCSyxHapGKD99854Vp5lQeFXOQgXl9Tx+aLEvlp/uTl+HVgwHdFf6MeFPhppMDWqsyz/3/gwTnqzfff8Xujpe93hGIBNEeB97OIcyIUPwodOuMOz/VlErO3Rakb/K4BITIpav+mOEdh8V+24ezoSuIyf4dqMKLsPH7KOY8LmjJEriCvb/1wPStleZsKl+lc4p9I5CnaT/NSaoDuf7ZKO+ZncuaScz4UzbwwtpGUzoF1fS9Mhkgmft/nV3s5ggZvRzZHUO1nfDoF/1PhPmTa3q2YvG0NWWns5nghf3EfecAiYMLq9IEC/GfmkzM7iNHgjhXSb0iAl+0T0jr9C7/qNBgY+ZsPhClCCtDMdHa41q4XMmlEmtq/4nw9l353f0ynZz9x60PG37/skKJhx8Pa/kn9UConkKL3HcsmrWl75mQpZYPkWhngl739t1taMnRFrWVDYwwW7Nx7+GjUz4tazON/8jlv/EGApRl12X+yObmJAqcz7wKvor1dFBoIUJm6V91QTRCWG/4gZ0bFlTm08bE35/Sr3/t4MJ0TLqxZSvOM7yL685o8N9uS9XdzLBf/OzJq8unC+1G9UDu5nwVZweZvqDCUcynm9vRRlt28rqe5hg3bN0rUkftq9dz/47ajR+zaDpFxOyc7//fjyA89FqJmflIBMW9K7ddnQUx+lRtK4TepzPOvEZG/tByZxejsbF3vimymFCqNxlLXvUoGzp35pxJpw9zshoQM3luq+sm2LCyEmBqi+oW1701qAZJvzlWb9q2RyONzCDorlYUEl9olqLikzM597JwwJmfOo4Ly8LYtrCpQTQbaG0ORE+FqgcEz09hEovPid6WJAF5/dVaD9CvQJfhFsxWfA0yD5tXIgFxflv+CfR1+yNSpnCLFBMd6FfFmFBZp7LyWeiLHApzdg9hF7c4RhxTIwFHkPP21dJsECq4nVL3kKsv6VdNx+9t9baecEiFhyx/6XaLMmCayPO2iwZFryybQ3tQTfeb4s+LsuC7AVx/tPo/OXrH61cygJR0dz6GnkWUL5OhhgrsOCJ6sW7+xRZ0JFM6VBXwue31xreW4HlUzu/tqEvXXyTzFRYMEvp5plbyYJ5IZWFK1exYKBFZ0hEnQU8TIufXqtZ4DozP2ocZUt6JNquwWuL4uNL17IgsWSJhhJae9vr8x1UOUrwwqAmlmNQ1XlAiwW6PoN1dRtYYLmJc68ebfwVG35nIwsOOegpDGuz4PTPxnm0TSz4abcKvFGTK1k31Aj2S8zaFXa6LPggvH1IQI8F5Ywfezv1WWD+eSWzG6W+ejvfzogFTS2psYfQpSe/hotsZUHNgqFyMXSFw1UbVzRRs6ZlizELPIULWzpNWLCkJB5Md+C4vijs+GfP3fE7vaiEcL9PixkLNEOVbrWiphuOKRw1Z0EaI+OItwULVqeqWyjsYkFFWv4LYyusf6mT6kHUJ2Wto/EeFjxO/iK+B60vUhzj2ceCnYdv9HXux37pbu1XtmNB6okU1VR0y83K9wcPsKC0pDIu7CC+V8XCDMHDOM5HypQFjrKga83AyVPoTNe+4d2OLPhFNxd75MSCWyFqjsrOLNwfLd8ELiwwrmp3oLmyYMzR2lAD3TZIW1KP9trVapq4scBqSO2hrDvW85WyzdKDBftO9g14oL9Xd4zKerLgRpNjUQkaVD667JU3jhfx2991mgX0U8I7uXxZUCW4rkwUFb59i6mNKsvw60T5YfkFDN9zZ/C90pWSlvnj+Dw8vIuNVp2iLh1HVbwlDysGsGAd+Sa9LpAFh9fc3HwqiAXHmm/OEwrGeRCampmI6u274M9zAeNQcufezSEscD9hFLs5lAXB13j3wEWcN41tc/vR6m11vYxLGHeNV6Ul0LCmqBs24Tiu2xavl4hkQY6zb+6ZKBbsHym8dhsN2jq7aAidUf6UZBKN4/m8Lf8PuslkycrsWBaoxW0p+n2ZBW+zkj9Rr+B42rM/HEYdr7Wc6EWNx3zvcMezQFz8fejlBMwPb5hF/zxYNZWWhUqH/z6anciC0cjhe+evYTvvFivcTWIBS0Kg1/E6zr9ImWfOyRjvTumBr9D9Ce71p26xIH19bl8I+kbT6UYtqja7avn6VBZoub6dfnwbx0et8lAfup666yZPGs6DiysmtqXjvDXj07RA34Vd5fRksKDZf+vqN3dYYDPQm9OMfj+9RyUkiwVrppk8d9E3Cz2SetHzqXxeT+9ifCUfDS29x4Lcj2l3uHNZIGgg5/EeFS85IcRBV3+L5h+8zwKxsMjTv/NZsEeTn2t3Id5ffy3NBT2fYyF0Ds0eWlOuVMSCH6/PrZIvxn62zfuSjq7XOhn2EIV28/yTD3Dcb1Wd3PyIBf3pCaQKbSu6sPV1CfZzj/rzi48xPt6pPbhfyoK7L4NOlqD5bP8wrqcsSLDkydBHz2+RlBUrw/fv3a6hjTrL1e5wQQdZ+SrZaPGry1dKnmE709+MPUFjzOTN5coxH/7dOPH5Oc6T5u5ByRcs+Hrx5SVTtGT76mrlCvx9tENY9SUL5j7GP4hDtS81C1hUsmDlSIiqCyryLfpcBfpFTmjNmyocv/kfn+XUsqAsRvB2HurwTuOQRh3G51D719x6FrT83PptewPO9/A2Y5n3LJBzFxUIQukHFeYt+ojx+eaLTxHKfcPttHALC6w7Fj6TQZWPqGszPrGAHOPjWfgZ8+5GFeu+NhZcjjyzKbAD43tJ8fZp1HjygaRnJ+bLkx9+qHzH+NLcaKTew8LzZmPs8z5cd5ZFcdL6WVBAbJsO/cJ6p5WH/+lYusvRDXV9pLPG/zfG6aPAgQ0DLLjEo6C3GV16z24BfRD7/7h3SNoQC8zawcRsmAUaq5fYmqMhZ4xOZaN13NqGcqMs+KzWzXTl4DwY26X+HU08MDATP4Z5ivA/T0D3O9j4N6PpVbtW842z4Prdu92x6KzdkTrJCRY4eSu38kxh/5Y6X74yjfU7/f3VgyZe0d5iOoN5MXg0LAdV97e+0IrWGgiEt82yYPmWfLdRVEmrcfWqv9h/tBCFftR9Py3nLLcQbNeymnyANsr29nWj2VfllbkoQnDDb2ZvC7rolbaAIo8QyOotSfinyusLuVtQOx7XR2dQtcVbhTOoQsAqvmH2gBfL6dgtVo2GDFqz6v9dF/fuH0XTk5O+PeITgqKdnUtG0L35oklMfiGYjDVT34ieqth7qgkNnFCSWiooBH0tsU3yqOV5vW9RaKrT0s1pqBxc/PkLbWtm9PTR/7WDMp7HFILjhwYV36KTv2R2/kD3rOmgFrGEYIqe+NZDSAgY8joX7IWFwGbT9dRwdGsZLSMAvb3bha8PLeP2Gl0yTwhqOTeueIgIwcru5h2B6AdLL9F/Ssnpx7mjz87TdnmjC57w1v1BowT7rD1EsT6vBeb+qPa+/sibaFnFgZoUVMlgLucimhwVcTYVtfHPD2aICcGqkLVTH/CaQh6viUDVZgS0HdGzN1QHXVB942UaOqjx/pfL72A9eeWlC8NQDZ8TxtLoYeF5pv9UWLi5fTFa+v0Xz3b0pVmHxnd8D5vg7qbV6KXT2VFZ+H4dRbVrG7EfsoSP5/9GpUqq1ing70PB48sBr93KdpjkYn9RlG79tUfvm6g+50U9uuM2B2J/qsk+mTVD83mF0r+gj0/ZRkyi5hMsnV2oVqORWT5DCI6t8XWwRN0qNV0HcZzepEfr7cZx+nD4atQJdHPVVadz6BdxEyKGhvN2q9XjuGsVR8jZoXVNyylyqPm8BxUDGDcDB0Lj4mjY7ujUBluMs6e7ip8UoqJ2N345YPxVL2Q5RWCcfjFb+KCESwh43Rf8n00ejcX/nPd+ttEJ4/r9Sju/Q2hWm3+VA8b/ua/hL0JxXjhc1XV8jG7KbE1wRUuXr0qUwHmV4xipIIJSTFRUzuC8UywVqLyNrm+3P9SDpjZsdX04ifnicnNJOqq/89KrW6gguRWvgvO17Mh9LwN0NOW6ryv+Hvgs4o0J+vuY7u6POM9HLknk30cN02D9QZznSXVyB/egBVObrxZifrgVBod0OP/2vz4XhdAdVuF3GajHoOEBUzYLolpKTqljXtkT8+517AjO58thizzQicnHCgWYjyrKos9UYn7ao8onfwltvT79Iw5d6JwWXYzmCrr3TGEee0MW1xuggaZbrQL+YD1RN1vl8PpW9Y5kK7weqH603Q/zYOK9ObtSlK0eybUf82Waaduh5ehLe7fd6zCPyjT//OmF+bYzQrTmLuZjw6J5j/kwP89zV3+p3cyC+57+AuVvWbA7sXYkFtcF09BrI9K43lTWNiisxvXNevqj7NYIFkRv3+Gng/uiO99XWyfgPkqZO39GzRTr6Wyaq9fA/Z32uW0maMXqzI1XZplQ2/Vhdf57PPeNPYuxw/PZeyP90/p5eK42TTSQw3NywP732hQPPN8kalt+cmJCx9/X8m/2MYE7KtkoVgLPe8d2DEvxMmG3lMta/Q8McFl0tSv9KQO02B0Z4WEMKDI+s7LmAgPPS5uMtc8w4GFgKO+S3QzYFqDi91yEAV9+Xppfh04sjAgqmKNDY2JQW+ITOgiYSHaIlNBBmp82vtODDumLHByHd9Hh8Y6lh9PKZojM889VjV0zRH+5gWbJ1VkS9tNjXqT0HKHaqQfsm5wjGxUsM+9Q/5JbEfulW3L/ko9eFZrV4lzgX3Qkr2EJF+i/z12isYcLXqT3T6je4YIOg4ANz3K5IPr0eEYPmwvOZCgo8y3ghr7dQX38h7jBRNeYLPTmhum7RR9WXeSG+ROub4rCuOGc67eSkHBu6PK8XjwTzQ2pNg4u2THcAOfXyJAEbmgJFjHovs0NFjCQeaCQG/Rh71B7BTdkPlsbsayFGx7Gqn2LQFNkf3qI/OaGhuYVj55TKVA7PnP/ggwFAvcVffopS4F+wx/B87Up8Map3/qMHgXE2hTL42wpQFt+1mVjEAUy3tm9yEmlgEVWi7d5BgUWp/1xPlhAgfmaRdu/11IgSF9o2ZuPFCiU+VbS1E6BBZYm4iY/KNBil1D+/jcFvgcmhJE/FFC+vkigFw294vQyZpICM0ZiW+fPUqBborafxs+D+wDmulImDzhWDcwdZPGApGytepcwD8i7+iSslOCBsv4jK9RRE0lb1r1FPHDPTHP6KXrsbfkZVVkeeN7jmHZQngdmw8v/KijwQMqhB6aty3igbo2udPZyHohj9y8eV+YB9hXGSMsKHhCmbV7yUIUH3hs8OGGJHlzJF6KL3tw/rsnE5zLvNTW4KfLA566Dr+5jeeYHQn/64rXUOl5F53//H3lUfAbtSUg8vFcJy9HRPjEnygM1F5rDW7D9nju47osPU2CRdPbmZ70UKNE5JtXcRQGpyyukND5RQMv62WT4G+xnY94Pj2socH42vmu2BO8PMgxTTcf+OaBlsCiZAk9kzYsexlDAT+VQjeV5CmwU/tWoj6ofXO8d5kcBuyeD7edRo6f9GwV8KSB+Z2VhmSsFiu1l5wrsKGDq8GnpnBkF+FrJ2WZjCpRKLPT21qBAw9s1Jvbog09B4qGqFBgWVLdIXUiB7MWh12UoFJCXod76NcANOd0v+iUxDktyg3cZSHLD05/qfAxubli0qcw98CMXbOnQ2ri1hgsWHLxuO+7NBbddjE2eeHLB3OATraJFXCDx+UiiDj8X1L4t8Jfn/CXDYqZV+5//JVJcSywFn/0lDFmlhltqf8mGpKUfDzyZIyu3T7zSfzRHfkcdk3u1c47ofXSKSn01S/bkz9iOP50l3reEWHnrZsmilonAUsVZknD+tPzxvBkyPnZT3fj2DPF+1lVkND1NrkmI3vjnzZqFp9TQ4DNdptd3TBPepdtWKRtOk1g9ryurcqbIVRiXDrw6RQTLo0RHVaeIral7dXHrJCm//UIksWKSJE64iuZbTJL4hy/CElQmiYzT3SHaskkS8XJLPO/0BIle1OA4XTdB/v7JoW9MnCA7K41Pbdk3Qaoj3uftMZ4g9TcVK/LXTpApw8eXNCbGCdVqbdDi/nGSEWdofKNpnOzvO5eVkjNObowdJavDx8mSY3b+X/eME+6dngV62uPEIiVjr5PcOEl09K6OGx4jHTzMXPg1RuZPva1P6xkjimNPXDQ+jxEtG7eThbVjxKIhPs6gZozw5PX9WV4wRpJHPj7rTh4jPpzI3PCkMaKzQLxHK2yMxN949LXNe4xs5M3OWOk8RiTK9L++OjxGktJXdn40HyPmtDsrbPTHyJLfx4YXaY8RFbme/S4bx8hlrgkf4Q1jpPhFgNLudWNEbEN+5h71MRI81Mp3RHGM3I6nudfIjpEBqzzOL6kx0jzisiRswRiRDzf5mCI+RrqknlzbMX+MnBWRP6jDHCPK8vLtcvQxwv1jj9plgTHyQmFJshf/GJlXlnjwBt8YORgjK36Pd4xsWpkJLyhjZMzw4hIrtNljB8Ofe4ws4ysvPY2+s0zqlER1T0m6zZvmkDHbqE+PORzyi/n8VCWbQ0bl2tq5BjjESofnQ0E/h5j2DW4s7OOQUJMFdWu7OSQ25XajWTOHiH5+uiCqlkMWBEyMhtRwiOuInI9GGYcc6nhpnviYQ3jXJ894FXKIotn5V19zOaSRufn30ywOqXdcrbAzk0MULqwQfHKLQ2hu+zwakzikZ4u8dko8Pu8s2hMdi/XuKA3yjuCQtscWQk/PcoiZiUn5MT8OOa/4SlnvNIdon6Vra7pwyLXcZLe3Dhwis/f7mnPovnXf1YUPcojO1kWFB6w4ZKbUYqbCkkMEZwdWaBqjt1ufGhEOsVT980BlM4dMagj/WajOISas3r898/F+2xw9I3TF6QjmoADeJ+0mKELlkGW3vkoPzbLJ2kmxZLdRNhFU+F3T0ccm5p8nz25oYpNHSrSl4zVsIiK/acfup2yyLvfbq5sP2eRBn90R7Rw2ebjlLVdQGptsHtGf1k1lk8fnE67vvcUmJ8WSBqOD2WTvkyWxX4LYRKPt/c5YPzZRzBkNP+HEJtL2ROPEQTb51hPW9c6WTdpVdRsrzNjkz2f/9o7NbKJb1aIaoI31ZwZVzq5nk/imvlt8mmzyrFQltGkBm1Tk2feFzmOT7TOndpnOjJLH1jOdp6ZGSbtEsPvxP6NE5VGEnNK7UaI7kcWj9maU9Op+mX+xYJRs8m31c742StZd3qEuHjpKZu88lJByHSVSqx5yvXMaJed/sXfH6o2SLTp7TvzWGSU+XR/LV8mNks83Fxo9ER4lktJrRS0E8XkrqVBZNLeoYZ5c3wjpN9jvw/kxQu5LBNQGN40Qg6RlqWvKRkh8o/OU5KMREv0w97LBgxES0zxxXf7+CClae8JA6foISZwN+Fx1bYT89Kp5Yuc3Qty57tbaovHV83U9fEeI3jGmc7vrCFHYvzX5A0p77u9zaM8IybRYQ00zHSG+hq1b120cIYsT+GqDlfH+w/EaUrIj5EBzhNVm8RHis9//Yy5zhAT/8XA0oI6QQxFDXKFjw+Q9T2G3e98wOZGaH0zpHiZ8m3RebG0ZJgaXBN79eDdMgoyMYnZWDhNtpvnQ5Ith8rPzvGrKs2HyOehxyvWCYZIxsfCacPYwqX1drH0kc5gYlimud701TCKyja+tu4nPf+Z+U3ZtmNzWIyu80eo9JvVLLg+TB38Wv9mEWnxzK49F3XesimhBVe0YSo/ih0lcwNakLXj/CdXjUkbXh4lyoKQb+/YwOb9d/bPTXSy/IvfP8L1hoiK9/oHGw2FSNi1Z1/Z4mKxqVxSTLR0mX6hL02Kw3bZ+thaa8Pw+S6r88Jsc5h0rT+QZ9Y1zDZ0Gnad61nmFRFHeezHMTftecZxnKwXZJyXd/nsD2bqle00kbIGUi1PyY4QiYDvblrRUbITZFL27fLjJA1P6ryheRGSHbN/MT5a0fIUQeD/Ckch20SZhEDJiOkI4XuYGw1QjRbFDJmbUeI2Iqk0289cVx3lmRfOztCpIYq4U3MCGn7dEptScIIUb4fLlFxC8tJv11Vj0YFpv/pyBwhXTY2XSVZI2TmeOcLl4IRMkv7W/y3fIScuDaygad+hHwaLtIqQD+si9djtoyQLUVPzq36hu2++918N8alx/vUtsjfIyRE0czjxuAICTpS2hIyOkIYvtXRJyijRINl9mgLdZRsPaczekJglFya+BX4SmiURJr+Xp4nNUqmfC/Hf5UeJVXm0i1usqNEvcPpNk1xlHz3740UQ68nfU7yVhsl5MQLvmiUvXz7fX6NUbKxoq5QVWuU3NH6QLNGhTcMGpbgfIr36WKkbRslV2z8l9YZjxJeecWEx2aj5C37Wu2AxSihl1UP0I+MktamGxXHT46S5jiJ3EcBo6RtngjjeiD+LirO339plFQeYC3RTxsla8wWXLmYju26YSBVWDpKLlxYsH5p+Si5US9ULvxxlKx6funrNcwHcZSzn3LGR8mhwLA5v9lRYulGXQ1/R8mxcZ9Lzahsk9aQABfmn+Onvkpxs8ndgtn4UVQy1aNSTpBNfqm8qluGvh+c5ZKhs8mWR/azZxhsEulYoFsthM9ls0NmFrGJ/oPfIypSbLItbp1mL5puoObXL4P50IwdFr6MTXJqJwYOLGeTngC6zkc1NuH5U3uhaA2bOP5JvB6Eee64bzX94gY2KQiMVNtgyCZNVvmz67eyiaoGn0WzCZv4qBzKLdrOJozGxaITu9kk6dSl3G972GQ2LsbM8wCb8G40MVjmzCYTEm4PQtAj/BbPHhzH96hwv67nyiYG54P0M1GHNHfn+afZxJLFZZmB2t1u/qDryyYxz8LS355hE5mdnnziF/D/xxfTalFHxTe/SASbGHv7h+ihEfydj7qjsR7XblphLJvoxDzauy6OTWbUTO69wutyU5vulkQ2eVez0kfpGpscPnOAsRR96VJ9xvc6m0y/z4jpTmGTn8cHauZQi/NDt/+tI4cE72zyz8D1SLx4JDCbTQ7k2K+IuMsmnYtqQph4LXaUe9vJXDbpd6KqmN/D8eiveLwC16EMk/3XdPLY5OtEcX0j/n/mka3hc1Sr1Hwr5z72w/J3O3YVson4Utq1vfm4jkld6XRCJQUcSvbg7z4bZpftRS9sphU9L2KT6JvCih4F+N6Fjne5H7HJ2fqLhv9UVVvl24nr3+9z8VX/TJgvIqGCv6dyLEs10MavB0LulLHJnh5K075yNnFKE6X9esEm9r2/HYdRJzPlKoFKNgkX0PyTWIvrYa3k5W9v2MQ0jr9A6RP2v6iyZiC6wdp5RX87judKvYAtnWwSYDB6fBGuy69HPO9OD+L7HGJLJ46xifzvJNVP6KdlfjO30Sj9t/3LUW/gKjmEPtnSHfpyhk3OXWE4++A6X7c9StwfdfC7fHg/Gv/4SMkJtDr2gacr6rTDuIjCwyETb9ojxVG3tUIxbXwcYlMNb3vpHNJnfPiKFoNDgpUmfujM45DdiYM9fv+kr/Has5BDIvZW7TKR5JCwbUn1hxdzyE6XL1WMpRzyLLX7fsUKDknen7A/RQP3Nd3SL3JQCyWG6y5NDrmxZPtDa9Tlmv/sxfUcUqTgaLdmE4dkvNpgkaDPIX+rR0sScJ9Tb/ml1dOcQxbOV7k3bcEhrZUZwft3cUjeC3nTa9YcIiYktbn5AIccvHDu5wDuo2o1u/eNoxJVH8LVPTnkc2vc8ntohUfw6YAzHNJf/cf5XgiHnHrR0/shlEPU87gEjl3nkOfKjWsfp2D7DTsrhvNx/9akcutaNYesFVMy9a3jEOHU3C0c9JdH7kW/d7iPa81lrWzhkMV9Ks12qO7Po0m6HRxyW/rU1Db0MtdeB7cerMdt+mPxFNavOmRUg562N/u1D/etqms+NjigP9fME2UIjpFwJWnmhmW4Hx4P3LkJ3XQ+ut1JeYwMUXm+OeP+3OqMzie1XXh/lLWS3gHcJ+sHLopAV65Yf6Tz/BixzvzlYhE8RjIH21ktl8aI3e+RCDaa9nTfqvQ7Y6QkZYl7C3p3p8cvi5IxMr0no8n28RhpjZPKYLwdI1vP9jbYoyujzl526R8jExt37r/8d4xUal5W+yk9TsRW8pc7m44TLnnJY4n7x8lvLZH6E3bjZM/8LA/2wXGiVffB7WLdOFm6hJ67+O04eVu9qV6XNkFa7d2vuqAxaslb56lOkKV/Gu3k0FGrpRPPN08QUy/1PwGpEyRqR2RBe/YEKWsbu8GVM0E+eQrnRuRPkLzkp5+dC/D5TfmJS79MkPeCyqcLhCfJverI8Z75k+TLATG1DPFJIrHAMLbafZL8nnGSUDg1SRZ7CZjZXZzE/eFILiN8kvh5/Co2qJwkOk/vV8ejigEGe7/8niSPctkH+9Hgv1ZdUVP4nOelyJXTk0T/bbfYIbEpkua4sykZNT13NHpCboq8klLl6Tk1RcQC42KivKZI5PCCW3mp+P8LrS3TQ1Pk1pmzocqsaaK+xLWFIjxNPodvsk1aM012yXpKb9gwTfxUOVcJ+mc4qDaWTBNixJdcj2ozf2c+Mp8mP2LNX7AspkmA22ji9NFpco+rXVPfcZqobRRs+HF/mhRcLPZi5U2TvTEaueUV0+TilmsvlV9PE5bwKbsENKd48IvGzDQxdL97z0dhhnxlJlJcDs2QWnuu70k+M2R96RX5BrSm7gs7NWeG+FsXizrfmyGORmNTfe9nyI+wghSBDzNEY+PHD++GZ4hQqnaD0OgMKV74rvgUmSXZ4Xquosmz5Bbr+IhR4ywZ2HBLrH/RHBk1GrzzzGiOaKxw3bpp6xwR2rP04OrQOXLujced/TVzJPTKt4F73+aI7lycucboHElWmejzZM8R+ldBk/Zlf8mDT7LRf1CHu4ENz9b9JeFrO9R2F/0lzYL6yzOK/xLf92kKAmJcME24o+pluGCNuf7cRRMuaBAWPiDuwAUdWbwD0fFcEGR/RfJMEnrVaOWbLC5okvob9g3NfRr4rPgJF0RP2Y/f4nDBAuPvrWd5uCG0Wfv5TXVu4B8tvy+5nhtmNYT5XE5xw7UKnc5vYdwQYK9idjeaG8r+9nTWXeGGGvqiFsl73JBYdYb16RU3REdI8oq2c4N+Z8W31dwUaK2JvFFOo8DR71UbXghSoFfigaetAQWyBi8yJ00o4HeBlGX4UEBSlcRMBlFAK6jyS1YcBb5Qjm8Ku0yBdSJLtp/OosBitSMhe/Ip0GBi3Zf9kQJWu7eG+vZR4HnJ6/KNf7Ec50NXlsnwAI0hVrBEjwduGDn7mxzhwX45ef2bIw+c+dlRJ3iaBxaU+ldyRfDAeT5+tZdXecDC/lpM4mse4N2/76LQGx4wJK1mzt08UGhTyfYc5YFrU5f4WbJUuMTeu4+6hAoRjbEf/HZQwVvFwHDYjAp+c8qTizypkLFUSCkkjwqPmkFPr54Ka69yddTOUcFm1YElJWtpcL6cV+/vVhrYuCx7Wx5Hg1bODt/3l2mQrhT83LydBm9KKN9Eu2ng5MXk5I/TYDwtQOLBIl5wUUxxzDPgheO8/VNZe3lBvdq9tjyFFyZ5Qu60p/JC9rwG8WM/eMHzw9Ob01N4X1OdvugiPvBd13D6tDofvJM9VOlgzgcRXiePae3jA+pbiTlFHz64szJrKN6XD8bzlecxM/jg5arL1Wn3+ECgqqLLt5YPLhs+Whz+hQ++lPdeEOjmA/Wv3HObRPlB+LmJjwU6L+bEwjdr+SFo5cXrZ7X5wXHkN1VKnx+Kn57K3n6EH/h9tZT0vPmBsTI7Qj2SHz768T19EcUPPFYHTXoe8oO4uErl4nJ+aJXawKx6zg89y+glte/5YeheznhkGz+kf9X5WMIlAH1//KWH0MbelWk5igKgZBjnvllPAMI2UNQGDQXgSrwH7ykHAXh/8/0KcV8BuHFioGclahZ7NOzeBQGo0XWg12cLwJugScN3jwWAspMWJdAgAMMu4tv1hgRgv4lWWv+YALAYdg79dEEISdX5cJIlCCuznpg/WyYIc9+06uirBCFU183HXFcQMkx4vgXYCUK0Y5rL3qOCIJcucm/1WUEoLPRwdU8WhLhjCzZ5puBzbF+Z4gxBABu+XKWHgqAp8VSPq1wQep8wRw3fCIJjsoXZ6wZBUC/bUHa1QxC2pKzXS+gVBOZo2brX/YJwT2WxSA0XHV65HZ1qo9DhrtVlnw4hOiSF6Jf/QMdUrwdOitOhV3QvvFhEB4M0jlujFB2qtIfJK1U6XD7odZihQwedLbxb56Me1NI5ewM66HH3JM83oQNT8nXOCQs6VJspKxSihuaelues6PB75sSD8+ieCSJQvo8Ok/Swff4n6ADK345bXKDDxrjTfI0hdDD1Txw8Gk6HfZvWXr4fR4czOmWVY1fp0Fm59dnCeDpc+rvozJLrdDimntE+fI8OZbGeL9tzsZ0HvgxXldKh56No8uhzOuQ5JPVFv8PytgUzWzvocNgsj8ygofd+fer7QYdpw7LL3X/o8CCgtYY2QIevHYdqNo1huYkSlHReBux6vrH3n3GXj7wI5mPA2iYJTSkJBjx3i4oNUGSAQ46Nj+9qBhhWv0iU0cH7F1+ZcdHD++xumN1Ed985VuZmxACZfVozjaYMqLfXj7Pey4DvhvTZINSqpfOCpC0DhBedzViMrucKtJc7xIBmmxrjw44MmG0b5nVC5Zxj7hcdw3rn9UxonML7VI/FpHgy4MPptMjAAAYEPb7y0y2QAf0q4heMLjHATGzFpyPxDOg6IhPci26NYT9ffZMBF/P9pzamMMDSJvFQbjYDzq2JbSzNYYD8/fdM43wGTCjs23Xuv+95//ue97/vef/7nve/73n/+573v+95//ue97/vef/7nve/73n/+573v+95//ue97/vef/7nve/73n/33/P+z9Jji1DAQAAAACAAACwVgAAEx4AAA==eJzt2/d/jff7B/Dr3OOcY4/YYm8RYhNba0ZRq2hRrb1KjKJKrCT2jNhBgtpbjfpYpVU1Y8QMjS1ij1J8X9d9vXPu+8fvH+CH5+N13e91v++RI3n09MxfbjoD/xzy0Mqdksx7wk1/HvbQsFUeX7KhjmR/lrKTdUedeZebbl8wrZqzG7IndFXawzfQWemktCnhoY7I6f+Y1Ba57qJJtUPdFD/QTcFYuylM3OG2RKI98xYP9YsT0zaLmfBRGQklNnnozEYPjUambBYZ1oteqAstEL9vEvm6mJQHye5BwF43zUTOh81Y56+5Hpq1x005cLwYsjuSNVIqbBRDNoguq0WJOiYdQHuL9eLIRjvZz8qs/m5Lpr12soroY7f2uH3JJqKegnXabpAcukEUTTQtLeZ7LLPXiDGYE/Kbm5pgDPtxnYhbL0ZHeWgqslM+05csCjWr7khWdZ2dbJQyda2d7OY6O9k1R7L7yg1HsgQlHmqDV6m8z+1L9gbPhx3BNZ7G/d6yRoQqVRbYyVZEY+5aW6EldrJp8+1k5daKbNdMGojzFDnothzFc01NVhfr/o7+Jb946D/kD4tFdbTvSjLpwX43zSlm0vZFHsvHVSLDPFFsjof8sP/nUWLYL+LfeWJJP5OKYz8n+HrQfr6uSZmjJSMwN9X7aPEnxvbFedL1MSl6sWQc1l+pHMM8tg597NQg01J8oPiAtdjhuXaymHliQrBJ8dhLp2qmxawqLjmSja1iUnOM/R6+w3H7apI9a5p0uLpJL2qIB8ojJQUeYs6/wXayR2AeNyxLypuW4Rg7Elb522J5HyVMSgqW+law6FNDlMK5T+P4LJyAc0qXGuI7JXM+8SqvScWwztPihpXtrhhUBOM7YMzXUPykQd+otOpw3VICY9hnuN6G0NiRjSCmuiiKNe+XwudfTTtzoL0M2t/FG7ShuMhU2qQzE3VKDBc7MY4Fom03snYtUcuRLF9V4VLaXjOoDNqrwWt4A4uHiCXwtpZ450g2Y4j4Y7CY87fhSzYJfew8znEI7VsTDepQ26S1kJr3Ak06iPZmwZIRGF8J1zr2tGFZhTGrgXNRafGolgjDM12F41IBooiyDAKKY/2SJhVAvfSsQfORMViPa84r5wxq+K0oh74f8Z6URXJdr7Sqy4hf+T4jT6l6J1RGf60y4oZyRamgFMGzYjsCTF+yY6VE4zLiKNZrhGwaII4q7cqIVqnQ1hp5AukKlPwbOpcR7ZXTAWJAcdGntJ3saAmRgNoP4vCOct0A95+9wBqsd2mxKa+drG1x0R7naHvGoAdIP+xnpnK+jEFtML6oblj1adzvEFzvjyelvoy1WXKAGIq27JiXS1lZSrTEmGFlJNlF7I1dDhAncQ72tKRoiTkDcLwZGYbxy7HHC7jvnKMq4J7j82H59wb9hJqllLOTnUV/FsheTlxBvbe8ZAj87UiWnKqcyIrrYEXOij1BpuX8ccOXzA97YzmCBJXFu4rrPgItsU40+gYHifo4/g3XwHr1EDG4RnYgSOyHbGkMOqiy/jlxnNfA2pzHlD+Vjl0Nyw9Ynw2CdHiWO3A/8peWmnMffj76lBFfKym49xeOGRTfDdeFfS3Bu3Qb78U4rBEC9f4y6HEJgwqN1Cl2g0FDcuO54ZwFcb7J6C+t6ow4F+tcCp/hyO9xL7tCdyUr2or2NCgDEM7ZrIedLMdQ3dK/l2F5hHH1YXkvyVPKtl7iiXKnp8jQ26CnuIZ7kTrlQnv6HlKf6oPPqd6Cuoun0zQ6hPbGaCvU085pUK+PeNtOtwyItJNFfGdYug3WqTnGRZUR9bHG/wrbktH3RV+8I72k3g05cFyws4t6Yp0fdakX9MfnGHBm7mdQtCPv9hOHUN9BHsX7fht7+Qr734I21guKXdOszIoxe/uLjUq3HwxqDt8oSwaIrrhXg0N0yjYQ5zfwbwYyapJOf4zRrTrTABE0ULyK0C09UY/DOrH9JSv1E79jX2wB5rCpg4Q/rvMk+ieo+gDW5/rGdc3SItSwjFJajdOpG3Id9sIWoJ4LCzFn1lSd7o7VrZpzJtbbPNigX7GPB+jj+hfMYRdQP8K957wII6foVj4YILZO0y1HhxiWgSG6L1nkUMOSFm2boSicgnyYw35CH+s4R6cMWHfnS53e4fgSxnDu7K9bMqMvzTCD8iiblJx4hsxAze5j3hNkfag+Xbcko+6Ods4xc3XqMcxNDVFvdolyGMPytcbPzjydTqNv/hw703ypU0nUpX5000bMfYwx61CzX2At3EL/gDwaRaJ/6Vy8UzgP6xelE+Hch4e7LYccybqNENUUb7RuWYL1WA6szR5inSj054aXeE/v/qzT9T6SvZHjcmu0OZdGV8d6rJqTXYM+mHMdmSXMQ+8ma9QQ64WNwzgYA6Mh8CdcG8ZthqJXNYrAOWpjf5HoO7zGZeUpXPfcUfg9HetwzkIbe4y1ZyJnwHS4gbWmq7pBhOZLVnW8x/JqnEarZmoUg76lsHGMh9bgvO+xVvGxbtqGc/yUQ6NY9EXP0qjReDtZE2gIEX01K3fgnGwv5rFH2TTallO8QP9reIxzvFM1Z3Fc61vkSty/HWF2sj9gepibtkZqtBd72AP7lNkThAfjZyH/DtYo63i35Q7m5caYRJWxhTWrDsV1NUH2xH3zw/ONGOO2ZBpvJ0sZpfmShU3E3+0FNUs46lRfh7stk1BHwoIBGgXjmhvjHKnJMvhpvmRh2AcLxXosYKLHlyz3BI0+4LxRqNn08XaySDzDxTguG+6hoHDJQFiAtpIOhfw1S0X0sQpQ3qE6hMw2rGQZxtvJxk3UaBeOU/0KfrjGVOdx39kBtB8Ml7yG86UqrxxOsJP9jr2xjeGi2UQ7WXiYqDdBxIzAev00as/rq2Tn0HcqwtYq3E42fbjmS9bvZztZAvpqTbA1n+D2JfsT/d9E2HZPtJP1xD1ITXYGfW3gjJLoSHbWkay5I9npSDvZP2EiBHtJTdb0Z9EHzy41WU28F6nJBmMNdifCTnYbnuDzMdUrXAPrHG4n257X60u2zZGsgSPZ4kGaL1kT1OxkHjvZzrwiPo+dLFF3+5K1zuv1JWvlSNYWYy85XHQka+dIdgHOcOb2UnOV7LIj2T+QBDcgEZoqbZW/colvlM6OZOdzeulr5MWcUl9UesA56K3kcSQbnAPjc0gOhUwqRyjZog1LOvwbkJrMD30sqyNZaBYvLcmOtXPYyWZlt5MtUgxlRjYvTVcm+nkpXzYxFfUEmAz5cTwFqWf1UiAyCoqonKHMUxYqi6AN9vSd4aa8mPc8o5f2Z8b1zjBoF477oC6E68qDOhhjK8/B77Kop2LOf+jLh3eVHUEbO6VEZxG9se5e/DtvoGaLlnppUGU8y2s6dZ6Nf9Oreqn7B5PWoe0p2jhzx+BeVfHSV1H4WcC7tgZzNqL9J9QvcS7Okph3A4LITTlxvB5jJuM4P/bKmQxVUO/FWguBc5TLTRmRLKWyeIX6NbxTnjuSDcLaL5DGMi9txn5cy6R2Ke/Rp3HbdZ08yHMYzy44kqVRfoA+0E95pYlkl1iE+rYm+QZeKivhCPqfObIqlIcpKtt8NKkUsgn6Of2QL0jSb7phKYh2Vm2KQblw/Bb9z6E/xiybaVDfaYYvm2G9FmjnDEByzTkE78AhPDNWAHNvTsXfglNEjVn4uwBjcqM94b1J85CJ6Occp0yAmi5xjsSsSPweiayLtkdYl3lxzIJdYjP2kao2js+gLw1J9oK4wRoFuSQfY355nofkOlipqNo4Wf73pi9ZT3L7khVQNSdr8UG8+mgnywVN/zPpMMY2+yB5GONTFUPbxfe2nB/EnPd2siWOZNnQ1x5rTcbak3EOzqgPtibvRYaPdrKl/9rJTr0TYZif8tHmP9lF0W9MWvxehAVqlocfTV+yPm9FxEc72fEeLksanIP1KqfR8CDNSjY7ymVJW14jL6RVPMrQIDuZXt5O9l05O9mgIDutGnsMVclKVtDoFvobV8TeVbIQR7LuzVy+ZO+CNF+y95Cvgm1GtMtyEDUbXMlOFupIdqiCnaxWRTvZaNz31GQ1Ue+uYMs91+VLtkvVnCyHI1lPnKu3Qw9HsswYU7WVy0r2J55J1bL4G2Cei3ZAugUuq8403UUvMX4bxpjQpZLkcDwDN/J7HHeFShfIqrdg3ztgG+yC7Sr3KoHK00riiSNZiiPZY0ey9sq7L1yWWz3Eb61ddBl/wxSZL3mnkrjtSJbkSHbLkWwM7in7XxXNst+RLD2kUckyQQaVoyvbxsCwNi5qUlWz7K0iGlYVlyoJF95Jf9yPDLjflytJdsK9/amyGOVINgKGQg