编辑注意：此文件被选为Matlab Central本周精选
Polarplot3d为三维极坐标数据生成表面、网格、线框和轮廓图。标记极轴以固定的高度绘制，或以最大半径沿表面轮廓绘制。极网格也可以绘制在表面的顶部。
本功能基于J De Freitas的polar3d，文件交换ID 7656。
输入参数是一个由大小、Zp和属性、值对组成的矩阵，这些属性和值对可以修改默认绘图行为。Zp的每一列包含沿单一半子午线的信息，每一行给出沿圆弧的高度值。默认情况下，假定Zp沿每一列递增半径，沿每一行递增角度(逆时针)。默认图绘制在一个单位半径的完整圆上。
“RadialRange”和“AngularRange”属性可以用来指定数据绘制的上、下角度和径向值。利用角距离向量和径向距离向量的相对顺序，确定了Zp的行、列的角方向和径向方向感。或者，可以指定一个向量，给出每一行或一列的位置。
在用户指定位置，极轴可以放置在最大半径处，在绘图盒的顶部或底部处的最小，最大值或平均值，或者在用户指定的位置处，或者它可以遵循数据周边的表面。可以使用“TickSpacing”属性调整极轴刻度标记间距。
默认的极轴方向是零度沿+X轴，递增的角度是逆时针方向。“极地方向”属性可以用来改变这是一个罗盘风格的绘图，沿+Y轴零度和增加角度顺时针围绕极点。
默认的表面着色是根据Zp中的值。这可以通过提供一个与Zp大小相同的矩阵作为'ColorData'属性的值来改变。
可以指定缩放参数来将数据插值到更细或更粗的网格上。输出的笛卡尔数据以三个矩阵的形式返回。
屏幕截图中的示例图是用以下命令生成的。
[t，r] = meshgrid（linspace（0,2 * pi，361），linspace（-4,4,101））;
(x, y) = pol2cart (t, r);
P =山峰(x, y);%峰值函数在极坐标上
％绘制3D极性图
图（'颜色'，'白色'，'numbertitle'，'off'，'name'，'polarplot3d v4.3'）;
polarplot3d（p，'plottype'，'surfn'，'pargid'，{4 24}，'tickspacing'，8，...
“AngularRange”,[270]*π/ 180,“RadialRange”,(。8 4],……
'Radleabels'，3，'Radlabellocation'，{180'Max'}，'Radlabelcolor'，'Red'）;
%设置绘图属性
集(gca、“DataAspectRatio”,[1 1 10],“视图”,(-12年38),…
'Xlim'，[-4.5 . 4]， 'Xtick'，[-4 -2 0 2 4]，…
'Ylim'，[-4.5 .5]，'Ytick'，[-4 -2 0 2 4]);
标题(“polarplot3d示例”);
这个zip文件包含了polarplot3d函数和一个m文件，其中包含示例绘图。
引用作为
肯Garrard(2021)。三维极坐标图（//www.tatmou.com/matlabcentral/fileexchange/13200-3d-polar-plot），matlab中央文件汇交流。检索到。
评论和评级（58）
你也可以从以下列表中选择一个网站:
来自XL的矩阵导入由Polarplot3d视为非数字'，'Sheet'，'Sheet1');
例如:matrix = readtable('
转换成数字
ZP = Table2Array（矩阵）;
应该在那之后工作
ive从Excel导入了19x25矩阵，我想要绘图，并输入它作为ZP，但我继续收到错误，颜色矩阵必须是数字和与ZP相同的大小，我如何解决这个问题？
非常漂亮!已经想象了这么长时间的数字!非常感谢。
有没有办法使整个表面是一种颜色？我正在使用此函数作为叠加覆盖以覆盖一些空值。
干杯
伟大的功能。我发现它有点奇怪，假设扇区从给定的θ值顺时针绘制数据扇区。在许多应用中（例如来自辐射扇区扫描的数据），绘制的扇区应以给定的θ值居中，而不是在它的一侧。无论如何，这很容易通过在弧度的半个部门的偏移量偏移中来修复，但我以为我会提到它。可能很好地添加为新选项。
你能举一个不使用峰峰函数的例子吗?例如，一个三角函数的例子，比如画一个球?
我认为如何使用它在更密切地审核脚本后用它与导入的数据一起使用它。只需将数据导入二维矩阵，并在PolarPlot3D功能中使用它。伟大的工作！
嗨，肯，这太壮观了。如果使用导入的数据，如何调整该命令。例如，对于一个给定的样品，我每旋转5度(即从0到355度)就有一个吸收光谱。
非常有用的脚本。有这么多的例子，使我能够立即理解参数选择的后果。干得好，肯!
你好
我有三个矩阵
r = [0:10:100]
PHI = 0：45：315
z =[11×8]矩阵
z有某个值(dBm)
我可以使用这个polarplot3d函数吗？
我想在任何程度上显示z值和任何半径。
真诚
请包括图形梯度的可视化。例如，心脏具有平均电梯度。拥有它是非常有用的。例如，这里的数学http://math.stackexchange.com/questions/586848/how-to-obtain-the-gradient-in-polar-coordinates.
令人敬畏的脚本！
非常感谢。
有可能在z轴上放一个标签吗?
你好肯恩,
感谢您提供的优秀代码
我有一个问题，我试图代表我的3D极性图中的数据。
对于所有值5，半径为常数5
Theta =[30 60 75 90]，每个Theta对应的z是[0.35 0.5 0.6 0.75]。所以这些峰一定代表z值。
提前谢谢你
嗨,肯,
感谢您在这里的所有优秀工作！
我有一个问题。在下面的14 × 16的数据矩阵中:
Data_Right_Arm= [0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0.8 1 1 0
0 0 0 0 0 0 0 0.5 0 0 0 0 0.9 1 1 0
0.95 0.5 0.5 0 0.5 1 1 0.95 0 0 0 0 0.6 1 1 1
1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1
1 1 1 1 0.7 0.95 0.8 0.8 0 0 0 0 0.5 1 0.95 0.95
1 1 0 0 0 0 0 0.55 1 1 1 1 1
1 1 0.85 0.85 0.85 0.8 0.9 0.85 0 0 0 0 0.55 1 0.95 1
1 0.85 0.95 1 0.8 0.95 0.95 0.95 0 0 0 0 0 0.55 0.95 1 1
0.9 0.8 1 1 0.5 1 0.95 1 0 0 0 0 0.75 1 0.95 0.95
0.85 0.95 0.85 1 0.8 0.95 0.9 0.95 0 0 0 0 0 1 1 0.75
0.9 1 1 0.7 1 0.9 0.8 0 0 0 0 0 0 0 0 0 0.75];
每个列的哪个角度对应关系如下：
Angles_deg = [022.5 45 67.5 90 112.5 135 157.5 180 202.5 225 247.5 270 292.5 315 337.5]
我需要创建一个3D极坐标网格，其中的单元格包含表示从0到1的数据矩阵值的颜色。我尝试使用“polarplot3d”函数如下:
T2 = [11.25 371.25]*pi/180;
R2 = [0 1.5];
%%使用默认参数绘制
图（'颜色'，'白色'）;
polarplot3d (Data_Right_Arm radialrange, r2,…
'极坐标'，{15 16}，'tickspacing'，22.5， 'GridScale'， [100 100]， 'angularrange'，t2);
视图(76 [-18]);
caxis（[0,1.0]）;
colorbar(垂直的);
但我没能得到正确的结果。任何建议都将非常感谢。
最好的
阿基
凯利,
您可以尝试使用返回矩阵xi，yi，zi进行内插功能。然后找到与所需任何x，y点对应的z值。
一个例子，使用散点interpolant，
[t，r] = meshgrid（linspace（0,2 * pi，361），...
linspace (4101));
(x, y) = pol2cart (t, r);
P =山峰(x, y);
ANGR = [20 350] * PI / 180;
radR =[。8 4];
图;
[xi,咦,子]=…
polarplot3d（p，...
'Angularrange'，Angr，......
RadialRange, radR);
zFun =…
散点interpolant（xi（:)，yi（:)，zi（:)，...
“天然”);
图;
冲浪(x, y, zFun (x, y)、“线型”,“没有一个”);
优秀的剧本;非常感谢你。
我想在图的中心(由第一行或输入向量构成)，以及输入向量的第一列和最后一列之间进行插值。
你能告诉我如何继续这件事吗？
Seb,
我相信这是由于OpenGL错误。例如，
图(“渲染器”、“opengl”);
polarplot3d(峰(101),“cartorigin”,[10 ^ 6 10 ^ 8]);
是支离破碎。但是,
图(“渲染器”、“zbuffer”);
polarplot3d(峰(101),“cartorigin”,[10 ^ 6 10 ^ 8]);
看起来不错。画家也有效，但网格里没有绘制得很好。
真正的好脚本，几乎所有我需要的。然而，当试图将CartOrigin移动到投影的UTM坐标(即x上的10^5和y上的10^6)时，网格完全被打乱了。知道发生在哪里吗?
干杯!
S.
做得好 ！
这段代码非常有助于关于人类细胞内部纳米粒子扩散的工作非常有用......可能不是您在创建此代码时所想的
安东尼
很好!编写良好，易于使用。
祝贺
做得好！
我想知道如何处理这种情况，如果数据点的分布不是均匀的，比如随着径向距离的增加，数据点更加密集，所以我们使用单元格而不是矩阵来存储数据。
做得好！
我想知道如何处理这种情况，如果数据点的分布不均匀，例如，数据点沿径向更密集。
如何下载
DD，
第589行使用~(波浪号)字符作为不需要的函数返回值的占位符。在介绍此特性之前，您必须使用Matlab的一个版本。将~替换为一个未使用的变量名，例如'dontcare'。
嗨,肯,
当运行演示代码时，我得到以下错误消息:
？??错误：文件：polarplot3d.m线：589列：3
表达式或语句不正确——可能不平衡(，{，或[。
你能帮我吗？
吉姆,
轮廓功能将笛卡尔轴添加到图中。我将这些轴的“可见”属性设置为“偏光”在Polarplot3d中。打开它们的类型
集(gca,“可见”,“上”);
在用polarplot3d绘制数据之后。
如果只需要轴标签和一个情节标题就可以输入
集(get (gca,“包含”),“可见”,“上”);
设置（get（gca，'ylabel'），'可见'，'开'）;
设置（get（gca，'标题'），'可见'，'开）;
这是非常有用的。,,,但我有“轮廓”的问题....当我运行与实例“surfn” canged为“轮廓”全部格式化，坐标轴标贴，标题都在前看不见！....什么是错误的....有人可以发给我一个如何使用轮廓选项的示例吗？
@dimaro - 我在这里找到了“Parse_pv_pairs.m”的版本：
http://www.kemt.fei.tuke.sk/predmety/kemt515_mt/_materialy/signal/signal/private/parse_pv_pairs.m.
话虽如此，我还没有机会尝试这个脚本(即3D Polar情节)，所以我不知道它是否可行
我找不到parse_pv_pairs函数？有人可以帮帮我吗？
谢谢。
这真是个好工作。救我很多时间。谢谢。
你好,
我正在处理地理数据，所以我的零方向(北方)是向上的，而度数以时钟的方式增加(所以与你的脚本相反)。
我怎样才能修改你的程序才能使用它呢?
提前谢谢!
卡米尔。
有不错的情节。这正是我制作钢管人像所需要的。谢谢!
很棒的工作! !我帮了我很多!谢谢!
书中,
您需要下载函数parse_pv_pabes（fex id 9082）。polarplot3d使用它来处理可选的输入参数。
(我使用Matlab b)版本，我不能运行这个函数(甚至在代码中给出的例子)。我得到以下错误:??使用==> polarplot3d时出错
解析变量列表出错
？??对于“struct”类型的输入参数，未定义函数或方法“parse_pv_pairs”。
谁能告诉我如何解决这个问题?
提交更新
谢谢你的评论，约翰。我也喜欢'属性'，值对时，参数列表很长，并已更新我的提交使用parse_pv_pairs。想象一下如果没有varmargin和p,v对，图形函数的句柄会是什么样子。
好的文档，许多例子，充满了特性。
如果我想要抱怨什么，我更喜欢在带有许多可选参数的代码中使用属性/值对，这些参数中的大多数将采用它们的默认值。