In this paper, a novel boundary reconstruction scheme for porous media flows is developed utilizing the reinitialization scheme of the level set method. The proposed scheme transfers usual digitally defined step-like boundary data to continuously smooth surface data. It is confirmed that the lattice Boltzmann simulation by the interpolated bounce-back boundary condition with the present scheme well reproduces the drag coefficients of a cylinder and a sphere. The results indicate that the present scheme improves predicition accuracy and it is more insensitive to grid resolutions than the standard scheme.