现在想编个程序，给定了x值、y值以及z值，现在现在想以x，y为坐标，得到一个关于z值的等高曲线图，二维、三维都可以。

x、y和z值之间没有函数关系。如下列数据：

x          y           z
-71       56  0
-70 5 1
74 5 1
75 57 9
0.5 18 15
0.5 53 9
-40 11 -16
-40 45 -9
-70 32 -4
35 11 10
39 45 5
-29 18 -5
-52 37 -3
52 25 6
-18 53 -5
15 15 10
-57 9 3
65 41 2
-56 52 3
15 5 9
19 33 10
-19 6 -6
-45 23 -12
-63 17 -2
57 17 12
matlab程序如下：d = [-71       56  0
-70 5 1
74 5 1
75 57 9
0.5 18 15
0.5 53 9
-40 11 -16
-40 45 -9
-70 32 -4
35 11 10
39 45 5
-29 18 -5
-52 37 -3
52 25 6
-18 53 -5
15 15 10
-57 9 3
65 41 2
-56 52 3
15 5 9
19 33 10
-19 6 -6
-45 23 -12
-63 17 -2
57 17 12];
[X Y] = meshgrid(-80:80, -80:80);
Z = griddata(d(:,1), d(:,2), d(:,3), X, Y);
C=contour(X, Y, Z)；clabel(C);
请问相同程序在R中如何实现