grdimage turns all black when region changes

[tjcrone]

I am using grdimage to plot data from a high-resolution NetCDF file. With smaller regions, the grid image shows just fine. With larger region sizes, the image turns completely black. Example here: https://github.com/ooicloud/axial-drilling/blob/master/maps/notebooks/grdimage_black.ipynb. Any thoughts on what could be going on? Suggestions on best ways to troubleshoot? Thank you!

[tjcrone]

I think I figured it out. Changing the region may lead to the incorporation of a bad data point, so we get this in the debug output: grdimage [INFORMATION]: Auto-stretching CPT file turbo to fit data range -1621.96 to inf. I should be able to fix this by either fixing the data file, or changing the range of the CPT.

[tjcrone]

Looks like the grid I have set masked areas to inf rather than NaN. I fixed the grid with:

$ gmt grdclip axial_auv_mbari09and11s_gcs_1m_v3.grd -Gaxial_auv_mbari09and11s_gcs_1m_v3_nan.grd -Srinf/NaN