  A possible way to study the reionization of cosmic hydrogen is by observing
the large ionized regions (bubbles) around bright individual sources, e.g.,
quasars, using the redshifted 21 cm signal. It has already been shown that
matched filter-based methods are not only able to detect the weak 21 cm signal
from these bubbles but also aid in constraining their properties. In this work,
we extend the previous studies to develop a rigorous Bayesian framework to
explore the possibility of constraining the parameters that characterize the
bubbles. To check the accuracy with which we can recover the bubble parameters,
we apply our method on mock observations appropriate for the upcoming SKA1-low.
For a region of size $\gtrsim 50$ cMpc around a typical quasar at redshift 7,
we find that $\approx 20$ h of integration with SKA1-low will be able to
constrain the size and location of the bubbles, as well as the difference in
the neutral hydrogen fraction inside and outside the bubble, with $\lesssim
10\%$ precision. The recovery of the parameters are more precise and the SNR of
the detected signal is higher when the bubble sizes are larger and their shapes
are close to spherical. Our method can be useful in identifying regions in the
observed field which contain large ionized regions and hence are interesting
for following up with deeper integration times.
