[Insight-users] watershed -- Repeated with different results (history dependent?)

Atwood, Robert C r.atwood at imperial.ac.uk
Mon, 19 Apr 2004 20:08:11 +0100


This is a multi-part message in MIME format.

------_=_NextPart_001_01C42641.A85EDC4C
Content-Type: text/plain;
	charset="us-ascii"
Content-Transfer-Encoding: quoted-printable


Here is the output of=20

cvs status itkWatershed* > wshd_status.txt

executed in the Code/Algorithms directory, the files seem to be up to
date. Are there other files I should examine? I'm sure it has been less
than seven weeks since I last built the package (Apr. 7 :-) =20


Thanks,
Robert











-----Original Message-----
From: Miller, James V (Research) [mailto:millerjv at crd.ge.com]=20
Sent: 19 April 2004 19:48
To: Atwood, Robert C; ITK Users (E-Mail)
Subject: RE: [Insight-users] watershed -- Repeated with different
results (history dependent?)


What versions of the Watershed code are you using?

There was a problem with the internal caching of results that I fixed 7
weeks ago.

Jim

-----Original Message-----
From: Atwood, Robert C [mailto:r.atwood at imperial.ac.uk]
Sent: Monday, April 19, 2004 1:06 PM
To: ITK Users (E-Mail)
Subject: [Insight-users] watershed -- Repeated with different results
(history dependent?)



Hi,

If I use the watershed in a loop, to get a series of segmentations ;
then choose one of the level values and run a single filter, I do not
seem to obtain the same results as I did within the loop for the same
parameters (level and threshold) (Also I get this behaviour if I do not
change the threshold at all i.e. it is related to the level setting)


I am guessing that this has something to do with the caching of the
initial watershed regions as described in the user guide p359:


"The relabeler object at the end of the
mini-pipeline uses the hierarchy and the initial segmentation to produce
an output image at any scale below the user-specified maximum. Data
objects are cached in the mini-pipeline so that changing watershed
depths only requires a (fast) relabeling of the basic segmentation."

In order to make efficient use this behaviour, I structured a loop so
that the "level" is progressively decreased, so the first update of the
watershed filter uses the highest value of the "level" desired.
However,runnign the routine again with a single value chosen from among
those tested does not result in the same output. I also tried updating
the watershed once with the highest 'level' value, then changing the
value to the target value. Any hints?

Thanks
Robert=20
_______________________________________________
Insight-users mailing list
Insight-users at itk.org http://www.itk.org/mailman/listinfo/insight-users

------_=_NextPart_001_01C42641.A85EDC4C
Content-Type: text/plain;
	name="wshd_status.txt"
Content-Transfer-Encoding: base64
Content-Description: wshd_status.txt
Content-Disposition: attachment;
	filename="wshd_status.txt"

PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PQ0KRmlsZTogaXRrV2F0ZXJzaGVkQm91bmRhcnkuaAlTdGF0dXM6IFVwLXRvLWRh
dGUNCg0KICAgV29ya2luZyByZXZpc2lvbjoJMS4xNQ0KICAgUmVwb3NpdG9yeSByZXZpc2lvbjoJ
MS4xNQkvY3Zzcm9vdC9JbnNpZ2h0L0luc2lnaHQvQ29kZS9BbGdvcml0aG1zL2l0a1dhdGVyc2hl
ZEJvdW5kYXJ5Lmgsdg0KICAgU3RpY2t5IFRhZzoJCShub25lKQ0KICAgU3RpY2t5IERhdGU6CQko
bm9uZSkNCiAgIFN0aWNreSBPcHRpb25zOgkobm9uZSkNCg0KPT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PQ0KRmlsZTogaXRr
V2F0ZXJzaGVkQm91bmRhcnlSZXNvbHZlci5oCVN0YXR1czogVXAtdG8tZGF0ZQ0KDQogICBXb3Jr
aW5nIHJldmlzaW9uOgkxLjkNCiAgIFJlcG9zaXRvcnkgcmV2aXNpb246CTEuOQkvY3Zzcm9vdC9J
bnNpZ2h0L0luc2lnaHQvQ29kZS9BbGdvcml0aG1zL2l0a1dhdGVyc2hlZEJvdW5kYXJ5UmVzb2x2
ZXIuaCx2DQogICBTdGlja3kgVGFnOgkJKG5vbmUpDQogICBTdGlja3kgRGF0ZToJCShub25lKQ0K
ICAgU3RpY2t5IE9wdGlvbnM6CShub25lKQ0KDQo9PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09DQpGaWxlOiBpdGtXYXRlcnNo
ZWRCb3VuZGFyeVJlc29sdmVyLnR4eAlTdGF0dXM6IFVwLXRvLWRhdGUNCg0KICAgV29ya2luZyBy
ZXZpc2lvbjoJMS44DQogICBSZXBvc2l0b3J5IHJldmlzaW9uOgkxLjgJL2N2c3Jvb3QvSW5zaWdo
dC9JbnNpZ2h0L0NvZGUvQWxnb3JpdGhtcy9pdGtXYXRlcnNoZWRCb3VuZGFyeVJlc29sdmVyLnR4
eCx2DQogICBTdGlja3kgVGFnOgkJKG5vbmUpDQogICBTdGlja3kgRGF0ZToJCShub25lKQ0KICAg
U3RpY2t5IE9wdGlvbnM6CShub25lKQ0KDQo9PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09DQpGaWxlOiBpdGtXYXRlcnNoZWRC
b3VuZGFyeS50eHgJU3RhdHVzOiBVcC10by1kYXRlDQoNCiAgIFdvcmtpbmcgcmV2aXNpb246CTEu
MTANCiAgIFJlcG9zaXRvcnkgcmV2aXNpb246CTEuMTAJL2N2c3Jvb3QvSW5zaWdodC9JbnNpZ2h0
L0NvZGUvQWxnb3JpdGhtcy9pdGtXYXRlcnNoZWRCb3VuZGFyeS50eHgsdg0KICAgU3RpY2t5IFRh
ZzoJCShub25lKQ0KICAgU3RpY2t5IERhdGU6CQkobm9uZSkNCiAgIFN0aWNreSBPcHRpb25zOgko
bm9uZSkNCg0KPT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PQ0KRmlsZTogaXRrV2F0ZXJzaGVkRXF1aXZhbGVuY2VSZWxhYmVs
ZXIuaAlTdGF0dXM6IFVwLXRvLWRhdGUNCg0KICAgV29ya2luZyByZXZpc2lvbjoJMS44DQogICBS
ZXBvc2l0b3J5IHJldmlzaW9uOgkxLjgJL2N2c3Jvb3QvSW5zaWdodC9JbnNpZ2h0L0NvZGUvQWxn
b3JpdGhtcy9pdGtXYXRlcnNoZWRFcXVpdmFsZW5jZVJlbGFiZWxlci5oLHYNCiAgIFN0aWNreSBU
YWc6CQkobm9uZSkNCiAgIFN0aWNreSBEYXRlOgkJKG5vbmUpDQogICBTdGlja3kgT3B0aW9uczoJ
KG5vbmUpDQoNCj09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT0NCkZpbGU6IGl0a1dhdGVyc2hlZEVxdWl2YWxlbmNlUmVsYWJl
bGVyLnR4eAlTdGF0dXM6IFVwLXRvLWRhdGUNCg0KICAgV29ya2luZyByZXZpc2lvbjoJMS44DQog
ICBSZXBvc2l0b3J5IHJldmlzaW9uOgkxLjgJL2N2c3Jvb3QvSW5zaWdodC9JbnNpZ2h0L0NvZGUv
QWxnb3JpdGhtcy9pdGtXYXRlcnNoZWRFcXVpdmFsZW5jZVJlbGFiZWxlci50eHgsdg0KICAgU3Rp
Y2t5IFRhZzoJCShub25lKQ0KICAgU3RpY2t5IERhdGU6CQkobm9uZSkNCiAgIFN0aWNreSBPcHRp
b25zOgkobm9uZSkNCg0KPT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PQ0KRmlsZTogaXRrV2F0ZXJzaGVkSW1hZ2VGaWx0ZXIu
aAlTdGF0dXM6IFVwLXRvLWRhdGUNCg0KICAgV29ya2luZyByZXZpc2lvbjoJMS40MQ0KICAgUmVw
b3NpdG9yeSByZXZpc2lvbjoJMS40MQkvY3Zzcm9vdC9JbnNpZ2h0L0luc2lnaHQvQ29kZS9BbGdv
cml0aG1zL2l0a1dhdGVyc2hlZEltYWdlRmlsdGVyLmgsdg0KICAgU3RpY2t5IFRhZzoJCShub25l
KQ0KICAgU3RpY2t5IERhdGU6CQkobm9uZSkNCiAgIFN0aWNreSBPcHRpb25zOgkobm9uZSkNCg0K
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PQ0KRmlsZTogaXRrV2F0ZXJzaGVkSW1hZ2VGaWx0ZXIudHh4CVN0YXR1czogVXAt
dG8tZGF0ZQ0KDQogICBXb3JraW5nIHJldmlzaW9uOgkxLjM2DQogICBSZXBvc2l0b3J5IHJldmlz
aW9uOgkxLjM2CS9jdnNyb290L0luc2lnaHQvSW5zaWdodC9Db2RlL0FsZ29yaXRobXMvaXRrV2F0
ZXJzaGVkSW1hZ2VGaWx0ZXIudHh4LHYNCiAgIFN0aWNreSBUYWc6CQkobm9uZSkNCiAgIFN0aWNr
eSBEYXRlOgkJKG5vbmUpDQogICBTdGlja3kgT3B0aW9uczoJKG5vbmUpDQoNCj09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT0N
CkZpbGU6IGl0a1dhdGVyc2hlZE1pbmlQaXBlbGluZVByb2dyZXNzQ29tbWFuZC5jeHgJU3RhdHVz
OiBVcC10by1kYXRlDQoNCiAgIFdvcmtpbmcgcmV2aXNpb246CTEuMw0KICAgUmVwb3NpdG9yeSBy
ZXZpc2lvbjoJMS4zCS9jdnNyb290L0luc2lnaHQvSW5zaWdodC9Db2RlL0FsZ29yaXRobXMvaXRr
V2F0ZXJzaGVkTWluaVBpcGVsaW5lUHJvZ3Jlc3NDb21tYW5kLmN4eCx2DQogICBTdGlja3kgVGFn
OgkJKG5vbmUpDQogICBTdGlja3kgRGF0ZToJCShub25lKQ0KICAgU3RpY2t5IE9wdGlvbnM6CShu
b25lKQ0KDQo9PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09DQpGaWxlOiBpdGtXYXRlcnNoZWRNaW5pUGlwZWxpbmVQcm9ncmVz
c0NvbW1hbmQuaAlTdGF0dXM6IFVwLXRvLWRhdGUNCg0KICAgV29ya2luZyByZXZpc2lvbjoJMS4z
DQogICBSZXBvc2l0b3J5IHJldmlzaW9uOgkxLjMJL2N2c3Jvb3QvSW5zaWdodC9JbnNpZ2h0L0Nv
ZGUvQWxnb3JpdGhtcy9pdGtXYXRlcnNoZWRNaW5pUGlwZWxpbmVQcm9ncmVzc0NvbW1hbmQuaCx2
DQogICBTdGlja3kgVGFnOgkJKG5vbmUpDQogICBTdGlja3kgRGF0ZToJCShub25lKQ0KICAgU3Rp
Y2t5IE9wdGlvbnM6CShub25lKQ0KDQo9PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09DQpGaWxlOiBpdGtXYXRlcnNoZWRSZWxh
YmVsZXIuaAlTdGF0dXM6IFVwLXRvLWRhdGUNCg0KICAgV29ya2luZyByZXZpc2lvbjoJMS45DQog
ICBSZXBvc2l0b3J5IHJldmlzaW9uOgkxLjkJL2N2c3Jvb3QvSW5zaWdodC9JbnNpZ2h0L0NvZGUv
QWxnb3JpdGhtcy9pdGtXYXRlcnNoZWRSZWxhYmVsZXIuaCx2DQogICBTdGlja3kgVGFnOgkJKG5v
bmUpDQogICBTdGlja3kgRGF0ZToJCShub25lKQ0KICAgU3RpY2t5IE9wdGlvbnM6CShub25lKQ0K
DQo9PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09DQpGaWxlOiBpdGtXYXRlcnNoZWRSZWxhYmVsZXIudHh4CVN0YXR1czogVXAt
dG8tZGF0ZQ0KDQogICBXb3JraW5nIHJldmlzaW9uOgkxLjExDQogICBSZXBvc2l0b3J5IHJldmlz
aW9uOgkxLjExCS9jdnNyb290L0luc2lnaHQvSW5zaWdodC9Db2RlL0FsZ29yaXRobXMvaXRrV2F0
ZXJzaGVkUmVsYWJlbGVyLnR4eCx2DQogICBTdGlja3kgVGFnOgkJKG5vbmUpDQogICBTdGlja3kg
RGF0ZToJCShub25lKQ0KICAgU3RpY2t5IE9wdGlvbnM6CShub25lKQ0KDQo9PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09DQpG
aWxlOiBpdGtXYXRlcnNoZWRTZWdtZW50ZXIuaAlTdGF0dXM6IFVwLXRvLWRhdGUNCg0KICAgV29y
a2luZyByZXZpc2lvbjoJMS4xMw0KICAgUmVwb3NpdG9yeSByZXZpc2lvbjoJMS4xMwkvY3Zzcm9v
dC9JbnNpZ2h0L0luc2lnaHQvQ29kZS9BbGdvcml0aG1zL2l0a1dhdGVyc2hlZFNlZ21lbnRlci5o
LHYNCiAgIFN0aWNreSBUYWc6CQkobm9uZSkNCiAgIFN0aWNreSBEYXRlOgkJKG5vbmUpDQogICBT
dGlja3kgT3B0aW9uczoJKG5vbmUpDQoNCj09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT0NCkZpbGU6IGl0a1dhdGVyc2hlZFNl
Z21lbnRlci50eHgJU3RhdHVzOiBVcC10by1kYXRlDQoNCiAgIFdvcmtpbmcgcmV2aXNpb246CTEu
MjkNCiAgIFJlcG9zaXRvcnkgcmV2aXNpb246CTEuMjkJL2N2c3Jvb3QvSW5zaWdodC9JbnNpZ2h0
L0NvZGUvQWxnb3JpdGhtcy9pdGtXYXRlcnNoZWRTZWdtZW50ZXIudHh4LHYNCiAgIFN0aWNreSBU
YWc6CQkobm9uZSkNCiAgIFN0aWNreSBEYXRlOgkJKG5vbmUpDQogICBTdGlja3kgT3B0aW9uczoJ
KG5vbmUpDQoNCj09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT0NCkZpbGU6IGl0a1dhdGVyc2hlZFNlZ21lbnRUYWJsZS5oCVN0
YXR1czogVXAtdG8tZGF0ZQ0KDQogICBXb3JraW5nIHJldmlzaW9uOgkxLjkNCiAgIFJlcG9zaXRv
cnkgcmV2aXNpb246CTEuOQkvY3Zzcm9vdC9JbnNpZ2h0L0luc2lnaHQvQ29kZS9BbGdvcml0aG1z
L2l0a1dhdGVyc2hlZFNlZ21lbnRUYWJsZS5oLHYNCiAgIFN0aWNreSBUYWc6CQkobm9uZSkNCiAg
IFN0aWNreSBEYXRlOgkJKG5vbmUpDQogICBTdGlja3kgT3B0aW9uczoJKG5vbmUpDQoNCj09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PT0NCkZpbGU6IGl0a1dhdGVyc2hlZFNlZ21lbnRUYWJsZS50eHgJU3RhdHVzOiBVcC10by1k
YXRlDQoNCiAgIFdvcmtpbmcgcmV2aXNpb246CTEuNw0KICAgUmVwb3NpdG9yeSByZXZpc2lvbjoJ
MS43CS9jdnNyb290L0luc2lnaHQvSW5zaWdodC9Db2RlL0FsZ29yaXRobXMvaXRrV2F0ZXJzaGVk
U2VnbWVudFRhYmxlLnR4eCx2DQogICBTdGlja3kgVGFnOgkJKG5vbmUpDQogICBTdGlja3kgRGF0
ZToJCShub25lKQ0KICAgU3RpY2t5IE9wdGlvbnM6CShub25lKQ0KDQo9PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09DQpGaWxl
OiBpdGtXYXRlcnNoZWRTZWdtZW50VHJlZUdlbmVyYXRvci5oCVN0YXR1czogVXAtdG8tZGF0ZQ0K
DQogICBXb3JraW5nIHJldmlzaW9uOgkxLjEyDQogICBSZXBvc2l0b3J5IHJldmlzaW9uOgkxLjEy
CS9jdnNyb290L0luc2lnaHQvSW5zaWdodC9Db2RlL0FsZ29yaXRobXMvaXRrV2F0ZXJzaGVkU2Vn
bWVudFRyZWVHZW5lcmF0b3IuaCx2DQogICBTdGlja3kgVGFnOgkJKG5vbmUpDQogICBTdGlja3kg
RGF0ZToJCShub25lKQ0KICAgU3RpY2t5IE9wdGlvbnM6CShub25lKQ0KDQo9PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09DQpG
aWxlOiBpdGtXYXRlcnNoZWRTZWdtZW50VHJlZUdlbmVyYXRvci50eHgJU3RhdHVzOiBVcC10by1k
YXRlDQoNCiAgIFdvcmtpbmcgcmV2aXNpb246CTEuMTcNCiAgIFJlcG9zaXRvcnkgcmV2aXNpb246
CTEuMTcJL2N2c3Jvb3QvSW5zaWdodC9JbnNpZ2h0L0NvZGUvQWxnb3JpdGhtcy9pdGtXYXRlcnNo
ZWRTZWdtZW50VHJlZUdlbmVyYXRvci50eHgsdg0KICAgU3RpY2t5IFRhZzoJCShub25lKQ0KICAg
U3RpY2t5IERhdGU6CQkobm9uZSkNCiAgIFN0aWNreSBPcHRpb25zOgkobm9uZSkNCg0KPT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09
PT09PQ0KRmlsZTogaXRrV2F0ZXJzaGVkU2VnbWVudFRyZWUuaAlTdGF0dXM6IFVwLXRvLWRhdGUN
Cg0KICAgV29ya2luZyByZXZpc2lvbjoJMS4xMA0KICAgUmVwb3NpdG9yeSByZXZpc2lvbjoJMS4x
MAkvY3Zzcm9vdC9JbnNpZ2h0L0luc2lnaHQvQ29kZS9BbGdvcml0aG1zL2l0a1dhdGVyc2hlZFNl
Z21lbnRUcmVlLmgsdg0KICAgU3RpY2t5IFRhZzoJCShub25lKQ0KICAgU3RpY2t5IERhdGU6CQko
bm9uZSkNCiAgIFN0aWNreSBPcHRpb25zOgkobm9uZSkNCg0KPT09PT09PT09PT09PT09PT09PT09
PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PQ0KRmlsZTogaXRr
V2F0ZXJzaGVkU2VnbWVudFRyZWUudHh4CVN0YXR1czogVXAtdG8tZGF0ZQ0KDQogICBXb3JraW5n
IHJldmlzaW9uOgkxLjcNCiAgIFJlcG9zaXRvcnkgcmV2aXNpb246CTEuNwkvY3Zzcm9vdC9JbnNp
Z2h0L0luc2lnaHQvQ29kZS9BbGdvcml0aG1zL2l0a1dhdGVyc2hlZFNlZ21lbnRUcmVlLnR4eCx2
DQogICBTdGlja3kgVGFnOgkJKG5vbmUpDQogICBTdGlja3kgRGF0ZToJCShub25lKQ0KICAgU3Rp
Y2t5IE9wdGlvbnM6CShub25lKQ0KDQo=

------_=_NextPart_001_01C42641.A85EDC4C--