A data driven method for OD matrix estimation