[Insight-developers] Re: [Insight-users] Using itkDeformableMesh3DFilter

[Thomas Boettger]

This is a multi-part message in MIME format.
--------------090607090706090409090604
Content-Type: text/plain; charset=us-ascii; format=flowed
Content-Transfer-Encoding: 7bit

Hi Pingkun, hello everybody,

the mesh source is in fact not producing a very regular mesh. I faced 
the same problems with the cell density around the poles. I can offer 
you a itkRegularSphereMeshSource.

If it was of general interest I could contribute this class to the 
toolkit. Right now it is not documented very well and the cells are 
produced in an inefficient way.

I am sending it to you anyway. I hope this can solve at least one of 
your problems.


Regards,

Thomas


To the ITK DEVELOPERS: Tell me if there is an interest/need for the 
class and what needs to be changed in the code/documentation.



Yan Pingkun wrote:
> Hi all,
> 
> I'm using itkDeformableMesh3DFilter to segment 3D object from volumetric
> images. To test the algorithm, I prepare some synthetic volumetric
> images with very simple object in the center. I followed the steps in
> the example DeformableModel2.cxx.
> However, I found that the initial mesh is hardly to snap on the surface
> of the object. In fact, the result changes very slightly from the
> initialization. I have changed the parameters within a large range, but
> the situation doesn't improve.
> 
> Now I have several questions.
> 1) Which paper is the algorithm based on? Some code segments look
> strange as indicated by Thomas Boettger in the mailing list.
> 
> 2) Is there any better way to initialize the mesh? Because I found that
> the mesh produced by using the SphereSource has too many cells near the
> poles. This may hinder the mesh to deform because it will produce large
> internal force.
> 
> 3) Can any one with successful experience on the filter share his
> experience?
> 
> Thanks a lot!
> 
> Regards,
> Pingkun
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users
> 


-- 
Dipl.-Inform. Thomas Boettger
Deutsches Krebsforschungszentrum         (German Cancer Research Center)
Div. Medical and Biological Informatics B010    Tel: (+49) 6221-42 2328
Im Neuenheimer Feld 280                          Fax: (+49) 6221-42 2345
D-69120 Heidelberg                            e-mail: t.boettger at dkfz.de
Germany                      http://www.dkfz.de/mbi/people/thomasb.shtml


--------------090607090706090409090604
Content-Type: application/x-zip-compressed;
 name="code_RegularSphereMeshSource.zip"
Content-Transfer-Encoding: base64
Content-Disposition: inline;
 filename="code_RegularSphereMeshSource.zip"

UEsDBBQAAAAIAFtrejBNfuv57gEAAA0FAAAiAAAAaXRrUmVndWxhclNwaGVyZU1lc2hTb3Vy
Y2VUZXN0LmNwcH2U32vbQAzH3wv9H4QD49JmJaVlHY6bl64be0g74qwvIZSrLSfHzj846ygh
9H/f3dlxm9guxhikj6SvpDsPRBJjAs+z8O756X5+ejIoFF+nHF65ykS2BgaxKPmLRPDh+ub7
NxgaBrNYJKcn9hmILJI6RvAE/ZthubnYeEfWOa615CosNqjQImGuVYRt8AcmXEsKiZOILLhQ
XFDpwA9oIPKSFPJ0as0iI+gvssCSmEFGEG24OoPlyjaws4EAtC3Qdm/Cff9PbrAgkTmn0dUU
IHGGhUEmLdbmr1EwrJPaCfaoCvYRtk5L8j5Tp6cWigog3R4DcNsX9ICvbFhlbfqCCG2iSf29
+CmkZONjyvefuNQuDZQRl/g7E7S8Wplau8uRed6O0zrKuBvaAfZtK/46DZHuXH1Wyajq95Bz
LHOpSeQZu/wUDG1t5hRYro+c5eYgC4xZf7K/RcwJ2T5LSbHvR7kmCALw2gE+eNbTdjRD2O/+
fZPp1trMxLrq/0J61FRoaiT0HMVpndDtAODw+H7wFNTZir0p9sJ/dm7BM6QLM38AWe09V8xd
wdvxBERQteJUP+j0BdVjdTBKox7E+fkQdtWgG8y5mRjBl4KqLRzqcv6lG6pwhlU94sL5D9W8
dfcF91ncoV0haZXB2M3DhP4HUEsDBBQAAAAIACJrejDob3+J2QUAAOgPAAAcAAAAaXRrUmVn
dWxhclNwaGVyZU1lc2hTb3VyY2UuaK1XbW/aSBD+jsR/GKVVD1IO0qhXnQitxJH2gi5vwm6r
SpWixR7wXm2vtbuGoKr3229mjbGBkLanIx+y9s48Oy/PzI57x6//r1+z0WwA3Go11yLpA8A4
NXIeWfBwnmBqhZUqhWcwwbk0VhePvlLxF2lZ80qFeYysCE8nI28m+UHaLySfx0J7WYQar9BE
nsp1gF17f99ZwFNWvRTpPBdzVh49f85vzoUtoAisWJ+enLzsnZz2XryCk9/7p6/6L08K5Q+o
DZnSd+fiQhYPL7pumwVGKltp50kraFdeqZldCo20nRqlrcyTLgzjGJyoAY0G9QLDLkN4iDD2
/9ogkfEWlIbI2qzf6y2Xyy552lV63rvwry57lWBkE5iRZIhWyNh0C5Po50fSgCmNoHXIUZXT
3GIIH8f+xc17H4bXn+DjcDIZXvufzmApbaRyC7jAFNYwNiLlJIslaRGSFqldgZrB1dvJ6ILU
hn+ML8f+Jzb23di/fut58O5mAkO4HU788ej95XBSQt2+n9zeeG+7hbsMLKZqgRBs4pcqKwM0
zqFEsdkpLRNHBefZ91n2g7/jXrPxRM7SEGdwd3eYRXcRyZGQTPG7coyYBnEeIhyRJG92o6MH
3q4Jynu93tbuBwys0vtavpZE4RhHGMcPYIqMSEb5T7HQbTZSkaDJRIBcIc3GV37XOz6Gz0Es
jIEDXlCiSGSqJUXFrcdplhNXmUGcr4DKFDWnn5/ohKjjViYQMUJLi1DmhpIGKIKICEdlzdXS
rms4XJGGILgEVJy7Os8ENQUk8A4sI0nKgdK0nak0NGCV09YY5K4UC9tCzGzE0sRNCDQSS9I5
oRrnlCSTwK7D5g7mQjt2mucq/cVCbuqQDs0ABWVOHtpIpPBbB6YYCJZjfZgjxbdoS3OkqCxQ
r8DEatmFDfRnkVMNaao+lQgDfyi0lgDLfWKdRaom6jgwKFLh3+SWgsxJeNNsPJoe6EOWT2MZ
QPVusK1PiS5E+lx0nHHPUrCFDuHIw3h2BHaVIRG664yB8vHgkeWPlc/qCsSrfv+AHSSeZ6id
M/tKXiK0vVWSuTRgXJJfP35HOKBWamGtwn3VVnqlv1dICQhdCylIQemykVb5PHIsUtO/qchg
JrjSVmUY6KxrXF6JQKsW47fPACrMSZ7+aonNzrR6U4IWM1kjJzQklvDRpl0D9UmhQD0Q4E4t
l+2aGxcqDoFNL7i/d7DJMJAz7srTVSFR8mpTSmY3x7UUQbVkE7fizv+5f+zIFNn2tZDEfYCb
6unHtF2qeMl0utk8HVCu0ElT3mNcajrl8sVWljyV8GWS0u0lqVGVoHtheNxCqv7q/UO0PKjO
zXkdnmr5E7ExmyZeGvLw6/8ECbCL5mDK2J0XV5wjkmuazFXq9zSv8Alm3ZWZf0XHz/WMb5ed
0LqSZefHbChLDKpkdWphofqFLbk9Kjio+r032JN/s+lOJOewD/KpJkA9hQp8mFtVhnkLo7bx
A1A7GmuofZwHkYwN+/1MSD3IeWhMqZhjlc47W09v6BIO8f6WxA67VzTj2hgw2FLqwA6kM4zk
DyNWBdrvF1MJr6FaOrV6+fFsS9QwtQ4Ehq4/vpVneRq4tsWTxPoO44altNAyXnGXS8QXxz8H
e0y71OlRGEkom27qoS2baTk51Fwje8sOWpcduaGlU/OofbYr4/EA04EdGbpLtbLkMYbuOj3Q
wlsF3j8Ht+HrNxZYKBnS9w8dwARsuewrmslRJM9AmU6RRU4cNS/p/rXBXXlrrxzAn8UUgvTV
Ilqlv5XdcC4XMsRW7U32ou48ZKftOt4wDJmurccbCWWhGPe2iER5ynh3HJrdHRne16+zRIUY
l/Njkc+aScldkaQz2NEw9RmxUKunmxQrJuyd5qZS89BpLt3r/sdJlgv+BHwsxdXk8YxGg14v
y3WmDBJ16ZvFfSAhf8NiuImryjhNSr/+Gd1vzgn++0aCgDRc7EzyO6P/Ix+/7jvgCUHIGWE9
9u1y0Wz8C1BLAwQUAAAACAAYa3owTofqIDgIAAA9IgAAHgAAAGl0a1JlZ3VsYXJTcGhlcmVN
ZXNoU291cmNlLnR4eM1aa2/aSBf+Hin/4exFkQkEbNN2VyaJRNN0izY3Ae02Xa0iFw8wb43N
a48h2Sr/fc+Z8dgGTEIaql0iLj7znMuceWZ8ZpzG/tG2Xrs7uzsAV1E4ityJAwCdIOajsYAe
G01YIFzBwwD2oMtGPBaRuuyHof+FC9I8D73EZ6QIP3dPekNOF1x8QXziu1FvOmYRO2fxuBcm
0YDVxe1tbQY/k+qZG4wSd0TKJ9UqSd64QplCY+q3bZovGqbdsF6B+atjv3JemEr5A4tiDMWR
ftmMqwurblEzAU7C6V0ke2IMKnmvwqGYuxHD5iAOI8GTSR3avg8SGkPEYhbNmFcnEz3GoNP/
PbOEwQsIIxgLMXUajfl8Xsee1sNo1HjXPz9r5MCxmMAQkR4TLvfjugoJX/0xjyHWQeBvj7LK
PyeCefBHp//u8n0f2hfX8Ee7221f9K9bMOdiHCYC2IwFkJoRY1SeTH2OWmgpcgNxB+EQzk+7
J+9Qrf26c9bpX1Owbzv9i9NeD95edqENV+1uv3Py/qzd1aau3nevLnunddVdMux+DmcMBln+
glDwAYtlhyYhhR3gz4mkguzZ4yzb8LXf2N35iQ8Djw3h5mYNhW6QQghDDA/YYzD6C9wJi6fu
gBEvd3e+kqyxj92Xb3QpGKYS6XY48N04hv5lIqaJIFPHuztrzB8uohxnDc6oSI8A0iPAPpxE
DH3JTIfShJI36EvcTRmFW4zBca5CHggWpXA4Wmy9YHOj0iJt7B6CoxAHK778/D82EI7TY+Ii
mXxm0eWwy/6f8Ih5Sjk2rIe1xFgBDbOWeq7/xkQai1FRypObE0bX9bfc9w1TC3sD12dKZmlZ
l8Whn8j14whsFN4/bSBmIfc2H47fWMBwtWK4jLjZGCS0CgQ4Z/wwGAHOuyn1Jv6z+ReG9NWs
WTX7vrUwDrnVPoqWh4LkqClwUh8cY3LShFVaAGrG57CDY8zpCfP9GFebcCAnzznDme0ZK06K
MOa9ucNIOCbUv6OG1/ITKKty5ZY9wNUMFxrssY4PJneqBcMrBqFHME4pQ+9LKSHXMLVaq3ni
3i1aMVtqNgGC/jQpYRaSOR1qElQzMuBVKwVaBDQLQGsBaGVAexloLwBtBVzKpwzcwPiq1Roa
0TnRAR785yM0Nw3Q+q8HeLB5hM8OUNE2C3HTMd4wPmu7CdxydAfbCK+xDz1+C3Q/oPoDqwdQ
KyH4dxxnPC7RdHtKAi4gliutXJ71Pffj1dn7HpgFwXnnAsDKBdcSYRcEEtHMBZ8k4kVBIBEv
9SKTL84yjdJgqyCVSfu0LJUp+qikZKXtebRaGnlSahm4lt6tNvJEPVxx9GlDR1apo2WT0s/1
c/zYpX5WkrTe0aaZa5Y7Ku3Qs4boxeZDJPmz6kjG9Kifl0/oUGniNnPzanM3H0vzdr2Zn1/0
KkQlu5HdzbENOK01Lfw6hGJJhpJqtSItf6WPB0ogWZ2sFBwTWZaUlBtSnFYbG1t1nI6g8g3D
H6RmUwcHx69xJxpkBh8xq8PDfV3iU+m8DMir54VCaGbVYGbju7nSto+NN1NRJrfXyJu5XOpi
HHszC6tEqSGvbLpq6ivtdYC7VLFUju2Ladq6KJdj3/Fuj8wa5Uz+ehB4ORzGTKyOmN4vyB7k
Q6fVUGHRQqmTgM0RsFpZK6w0fe5OFwdp7Aaez7xTb8RoxBdBeqCUgVh4jjPAwOHwEH6UHKnH
/G/cbjnwI8kyuvSklERSh6GLzMp8zH1mpBT7ocCx08BDnUoabTojkGZTDEs5++D6CdrNK+qO
l6cKllMqEYaY4oSvSQY8jrQQaW+EtBHZTJEaq/heqACUb6Rfbq8MYxHGfhhjE6aZecud8iEY
PxQH8eC4g3v529NbLDJiQ+Z/4n5hN1OXRzoi6bRSgb09bUa/nmTKUqZMNFXRhr7mFlM6mn8d
ae628sZlPzGLxKnP6PjrwaBr2USoFKytpEyDqAJ7w2fcYwYuL5jlXOte/2B+zErCp8x+Y2Ir
RTuFRCCTFy0ipTbpdLGr9/lPCnwrngoDWebpXheI9Go01AhQpakkT6Zg6s7eAgXtQtLXU1De
3Z9NwjzsZ5DQxnXjO5KwkNpyElrfRg376SR8qif7YbqXk7D5zSS0c84/l4RmIenrSWhvhYR5
2M8gYRNviN+RhIXUlpPQ/jZqrFufHiDhUz2ZD9P9Pqfg0oZC6rZgtVVOg2xdLkPYBYRdQOjt
hhrSwlYjLTXzCFNBtdpaG18ewZoQ5dQraSxGVwB8j+BK/evg7MeC+/dSl4f/3Uc3l8TZdcrO
YqnaCbjgri/L8NZqs14cYkMvF4Uz6jUKak9ZwKebzAfg8jnAkoZ6NqCU7tPnEVt4GuQ4+eZv
8wcW6Xq4cBhfWzibt+WCnq4sxYbaMO33BzbALbMUerVBxiIPRx3V4UCf7+MVrf8WVMHwYB/M
+kudugE1IFIfK6bioaJffpTZGNDToIswmmSjOlQMzI8xSzF2EWMvY9IgyNv+kfSq4yLjUmTl
IjsVZfNRfUZMJFEAhu4DdnNAxu8fHd6nPWPS0+bR04cJHc8sb+P1LTNebuE44QojTS7aiQi1
McFiyd60y/qy3sd7x+U8YFE85lODpjn0I05NkhCVJXw+82jnqmPRqElx7tDxcS1ThFKMnEsS
hxTmt0w6rTTrZja9tpz8qwgj7jF/aKiNfRiLiLmTPQgxoVQT0GGX+qqoc5T0IWAvmbJI+i8a
QS2F1iwMYzozSC3QIYPikj5e0NRaOlco0ZNMz9Tk1QZa+dFcplp4gLqkn2b3HmtRlNA/Aiw9
86a/n7CJD3d3/gFQSwECFAAUAAAACABba3owTX7r+e4BAAANBQAAIgAAAAAAAAABACAAtoEA
AAAAaXRrUmVndWxhclNwaGVyZU1lc2hTb3VyY2VUZXN0LmNwcFBLAQIUABQAAAAIACJrejDo
b3+J2QUAAOgPAAAcAAAAAAAAAAEAIAC2gS4CAABpdGtSZWd1bGFyU3BoZXJlTWVzaFNvdXJj
ZS5oUEsBAhQAFAAAAAgAGGt6ME6H6iA4CAAAPSIAAB4AAAAAAAAAAQAgALaBQQgAAGl0a1Jl
Z3VsYXJTcGhlcmVNZXNoU291cmNlLnR4eFBLBQYAAAAAAwADAOYAAAC1EAAAAAA=
--------------090607090706090409090604--