matlab - 将图形转换为 2D 图表
问题描述
我想创建一个如下所示的二维图,
上面的图像是Graph
在 MATLAB ( ref )中创建的
s = [1 1 1 1 2 2 3 4 4 5 6];
t = [2 3 4 5 3 6 6 5 7 7 7];
weights = [50 10 20 80 90 90 30 20 100 40 60];
G = graph(s,t,weights)
plot(G,'EdgeLabel',G.Edges.Weight)
信息存储为 a 的节点、边、边权重graph
。我想使用此信息创建 2D CAD 绘图。可以使用边权重指定线的长度。但是,我不确定如何从图中检索角度。据我了解,边缘的方向因选择用于创建图形对象的布局而异。我想创建一个 [x,y] 坐标文件并导入 Autocad。
编辑:从下面解释的答案中,我知道将边缘权重分配为长度并不简单。作为替代方案,我想从图像中获取节点的坐标,计算节点之间的距离并将距离分配为边缘权重(忽略上面提供的权重)。使用一组坐标、节点-节点连接和节点-节点距离,我想以编程方式生成一维 CAD 图。
EDIT2:由于节点的坐标不能直接从 MATLAB 输出中获取,并且边权重(在原始输入中)不能指定为边长,我想尝试另一种方法。例如,如果这些是节点的坐标((75 25) (115 45) (90 60) (10 5) (45 0) (45 55) (0 25))
,我想计算坐标之间的欧几里德距离并将距离分配为边权重。据我了解,dimension
AutoCAD 中的选项卡计算欧几里得距离。但是,我不确定如何将此输出分配为边权重。
任何有关如何进行的建议将不胜感激。
解决方案
首先,对于您的特定示例,无法生成边缘权重为线长的图形。
例如,如果节点之间的距离1
,2
和3
是每个数组:
1 → 2 = 50
1 → 3 = 10
那么距离2 → 3
必须在40
&之间60
,否则三角形不存在。而您的数组将此距离指定为90
.
为了直观地证明这一点,如果您要描绘一条长度50
跨越节点1
&的线2
,如果您在这条线的任一端构造圆,半径等于节点之间的距离1 → 3
和2 → 3
,那么这些圆必须相交才能存在三角形。
目前,使用您指定的权重,没有这样的交集:
因此,假设每个节点具有任意位置,并将节点坐标作为参数提供给函数,您可以使用如下 AutoLISP 示例之类的函数构造所需的图形:
(defun graph ( pts sls tls wgt )
( (lambda ( l )
(foreach x l (text (cdr x) (itoa (car x)) 0.0 1))
(mapcar
'(lambda ( a b c / p q r )
(setq p (cdr (assoc a l))
q (cdr (assoc b l))
r (angle p q)
)
(entmake (list '(0 . "LINE") (cons 10 p) (cons 11 q) '(62 . 8)))
(text
(mapcar '(lambda ( x y ) (/ (+ x y) 2.0)) p q)
(itoa c)
(if (and (< (* pi 0.5) r) (<= r (* pi 1.5))) (+ r pi) r)
2
)
)
sls tls wgt
)
)
(mapcar 'cons (vl-sort (append sls tls) '<) pts)
)
)
(defun text ( p s a c )
(entmake
(list
'(0 . "TEXT")
(cons 10 p)
(cons 11 p)
(cons 50 a)
(cons 01 s)
(cons 62 c)
'(40 . 2)
'(72 . 1)
'(73 . 2)
)
)
)
当使用以下参数评估上述函数时(其中第一个参数表示七个节点的坐标):
(graph
'((75 25) (115 45) (90 60) (10 5) (45 0) (45 55) (0 25))
'( 1 1 1 1 2 2 3 4 4 5 6)
'( 2 3 4 5 3 6 6 5 7 7 7)
'(50 10 20 80 90 90 30 20 100 40 60)
)
它将在 AutoCAD 中产生以下结果:
但是,如果您希望权重由每个提供的节点坐标之间的 2D 距离确定,则可能需要考虑以下 AutoLISP 函数:
(defun graph ( pts sls tls )
( (lambda ( l )
(foreach x l (text (cdr x) (itoa (car x)) 0.0 1))
(mapcar
'(lambda ( a b / p q r )
(setq p (cdr (assoc a l))
q (cdr (assoc b l))
r (angle p q)
)
(entmake (list '(0 . "LINE") (cons 10 p) (cons 11 q) '(62 . 8)))
(text
(mapcar '(lambda ( x y ) (/ (+ x y) 2.0)) p q)
(rtos (distance p q) 2)
(if (and (< (* pi 0.5) r) (<= r (* pi 1.5))) (+ r pi) r)
2
)
)
sls tls
)
)
(mapcar 'cons (vl-sort (append sls tls) '<) pts)
)
)
(defun text ( p s a c )
(entmake
(list
'(0 . "TEXT")
(cons 10 p)
(cons 11 p)
(cons 50 a)
(cons 01 s)
(cons 62 c)
'(40 . 2)
'(72 . 1)
'(73 . 2)
)
)
)
提供节点坐标列表和两个连接节点列表:
(graph
'((75 25) (115 45) (90 60) (10 5) (45 0) (45 55) (0 25))
'(1 1 1 1 2 2 3 4 4 5 6)
'(2 3 4 5 3 6 6 5 7 7 7)
)
此函数将产生以下结果:
在这里,权重的准确性将由LUPREC
AutoCAD 中系统变量的值(4
在上例中设置为)确定。您也可以通过在我的代码中为函数提供精度参数来覆盖它rtos
,例如,对于3
小数位精度,表达式为:
(rtos (distance p q) 2 3)
推荐阅读
- solr - Solr 中的原始查询参数被忽略
- javascript - 如何通过多个页面的异步获取请求将数据存储在数组中?
- sharepoint-online - 组合框选定的值不保留
- php - 使用数组数据格式化 Excel Laravel
- python - 为什么这个 SQL INSERT INTO 语句会抛出错误:TypeError: not all arguments convert during string formatting
- ios - currentContext 表示样式改变子视图控制器视图的大小
- mysql - 无法使用 java.util.UUID 从主键为二进制 (128) 的 Mysql db 检索数据
- firebase - Firebase 即时用户禁用?
- calayer - 通过 AVAssetExportSession 导出时,CABasicAnimation 不更新自定义 CALayer 属性
- scala - Scala Recursive For Comprehension 只在空列表前面添加一次,为什么?