Increase the dispersion range of the points plotted within the sphere
3 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Nnamdi Chukwunenye
le 31 Août 2020
Modifié(e) : Bruno Luong
le 31 Août 2020
I am working with the FAQ Solution to 'How do I create a set of random locations on the surface of a sphere?' and I am trying to increase the dispersion of the points within the hemisphere I have setup. Currently, the points are a bit too clustered around the center of the sphere and I haven't been able to get more points farther away from the center. I know I have to edit line 55 r = bsxfun(@rdivide, r, sqrt(sum(r.^2,1))); But I haven't been able to change it properly to get the result Im looking for. Here is what I have below any help would be greatly appreciated.
clear
clc
% Create Hemisphere Domain
[x_dom,y_dom,z_dom] = sphere(80); %Create Sphere
x_dom = x_dom(41:end,:); % Keep top 41 x points
y_dom = y_dom(41:end,:); % Keep top 41 y points
z_dom = z_dom(41:end,:); % Keep top 41 z points
hemisphere_radius = 80;
figure();
Hemi_sf = surf(hemisphere_radius.*x_dom,hemisphere_radius.*y_dom,hemisphere_radius.*z_dom, 'FaceColor','#4DBEEE','EdgeColor', 'none');
alpha 0.2
x_ax_lab = xlabel('x axis', 'Color', '#4DBEEE');
y_ax_lab = ylabel('y axis', 'Color', '#4DBEEE');
z_ax_lab = zlabel('z axis', 'Color', '#4DBEEE');
% Plot Outerboundary Circular Plane
x_c = 0;
y_c = 0;
z_c = 0;
radii_plane = 80;
radii_vein = 4;
radii_node = 4;
center_plane = [x_c, y_c]; % center point of circular plane
viscircles(center_plane, radii_plane, 'color', '#77AC30');
hold on
%%
% Place Starting Center Sphere
[x_c,y_c,z_c] = sphere();
x_c = x_c(11:end,:); % Keep top 11 x points
y_c = y_c(11:end,:); % Keep top 11 y points
z_c = z_c(11:end,:); % Keep top 11 z points
Center_Root = surf(radii_vein.*x_c,radii_vein.*y_c,radii_vein.*z_c, 'FaceColor', 'k');
%%
% Create Set of Random Points to be Plotted inside Hemisphere
n = 200;
A = zeros(n,3);
radius = 76;
% Generate random points
r = randn(n,3);
r = bsxfun(@rdivide, r, sqrt(sum(r.^2,1)));
r = radius*r;
% Extract the x, y, and z coordinates from the array.
x = r(:,1); % Extract x from column #1.
y = r(:,2); % Extract y from column #2.
z = abs(r(:,3)); % Extract z from column #3.
A = [x, y, z]; % stores random points generated
BioM_3Dgraph1 = scatter3(x, y, z, 'filled', 'r');
4 commentaires
Bruno Luong
le 31 Août 2020
Modifié(e) : Bruno Luong
le 31 Août 2020
Yes you are summing along a wrong dimension. It should be along the dimension of the 3-coordinates, meaning dimension 2 in YOUR case.
Réponse acceptée
Bruno Luong
le 31 Août 2020
Modifié(e) : Bruno Luong
le 31 Août 2020
This code generates randomly n 3D points unifformly distributed inside the hemi sphere r < radius and z > 0.
radius = 76;
n = 2000;
s = randn(3,n); % PLEASE DO NOT TRANSPOSE
r = (rand(1,n)).^(1/3);
c = r./sqrt(sum(s.^2,1));
r = radius * (s .* c);
x = r(1,:);
y = r(2,:);
z = abs(r(3,:));
scatter3(x, y, z, 'filled', 'r')
axis equal
0 commentaires
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Surface and Mesh Plots dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!