<div dir="ltr"><div>Hi,</div><div>This filter is indeed not wrapped currently due to the notwrap extension:</div><div><a href="https://github.com/SimonRit/RTK/blob/master/wrapping/rtkReg1DExtractShroudSignalImageFilter.notwrap">https://github.com/SimonRit/RTK/blob/master/wrapping/rtkReg1DExtractShroudSignalImageFilter.notwrap</a></div><div>I have opened an issue to track this down in the future:</div><div><a href="https://github.com/SimonRit/RTK/issues/506">https://github.com/SimonRit/RTK/issues/506</a></div><div>Simon<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jul 29, 2022 at 8:01 AM dongping.wang--- via Rtk-users <<a href="mailto:rtk-users@public.kitware.com">rtk-users@public.kitware.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;" lang="ZH-CN"><div><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">Hi all<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">I am trying to use Reg1DExtractShroudSignalImageFilter with python. But I met  one error as the following:<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><img style="width: 8.4416in; height: 5.3583in;" id="gmail-m_-8825305984279446788图片_x0020_1" src="cid:18248a419af692e331" width="810" height="514"></span><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">Please help me check what this problem may be?<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">Regards<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">Wang Dongping<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><div><div style="border-color:rgb(225,225,225) currentcolor currentcolor;border-style:solid none none;border-width:1pt medium medium;padding:3pt 0cm 0cm"><p class="MsoNormal"><b><span style="font-size:11pt;font-family:等线">发件人<span lang="EN-US">:</span></span></b><span style="font-size:11pt;font-family:等线" lang="EN-US"> <a href="mailto:dongping.wang@simmhealthcare.com" target="_blank">dongping.wang@simmhealthcare.com</a> <<a href="mailto:dongping.wang@simmhealthcare.com" target="_blank">dongping.wang@simmhealthcare.com</a>> <br></span><b><span style="font-size:11pt;font-family:等线">发送时间<span lang="EN-US">:</span></span></b><span style="font-size:11pt;font-family:等线" lang="EN-US"> 2022</span><span style="font-size:11pt;font-family:等线">年<span lang="EN-US">7</span>月<span lang="EN-US">25</span>日<span lang="EN-US"> 12:33<br></span><b>收件人<span lang="EN-US">:</span></b><span lang="EN-US"> <a href="mailto:rtk-users@openrtk.org" target="_blank">rtk-users@openrtk.org</a><br></span><b>主题<span lang="EN-US">:</span></b><span lang="EN-US"> </span>回复<span lang="EN-US">: scatterglarecorrection <u></u><u></u></span></span></p></div></div><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">Hi all<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">I have a 3d raw data which has a size of 512*512*1012.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">When I use the following code:<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">inputImage=ROIRegionFilter.GetOutput()<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">SGFilterType=rtk.ScatterGlareCorrectionImageFilter[ImageType,ImageType,itk.F]<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">SGFilter=SGFilterType.New()<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">SGFilter.SetInput(ROIRegionFilter.GetOutput())<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">SGFilter.SetCoefficients([0.0787,10.6244])<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">SGFilter.SetTruncationCorrection(0.5)<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">And Extract a full size:<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">I got one error said:<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal" style="background:white none repeat scroll 0% 0%;vertical-align:baseline;word-break:break-all"><b><span style="font-size:10.5pt;font-family:"Courier New";color:rgb(178,43,49)" lang="EN-US">RuntimeError</span></b><span style="font-size:10.5pt;font-family:"Courier New";color:black" lang="EN-US">: C:\P\IPP\ITK-source\ITK\Modules\Filtering\FFT\include\itkVnlRealToHalfHermitianForwardFFTImageFilter.hxx:58:<u></u><u></u></span></p><p class="MsoNormal" style="background:white none repeat scroll 0% 0%;vertical-align:baseline;word-break:break-all"><span style="font-size:10.5pt;font-family:"Courier New";color:black" lang="EN-US">ITK ERROR: VnlRealToHalfHermitianForwardFFTImageFilter(000002AD61E8C4E0): Cannot compute FFT of image with size [1024, 1024, 1012]. VnlRealToHalfHermitianForwardFFTImageFilter operates only on images whose size in each dimension has a prime factorization consisting of only 2s, 3s, or 5s.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">When I extract the slices upto 512 or 2,  then the python hung up.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">If I just extracted one of the slices, it is ok.<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">So my question is:<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">ScatterGlareCorrectionImageFilter is only working with a slice (for example 512X512X1)</span><span style="font-size:10.5pt;font-family:等线">？<span lang="EN-US"><u></u><u></u></span></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">If I wanted to use this filter, I have to walk through the raw data slice by slice? Thanks!<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">Regards<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">Wang Dongping<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><div><div style="border-color:rgb(225,225,225) currentcolor currentcolor;border-style:solid none none;border-width:1pt medium medium;padding:3pt 0cm 0cm"><p class="MsoNormal"><b><span style="font-size:11pt;font-family:等线">发件人<span lang="EN-US">:</span></span></b><span style="font-size:11pt;font-family:等线" lang="EN-US"> <a href="mailto:dongping.wang@simmhealthcare.com" target="_blank">dongping.wang@simmhealthcare.com</a> <<a href="mailto:dongping.wang@simmhealthcare.com" target="_blank">dongping.wang@simmhealthcare.com</a>> <br></span><b><span style="font-size:11pt;font-family:等线">发送时间<span lang="EN-US">:</span></span></b><span style="font-size:11pt;font-family:等线" lang="EN-US"> 2022</span><span style="font-size:11pt;font-family:等线">年<span lang="EN-US">7</span>月<span lang="EN-US">19</span>日<span lang="EN-US"> 11:29<br></span><b>收件人<span lang="EN-US">:</span></b><span lang="EN-US"> '<a href="mailto:rtk-users@openrtk.org" target="_blank">rtk-users@openrtk.org</a>' <<a href="mailto:rtk-users@openrtk.org" target="_blank">rtk-users@openrtk.org</a>><br></span><b>主题<span lang="EN-US">:</span></b><span lang="EN-US"> scatterglarecorrection <u></u><u></u></span></span></p></div></div><p class="MsoNormal"><span lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">Hi all<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US">I am trying to use ScatterGlareCorrectionImageFilter to do some scattering correction. From the test, we can see we just need to set the coef, then call the update().<u></u><u></u></span></p><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><table style="background:white none repeat scroll 0% 0%;border-collapse:collapse" cellspacing="0" cellpadding="0" border="0"><tbody><tr><td style="padding:0cm 7.5pt" valign="top"><p class="MsoNormal" style="line-height:15pt"><span style="font-size:9pt;font-family:Consolas;color:rgb(36,41,47)" lang="EN-US">ImageType::Pointer testImage = createInputImage(coef);<u></u><u></u></span></p></td></tr><tr><td style="width:37.5pt;background:transparent none repeat scroll 0% 0%;padding:0cm 7.5pt" width="50" valign="top" nowrap></td><td style="background:transparent none repeat scroll 0% 0%;padding:0cm 7.5pt" valign="top"><p class="MsoNormal" style="line-height:15pt"><span style="font-size:9pt;font-family:Consolas;color:rgb(36,41,47)" lang="EN-US">  SFilter->SetInput(testImage);<u></u><u></u></span></p></td></tr><tr><td style="width:37.5pt;padding:0cm 7.5pt" width="50" valign="top" nowrap></td><td style="padding:0cm 7.5pt" valign="top"><p class="MsoNormal" style="line-height:15pt"><span style="font-size:9pt;font-family:Consolas;color:rgb(36,41,47)" lang="EN-US">  TRY_AND_EXIT_ON_ITK_EXCEPTION(SFilter->Update())<u></u><u></u></span></p></td></tr></tbody></table><p class="MsoNormal"><span style="font-size:10.5pt;font-family:等线" lang="EN-US"><u></u> <u></u></span></p><div style="border-color:rgb(168,184,217) rgb(168,184,217) currentcolor;border-style:solid solid none;border-width:1pt 1pt medium;padding:6pt 6pt 0cm"><h2 style="margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:rgb(226,232,242) none repeat scroll 0% 0%"><span style="font-size:10.5pt;font-family:等线;color:black;font-weight:normal" lang="EN-US">The question is if we need to call the UpdateFFTProjectionsConvolutionKernel() to build the filter before we call update()?</span><span style="font-size:10.5pt;font-family:等线;font-weight:normal" lang="EN-US"><u></u><u></u></span></h2><h2 style="margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:rgb(226,232,242) none repeat scroll 0% 0%"><span style="font-size:10.5pt;font-family:等线;color:black;font-weight:normal" lang="EN-US">If we wanted to correct one 3D raw project data, how we do the correction slice by slice?  Any example for it?</span><span style="font-size:10.5pt;font-family:等线;font-weight:normal" lang="EN-US"><u></u><u></u></span></h2><h2 style="margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:rgb(226,232,242) none repeat scroll 0% 0%"><span style="font-size:10.5pt;font-family:等线;font-weight:normal" lang="EN-US"><u></u> <u></u></span></h2><h2 style="margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:rgb(226,232,242) none repeat scroll 0% 0%"><span style="font-size:10.5pt;font-family:等线;font-weight:normal" lang="EN-US"><u></u> <u></u></span></h2><h2 style="margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:rgb(226,232,242) none repeat scroll 0% 0%"><span style="font-size:10.5pt;font-family:等线;color:black;font-weight:normal" lang="EN-US">Regards</span><span style="font-size:10.5pt;font-family:等线;font-weight:normal" lang="EN-US"><u></u><u></u></span></h2><h2 style="margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:rgb(226,232,242) none repeat scroll 0% 0%"><span style="font-size:10.5pt;font-family:等线;color:black;font-weight:normal" lang="EN-US">Wang Dongping</span><span style="font-size:10.5pt;font-family:等线;font-weight:normal" lang="EN-US"><u></u><u></u></span></h2><h2 style="margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:rgb(226,232,242) none repeat scroll 0% 0%"><span style="font-size:10.5pt;font-family:等线;font-weight:normal" lang="EN-US"><u></u> <u></u></span></h2></div></div></div>_______________________________________________<br>
Rtk-users mailing list<br>
<a href="mailto:Rtk-users@public.kitware.com" target="_blank">Rtk-users@public.kitware.com</a><br>
<a href="https://public.kitware.com/cgi-bin/mailman/listinfo/rtk-users" rel="noreferrer" target="_blank">https://public.kitware.com/cgi-bin/mailman/listinfo/rtk-users</a><br>
</blockquote></div>