Humm, yes, I think you are right. Mea Culpa
It should works. I will try to review your script tomorrow.
Get Outlook for iOS<https://aka.ms/o0ukef>
________________________________
From: Andrew E. Slaughter <***@inl.gov>
Sent: Tuesday, November 27, 2018 12:39:35 AM
To: Jonathan Borduas; vtkusers
Subject: Re: [vtkusers] Mesh Interpolation
Thanks for the information, but I dont see where the vtkInterpolateDataSetAttribute documentation discusses implicit functions. My reason for going down this path is as follows:
The documentation states vtkInterpolateDataSetAttributes is a filter that interpolates data set attribute values between input data sets. This is exactly my situation. I have two vtkUnstructuredGrid objects (which is a vtkDataSet). Each of these objects contain vtkPointData (which is a vtkDataSetAttributes). These two objects are exactly the same with respect to the number of nodes, cells, etc.
If you can explain why my understanding of this class is wrong I would appreciate it.
I used the example you provided as a guide to setting up the vtkInterpolateDataSetAttribute object, but a problem I am seeing with the interpolation is that the number of input ports is always reporting 1 in the vtkInterpolateDataSetAttribute object, regardless of how many items I add. I cant set the number of points in python, since the SetNumberOfInputPorts is a protected member, thus not part of the python API. I will keep investigating this, thanks for your help.
From: Jonathan Borduas <***@caboma.com>
Date: Monday, November 26, 2018 at 10:00 PM
To: "Andrew E. Slaughter" <***@inl.gov>, vtkusers <***@vtk.org>
Subject: RE: [vtkusers] Mesh Interpolation
See the test of the vtkInterpolateDataSetAttributes:
https://www.vtk.org/gitweb?p=VTK.git;a=blob;f=Filters/Hybrid/Testing/Python/3DMorph.py<https://urldefense.proofpoint.com/v2/url?u=https-3A__www.vtk.org_gitweb-3Fp-3DVTK.git-3Ba-3Dblob-3Bf-3DFilters_Hybrid_Testing_Python_3DMorph.py&d=DwMGaQ&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=h7heP8xwI1i_HikChvhFbEBurKirgfOCdwgBxB9lM8c&m=axwcLPpoiBAwLGukOnCod5sBNR8H6iwjHEZzqOt42f0&s=E0pqg65BmyPlyhMeb9i9kudpH5yoQmJsE7DtjwAH28w&e=>
Didnt understand why you would need to use that class though.
I think the name is misleading you. If you look at the example, this isnt for data interpolation ! It is for the interpolation of implicit functions (as seen in the test).
More on implicit functions: https://www.vtk.org/doc/nightly/html/classvtkImplicitFunction.html<https://urldefense.proofpoint.com/v2/url?u=https-3A__www.vtk.org_doc_nightly_html_classvtkImplicitFunction.html&d=DwMGaQ&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=h7heP8xwI1i_HikChvhFbEBurKirgfOCdwgBxB9lM8c&m=axwcLPpoiBAwLGukOnCod5sBNR8H6iwjHEZzqOt42f0&s=RSejo5gmFIa-9ILQrG7UZD5NgvcpxnwvBh9-mXzC0fo&e=>
vtkInterpolateDataSetAttributes actually looks more like an Optimal Transport Interpolator: https://www.youtube.com/watch?v=6nhfArEM4Sg<https://urldefense.proofpoint.com/v2/url?u=https-3A__www.youtube.com_watch-3Fv-3D6nhfArEM4Sg&d=DwMGaQ&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=h7heP8xwI1i_HikChvhFbEBurKirgfOCdwgBxB9lM8c&m=axwcLPpoiBAwLGukOnCod5sBNR8H6iwjHEZzqOt42f0&s=GR9UMU1eEHx7lhmZz0mJ_cSBRIGz2A56NzXupdmRRLk&e=>
Jonathan
From: Andrew E. Slaughter <***@inl.gov>
Sent: Monday, November 26, 2018 11:43 PM
To: Jonathan Borduas <***@caboma.com>; vtkusers <***@vtk.org>
Subject: Re: [vtkusers] Mesh Interpolation
I nearly have the problem figured out, but I cannot get vtkInterpolateDataSetAttributes to run. If anyone has experience with this class I would appreciate some insight into it use. I am using the object on line 114 of the attached script.
Thanks.
From: Jonathan Borduas <***@caboma.com<mailto:***@caboma.com>>
Date: Monday, November 26, 2018 at 11:17 AM
To: "Andrew E. Slaughter" <***@inl.gov<mailto:***@inl.gov>>, vtkusers <***@vtk.org<mailto:***@vtk.org>>
Subject: Re: [vtkusers] Mesh Interpolation
Hi Andrew,
One way to interpolate point and cell data between two mesh: you can use the vtkPointInterpolator with a Voronoi Kernel (or other kernel), do a itk::DiscreteGaussianImageFilter over the 3D image to smooth the interpolation and then do a vtkProbeFilter.
Best,
Jonathan
Get Outlook for iOS<https://urldefense.proofpoint.com/v2/url?u=https-3A__aka.ms_o0ukef&d=DwMF-g&c=54IZrppPQZKX9mLzcGdPfFD1hxrcB__aEkJFOKJFd00&r=h7heP8xwI1i_HikChvhFbEBurKirgfOCdwgBxB9lM8c&m=aAOXM3i-LCmwbKnx-O9gSjuXokpbn-5eRK65fXdzvCg&s=fr7AzqyWyyeMEM2RozZhbe697LCCwggGrGZtJhKSQAY&e=>
________________________________
From: vtkusers <vtkusers-***@public.kitware.com<mailto:vtkusers-***@public.kitware.com>> on behalf of Andrew E. Slaughter via vtkusers <***@public.kitware.com<mailto:***@public.kitware.com>>
Sent: Monday, November 26, 2018 1:09:55 PM
To: vtkusers
Subject: [vtkusers] Mesh Interpolation
I need to perform interpolation of data between two ExodusII files that differ because of refinement. I was hoping someone could point me in the
correct direction for performing this action. I have attached the Exodus files that I am working with as well as a script that reads and displays the two files.
I am using VTK7.1 (I can use 8 if needed) on MacOS.
Here is the steps that I envision, but I am open to all options:
1. Interpolate the result from the coarse mesh to the fine mesh (I dont know how to do this step).
2. Use vtkInterpolateDataSetAttributes to interpolate between the result from (1) and the fine mesh.
Thanks for the help,
Andrew