1、给出一幅单通道图片。
2、为了节省计算时间，这里把图片缩小。
img = ImageResize[img, 30]
3、获取图片数据：
data=ImageData[img]
4、把每一个像素的坐标转化为点坐标，并把像素值转化为数据集，得到一个三元数据集：
shuju = Flatten[
Table[{m, n, data[[m, n]]}, {m, 1, Length[data], 1}, {n, 1,
Length[data[[1]]], 1}], 1]
5、画出数据的散点图：
sandiantu = ListPlot3D[shuju, PlotRange -> All]
6、用一种方案拟合数据：
nihe = Fit[shuju, {1, x, y, x y}, {x, y}]
并画出曲面图：
Plot3D[nihe, {x, 1, 30}, {y, 1, 30}]
与上面的散点图相差甚远。
7、改变拟合方案，增加复杂度：
nihe = Fit[
shuju, {x, y, x y, x^2, y^2, x^3, y^2, x y^2, x^2 y, x^4, y^4,
x y^3, x^3 y, x^2 y^2}, {x, y}]
看样子，复杂度还是不够。
8、引进正弦函数：
nihe = Fit[shuju,
Sin /@ {x, y, x y, x^2, y^2, x^3, y^2, x y^2, x^2 y}, {x, y}]
9、再增加复杂度：
nihe = Fit[shuju,
Flatten[{a0 =
Flatten[{1, Table[x^m y^n, {m, 1, 9, 1}, {n, 1, 9, 1}]}],
Sin /@ a0}], {x, y}]
10、出图效果。
11、用这个拟合出来的曲面，来还原原来的img：
Image[Table[nihe, {x, 1, 30, 1}, {y, 1, 30, 1}]]
把图片放大到500像素，结果得到下图，与原图相差太远。
