Is there another way to plot 3d intensity levels in an image? Using matrixplot takes a long time.
restart;gc():
with(ImageTools):
with(plots):
a:=Read("c:/lighthouse1.jpg"):
b:=convert(ToGrayscale(a),Array):
c:=convert(b,Matrix):
matrixplot(c,style=point,shading=zhue)

martixplot(c,style=point,shading=zhue,orientation=[0,0])