【正文】 a(81:146,10)。 ti=0:10:28654。 tu=0:10:18397。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %3區(qū) Ni分布圖 x=data(81:146,2)。 y=data(81:146,3)。 z=data(81:146,11)。 ti=0:10:28654。 tu=0:10:18397。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %3區(qū) Pb分布圖 x=data(81:146,2)。 y=data(81:146,3)。 z=data(81:146,12)。 ti=0:10:28654。 35 tu=0:10:18397。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %3區(qū) Zn分布圖 x=data(81:146,2)。 y=data(81:146,3)。 z=data(81:146,13)。 ti=0:10:28654。 tu=0:10:18397。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %4區(qū) As分布圖 x=data(147:284,2)。 y=data(147:284,3)。 z=data(147:284,6)。 ti=0:10:26424。 tu=0:10:17949。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %4區(qū) Cd分布圖 x=data(147:284,2)。 y=data(147:284,3)。 z=data(147:284,7)。 ti=0:10:26424。 tu=0:10:17949。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) 36 %4區(qū) Cr分布圖 x=data(147:284,2)。 y=data(147:284,3)。 z=data(147:284,8)。 ti=0:10:26424。 tu=0:10:17949。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %4區(qū) Cu分布圖 x=data(147:284,2)。 y=data(147:284,3)。 z=data(147:284,9)。 ti=0:10:26424。 tu=0:10:17949。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %4區(qū) Hg分布圖 x=data(147:284,2)。 y=data(147:284,3)。 z=data(147:284,10)。 ti=0:10:26424。 tu=0:10:17949。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %4區(qū) Ni分布圖 x=data(147:284,2)。 y=data(147:284,3)。 z=data(147:284,11)。 37 ti=0:10:26424。 tu=0:10:17949。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %4區(qū) Pb分布圖 x=data(147:284,2)。 y=data(147:284,3)。 z=data(147:284,12)。 ti=0:10:26424。 tu=0:10:17949。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %4區(qū) Zn分布圖 x=data(147:284,2)。 y=data(147:284,3)。 z=data(147:284,13)。 ti=0:10:26424。 tu=0:10:17949。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %5區(qū) As分布圖 x=data(285:319,2)。 y=data(285:319,3)。 z=data(285:319,6)。 ti=0:10:20582。 tu=0:10:17133。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) 38 hold on plot3(x,y,z,39。o39。) %5區(qū) Cd分布圖 x=data(285:319,2)。 y=data(285:319,3)。 z=data(285:319,7)。 ti=0:10:20582。 tu=0:10:17133。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %5區(qū) Cr分布圖 x=data(285:319,2)。 y=data(285:319,3)。 z=data(285:319,8)。 ti=0:10:20582。 tu=0:10:17133。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %5區(qū) Cu分布圖 x=data(285:319,2)。 y=data(285:319,3)。 z=data(285:319,9)。 ti=0:10:20582。 tu=0:10:17133。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %5區(qū) Hg分布圖 39 x=data(285:319,2)。 y=data(285:319,3)。 z=data(285:319,10)。 ti=0:10:20582。 tu=0:10:17133。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %5區(qū) Ni分布圖 x=data(285:319,2)。 y=data(285:319,3)。 z=data(285:319,11)。 ti=0:10:20582。 tu=0:10:17133。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %5區(qū) Pb分布圖 x=data(285:319,2)。 y=data(285:319,3)。 z=data(285:319,12)。 ti=0:10:20582。 tu=0:10:17133。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %5區(qū) Zn分布圖 x=data(285:319,2)。 y=data(285:319,3)。 z=data(285:319,13)。 ti=0:10:20582。 tu=0:10:17133。 40 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。) %城區(qū)地形圖 x=data(1:319,2)。 y=data(1:319,3)。 z=data(1:319,4)。 ti=0:10:28654。 tu=0:10:18449。 [xi,yi]=meshgrid(ti,tu)。 zi=griddata(x,y,z,xi,yi,39。cubic39。)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。o39。