一个IDL画等值线图问题

墨小代MD

主程序完全照着书上的例子z=randomu(-100L,50,50)for i=0,4 do z=smooth(z,15)
z=(z-min(z))*15000.0+100.0
x=findgen(50)-100.0
y=findgen(50)+10.0

map_set,35.0,-75.0,/mercator,scale=50e6
map_continents
map_grid

levels=[150,200,250,300,350,400,450,500]
c_labels=[0,1,0,1,0,1,0,1]
contour,z,x,y,levels=levels,c_labels=c_labels,/overplot

为啥输出来是这样的图？等值线呢？

墨小代MD

图不知怎么传不上来，大概就是一个地图投影，没有等值线，为啥？

fatcat916

for i=0,4 do z=smooth(z,15, /EDGE_TRUNCATE)
平滑要加边界处理参数。

fatcat916

此外，建议把
z=(z-min(z))*15000.0+100.0
改成
z=(z-min(z))/(max(z)-min(z))*500.0+100.0
看着要舒服一点。

墨小代MD

fatcat916 发表于 2015-5-22 01:06
for i=0,4 do z=smooth(z,15, /EDGE_TRUNCATE)
平滑要加边界处理参数。

的确是这样！谢谢！还要请教简单说下edge_mirror和edge_truncate区别在哪么？

fatcat916

墨小代MD 发表于 2015-5-22 10:08
的确是这样！谢谢！还要请教简单说下edge_mirror和edge_truncate区别在哪么？

字面上看，应该是平滑算法用到的超过边界的数据采用镜像（以边界为对称轴）还是复制边界（重复边界）的区别吧？

fatcat916

墨小代MD 发表于 2015-5-22 10:08
的确是这样！谢谢！还要请教简单说下edge_mirror和edge_truncate区别在哪么？

字面上看，应该是平滑算法用到的超过边界的数据是采用镜像（以边界为对称轴）还是复制边界（重复边界）的区别吧？