Simple example how to make a pcolor plot with missing values (equal to NaN) in Julia

pcolor_NaN.jl

# pcolor plot with missing values (equal to NaN) in Julia
using PyCall
@pyimport numpy.ma as ma

# some random data
A = randn(10,10); 
A[2,2] = NaN;

# created a Masked Array
# see also 
# http://stackoverflow.com/questions/35709042/how-to-create-a-masked-array-using-numpy-ma-imported-by-pycall-in-julia/35710432
x = pycall(ma.array, Any, A, mask=isnan(A))
pcolor(x)