<div dir="ltr">I think you should interpolate the dose to a finer resolution ( by changing spacing). <br><br>Jothy<br><br><div class="gmail_quote">On Tue, Jun 29, 2010 at 4:37 PM, Xiaofeng Z <span dir="ltr"><<a href="mailto:xf10036@hotmail.com">xf10036@hotmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
<div>
Hi vtkusers,<br>
<br>
I need to overlay a iso-dose on top of a CT image. The dose are computed in a grid much coarser (2.5x2.5 pixel size) than the CT image (0.98x0.98 pixel size). Consequently, the boundary of the overlay appeared not smooth (see the attached image). I would like to know what filter to use to smooth the boundary.<br>
<br>
The code that places the overlay are as following:<br>
<br><font face="Consolas" size="2"><font face="Consolas" size="2">
vtkLookupTable* table = vtkLookupTable::New();<br>
table->SetRange(min, max); </font></font><font color="#008000" face="Consolas" size="2"><font color="#008000" face="Consolas" size="2"><font color="#008000" face="Consolas" size="2">// image intensity range<br></font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2">
table->SetValueRange(0.5, 1.0); </font></font><font color="#008000" face="Consolas" size="2"><font color="#008000" face="Consolas" size="2"><font color="#008000" face="Consolas" size="2">// from black to white<br></font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2">
table->SetSaturationRange(1.0, 1.0);<br>
table->SetHueRange(0.666, 0.0);<br>
table->SetAlpha(0.4);<br>
table->SetRampToLinear();<br>
table->Build();<br>
table->SetTableValue(0,0,0,0,0);<br></font></font>
<br><font face="Consolas" size="2"><font face="Consolas" size="2">
doseReslice = vtkImageReslice::New();<br>
doseReslice->SetOutputDimensionality(2);<br>
doseReslice->SetResliceAxes(reslice->GetResliceAxes());<br>
doseReslice->AddInput(dose);<br>
</font></font><font color="#008000" face="Consolas" size="2"><font color="#008000" face="Consolas" size="2"><font color="#008000" face="Consolas" size="2">// Map the image through the lookup table<br></font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2">
vtkImageMapToColors *color = vtkImageMapToColors::New();<br>
color->SetLookupTable(lut);<br>
color->SetInputConnection(doseReslice->GetOutputPort());<br>
color->Update();<br>
</font></font><font color="#008000" face="Consolas" size="2"><font color="#008000" face="Consolas" size="2"><font color="#008000" face="Consolas" size="2">// create an actor<br></font></font></font><font face="Consolas" size="2"><font face="Consolas" size="2">
vtkImageActor *actor = vtkImageActor::New();<br>
actor->SetInput(color->GetOutput());<br></font></font><br><br>Xiaofeng<br><br> <br><hr>The New Busy is not the old busy. Search, chat and e-mail from your inbox. <a href="http://www.windowslive.com/campaign/thenewbusy?ocid=PID28326::T:WLMTAGL:ON:WL:en-US:WM_HMP:042010_3" target="_blank">Get started.</a></div>
<br>_______________________________________________<br>
Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
<br>
Please keep messages on-topic and check the VTK FAQ at: <a href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">http://www.vtk.org/Wiki/VTK_FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>
<br></blockquote></div><br></div>