Construction of initial data for 3+1 numerical relativity