如何情节上的经度和纬度地球不使用映射工具箱?
58视图(30天)
显示旧的评论
有一个函数显示世界的国家在全球的阴谋,我需要知道如何情节甚至增量的纬度和经度10度到全球阴谋没有使用geoshow命令映射工具箱中找到。
接受的答案
乍得格林
2017年8月2日
编辑:MathWorks支万博1manbetx持团队
2022年6月8日
* * * * * * *
你可以轻易与标准的matlab函数
sph2cart
。治疗的纬度和经度分别如下:
R = 6371;%在千米地球半径
latspacing = 10;
lonspacing = 20;
%行经度:
[lon1, lat1] = meshgrid (-180: lonspacing: 180年,linspace (-90、90300));
(x1, y1, z1] = sph2cart (lon1 *π/ 180,lat1 *π/ 180,R);
plot3 (x1, y1, z1,“- - -”,“颜色”,0.5 * (1 1 1))
持有在
%的纬度线:
[lat2, lon2] = meshgrid (-90: latspacing: 90年,linspace (-180180300);
(x2, y2、z2) = sph2cart (lon2 *π/ 180,lat2 *π/ 180,R);
plot3 (x2, y2、z2,“- - -”,“颜色”,0.5 * (1 1 1))
轴平等的紧张了
因为你可以看到在世界各地,也许你想把一个不透明的球体在世界各地。做到这样的。我让球0.99倍地球的大小来确保它不重叠的纬度和经度:
[X, Y, Z] =球(100);
冲浪(X * * 0, Y * * 0, Z。* R *,“facecolor”,' w ',“edgecolor”,“没有”)
C =负载(“borderdata.mat”);
为k = 1:246
[xtmp, ytmp ztmp] = sph2cart(函数(C.lon {k}),函数(C.lat {k}), R);
plot3 (xtmp ytmp ztmp,“k”)
结束
6个评论
山姆·波尔克
2021年6月10日
编辑:山姆·波尔克
2021年6月10日
%经度的情节:
[lon_lon, lat_lon] = meshgrid (-180: lonspacing: 180年,linspace (-90、90300));%的经度坐标坐标网格
[xlon, ylon zlon] = geodetic2ecef (wgs84Ellipsoid (“米”),lon_lon lat_lon 0);%皈依(x, y, z)在地球的表面
plot3 (xlon ylon zlon,“- - -”,“颜色”,0.5 * (1 1 1))
持有在
%纬度的情节:
[lat_lat, lon_lat] = meshgrid (-90: latspacing: 90年,linspace (-180180300);%坐标lattitude网格坐标
[xlat, ylat zlat] = geodetic2ecef (wgs84Ellipsoid (“米”),lon_lat lat_lat 0);%皈依(x, y, z)在地球的表面
plot3 (xlat ylat zlat,“- - -”,“颜色”,0.5 * (1 1 1))
我遇到了一些麻烦添加地图的情节,然而。让我知道如果你有任何建议!