- 积分
- 8082
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-12-27
- 最后登录
- 1970-1-1
|
发表于 2020-11-28 20:24:53
|
显示全部楼层
按照你的思路画了立体包络图,但是效果出来怎么不对呀?数据应该没问题,因为我画slice图效果是没问题的
楼主能帮我看看哪里不对么?
立体包络图效果图:
垂直剖面图效果:
主要代码如下:- Point=rdBz_xyz(:,1:3);%--x,y,z坐标
- V=rdBz_xyz(:,4);%--反射率
- F = scatteredInterpolant(Point,V);
- % F = scatteredInterpolant(P,v) 以数组形式指定样本点坐标。P 的行包含 v 中值的 (x, y) 或 (x, y, z) 坐标。
- xx=-100:1:100;
- yy=-100:1:100;
- zz=0:0.5:15;
- [xq,yq,zq] = meshgrid(xx,yy,zz);
- vq = F(xq,yq,zq);
- vq=smooth3(vq);
- fw=5;
- % --立体包络数值fw
- set(h1,'name','dBz定值包络立体图','MenuBar','none','ToolBar','none');
- title(['最外层表面的值为: ' , num2str(fw)]);
- fv=isosurface(xq,yq,zq,vq,fw);
- hiso=patch(fv);
- set(hiso,'facecolor','b','edgecolor','none');
- patch(isocaps(xq,yq,zq,vq,fw), 'FaceColor', 'interp', 'EdgeColor', 'none');
复制代码
|
|