nH is calculated using all points (from https://vizier.u-strasbg.fr site.) within r_arcmin
nh(ra = NULL, dec = NULL, r_arcmin = 10, showInfo = FALSE)
RA in degrees
DEC in degrees
search radius in arcmin
show additional information ? Default=false
weighted nH
if (FALSE) {
xrayr::nh(213.2, -65.4, r_arcmin = 30, showInfo = T)
}