I had a little time on a flight today to think about your problem, and so I applied the standard integration method to see whether or not your equation could be explicitly integrated (in the sense that the Monge-Ampère equation $u_{xx}u_{yy}-{u_{xy}}^2=1$ can be integrated by transforming it to Laplace's equation). The answer is that it cannot be integrated that way. However, one can make the problem equivalent to a linear one, at least locally, that has an explicit solution in series, and this may or may not help you. I'll record the results here, just in case you find it useful.
Suppose that we have a solution $B(x,y)$ on some simply connected domain $D$ in the $xy$-plane and suppose that it satisfies $B_{xx} = -2p^2 < 0$, $B_{yy} = -q^2 < 0$, and $B_{xy} = 2pq$ for some positive functions $p$ and $q$. (One can also treat the case when one of $p$ or $q$ is negative, but that's a minor variation that I'll leave to you.) Then, we have
$$
\mathrm{d} B_x = -2p^2\,\mathrm{d} x + 2pq\,\mathrm{d} y
\quad\text{and}\quad
\mathrm{d} B_y = 2pq\,\mathrm{d} x - q^2\,\mathrm{d} y
$$
Thus, the functions $p$ and $q$ must satisfy
$$
\mathrm{d}\bigl(-2p^2\,\mathrm{d} x + 2pq\,\mathrm{d} y\bigr)
=\mathrm{d}\bigl(2pq\,\mathrm{d} x - q^2\,\mathrm{d} y\bigr) = 0,
$$
and, conversely, if $p$ and $q$ satisfy these two conditions (which are first order PDE), then the above equations determine $B_x$ and $B_y$ (up to an additive constant) and then $\mathrm{d}B = B_x\,\mathrm{d}x + B_x\,\mathrm{d}y$ determines $B$ up to an additive constant, so the two PDE on $p$ and $q$ are essentially equivalent to the original equations.
Next, make a change of variables: Set $s = px$ and $t = qy - px$ and then $p= e^u$ and $q = e^v$. Then set $w = s + it$ and $z = -\tfrac12(v-iu)$. Then the above two PDE relating $(x,y,p,q)$ simply become the real and imaginary parts of the linear elliptic complex equation
$$
\frac{\partial w}{\partial \bar z} = \bar w
$$
Now, this equation is known not to be integrable by the method of Darboux, so, in particular, you cannot write down its general solution in terms of a single holomorphic function of $z$. However, all the $C^1$ solutions are real analytic (because the equation is elliptic), and there is an explicit representation of the analytic solution in terms of power series:
$$
w(z)=\sum_{k=0}^\infty c_k\ f^{(k)}(z\bar z) z^k+\overline{c_k}\ f^{(k+1)}(z\bar z) \bar z^{k+1}
$$
where $f^{(k)}$ is the $k$-th derivative of the modified Bessel $I$-function whose series representation is
$$
f(r) = 1 + \sum_{j=1}^\infty \frac{r^j}{(j!)^2}
$$
Using this representation, you can trace back through and integrate by parts to get series representations of $B(x,y)$, $x$, and $y$ in terms of $u$ and $v$, which gives you a graphical representation of the local solutions in this case.
Whether this will help you with the global solutions, I don't know. However, it does give you the 'general' solution.