If I have two spheres plotted (one inside the other), how can I rotate only the interior one?
3 次查看(过去 30 天)
显示 更早的评论
Hi guys, so I'm currently working on a code which generates a model of the earth, plots your position on the earth, and then a sphere of stars around it. I am able to make this plot, and I am also able to get the earth to rotate when plotted alone (using the rotate function). However, whenever I try to plot both and then rotate the earth either the outer sphere doesn't get plotted, or both are plotted but I get the following error:
Error using rotate (line 29)
H must contain axes children only.
I'm really not sure what this error means, and I don't know how to fix it... if you need to look at my code, here it is:
%%Request the user for GPS Coordinates:
Lat = input('Please input the Lattitude of your desired Location: ')
Long = input('Please input the Longitude of your desired Location: ')
%%Generate the sphere that will later become Earth
alt = 1;
skyDome = figure;
set(skyDome,'ToolBar','none','MenuBar','none','name','Model of Sky','numbertitle','off',...
'units','normalized','outerposition',[0 .05 1 .95]);
subplot(1,2,1);
sphere;
positionGlobe = plot3(0,0,1,'.','markersize',20,'color','red'); %Place point on earth's surface
[x y z] = sphere(128);
r = 5;
stars = surfl(r*x, r*y, r*z);
set(stars, 'FaceAlpha', .25)
shading interp
%%Import an image of the Earth and map it to the sphere
img = imread('Earth Texture.jpg');
[imgInd,map] = rgb2ind(img,256);
[imgIndRows,imgIndCols] = size(imgInd);
[X,Y,Z] = sphere(imgIndRows,imgIndCols);
surface(X,Y,Z,flipud(imgInd),...
'FaceColor','texturemap',...
'EdgeColor','none',...
'CDataMapping','direct');
colormap(map);
view(-35,45);
%%Reproduce Earth's axial tilt
earth = findobj(gca,'type','surface');
%Define the prime meridian/equator intersection
rotate(positionGlobe,[1 0 0],90); %Align with prime meridian
rotate(positionGlobe,[0 0 1],90); %Align with equator
%Input GPS coordinates:
rotate(positionGlobe,[0 1 0],-Lat);
rotate(positionGlobe,[0 0 1],Long);
rotate(earth,[0 1 0],23.6);
rotate(positionGlobe,[0 1 0],23.6); %Orient coordinates along axial North Pole
%%Produce a standard map to show position location
subplot(1,2,2);
imshow(img);
hold on
positionMap = plot(800+4.444444*(Long),400-4.4444444*(Lat),'.','markersize',20,'color','red'); %Display point on map
positionData = text(0, 820,['Latitude: ', num2str(Lat), char(176), ' | Longitude: ', num2str(Long), char(176),...
' | Time of Day: '],'Fontsize',12)
%%Spin the model of the earth
for i = 1:1000
hold on
deg = 1;
rotate(earth,[.22 0 .5],deg);
rotate(positionGlobe,[.22 0 .5],deg);
M(i) = getframe(gca);
set(positionData,'string',['Latitude: ', num2str(Lat), char(176), ' | Longitude: ', num2str(Long), char(176),...
' | Time of Day: ',num2str(((deg/360)*24)*i)])
end
Any help would be greatly appreciated!
0 个评论
采纳的回答
Walter Roberson
2016-5-28
You have
earth = findobj(gca,'type','surface');
That is risky, as the object might be part of a different axes. It would be safer to assign the result of the surface() call to earth and not do the findobj(). As it is now, there is the possibility that earth will come out as the empty graphics placeholder.
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Map Display 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!