The data now contains negative values. This occurs regardless of which random gridding method I use. Am I doing something wrong? Or is it a bug?
I know a similar problem was discussed here: https://www.originlab.com/forum/topic.asp?TOPIC_ID=9859, but I need a heat map with a spectrum bar for the z-values and it doesn't seem like that can be done with the contour plots.
Gridding routines typically make use of cubic polynomials to maintain smooth continuance in the interpolated curves and their first derivatives. If there are large gaps, the powerful cubic term can drive the interpolated results far from the data's neighborhood.
You could try using linear interpolation if you do not require smooth first-derivative curves.
Or, based on the large range of your data, I would recommend gridding in the log domain. This will guarantee no negative values.