Oh right - in the x-y plane the lower limits would be 0, and the upper limit would be defined by the curve where z(x,y) intersects the x-y plane. The trouble is that you need to find, say, limits of x from 0 to f(y) ... where f is a function of y alone ... which is quite difficult from the expression. This is why you change the coordinates around: so the limits in the u-v axes are findable.
So:
Draw axis u vs v.
You know x>0 ... that should give you some relationships for u and v ... it may tell you that u<a: a is some number ... so you draw a line u=a ... all the u values are less than that. You may find that u<f(v) or v>f(u) ... so sketch the curve for that function etc.
Do it again for the limit y>0.
And again for the curve where z(u,v) intersects the u-v plane.
Once you have all those lines, you should be able to shade in the region you are doing the integration over.
Now you should be back on familiar territory.
It may be that your choice for u and v transformations is not good for this - if so, you need to find another transformation that works better.