julia中等面积变换系统的等值线

假设我有一个系统，在2d系统a，b，c中有3个粒子，位置如下

a = [0.1,0.0]
b = [-0.1,0.0]
c = [0.0,0.1]

每个粒子由高斯函数描述，因此系统函数如下：

σ = 0.08
f(x,y) = exp(-((x-a[1])^2+(y-a[2])^2)/(2(σ^2)))+exp(-((x-b[1])^2+(y-b[2])^2)/(2(σ^2)))+exp(-((x-c[1])^2+(y-c[2])^2)/(2(σ^2)))

我画出了系统的初始形状

x = range(-0.5,stop=0.5,length=50)
y = range(-0.5,stop=0.5,length=50)
contour(x,y,f,levels=[1.0],aspect_ratio=:equal)

使用Plots包。一段时间后，粒子移动到其他位置，例如：

a = [0.1,-0.05]
b = [-0.1,0.0]
c = [0.05,0.1]

它根据相同的f（x，y）函数（具有相同的sigma）改变轮廓的形状。 我需要生成一个闭合轮廓，该轮廓对于新粒子位置保持相同的面积。 我该怎么做？或者我怎样才能找到达到相同区域的“水平”？ 这是一个例子，事实上我的函数中有很多粒子