I'm working through converting the ImageJ/Fiji script templates to use SciJava Ops. The Richardson-Lucy Total Variation deconvolution imagej-ops Groovy script uses the equation function in the image namespace to create a grid image:

grid_img = ops.run("create.img", [x_dim, y_dim], new FloatType())
formula = "p[0]^2 * p[1]"
ops.image().equation(grid_img, formula)

If I remember correctly the equation function depends on JavaScript (Nashorn) which may complicate things. I'm not necessarily advocating that we add this to SciJava Ops, and as I understand it this is not a highly used function.