<br><br><blockquote>Hi Bill,<br><br>Thanks for your reply.. I implemented that technique and it works fine. As i said in my previous post, outsid range values are not interpolated.<br><br>For example, pls see the attached image file for reference.. In the figure, the mesh which is generated, doesnt completely fall in the area of the scatter data.. <br>In the rwo white patches, the values are not interpolated. I mean i dont get the actual result there in those areas (as there is no scatter points avail there)<br><br>So, what can i do to get values for my entire mesh instead of getitng for partial mesh.??Either interolation or extrapolation.. Any means is available in VTK <br>to implement this..??<br><br>Thanks<br><br>Regards<br>Rakesh Patil<br>---------- Original message ----------<br>From:Bill Lorensen&lt; bill.lorensen@gmail.com &gt;<br>Date: 24 Apr 10 18:15:33<br>Subject:  Re: [vtkusers] Interpolation of data..<br>To: Rakesh Patil <rakeshthp@in.com><br><br>The last part of
  this example:<br>http://www.vtk.org/Wiki/VTK/Examples/InterpolateTerrain<br>does what you want.<br><br>Basically,<br>1) Add the lat/log and elevations to a polydata<br>2) Create a surface with Delaunay2D<br>3) Add the points with unknown elevation to a vtkPoints<br>4) Use a cell locator to cast infinite rays into the surface to find<br>the elevation.<br><br><br>On Fri, Apr 23, 2010 at 9:50 PM, Rakesh Patil <rakeshthp@in.com> wrote:<br>&gt;<br>&gt;<br>&gt; Well,<br>&gt;<br>&gt; I was trying with this example:<br>&gt;<br>&gt; <a target=\"_blank\" href="http://www.vtk.org/Wiki/VTK/Examples/InterpolateTerrain" target="_blank">http://www.vtk.org/Wiki/VTK/Examples/InterpolateTerrain</a><br>&gt;<br>&gt; In this example, there are two ways specified for interpolation.. One is<br>&gt; using<br>&gt; vtkProbeFilter and another is using vtkCellLocator<br>&gt;<br>&gt; First i tried using vtkProbeFilter. But failed getting the output.. But<br>&gt; using the<br>&gt; vtkCellLocator, i was 
 able to produce interpolated values for the third<br>&gt; column of MeshG.<br>&gt;<br>&gt; After observation, i foound that only the region which is common in both the<br>&gt; grids,<br>&gt; was interpolated. So i was fully stuck wondering how to find the remaining<br>&gt; values..<br>&gt;<br>&gt; any comments..??<br>&gt;<br>&gt; Thanks<br>&gt;<br>&gt; Rakesh Patil<br>&gt;<br>&gt; ---------- Original message ----------<br>&gt; From:"Dean Inglis"&lt; dean.inglis@camris.ca &gt;<br>&gt; Date: 23 Apr 10 23:59:48<br>&gt; Subject: Re: Re:[vtkusers]Interpolation of data..<br>&gt; To: "Rakesh Patil"<br>&gt;<br>&gt; How can you expect your mesh to get interpolated<br>&gt; when its z coordinate values are all 0<br>&gt; but the scatter data range in z from -700 to -900?<br>&gt; The bounds of each unstructured grid do not even<br>&gt; overlap / intersect!<br>&gt;<br>&gt; Dean<br>&gt;<br>&gt;&gt; Hi,My both the meshes are unstructured grid. My application is based on<br>&gt;&gt; FEA meth
 od.. ScatterG does contains scalars and the depth values are<br>&gt;&gt; stored as a scalar array for scatterG..I dont want coordinate values to be<br>&gt;&gt; interpolated.. As an example,scatterG contains data something like<br>&gt;&gt; this,Points: 87.23 16.2690086.24 17.8753087.18 17.27700......and the<br>&gt;&gt; scalar for scatterG is defined as 900530700......meshG contains<br>&gt;&gt; points84.1217.83 081.3315.28 083.2418.28 0......So,, what i want to do is,<br>&gt;&gt; using the third column of scatterG, i need to estimate the third column of<br>&gt;&gt; meshG..Any ideas how it can be done.??ThanksRakesh Patil Original message<br>&gt;&gt; From:"De an Inglis"&lt; dean.inglis@camris.ca &gt;Date: 23 Apr 10<br>&gt;&gt; 20:16:47Subject:Re: Re:[vtkusers]Interpolation of data..To: "Rakesh Patil"<br>&gt;&gt; vtkProbeFlter will probe the scalar values ateach spatial location of your<br>&gt;&gt; scattered points: youneed to have scalar, vector, or tensor data<br>&gt;&gt; asso
 ciatedwith the grid's points:scatterG&gt;GetPointData()&gt;SetScalars(<br>&gt;&gt; somevt<br>&gt;&gt; kDataArray );the vtkProbeFilter will not interpolate any coordinatevalues<br>&gt;&gt; like x , y,z locations.If you have depth valuesin a vtkDataArray that<br>&gt;&gt; correspond to the x, y, z locationsin your scatterG grid, then this should<br>&gt;&gt; work.Does anyone know ifmaybe the input to the probe filter needs to be a<br>&gt;&gt; structured data set likevtkPolyData? By the sound of your grid's name<br>&gt;&gt; "mesh" , this sounds like structureddata.Can it be converted to a<br>&gt;&gt; vtkPolyData?Also, you allocate a smart pointer as:vtkSmartPointer probe =<br>&gt;&gt; vtkSmartPointer: :New()Dean Original Message From: "Rakesh Patil" To: Cc:<br>&gt;&gt; Sent: Friday, April 23, 2010 9:40 AMSubject: Re: Re:[vtkusers]<br>&gt;&gt; Interpolation of data..&gt; Hi..I can independently view those grids... This<br>&gt;&gt; is what I have done..vtkUnstructuredGrid *scatterG
  =<br>&gt;&gt; MyApp::getInstance()&gt;getScatterGrid();vtkUnstructuredGrid *meshG =<br>&gt;&gt; MyApp::getInstance()&gt;getMeshGrid();vtkSmartPointer probe =<br>&gt;&gt;<br>&gt;&gt; vtkSmartPointer::New();probe&gt;SetSource(scatterG);probe&gt;SetInput(meshG);probe&gt;Update();vtkDataArr<br>&gt;&gt; ay *data = probe&gt;GetOutput()&gt;GetPointData()&gt;GetScalar();vtkDoubleArray<br>&gt;&gt; *doubleArray =<br>&gt;&gt;<br>&gt;&gt; vtkDoubleArray::SafeDownCast(data);MyContour::getInstance()&gt;setDisplayArray(doubleArray);After<br>&gt;&gt; doing this, i click on contour display button, and the contour for<br>&gt;&gt; bathymetry gets displayed for meshG.Initially, meshG doesnt have<br>&gt;&gt; depth/elevation (z) value. z value is always 0 for every node point. x and<br>&gt;&gt; y corresponds to the lon and lat respectively.But the scatter data<br>&gt;&gt; contains lon, lat, and depth.. So, once the interpolation is done, i<br>&gt;&gt; should get the bathymetry (depth/elevation) 
 for meshG also..'MyContour' in<br>&gt;&gt; the above code is the class which displays the contours for wh&gt; atever<br>&gt;&gt; data is being passed.So, any suggestions or &gt; comments..??ThanksRegards<br>&gt;&gt; Rakesh Patil Original messageFrom:Dean Inglis&lt; &gt; dean.inglis@camris.ca &gt;<br>&gt;&gt; Date: 23 Apr 10 16:49:50 Subject:Re: &gt; Re:[vtkusers]Interpolation of<br>&gt;&gt; data.. To: "Rakesh Patil"&gt;<br>&gt;&gt;<br>&gt;<br>&gt;<br>&gt; _______________________________________________<br>&gt; Powered by <a target=\"_blank\" href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>&gt;<br>&gt; Visit other Kitware open-source projects at<br>&gt; <a target=\"_blank\" href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>&gt;<br>&gt; Please keep messages on-topic and check the VTK FAQ at:<br>&gt; <a target=\"_blank\" href="http://www.vtk.org/Wiki/VTK_FAQ" target="_blank">htt
 p://www.vtk.org/Wiki/VTK_FAQ</a><br>&gt;<br>&gt; Follow this link to subscribe/unsubscribe:<br>&gt; <a target=\"_blank\" href="http://www.vtk.org/mailman/listinfo/vtkusers" target="_blank">http://www.vtk.org/mailman/listinfo/vtkusers</a><br>&gt;<br>&gt;<br></rakeshthp@in.com></rakeshthp@in.com></blockquote>