edwardferdian / wss_mri_calculator Goto Github PK

Hello,
thanks for sharing this beautiful work. I have a soubt about the parabolic wss computation. At line

https://github.com/EdwardFerdian/wss_mri_calculator/blob/e206a83ccd91833436f1fad27f520678fced77d0/src/calculate_wss.py#LL101C17-L101C35

you perform the "sampling" of the velocity at specific points inside the domain, which do not coincide with voxel centers. This involves an interpolation that I think is linear (but i am not sure as I could't find details in the pyvista documentation), so the overall interpolation procedure would give results that are not parabolic.

If I understood well in this paper, you should pick the velocity values at voxel centers that are nearest to the locations of the surface wall, without interpolating them. At that point, you will perform the parabolic interpolation.

Am I wrong ?