LineIntegralConvolutionPlot[{{Cos[y] - Sin[x]^3, -.1 y - Sin[x]}, ExampleData[{"TestImage", "Lena"}]}, {x, -2, 4}, {y, -2, 4}, LineIntegralConvolutionScale -> 0.3, RasterSize -> 300, ImageSize -> Large]