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