爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 270622|回复: 320

[源代码] 利用scipy和sympy扩大三角网格形成的轮廓

  [复制链接]

新浪微博达人勋

发表于 2015-11-5 16:11:59 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
本帖最后由 非对称 于 2015-11-10 10:29 编辑


问题的引出:当我们要对一组离散点利用 构建三角网格算法 进行插值后绘制色斑图,会碰到一个问题,它插值后在构建的三角形网之外的格点的值都为null,如果我们直接绘图,然后利用边界切图,将会碰到三角网轮廓到切图边界之间没有任何填色的情况,如下图。
2015110515_01_mark.png 2015111010_01.png


问题的解决办法:
1、首先利用scipy.spatial.ConvexHull函数求得所有离散点的凸多边形轮廓。
figure_1.png
2、然后以构成这个轮廓的每条线段为一条直角边构建一个直角三角形,这些直角三角形的另外所求的那一点,全部连接起来即可扩大轮廓,然后利用线性插值,求得这些扩充点的值(或者直接用直角三角形另外2点的任意一点值)。
a.png
3、利用扩充后的点集来重新构成包含切图区域的离散点集。再插值绘图,即可得到完美的图像。
2015110515_01.png 2015111010_01.png

具体代码如下:见附件

游客,如果您要查看本帖隐藏内容请回复


评分

参与人数 5金钱 +57 贡献 +12 收起 理由
道狼52203 + 2 + 1 很给力!
云知道 + 10 + 2 很给力!
mofangbao + 15 + 5
风子 + 20 + 2 感谢分享
lovechang1314 + 10 + 2 很给力!

查看全部评分

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-11-5 18:11:41 | 显示全部楼层
nice啊,大神
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-11-5 20:39:54 | 显示全部楼层
又见大神。。。
{:eb502:}{:eb502:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2015-11-5 22:46:29 | 显示全部楼层
{:5_201:}{:5_201:}太厉害了
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2015-11-6 15:13:58 | 显示全部楼层
谢谢分享,学习一下
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-11-7 09:50:34 | 显示全部楼层
好厉害啊  哪位大神  用MATLAB  或者 它的  m_map  做出这种效果图的 教我下呗   求指教
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-11-7 10:38:25 | 显示全部楼层
太牛了太牛了
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-11-8 12:31:57 | 显示全部楼层
{:eb502:}
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

新浪微博达人勋

发表于 2015-11-8 20:42:58 | 显示全部楼层
谢谢分享 楼主
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2015-11-9 15:36:21 | 显示全部楼层
看看你】~~~~~~~~~~~~~~~
密码修改失败请联系微信:mofangbao
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表