<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
{font-family:宋体;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:等线;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"\@宋体";
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"\@等线";
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:12.0pt;
font-family:宋体;}
h2
{mso-style-priority:9;
mso-style-link:"标题 2 字符";
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:18.0pt;
font-family:宋体;}
span.2
{mso-style-name:"标题 2 字符";
mso-style-priority:9;
mso-style-link:"标题 2";
font-family:宋体;
font-weight:bold;}
span.EmailStyle20
{mso-style-type:personal-reply;
font-family:等线;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 90.0pt 72.0pt 90.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=ZH-CN link=blue vlink=purple style='word-wrap:break-word'><div class=WordSection1><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>Hi all<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>I have a 3d raw data which has a size of 512*512*1012.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>When I use the following code:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>inputImage=ROIRegionFilter.GetOutput()<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>SGFilterType=rtk.ScatterGlareCorrectionImageFilter[ImageType,ImageType,itk.F]<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>SGFilter=SGFilterType.New()<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>SGFilter.SetInput(ROIRegionFilter.GetOutput())<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>SGFilter.SetCoefficients([0.0787,10.6244])<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>SGFilter.SetTruncationCorrection(0.5)<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>And Extract a full size:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>I got one error said:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><p class=MsoNormal style='background:white;vertical-align:baseline;word-break:break-all'><b><span lang=EN-US style='font-size:10.5pt;font-family:"Courier New";color:#B22B31'>RuntimeError</span></b><span lang=EN-US style='font-size:10.5pt;font-family:"Courier New";color:black'>: C:\P\IPP\ITK-source\ITK\Modules\Filtering\FFT\include\itkVnlRealToHalfHermitianForwardFFTImageFilter.hxx:58:<o:p></o:p></span></p><p class=MsoNormal style='background:white;vertical-align:baseline;word-break:break-all'><span lang=EN-US style='font-size:10.5pt;font-family:"Courier New";color:black'>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.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>When I extract the slices upto 512 or 2, then the python hung up.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>If I just extracted one of the slices, it is ok.<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>So my question is:<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>ScatterGlareCorrectionImageFilter is only working with a slice (for example 512X512X1)</span><span style='font-size:10.5pt;font-family:等线'>？<span lang=EN-US><o:p></o:p></span></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>If I wanted to use this filter, I have to walk through the raw data slice by slice? Thanks!<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>Regards<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>Wang Dongping<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><div><div style='border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm'><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:等线'>发件人<span lang=EN-US>:</span></span></b><span lang=EN-US style='font-size:11.0pt;font-family:等线'> dongping.wang@simmhealthcare.com <dongping.wang@simmhealthcare.com> <br></span><b><span style='font-size:11.0pt;font-family:等线'>发送时间<span lang=EN-US>:</span></span></b><span lang=EN-US style='font-size:11.0pt;font-family:等线'> 2022</span><span style='font-size:11.0pt;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> 'rtk-users@openrtk.org' <rtk-users@openrtk.org><br></span><b>主题<span lang=EN-US>:</span></b><span lang=EN-US> scatterglarecorrection <o:p></o:p></span></span></p></div></div><p class=MsoNormal><span lang=EN-US><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>Hi all<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'>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().<o:p></o:p></span></p><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><table class=MsoNormalTable border=0 cellspacing=0 cellpadding=0 style='background:white;border-collapse:collapse'><tr><td valign=top style='padding:0cm 7.5pt 0cm 7.5pt'><p class=MsoNormal style='line-height:15.0pt'><span lang=EN-US style='font-size:9.0pt;font-family:Consolas;color:#24292F'>ImageType::Pointer testImage = createInputImage(coef);<o:p></o:p></span></p></td></tr><tr><td width=50 nowrap valign=top style='width:37.5pt;background:transparent;padding:0cm 7.5pt 0cm 7.5pt'></td><td valign=top style='background:transparent;padding:0cm 7.5pt 0cm 7.5pt'><p class=MsoNormal style='line-height:15.0pt'><span lang=EN-US style='font-size:9.0pt;font-family:Consolas;color:#24292F'> SFilter->SetInput(testImage);<o:p></o:p></span></p></td></tr><tr><td width=50 nowrap valign=top style='width:37.5pt;padding:0cm 7.5pt 0cm 7.5pt'></td><td valign=top style='padding:0cm 7.5pt 0cm 7.5pt'><p class=MsoNormal style='line-height:15.0pt'><span lang=EN-US style='font-size:9.0pt;font-family:Consolas;color:#24292F'> TRY_AND_EXIT_ON_ITK_EXCEPTION(SFilter->Update())<o:p></o:p></span></p></td></tr></table><p class=MsoNormal><span lang=EN-US style='font-size:10.5pt;font-family:等线'><o:p> </o:p></span></p><div style='border:solid #A8B8D9 1.0pt;border-bottom:none;padding:6.0pt 6.0pt 0cm 6.0pt'><h2 style='mso-margin-top-alt:5.0pt;margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:#E2E8F2'><span lang=EN-US style='font-size:10.5pt;font-family:等线;color:black;font-weight:normal'>The question is if we need to call the UpdateFFTProjectionsConvolutionKernel() to build the filter before we call update()?</span><span lang=EN-US style='font-size:10.5pt;font-family:等线;font-weight:normal'><o:p></o:p></span></h2><h2 style='mso-margin-top-alt:5.0pt;margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:#E2E8F2'><span lang=EN-US style='font-size:10.5pt;font-family:等线;color:black;font-weight:normal'>If we wanted to correct one 3D raw project data, how we do the correction slice by slice? Any example for it?</span><span lang=EN-US style='font-size:10.5pt;font-family:等线;font-weight:normal'><o:p></o:p></span></h2><h2 style='mso-margin-top-alt:5.0pt;margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:#E2E8F2'><span lang=EN-US style='font-size:10.5pt;font-family:等线;font-weight:normal'><o:p> </o:p></span></h2><h2 style='mso-margin-top-alt:5.0pt;margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:#E2E8F2'><span lang=EN-US style='font-size:10.5pt;font-family:等线;font-weight:normal'><o:p> </o:p></span></h2><h2 style='mso-margin-top-alt:5.0pt;margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:#E2E8F2'><span lang=EN-US style='font-size:10.5pt;font-family:等线;color:black;font-weight:normal'>Regards</span><span lang=EN-US style='font-size:10.5pt;font-family:等线;font-weight:normal'><o:p></o:p></span></h2><h2 style='mso-margin-top-alt:5.0pt;margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:#E2E8F2'><span lang=EN-US style='font-size:10.5pt;font-family:等线;color:black;font-weight:normal'>Wang Dongping</span><span lang=EN-US style='font-size:10.5pt;font-family:等线;font-weight:normal'><o:p></o:p></span></h2><h2 style='mso-margin-top-alt:5.0pt;margin-right:11.25pt;margin-bottom:0cm;margin-left:0cm;background:#E2E8F2'><span lang=EN-US style='font-size:10.5pt;font-family:等线;font-weight:normal'><o:p> </o:p></span></h2></div></div></body></html>