Hi David,
The process of generating this figure is a bit complicated including retrieve large amount of time-series data from NSC DR1. However, there is a much simple way to have a similar figure. You can just load the data from nsc_dr1.object table (my figure shows the data for the SMC region with 3<RA<25 and -75.5<Dec<-70; note that I have also deblened the catalog by excluding any targets having neighbors within 3 arcsec, but I don't think it has any major influence on this figure), choose targets with (u,g,r,i,z)rms<10 to avoid the 999999.0 value, and plot the spatial distribution with z-axis color coded with rms (you may need to adjust the range of z-axis to make it more obvious, e.g., 0~0.1 or 0~0.2). I have checked my data and found the patterns in g,r,i,z bands.