mountain =
；hut =
；
indexFunc[s_, w_] := {Quotient[s, w, 1] 1, Mod[s, w, 1]};
edgeWeight[_[s_, t_], data_, w_] :=
Block[{sindex, tindex}, sindex = indexFunc[s, w];
tindex = indexFunc[t, w];
Norm[Extract[data, sindex]] Norm[Extract[data, tindex]]];
weightImageGraph[hut_, mountain_] :=
Block[{g, dimx, dimy, data, elist, vcount, eweights}, {dimx, dimy} =
ImageDimensions[mountain];
g = GridGraph[{dimx, dimy}];
elist = EdgeList[g];
vcount = VertexCount[g];
data = ImageData[hut][[;; dimy]] - ImageData[mountain];
eweights = edgeWeight[#, data, dimx] & /@ elist;
Graph[Join[Table[UndirectedEdge["s", i], {i, 1, dimx}], elist,
Table[UndirectedEdge[vcount - i, "t"], {i, 0, dimx - 1}]],
EdgeCapacity -> (Join[#, eweights, #] &[
Table[Total[eweights], {dimx}]])]];
graph = weightImageGraph[hut, mountain];
edgecut = FindEdgeCut[graph, "s", "t"];
synthesis[hut_, mountain_, edgecut_, graph_] :=
Block[{dimx, dimy, comp, newpart, mountaindata,
hutdata}, {dimx, dimy} = ImageDimensions[mountain];
mountaindata = ImageData[hut];
hutdata = ImageData[mountain];
comp = Select[ConnectedComponents[EdgeDelete[graph, edgecut]],
MemberQ[#, "s"] &] // First;
newpart = indexFunc[#, dimx] & /@ DeleteCases[comp, "s"];
ReplacePart[
mountaindata, ({#1, #2} -> hutdata[[#1, #2]]) & @@@ newpart]]
newdata = synthesis[hut, mountain, edgecut, graph];
Image[newdata, ColorSpace -> "RGB"]
--------By 官网（原谅小编不放链接，因为俺觉得不去官网混一圈的人，是很难理解Mathematica 与 Wolfram 本人是多么的伟大）
