diff --git a/src/image.jl b/src/image.jl
index 8a84e93a78501b4fc1fffd8004d955220e78d494..e3629c4208f1d7824c02ec963b9feaa975009ca4 100644
--- a/src/image.jl
+++ b/src/image.jl
@@ -120,7 +120,8 @@ function _sample(f::FeFunction)
                 (x .- A[:, 1])
             ε = eps()
             if all(xloc .>= -ε) && sum(xloc) <= 1 + ε
-                out[fcolon..., I] .= evaluate(f, xloc)
+                # use view here to handle scalar data assigment as well
+                view(out, fcolon..., I) .= evaluate(f, xloc)
             end
         end
     end