- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
亲爱的 stackoverflow 用户,
几年前,我用了几个月的 Mathematica。在几年没有编程之后,我现在做一个研究项目,作为一名学生,我在其中使用 Matlab。我在 stackoverflow 上找到了很多很好的帮助,但现在我遇到了以下问题:
我有一个矩形网格上节点之间连接的数据集,每个节点都可能与其 8 个邻居有连接。我的测量采用 3 x n 矩阵的形式,其中前两个值指定一个节点,第三个值指定它们是否连接,网格的大小是预先确定的。通常有大约十条线路来自至少彼此相邻的两个或三个节点。我的研究项目的目标是计算这组线周围距离 r 处的面积。
到目前为止,我已经能够使用下面的代码绘制线条,为此我使用了 stackoverflow 上的一些代码,这非常有用。但是我无法在一定距离处获得围绕它的等高线(我希望用它来计算该等高线内的面积)。 gplot 函数返回两个向量,每行有两个坐标,我发现很难将其转换为更有用的东西。我尝试在离线一定距离处定义一个值 Z,随着离线的距离而下降,所以我从这些线得到一个斜率。从这个斜率我可以计算等高线。但是,因为这些线只是坐标,所以我不知道如何计算到那条线的距离,而不是它们何时会成为函数。
我真的很茫然。我希望我在这里清楚地发布了我的问题。这是我第二次发布这个问题,我现在已经在代码和图片中添加了注释以更好地解释自己。感谢您提供的任何建议!
到目前为止,xls 文件是我上面提到的 3 x n 矩阵,我还在下面的代码中以矩阵形式编写了它的内容,所以我的问题更容易理解:
%# here i set my data/constants
filename='random.xls'; file=xlsread(filename); y=width; x=length;
%# random.xls looks approximately like this, after xlsread(filename) you get
file=[21 22 1;
21 20 1;
15 16 1;
15 14 1;
15 23 1;
14 22 1;
14 21 1;
22 15 1;
23 14 1;
24 15 1;
6 15 1;
5 14 1;
7 14 1;
8 15 1];
%# predefined width and length, i usually get this from the file
width=8; length=4;
%# here i create my adjaceny matrix in a elegant way user amro posted on stackoverflow
%# however i immediately multiply it by 0, creating a y*x by y*x matrix with all zeroes
[X Y] = meshgrid(1:x,1:y); X = X(:); Y = Y(:);
adjacency = squareform( pdist([X Y], 'chebychev') == 1 ); adjacency=adjacency*0;
%# here i take the matrix "file" for which the first two values are node numbers and
%# the third value designates whether there is a connection between the two nodes to
%# fill in the connections in the adjacencymatrix
[nrows,ncols]=size(file);
for r = 1:nrows
if file(r,3)==1
adjacency(file(r,1),file(r,2))=1;
end
end
adjacency=(adjacency+adjacency.');
%# plots the adjacencymatrix
subplot(121), spy(adjacency)
%# plots the connections and designates the nodes, note that the numbers designating
%# the nodes do not match original data, this is a separate problem i have not solved
[xx yy] = gplot(adjacency, [X Y]);
subplot(122), plot(xx, yy, 'ks-', 'MarkerFaceColor','b')
%# these last lines of code for plotting the numbers of the grid i do not fully
%# understand, in here is the cause for the numbers not matching the original data
axis([0 x+1 0 y+1])
[X Y] = meshgrid(1:x,1:y);
X = reshape(X',[],1) + 0.1; Y = reshape(Y',[],1) + 0.1;
text(X, Y(end:-1:1), cellstr(num2str((1:x*y)')) )
xlabel('length')
ylabel('width')
title(filename)
为了澄清我的问题,我添加了这两张图片:当前地 block http://imgur.com/5uPd4我想知道的区域http://imgur.com/WsIbg
最佳答案
从 Matlab 中的线集合中找到距离为 r 的等值线或等高线内表面积的解决方案,通过图形处理(扩张)进行近似,而不是精确或有效的 awnser!
我做了一个近似值,所以这不是一个精确的 awnser,也不是有效的编码。我的一位研究机器视觉的 friend 建议将线条转换为像素,然后用圆盘对图像进行膨胀,之后像素数就是表面积的度量:
%# constants and variables
minx = min(xx);
miny = min(yy);
maxx = max(xx);
maxy = max(yy);
rangex = maxx - minx;
rangey = maxy - miny;
borderRelNum = sqrt(2);
electrodeToImageScaleFactor = 100;
imsizex = 2*(maxx+borderRelNum)*electrodeToImageScaleFactor+2;
imsizey = 2*(maxy+borderRelNum)*electrodeToImageScaleFactor+2;
im = zeros(imsizex, imsizey);
grayscalevalue = 255;
disksize = round(borderRelNum*electrodeToImageScaleFactor);
%# transformation matrices
centerElectrodeSpace = [1, 0, -(minx + maxx) / 2;
0, 1, -(miny + maxy) / 2;
0, 0, 1 ];
scaleElectrodeToImage = [electrodeToImageScaleFactor , 0, 0;
0, electrodeToImageScaleFactor , 0;
0, 0, 1 ];
centerImageSpace = [ 1, 0, imsizex / 2;
0, 1, imsizey / 2;
0, 0, 1 ];
electrodeToImage = centerImageSpace * scaleElectrodeToImage * centerElectrodeSpace;
%# transformation
for i = 0:(size(xx,1) / 3 - 1)
p1 = [xx(i*3 + 1); yy(i*3 + 1); 1];
p2 = [xx(i*3 + 2); yy(i*3 + 2); 1];
p1im = electrodeToImage * p1
p2im = electrodeToImage * p2
lx = linspace( min( p1im(1), p2im(1) ), max( p1im(1), p2im(1) ), borderRelNum*electrodeToImageScaleFactor )
ly = linspace( min( p1im(2), p2im(2) ), max( p1im(2), p2im(2) ), borderRelNum*electrodeToImageScaleFactor )
index = sub2ind(size(im),round(lx),round(ly));
im(index) = grayscalevalue;
end
%# Now dilate and count pixels
se = strel('disk', disksize, 0);
im = imdilate(im, se);
image(im)
colormap(gray)
sum(sum(im/grayscalevalue))*(1/electrodeToImageScaleFactor^2)
如果有人能够更优雅、高效或更精确地解决我的问题,我仍然会非常感激。但现在就可以了。
-edit- 好吧,这确实非常低效,我的电脑已经在我的数据集(10 个 xls 文件,不多)上处理数字 30 分钟了,现在仍然在文件 1,如果我查看工作空间值,似乎
关于arrays - 距Matlab中的线集合一定距离的表面积,类似于轮廓?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/12138160/
A是不同元素的序列,B是A的子序列,A-B是A中的所有元素,但不是B中的所有元素距离(A) = 总和|a(i)-a(i+1)|从 i=1 到 n-1找到一个子序列 B 使得 Dist(B)+Dist(
我想通过计算每对中所有(多维)点集之间距离的平均值来量化组相似性。 我可以很容易地手动为每对组手动完成此操作,如下所示: library(dplyr) library(tibble) library(
在 OpenXML 中用于指定大小或 X、Y 坐标的度量单位是什么? (介绍)。 将那些与像素匹配是否有意义,如果是这样,那些如何转换为像素? graphicFrame.Transform = new
我想知道是否有人可以帮助我替换过渡层中的值。 如果我尝试: transitionlayer[transitionlayer >= 0.14] = 0.14 : comparison (5) is
我在 firebase 中有一个列表,其中包括地理位置(经度和纬度),并且我想获得距给定坐标最近的 10 个位置。 我正在从 MySQL 过渡,在那里我将计算 SELECT 中的距离, 并在 ORDE
如何在 Python 中根据 2 个 GPS 坐标计算速度、距离和方向(度)?每个点都有纬度、经度和时间。 我在这篇文章中找到了半正矢距离计算: Calculate distance between
关闭。这个问题需要多问focused 。目前不接受答案。 想要改进此问题吗?更新问题,使其仅关注一个问题 editing this post . 已关闭 6 年前。 Improve this ques
我只想使用 matplotlib 标记两条曲线之间发生最大偏差的位置。请帮助我。 垂直距离适用于 Kolmogorov–Smirnov test import numpy as np %matplot
我有一个包含数万行重复项的文件。我想根据行号找到重复项之间的平均时间/距离。 例如:(其中第一列是行号) 1 string1 2 string2 3 string2 4 string1 5 strin
用公式speed=distance/time计算时间 但时间总是0我的输入是 distance=10 和 speed=5 我的输出必须 = 2 #include int main() { in
我正在使用 Levenshtein 算法来查找两个字符串之间的相似性。这是我正在制作的程序的一个非常重要的部分,因此它需要有效。问题是该算法没有发现以下示例相似: CONAIR AIRCON 算法给出
对于一个房地产网站,我需要实现一个允许搜索文本和距离的搜索机制。 当 lat 和 lon 记录在单独的列中时,在 MySQL 表上进行距离计算很容易,但房子往往有 LOT true/false 属性。
是否可以在触发前更改 UIPanGestureRecognizer 的距离?目前的实现似乎在触发前有 5-10 像素的距离余量,我想降低它如果可能的话。 原因是我将 UIPanGestureRecog
我试图找到两个网格之间的偏差。例如在 3d 空间中定义的两组点之间的差异,我计划使用一些 3d 可视化工具来可视化距离,例如QT3d 或一些基于开放式 gl 的库。 我有两组网格,基本上是两个 .ST
所以,我有这个函数可以快速返回两个字符串之间的 Levenshtein 距离: Function Levenshtein(ByVal string1 As String, ByVal string2
我正在尝试用字典创建一个光学字符识别系统。 事实上,我还没有实现字典=) 我听说有一些基于 Levenstein 距离的简单指标,这些指标考虑了不同符号之间的不同距离。例如。 'N' 和 'H' 彼此
我在PostGIS数据库(-4326)中使用经纬度/经度SRID。我想以一种有效的方式找到最接近给定点的点。我试图做一个 ORDER BY ST_Distance(point, ST_GeomF
我想从线串的一端开始提取沿线串已知距离处的点的坐标。 例如: library(sf) path % group_by(L1) %>% summarise(do_union =
我已经编写了这些用于聚类基于序列的数据的函数: library(TraMineR) library(cluster) clustering <- function(data){ data <- s
是否可以设置 UILabel 的行之间的距离,因为我有一个 UILabel 包含 3 行,并且换行模式是自动换行? 最佳答案 如果您指的是“前导”,它指的是类型行之间的间隙 - 您无法在 UILabe
我是一名优秀的程序员,十分优秀!