自己画个球体玩玩

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

1. library(rgl)  
   
2. library(MASS)  
   
3. library(evd)
   
4. N1=100;N2=100
   
5. A=seq(0,2*pi,length.out=N1)
   
6. B=seq(0,pi,length.out=N2)
   
7. r=2
   
8. x=outer(r*sin(B),cos(A),"*")
   
9. y=outer(r*sin(B),sin(A),"*")
   
10. z=r-cos(B)
    
11. #col=sample(rainbow(N2),N2)
    
12. plot3d(x,y,z,col=rainbow(100),axes=F,xlab="",ylab="",zlab="")

复制代码

想画个地球仪来着的

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：Library Rainbow length plot3D Sample rainbow library 地球仪

不错，如果有世界地图在上面就更好了。

1. lat <- matrix(seq(90,-90, len=50)*pi/180, 50, 50, byrow=TRUE)
   
2. long <- matrix(seq(-180, 180, len=50)*pi/180, 50, 50)
   
3. r <- 6378.1 # radius of Earth in km
   
4. x <- r*cos(lat)*cos(long)
   
5. y <- r*cos(lat)*sin(long)
   
6. z <- r*sin(lat)
   
7. open3d()
   
8. persp3d(x, y, z, col="white",
   
9.        texture=system.file("textures/worldsmall.png",package="rgl"),
   
10.        specular="black", axes=FALSE, box=FALSE, xlab="", ylab="", zlab="",
    
11.        normal_x=x, normal_y=y, normal_z=z)

复制代码

求证1加1 发表于 2013-11-8 21:59
不错，如果有世界地图在上面就更好了。

这个只能显示西半球么？如果能够全部显示的话，还能让它旋转哈

yywan0913 发表于 2013-11-8 22:13
这个只能显示西半球么？如果能够全部显示的话，还能让它旋转哈

全部都有的哈，转一转就有了

求证1加1 发表于 2013-11-8 22:15
全部都有的哈，转一转就有了

刚在你那基础上把X,Y角度repeat了一下，的确都能看到，可惜关不掉了，折腾了好久，被迫关机了。