'pcolor' equavalent plotting method

Hi,

I am wondering if there could be a more convenient way to make a plot like 'pcolor' (flat shading) function in MATLAB/matplotlib, which makes a color plot for a non-gridded XYZ data without interpolating.

This type of plot is used a lot in our data analysis that the parameters (such as temperature, magnetic field) cannot be perfectly controlled with ideally same intervals. Given a case that a resistance is measured by varying temperature and magnetic field, although the 'set' parameter points are ideally gridded, while due to fluctuation and PID controls, the 'real measured'  parameter points are not gridded. In these measurements, the final results will be returned as three 2-D waves X(u,v), Y(u,v), Z(u,v), where u,v are the 'set point' of parameters, X, Y are the 'measured' values of parameters, and Z is the quantity of interest, and the goal is to visualize the Z(X, Y). 

I found two viable approach in current version of Igor pro, but both of them cannot provide a figure with satisfying quality as MATLAB or matplotlib provides.

1) 'Plotting 1D X, Y and Z Waves With Non-Gridded XY Data' - The interpolation is not preferred, with dense parameter points, 'flat shading' option in MATLAB can visualize data in a smooth and accurate representation. In addition, the data has to convert into 1D wave first, during which it lost the information of the 'mesh unit' and make the operating time longer. 

2) 'Parametric Surface Data Formats' for Gizmo figure - It provides a figure that is exactly same mechanism as 'pcolor' function do, however, the detailed modification is difficult in Gizmo panel. Especially the axis properties cannot be easily changed for X, Y axis, and it frequently outputs the axis with unclear fonts or labels.

In all, I believe it will be greatly helpful with a separate function that visualize Gizmo figure with 'Parametric Surface Data Formats' input and fixed top view direction, but with more handy axis help as in image plot panels.