<html><div style='background-color:'><P>Louis,</P>
<P>&nbsp; vtkImplicitDataSet interpolates the scalar data associated with its input (PolyDataA) at any given geometrical position (points in PolyDataB). This will not give you the geometrical intersection of the two sets of PolyData. I am not sure what you want mean by intersection: is it the set of lines (and points) that lie at the intersection of two polygonal surfaces? If so, I do not think there is any way to do this (generally) in vtk. </P>
<P>The class I have created, vtkImplicitPolyData, can be used to clip one PolyData with another, although with some error. The error will be dependant on the relative resolution (point density) of the two data sets, and the limitations of linear interpolation generally used in vtk. The output will generally be&nbsp;polygons, not lines. Here is a simple example: PolyData A (1) is overlapped by a similarly shaped PolyData B (2). The class vtkImplicitPolyData is used to clip PolyData A, using PolyData B as an implicit surface (3).</P><PRE><BR>&nbsp;A&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; A B&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; A'<BR>&nbsp;/\&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /\/\&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; /\<BR>/&nbsp; \&nbsp;&nbsp;&nbsp; / /\ \&nbsp;&nbsp;&nbsp; /<BR>\&nbsp; /&nbsp;&nbsp;&nbsp; \ \/ /&nbsp;&nbsp;&nbsp; \<BR>&nbsp;\/&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; \/\/&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; \/<BR>1.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.<BR></PRE><PRE>If you think this might be useful to you I can send you the code for vtkImplicitPolyData.</PRE><PRE>  regards</PRE><PRE>      Dave Pont</PRE><PRE>&nbsp;</PRE>&gt;From: louis@itri.org.tw 
<DIV></DIV>&gt;To: david_pont@hotmail.com 
<DIV></DIV>&gt;Subject: About Clipping PolyData with PolyData!! 
<DIV></DIV>&gt;Date: Wed, 5 Dec 2001 11:21:37 +0800 
<DIV></DIV>&gt; 
<DIV></DIV>&gt;Hi Dave Pont: 
<DIV></DIV>&gt;May I ask you about vtk problem? I had seen your mail about "Clipping 
<DIV></DIV>&gt;PolyData with PolyData" in vtk mailing list. 
<DIV></DIV>&gt; 
<DIV></DIV>&gt;My problem: 
<DIV></DIV>&gt;I have two PolyData surfaces(PolyDataA and PolyDataB),I want to find the 
<DIV></DIV>&gt;intersection of two PolyData surfaces. This problem like Chris's problem in 
<DIV></DIV>&gt;vtk mailing list. I had written some code in C++. Coulde you give me some 
<DIV></DIV>&gt;ideas? 
<DIV></DIV>&gt; 
<DIV></DIV>&gt;vtkImplicitDataSet* ImplicitDataSet=vtkImplicitDataSet::New(); 
<DIV></DIV>&gt;ImplicitDataSet-&gt;SetDataSet( PolyDataA ); 
<DIV></DIV>&gt;vtkClipPolyData* ClipPolyData=vtkClipPolyData::New(); 
<DIV></DIV>&gt;ClipPolyData-&gt;SetInput( PolyDataB ); 
<DIV></DIV>&gt;ClipPolyData-&gt;SetClipFunction( ImplicitDataSet ); 
<DIV></DIV>&gt;ClipPolyData-&gt;InsideOutOff(); 
<DIV></DIV>&gt; 
<DIV></DIV>&gt;My program will make mistakes in "ImplicitDataSet-&gt;SetDataSet( 
<DIV></DIV>&gt;PolyDataA-&gt;GetOutput())". Hope you can know my problem. My English is very 
<DIV></DIV>&gt;poor. ^_^ 
<DIV></DIV>&gt; 
<DIV></DIV>&gt;thanks with regards 
<DIV></DIV>&gt; 
<DIV></DIV>&gt; 
<DIV></DIV>&gt;Louis Lee 
<DIV></DIV>&gt;Color Technology Section 
<DIV></DIV>&gt;Industrial Technology Research Institute 
<DIV></DIV>&gt;louis@itri.org.tw 
<DIV></DIV>&gt;03-5913714 
<DIV></DIV>&gt; 
<DIV></DIV></div><br clear=all><hr>Get your FREE download of MSN Explorer at <a href='http://go.msn.com/bql/hmtag_itl_EN.asp'>http://explorer.msn.com</a><br></html>