headers
arpitagrawal | 1 Jun 2010 09:51
Picon
Favicon

Difference of Gaussian Filter

Hey,
I have been working on segmentation part extensively with itk since a
month...I wanted to know how can we implement a difference of gaussian
filter("NOT" DifferenceOfGaussianGradientFilter) on a grayscale image.

Regards,
Arpit

_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.html

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users
Permalink | Reply | 2 comments |
headers
Sergio Vera | 1 Jun 2010 16:00
Gravatar

Re: Difference of Gaussian Filter

Hello Arpit

If you are referring to this method: http://en.wikipedia.org/wiki/Difference_of_Gaussians
you can try to simply use two RecursiveGaussianImageFilter with two different sigmas, and substract the smoothed images with DifferenceImageFilter.

Regards.


On Tue, Jun 1, 2010 at 9:51 AM, <arpitagrawal-dOeZK83GKW6f0DUV/oxz1A@public.gmane.org> wrote:
Hey,
I have been working on segmentation part extensively with itk since a
month...I wanted to know how can we implement a difference of gaussian
filter("NOT" DifferenceOfGaussianGradientFilter) on a grayscale image.

Regards,
Arpit

_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.html

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users



--
Sergio Vera

Alma IT Systems
C/ Vilana, 4B, 4º 1ª
08022 Barcelona
T. (+34) 932 380 592
www.alma3d.com
_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.html

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users
Permalink | Reply | 1 comment |
headers
Matt McCormick | 1 Jun 2010 17:47
Favicon
Gravatar

Re: Difference of Gaussian Filter

Sergio Vera <sergio.vera <at> ...> writes:

> 
> 
> Hello ArpitIf you are referring to this method: 
http://en.wikipedia.org/wiki/Difference_of_Gaussiansyou can try to simply use two

RecursiveGaussianImageFilter with two different sigmas, and substract the 
smoothed images with DifferenceImageFilter.Regards.

DifferenceImageFilter is specialized for making comparisons between two images to
 tell if they are the same -- SubtractImageFilter may be more appropriate.

There is a good example at the end of the ITK Software Guide on creating 
composite filters.

Regards,
Matt

_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.html

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users
Permalink | Reply | 0 comments |
headers
Charlotte Curtis | 1 Jun 2010 18:08
Picon

TransformPoint with non-invertible transform

Hi all,

I'm transforming an image using an elasticbodysplinekerneltransform,
followed by a bsplinedeformabletransform.  I've got a set of points
defined on the moving image that I would like to transform as well.
Looking at the code (and the behaviour when I try to GetInverse()), it
seems that the inverse is not implemented (or my transform is not
invertible).

How can I transform my point set so that it matches the transformed
moving image?  Thanks,

Charlotte
_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.html

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users
Permalink | Reply | 3 comments |
headers
Ken Urish | 1 Jun 2010 18:48
Picon

Re: registration question

> What I'm reading from your email is that the image
> origin and spacing got lost during the process of
> writing out the data to a DICOM file.
> Is that right ?

Yes. When I write the image using the series writer, the spacing and
origin for the z-axis were being thrown out. I found a way to fix the
problem although I do not know why it works. In my series writer (code
below), when I comment out:

//seriesWriter->SetMetaDataDictionaryArray(fixedReader->GetMetaDataDictionaryArray()
);

everything gets fixed and the origin and spacing are preserved. I
noticed this line of code is also in the registration/example7.cxx
file.
I am assuming that if it was in the example that it should be
important. Any thought on why it now is getting me in trouble? I am
running a cvs itk build from 04/08.

SeriesWriterType::Pointer seriesWriter = SeriesWriterType::New();
  fixedNamesGenerator->SetOutputDirectory(outputDirectory.c_str());
  seriesWriter->SetInput(caster->GetOutput() );
  seriesWriter->SetImageIO( gdcmIO );
  seriesWriter->SetFileNames( fixedNamesGenerator->GetOutputFileNames() );
  //seriesWriter->SetMetaDataDictionaryArray(fixedReader->GetMetaDataDictionaryArray()
);

I have one last question.

I have not worked with the image direction matrix before. I notice
that my output registered image has a new direction matrix as compared
to the fixed image:
1 0 0
0 1 0
0 0 1

The spacing and origin are otherwise identical and the images line up
perfectly when I put them through the checkerboard grid. I wanted to
make sure that its ok that the direction matrix is not the same and
that I am not putting my head in the sand on another big problem. If I
do use the resampler to make the output image have the same direction
matrix as the fixed image then it starts to screw w my spacing and
origin on my z-axis again.

I really appreciate your help.

Thanks
--Ken--

>
>
> You could make an additional verification by
> replacing the ImageSeriesWriter with an
> ImageFileWriter and saving the resampled moving
> image as a 3D single file, by using the MetaImage
> file format (.mhd), or Nrrd (.nrrd), or Nifti (nii).
>
> All of these file formats should preserve origin,
> spacing and orientation.   (DICOM should preserve
> it too... so.. once you run this verfication, we should
> go back an take a closer look at how you are writing
> out the DICOM series).
>
>
>      Regards,
>
>
>          Luis
>
>
> -------------------------------------------------------------------------
> On Fri, May 28, 2010 at 7:33 AM, Ken Urish <ken.urish@...> wrote:
>>
>> Luis
>>
>> Thanks for your response. Ive always been impressed by the excellent
>> support in the itk community.
>>
>> Indeed, you are right. My z-axis spacing and origin is getting lost.......
>> Im using dicom so I am surprised this is happening.
>>
>> When I look at the output of the registered and resampled output of
>> the moving image - the spacing and origin is identical to the input of
>> he fixed image. This is expected. I am doing a mattes registration so
>> my code is similar to the example itk code. I am using the verser
>> transform.
>>
>> When I take that identical output registered image and use it as input
>> for the checkerboard filter program the z-axis origin and spacing are
>> changed to 0 and 1 respectively.
>>
>> My registration output image info is collected from the resampled
>> image just before it is sent to the image series writer using gdcmIO.
>> The checkerboard input image info is collected just after the input
>> reader. Thus my registration output writer is somehow loosing this
>> z-axis info. Im keeping it in dicom so Im stuck on how Im loosing it.
>>
>> One solution would be to change the spacing and origin on my
>> checkerboard program. But I would be interested in having this fixed
>> on the initial registration program. Any thoughts on my error? I
>> appreciate the help. I included the specific output below.
>>
>> Thanks
>> --Ken--
>>
>>
>> Registration Fixed input image:
>>  Spacing: [0.3125, 0.3125, 3.48]
>>  Origin: [-18.1241, -60.8059, 58.0112]
>>  Direction:
>> -0.0348995 0 -0.999391
>> 0.999391 0 -0.0348995
>> 0 -1 0
>>
>> Registration Moving input image:
>>  Spacing: [0.364583, 0.364583, 0.7]
>>  Origin: [-3.39372, -70.4365, 70]
>>  Direction:
>> -0.0348995 0 -0.999391
>> 0.999391 0 -0.0348995
>> 0 -1 0
>>
>> Registration OUTPUT Image of the resampled moving image:
>>  Spacing: [0.3125, 0.3125, 3.48]
>>  Origin: [-18.1241, -60.8059, 58.0112]
>>  Direction:
>> 1 0 0
>> 0 1 0
>> 0 0 1
>>
>>
>> =====================================================================
>> Output of when the registered images are put straight through a checker
>> filter:
>>
>> Checker Input image of the registrated output moving image
>>  Spacing: [0.3125, 0.3125, 1]
>>  Origin: [-18.1241, -60.8059, 0]
>>  Direction:
>> 1 0 0
>> 0 1 0
>> 0 0 1
>>
>> Checker Input image of the fixed
>>  Spacing: [0.3125, 0.3125, 3.48]
>>  Origin: [-18.1241, -60.8059, 58.0112]
>>  Direction:
>> -0.0348995 0 -0.999391
>> 0.999391 0 -0.0348995
>> 0 -1 0
>>
>>
>> Checker Output with no resampling
>>  Spacing: [0.3125, 0.3125, 1]
>>  Origin: [-18.1241, -60.8059, 0]
>>  Direction:
>> 1 0 0
>> 0 1 0
>> 0 0 1
>
>
_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.html

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users
Permalink | Reply | 0 comments |
headers
Roger Bramon Feixas | 1 Jun 2010 21:40
Picon

Re: Performance regression in ITK + gdcm 2.0 reading DICOM JPEG files

Hi Mathieu,


Thanks for you time! I just tested it on Mac OS X (I will also test it on Windows soon) and there is a really nice performance improvement, about 2x. 

However, I have two code comments: First, there is an abort() at line 350 of gdcmBitmap.cxx which causes a crash, I commented it to be able to test it. The other problem is related to the use of gdcm::PixelFormat::UINT8 at line 342 of gdcmBitmap.cxx. If you're working with JPEG files and its pixel format is more than 8 bits, a downscale has to be done. Perhaps, it could be replaced for GetPixelFormat().

On the other hand, we want to include this commit in our GDCM compilation but I would like to know which is the best way:
- Using GDCM 2.0.14 patched. It's the last stable version but I don't know if patch it is a good idea because some changes have been done.
- Using the trunk version even though it isn't stable, or
- waiting for GDCM 2.1 if it will be released soon.

I will try to test it on Windows tomorrow and send you the time results.

Many thanks!

Roger


On Mon, May 31, 2010 at 11:21 PM, Mathieu Malaterre <mathieu.malaterre-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:
Hi Roger,

 Yes ! I finally found some time, Snakes on a Plane gave this opportunity ;)

 I think I made some progress with the JPEG DICOM file. JPEG 2000
files will be slightly more difficult. See log at:
http://sourceforge.net/tracker/index.php?func=detail&aid=3009661&group_id=137895&atid=739587

 For reference, this issue entered GDCM 2.x, when handling the
lossy/lossless duality of J2K and JPEG-LS in DICOM. I am required to
know after ExecuteInformation whether the encapsulated JPEG fragment
was compressed with lossy algorithm (this is not indicated in the
DICOM dataset).

 Could you please confirm this on your side (again JPEG 2000 will
take me some more time) ?

Thanks again for your report

On Fri, May 28, 2010 at 12:00 AM, Roger Bramon Feixas
<rogerbramon-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:
> Hi,
> Have you had time to think about the causes of the performance regression
> problem? I can work on it, but I would like to have, if you have it, some
> information which helps me to know where I need to focus the effort.
>
> Thanks,
> Roger
>
> On Tue, May 18, 2010 at 9:14 AM, Roger Bramon Feixas <rogerbramon-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org>
> wrote:
>>
>> Hi Mathieu,
>> Perhaps, the answer of question 3 wasn't so clear:
>> ITK versions:
>>   - CMAKE_BUILD_TYPE:RelWithDebInfo
>>   - BUILD_SHARED_LIBS: OFF
>> GDCM: Release from  http://gdcm.sourceforge.net
>> Test program: Release.
>> Compiling the test program in Debug mode, it takes more time but we also
>> have large differences between GDCM 1.2.x and GDCM 2.0.
>>
>> Roger
>> On Mon, May 17, 2010 at 11:50 PM, Roger Bramon Feixas
>> <rogerbramon-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:
>>>
>>> Hi Mathieu,
>>> Thanks for your attention. I answer your questions:
>>> 1. What is your platform ?
>>> I used a Windows XP 32bits to do the test. However, this behaviour also
>>> occurs in Mac OS X platform.
>>> 2. What is your compiler ?
>>> All ITK versions were compiled with Visual Studio 2008 and we downloaded
>>> GDCM 2.0.14 binary from http://gdcm.sourceforge.net
>>> 3. What are your compiler option (eg. I hope you do not run bench
>>> without CMAKE_BUILD_TYPE:Release) ?
>>> RelWithDebInfo
>>> 4. Where is the dataset used ?
>>> Dataset 5 is from OSIRIX: http://pubimage.hcuge.ch:8080/DATA/CALIX.zip
>>> -> CALIX/CT1 abdomen/D30MN BILISCOPIN
>>> We are really interested to improve the time needed to load jpeg lossless
>>> datasets. The jpeg lossless datasets used are non-anonymized datasets, I
>>> will try to anonymize them if you want them.
>>> 5. Where is the source code used ?
>>> I attach it in this mail.
>>>
>>> Thanks!
>>> Roger
>>> On Mon, May 17, 2010 at 10:44 PM, Mathieu Malaterre
>>> <mathieu.malaterre-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:
>>>>
>>>> Hi Roger,
>>>>
>>>>  This is a very interesting post !
>>>>
>>>> On Thu, May 13, 2010 at 11:00 PM, Roger Bramon Feixas
>>>> <rogerbramon <at> gmail.com> wrote:
>>>> > Hi,
>>>> > Recently we updated our ITK version from ITK 3.16 to ITK 3.18 and we
>>>> > decided
>>>> > to relink ITK with system GDCM 2.0. After some tests, we concluded ITK
>>>> > 3.18+GDCM 2.0 is faster than ITK 3.18 reading DICOM Little Endian
>>>> > files,
>>>> > however ITK 3.18+GDCM 2.0 is rather slower reading DICOM JPEG files. I
>>>> > attach a PDF file which is a time comparison of ITK 2.8, 3.16, 3.18
>>>> > and
>>>> > 3.18+GDCM2.0 versions.
>>>> > I don't know if it's just a GDCM problem or if it depends on how ITK
>>>> > uses
>>>> > GDCM.
>>>>
>>>>                       GDCM             Reading
>>>> UpdateOutput-        UpdateLargest-
>>>> Data       ITK version version          directory         Information
>>>>        PossibleRegion
>>>>
>>>>
>>>>       5           2,8            1,2           642
>>>>  6                  20153
>>>>                  3,16            1,2           547
>>>>  0                  25121
>>>>                  3,18            1,2           547
>>>>  0                  25124
>>>>                  3,18              2         25297
>>>> 187                  71437
>>>>
>>>>
>>>>
>>>> Data     Transfer syntax                Modality          Files
>>>>        Size (MB)
>>>>
>>>>       5 JPEG 2000                      CT
>>>> 243                     23
>>>>
>>>>
>>>>
>>>>  Could you please post a little bit more on :
>>>> 1. What is your platform ?
>>>> 2. What is your compiler ?
>>>> 3. What are your compiler option (eg. I hope you do not run bench
>>>> without CMAKE_BUILD_TYPE:Release) ?
>>>> 4. Where is the dataset used ?
>>>> 5. Where is the source code used ?
>>>>
>>>> Thanks a bunch !
>>>> --
>>>> Mathieu
>>>
>>
>
>



--
Mathieu

_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.html

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users
Permalink | Reply | 1 comment |
headers
Ganesh Narayanasamy | 2 Jun 2010 01:34
Picon
Favicon

Deformation field in Demons registration


When I ran DiffeomorphicDemons fast symmetric algorithm between a pair of 3D data sets, I saved the deformation field which is attached here. Can you tell me what is wrong with this color image attached? How come this file is a "vector volume" while both the input images are "volume" and why cannot this deformation field be saved as a tiff file in Slicer?
-Ganesh
The New Busy think 9 to 5 is a cute idea. Combine multiple calendars with Hotmail. Get busy. 
_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.html

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users
Permalink | Reply | 1 comment |
headers
Tom Vercauteren | 2 Jun 2010 09:23

Re: Deformation field in Demons registration

Hi Ganesh,

Not sure I understand your question. A deformation field is a vector
field not an scalar field. If you are looking for the deformed image
you need to warp your input image with the deformation field you have.
This will give you a "volume" and not a "vector volume"...

As for slicer IO problems, you may want to ask them directly.

Cheers,
Tom

On Wed, Jun 2, 2010 at 01:34, Ganesh Narayanasamy
<nganesh76@...> wrote:
>
> When I ran DiffeomorphicDemons fast symmetric algorithm between a pair of 3D
> data sets, I saved the deformation field which is attached here. Can you
> tell me what is wrong with this color image attached? How come this file is
> a "vector volume" while both the input images are "volume" and why cannot
> this deformation field be saved as a tiff file in Slicer?
> -Ganesh
>
> ________________________________
> The New Busy think 9 to 5 is a cute idea. Combine multiple calendars with
> Hotmail. Get busy.
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Kitware offers ITK Training Courses, for more information visit:
> http://www.kitware.com/products/protraining.html
>
> Please keep messages on-topic and check the ITK FAQ at:
> http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-users
>
>
_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.html

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users
Permalink | Reply | 0 comments |
headers
Roger Bramon Feixas | 2 Jun 2010 14:29
Picon

Re: Performance regression in ITK + gdcm 2.0 reading DICOM JPEG files

Hi,


I attach the test results on Windows. We can observe that trunk version is about 2x faster compared to gdcm2.0.14, however gdcm1.2 is still 2x faster than trunk version. 

I'm looking forward your answers and let's me know if you need some help.

Thanks,

Roger

On Tue, Jun 1, 2010 at 9:40 PM, Roger Bramon Feixas <rogerbramon <at> gmail.com> wrote:
Hi Mathieu,

Thanks for you time! I just tested it on Mac OS X (I will also test it on Windows soon) and there is a really nice performance improvement, about 2x. 

However, I have two code comments: First, there is an abort() at line 350 of gdcmBitmap.cxx which causes a crash, I commented it to be able to test it. The other problem is related to the use of gdcm::PixelFormat::UINT8 at line 342 of gdcmBitmap.cxx. If you're working with JPEG files and its pixel format is more than 8 bits, a downscale has to be done. Perhaps, it could be replaced for GetPixelFormat().

On the other hand, we want to include this commit in our GDCM compilation but I would like to know which is the best way:
- Using GDCM 2.0.14 patched. It's the last stable version but I don't know if patch it is a good idea because some changes have been done.
- Using the trunk version even though it isn't stable, or
- waiting for GDCM 2.1 if it will be released soon.

I will try to test it on Windows tomorrow and send you the time results.

Many thanks!

Roger


On Mon, May 31, 2010 at 11:21 PM, Mathieu Malaterre <mathieu.malaterre-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:
Hi Roger,

 Yes ! I finally found some time, Snakes on a Plane gave this opportunity ;)

 I think I made some progress with the JPEG DICOM file. JPEG 2000
files will be slightly more difficult. See log at:
http://sourceforge.net/tracker/index.php?func=detail&aid=3009661&group_id=137895&atid=739587

 For reference, this issue entered GDCM 2.x, when handling the
lossy/lossless duality of J2K and JPEG-LS in DICOM. I am required to
know after ExecuteInformation whether the encapsulated JPEG fragment
was compressed with lossy algorithm (this is not indicated in the
DICOM dataset).

 Could you please confirm this on your side (again JPEG 2000 will
take me some more time) ?

Thanks again for your report

On Fri, May 28, 2010 at 12:00 AM, Roger Bramon Feixas
<rogerbramon-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:
> Hi,
> Have you had time to think about the causes of the performance regression
> problem? I can work on it, but I would like to have, if you have it, some
> information which helps me to know where I need to focus the effort.
>
> Thanks,
> Roger
>
> On Tue, May 18, 2010 at 9:14 AM, Roger Bramon Feixas <rogerbramon-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org>
> wrote:
>>
>> Hi Mathieu,
>> Perhaps, the answer of question 3 wasn't so clear:
>> ITK versions:
>>   - CMAKE_BUILD_TYPE:RelWithDebInfo
>>   - BUILD_SHARED_LIBS: OFF
>> GDCM: Release from  http://gdcm.sourceforge.net
>> Test program: Release.
>> Compiling the test program in Debug mode, it takes more time but we also
>> have large differences between GDCM 1.2.x and GDCM 2.0.
>>
>> Roger
>> On Mon, May 17, 2010 at 11:50 PM, Roger Bramon Feixas
>> <rogerbramon-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:
>>>
>>> Hi Mathieu,
>>> Thanks for your attention. I answer your questions:
>>> 1. What is your platform ?
>>> I used a Windows XP 32bits to do the test. However, this behaviour also
>>> occurs in Mac OS X platform.
>>> 2. What is your compiler ?
>>> All ITK versions were compiled with Visual Studio 2008 and we downloaded
>>> GDCM 2.0.14 binary from http://gdcm.sourceforge.net
>>> 3. What are your compiler option (eg. I hope you do not run bench
>>> without CMAKE_BUILD_TYPE:Release) ?
>>> RelWithDebInfo
>>> 4. Where is the dataset used ?
>>> Dataset 5 is from OSIRIX: http://pubimage.hcuge.ch:8080/DATA/CALIX.zip
>>> -> CALIX/CT1 abdomen/D30MN BILISCOPIN
>>> We are really interested to improve the time needed to load jpeg lossless
>>> datasets. The jpeg lossless datasets used are non-anonymized datasets, I
>>> will try to anonymize them if you want them.
>>> 5. Where is the source code used ?
>>> I attach it in this mail.
>>>
>>> Thanks!
>>> Roger
>>> On Mon, May 17, 2010 at 10:44 PM, Mathieu Malaterre
>>> <mathieu.malaterre-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:
>>>>
>>>> Hi Roger,
>>>>
>>>>  This is a very interesting post !
>>>>
>>>> On Thu, May 13, 2010 at 11:00 PM, Roger Bramon Feixas
>>>> <rogerbramon-Re5JQEeQqe8AvxtiuMwx3w@public.gmane.org> wrote:
>>>> > Hi,
>>>> > Recently we updated our ITK version from ITK 3.16 to ITK 3.18 and we
>>>> > decided
>>>> > to relink ITK with system GDCM 2.0. After some tests, we concluded ITK
>>>> > 3.18+GDCM 2.0 is faster than ITK 3.18 reading DICOM Little Endian
>>>> > files,
>>>> > however ITK 3.18+GDCM 2.0 is rather slower reading DICOM JPEG files. I
>>>> > attach a PDF file which is a time comparison of ITK 2.8, 3.16, 3.18
>>>> > and
>>>> > 3.18+GDCM2.0 versions.
>>>> > I don't know if it's just a GDCM problem or if it depends on how ITK
>>>> > uses
>>>> > GDCM.
>>>>
>>>>                       GDCM             Reading
>>>> UpdateOutput-        UpdateLargest-
>>>> Data       ITK version version          directory         Information
>>>>        PossibleRegion
>>>>
>>>>
>>>>       5           2,8            1,2           642
>>>>  6                  20153
>>>>                  3,16            1,2           547
>>>>  0                  25121
>>>>                  3,18            1,2           547
>>>>  0                  25124
>>>>                  3,18              2         25297
>>>> 187                  71437
>>>>
>>>>
>>>>
>>>> Data     Transfer syntax                Modality          Files
>>>>        Size (MB)
>>>>
>>>>       5 JPEG 2000                      CT
>>>> 243                     23
>>>>
>>>>
>>>>
>>>>  Could you please post a little bit more on :
>>>> 1. What is your platform ?
>>>> 2. What is your compiler ?
>>>> 3. What are your compiler option (eg. I hope you do not run bench
>>>> without CMAKE_BUILD_TYPE:Release) ?
>>>> 4. Where is the dataset used ?
>>>> 5. Where is the source code used ?
>>>>
>>>> Thanks a bunch !
>>>> --
>>>> Mathieu
>>>
>>
>
>



--
Mathieu


Attachment (Performance-itk-gdcm.pdf): application/pdf, 25 KiB
_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.html

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users
Permalink | Reply | 0 comments |
headers
Luca Antiga | 2 Jun 2010 16:37
Picon
Gravatar

Re: MultiScaleHessianBasedMeasureImageFilter

Hi Juan Carlos,
  I ran the test for the filter in Testing/Code/Review on one of my  
images and it worked correctly, the scales image was filled with scale  
values. Find the code I ran attached.
Can you try it out on your image?

Luca
Attachment (foo.cxx): application/octet-stream, 8 KiB
cmake_minimum_required(VERSION 2.4)

PROJECT(FOO)

SET(FOO_SRCS
  foo.cxx
  )

INCLUDE (${CMAKE_ROOT}/Modules/FindITK.cmake)

IF (USE_ITK_FILE)
  INCLUDE (${USE_ITK_FILE})
ENDIF (USE_ITK_FILE)

ADD_EXECUTABLE(foo ${FOO_SRCS})
TARGET_LINK_LIBRARIES(foo ITKCommon ITKBasicFilters ITKAlgorithms ITKIO ITKStatistics)



On May 31, 2010, at 8:18 PM, Juan Carlos Prieto wrote:

> Hello Luis,
>
> I'm trying to use the MultiScaleHessianBasedMeasureImageFilter.
> My problem is that I don't know how to get the "GetScalesOutput()"  
> to work, I have both methods to On "GenerateHessianOutputOn()" and
>         "GenerateScalesOutputOn()", but the ScalesOutput produces  
> nothing.
>
> This is my code
>
> _hessianvessel->SetAlpha1 (0.5);
>         _hessianvessel->SetAlpha2 (2);
>
>         _hessianmultiscale->SetInput(image);
>         _hessianmultiscale->SetHessianToMeasureFilter(_hessianvessel);
>         _hessianmultiscale->SetSigmaMinimum(0.5);
>         _hessianmultiscale->SetSigmaMaximum(4.0);
>         _hessianmultiscale->SetNumberOfSigmaSteps(5);
>         _hessianmultiscale->GenerateHessianOutputOn();
>         _hessianmultiscale->GenerateScalesOutputOn();
>         _hessianmultiscale->Update();
>
>         _outputimagefloat = ImageTypeFloat::New();
>         _outputimagefloat->SetRegions(_hessianmultiscale- 
> >GetOutput()->GetRequestedRegion());
>         _outputimagefloat->Allocate();
>
>         TypeConstRegionIterator constit( _hessianmultiscale- 
> >GetScalesOutput(),
>                          _hessianmultiscale->GetScalesOutput()- 
> >GetRequestedRegion() );
>
>         TypeRegionIterator outputit( _outputimagefloat,  
> _outputimagefloat->GetRequestedRegion() );
>
>         for ( constit.GoToBegin(); !constit.IsAtEnd(); ++constit);
>         {
>             cout<<constit.Get()<< " "; //this output produces "0" only
>         }
>         /*for ( constit.GoToBegin(), outputit.GoToBegin(); ! 
> constit.IsAtEnd(); ++constit, ++outputit)
>         {
>             outputit.Set(constit.Get()); //does not work
>         }*/
>
>
> I hope this is enough information
>
> Thank you for your time
>
> Juan Carlos Prieto
>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Kitware offers ITK Training Courses, for more information visit:
> http://www.kitware.com/products/protraining.html
>
> Please keep messages on-topic and check the ITK FAQ at:
> http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-users


_____________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at
http://www.kitware.com/opensource/opensource.html

Kitware offers ITK Training Courses, for more information visit:
http://www.kitware.com/products/protraining.html

Please keep messages on-topic and check the ITK FAQ at:
http://www.itk.org/Wiki/ITK_FAQ

Follow this link to subscribe/unsubscribe:
http://www.itk.org/mailman/listinfo/insight-users
Permalink | Reply | 2 comments |

Gmane