code: def zoomedHealpixMap(title, map, lon, lat, radius,
                       xsize=1000, **kwargs):
      
      reso = 60. * 2. * radius / xsize # Deg to arcmin
      hp.gnomview(map=map, rot=[lon, lat, 0], title=title, xsize=xsize, reso=reso, degree=False, **kwargs)
docs: Inputs: lon (deg), lat (deg), radius (deg)