简介

站心坐标系，是以接收机位置作为zh，正北方向为方位角零度方向，天顶方向为仰角90°方向的坐标系。计算不同卫星相对于接收机位置的仰角与方位角，确定卫星在坐标系中的位置，即可绘制站心星空图。中心

读取文件

上图为不同编号卫星的仰角与方位角信息。首先读取文件：

绘制星空图

function hpol = skyPlot(az, el, prn, ss)
%% 准备坐标轴 ===========================================================
hAxis = newplot([]);
hold(hAxis, 'on');
%--- 绘制圆环 --------------------------------------------------------
rectangle('position', [-90, -90, 180, 180], ...'Curvature', [1 1], ...'facecolor', 'white', ...'edgecolor', 'k');
%--- 绘制辐条 ----------------------------------------------------
th = (1:6) * 2*pi / 12;
cst = cos(th); snt = sin(th);
cs = [cst; -cst];
sn = [snt; -snt];
line(90*sn, 90*cs, 'linestyle', ':', 'color', 'k', 'linewidth', 0.5, ...'handlevisibility', 'off');
rt = 1.1 * 90;
for i = 1:max(size(th))%--- 在绘图的右半部分标注方位角 -------------------------text(rt*snt(i), rt*cst(i), int2str(i*30), ...'horizontalalignment', 'center', 'handlevisibility', 'off');if i == max(size(th))loc = int2str(0);elseloc = int2str(180 + i*30);end%--- 在绘图的另一半标注方位角 ----------------------text(-rt*snt(i), -rt*cst(i), loc, ...'handlevisibility', 'off', 'horizontalalignment', 'center');
end%--- 绘制仰角网格 ----------------------------------------
th = 0 : pi/50 : 2*pi;
xunit = cos(th);
yunit = sin(th);
for elevation = 0 : 15 : 90elevationSpherical = 90*cos((pi/180) * elevation);line(yunit * elevationSpherical, xunit * elevationSpherical, ...'lineStyle', ':', 'color', 'k', 'linewidth', 0.5, ...'handlevisibility', 'off');text(0, elevationSpherical, num2str(elevation), ...'BackgroundColor', 'white', 'horizontalalignment','center', ...'handlevisibility', 'off');
end
%--- 切换为二维视角 ------------------------------------------------------
view(0, 90);
%--- 设定坐标轴范围 ------------------------------------------------------
axis([-95 95 -90 101]);
%% 将仰角转换为到中心的距离 ------
elSpherical = 90*cos(el * pi/180);
%--- 转换为笛卡尔坐标 ------------------------------
yy = elSpherical .* cos(az * pi/180);
xx = elSpherical .* sin(az * pi/180);
%% 在网格顶部绘制数据 ===========================================
hpol = plot(hAxis, xx', yy', '.-');
%--- 标注卫星最后位置 ------------------------------
plot(hAxis, xx(:,end)', yy(:,end)', 'o', 'MarkerSize', 7);%--- 标注卫星编号 -------------------
for i = 1:length(prn)if(prn(i) ~= 0)text(xx(i, end), yy(i, end), ['  ',ss, int2str(prn(i))], 'color', 'b');end
end
%--- 设定xy轴具有相同的数据纵横比 ------------------
axis(hAxis, 'equal');
%--- 关闭笛卡尔坐标系坐标轴 -------------------------------
axis(hAxis, 'off');

利用不同卫星的仰角方位角绘制站心星空图相关推荐

计算卫星高度角、方位角
最小二乘定权.电离层对流层改正,都需要卫星的高度角.方位角.本章将介绍求解完卫星的地固坐标系的位置后,如何求解卫星的高度角.方位角. 卫星位置求解请参考之前的博客:卫星位置解算原理与程序设计参考书籍 ...
Python之pyecharts：利用pyecharts(十五种风格)绘制2018年~2020年中国互联网Top巨头近三年市值(美元)柱状图
Python之pyecharts:利用pyecharts(十五种风格)绘制2018年~2020年中国互联网Top巨头近三年市值(美元)柱状图目录 pyecharts(十五种风格)绘制2018年~20 ...
Python：利用python编程实现三维图像绘制展示(六面体旋转、三维球柱状体、下雪场景等)
Python:利用python编程实现三维图像绘制展示(六面体旋转.三维球柱状体.下雪场景等) 目录利用python编程实现三维图像绘制展示(六面体旋转.三维球柱状体.下雪场景等) 1.实现六面体旋 ...
利用人工智能和进化分析，绘制出真核生物的蛋白质之间相互作用的3D模型
美国科学家主导的国际科研团队在最新一期<科学>杂志撰文指出,他们利用人工智能和进化分析,绘制出了真核生物的蛋白质之间相互作用的3D模型,首次确定了100多个可能的蛋白质复合物,并为700多 ...
python之利用Pygal+2010年世界人口=绘制世界地图
python之利用Pygal+2010年世界人口=绘制世界地图工具: 博主是在Pycharm下进行的. 目的: Pygal+2010年世界人口=绘制世界地图科普: JSON格式 JSON格式介绍 ...
基于gan和孪生网络框架，利用双时间光学遥感图像绘制滑坡清单Landslide Inventory Mapping Using Bi-Temporal Optical Remote Sensin
Landslide Inventory Mapping Using Bi-Temporal Optical Remote Sensing Images Remote Sensing Images 基于 ...
利用Scrapy框架爬取LOL皮肤站高清壁纸
利用Scrapy框架爬取LOL皮肤站高清壁纸 Lan 2020-03-06 21:22 81 人阅读 0 条评论成品打包:点击进入代码: 爬虫文件 # -*- coding: utf- ...
【matlab】动态绘制三维心形图
[转载]matlab动态绘制三维心形图 matlab动态绘制三维心形图程序如下,核心是isosurface函数,verts存储了心形曲面各点的坐标,而faces存储了各点的连接顺序. clear fi ...
Unity中的坐标与绘制准心
unity3d中的屏幕坐标系是以屏幕左下角为(0,0)点右上角为(Screen.Width,Screen.Height) 鼠标位置坐标与屏幕坐标系一致视口坐标是以摄像机为准以屏幕的左下 ...

利用不同卫星的仰角方位角绘制站心星空图

简介

读取文件

绘制星空图

利用不同卫星的仰角方位角绘制站心星空图相关推荐

最新文章

热门文章