During binary-black-hole (BBH) mergers, energy and momenta are carried away from the binary system as gravitational radiation. Access to the radiated energy and momenta allows us to accurately predict the properties of the remnant black hole. We develop a python package gw_remnant to efficiently extract the remnant mass, remnant spin, peak luminosity and the final kick imparted on the remnant black hole directly from the gravitational radiation. We then compute the remnant properties of the final black hole in case of non-spinning BBH mergers with mass ratios ranging from q = 2.5 to q = 1000 using waveforms generated from BHPTNRSur1dq1e4, a recently developed numerical relativity informed surrogate model based on black-hole perturbation theory framework. We validate our results against the remnant properties estimated from numerical relativity (NR) surrogate models in the comparable mass ratio regime and against recently available high-mass ratio RIT NR simulations at q = [15,32,64]. We find that our remnant property estimates match very closely to the estimates obtained from NR surrogate model and the NR data respectively in both the regimes. We then present BHPTNR_Remnant, a surrogate model for the properties of the remnant black hole in BBH mergers with q = 2.5 to q = 1000, using Gaussian process regression fitting methods. Finally, we comment on the possible implication of remnant information in gravitational waveform modelling. We make both the gw_remnant and BHPTNR_Remnant packages publicly available.