绘制3 d穹顶网在全球范围内
这个例子展示了如何从一个3 d功能的系统本地east-north-up (ENU表示)坐标,然后变换,把它与一个全球显示在地球,地球定点(ECEF)坐标。
步骤1:设置定义参数
使用大地参考系统1980公里(GRS80)和工作单位。华盛顿附近的本地系统的起源,华盛顿,美国。
grs80 = referenceEllipsoid (“grs80”,“公里”);domeRadius = 3000;%公里domeLat = 39;%度domeLon = -77;%度domeAlt = 0;%公里
步骤2:在当地建造圆顶East-North-Up坐标
定义本地ENU表示系统对大地测量参考点,在这种情况下,指定(domeLat,domeLon,domeAlt)。它是一个三维笛卡尔坐标系统的正轴指向东方,正轴指向北方,z轴是正常的参考椭球体和向上。
在这个例子中,3 d功能是半空间z > = 0的半球的半径3000公里。这个半球可以附上,假设,体积的空间范围内的一种理想化的雷达系统有统一的覆盖范围从地平线到天顶,方位角》。这样的空间,当代表区域有效的监测范围,有时非正式称为“雷达穹顶。”
快速构建坐标数组列出一个封闭的半球圆顶与单位球开始,扩展半径、下半球崩溃。很容易想象如果你让它半透明——设置FaceAlpha在这种情况下,至0.5。
[x, y, z] =球体(20);xEast = domeRadius * x;yNorth = domeRadius * y;z上= domeRadius * z;z上(z上< 0)= 0;冲浪(xEast yNorth z上,“FaceColor”,“黄色”,“FaceAlpha”,0.5)轴平等的
步骤3:将圆顶的地球地球定点(ECEF)系统
使用enu2ecef函数将从当地ENU表示一个圆顶ECEF系统,基于GRS 80参考椭球体。它适用于三维平移和旋转。注意西半球变得倾斜以及其中心数千公里从原点。
[xECEF, yECEF zECEF]…= enu2ecef (xEast yNorth、z上domeLat, domeLon, domeAlt, grs80);冲浪(xECEF yECEF zECEF,“FaceColor”,“黄色”,“FaceAlpha”,0.5)轴平等的
第四步:构建一个全球范围内显示
构建一个基本的全球显示使用axesm和全球。
图ax = axesm (“全球”,“大地水准面”grs80,“网格”,“上”,…“GLineWidth”,1“GLineStyle”,“-”,…“Gcolor”(0.9 0.9 0.1),“Galtitude”,100);斧子。位置= (0 0 1 1);轴平等的从视图(3)
第五步:添加各种全球地图数据
添加低分辨率的全球地形、海岸线和河流。
负载topo60cgeoshow (topo60c topo60cR,“DisplayType”,“texturemap”)demcmap (topo60c)负载海岸线geoshow (coastlat coastlon,“颜色”,“黑色”)河流= readgeotable (“worldrivers.shp”);geoshow(河流、“颜色”,“蓝色”)
步骤6:添加穹顶全球显示
圆顶的ECEF版本添加到全球轴半透明网状。
冲浪(xECEF yECEF zECEF,“FaceColor”,“黄色”,“FaceAlpha”,0.5)
您可以查看圆顶和全球交互式地从不同角度的旋转坐标轴。
信贷
由于爱德华·j·梅休jr .)提供技术背景对“雷达穹顶”,把我们的注意力问题的可视化映射工具箱™。
