“Computation” of IJK (voxel coordinates) from RAS or LPS (physical coordinates) is multiplication of the coordinates by a 4x4 matrix, so it cannot really be any simpler than that. The matrix you need is the RAS to IJK matrix, available in the volume node.