如何情节上的经度和纬度地球不使用映射工具箱?

58视图(30天)
有一个函数显示世界的国家在全球的阴谋,我需要知道如何情节甚至增量的纬度和经度10度到全球阴谋没有使用geoshow命令映射工具箱中找到。
3评论

登录置评。

接受的答案

乍得格林
乍得格林 2017年8月2日
编辑:MathWorks支万博1manbetx持团队 2022年6月8日
更新——至于R2018b, geoplot 函数可以使用MATLAB作为一种阴谋经度和纬度数据二维全球地图。
* * * * * * *
你可以轻易与标准的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日
@Chad格林 这是伟大的!我修改你的代码使用 wgs84Ellipsoid geodetic2ecef 函数来做这些情节在一个球体(而不是一个球体)。下面这段代码:
%经度的情节:
[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))
我遇到了一些麻烦添加地图的情节,然而。让我知道如果你有任何建议!

登录置评。

答案(1)

乍得格林
乍得格林 2019年3月8日
我把我的答案变成了要调用的函数的集合 globeplot 、globesurf globegraticule globeborders等等,你可以找到的 气候数据为Matlab工具箱
2的评论

登录置评。

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!