如何计算遮蔽角

这篇具有很好参考价值的文章主要介绍了如何计算遮蔽角。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

1. 问题描述

对于指定地点(E94.36488,N29.29815)安置的雷达系统,假设已知雷达天线的安装高度,如10米,根据DEM数据,计算某一方位的遮蔽角随距离的变化情况。

遮蔽角是指由于地物影响,导致相应地物远方特定俯仰角以下的区域无法被雷达探测,该角度被称为遮蔽角,如下图所示。

雷达阵地遮蔽角计算,电磁波覆盖范围分析,算法,c++,开发语言

2. 问题分析

求解遮蔽角,首先,需要获得雷达站的站心坐标系中所求方位上的高程在该坐标系中的高程随距离的变化。然后,计算各点相对雷达天线的俯仰角。由于遮蔽的存在,远处较小俯仰角将由近处的较大俯仰角所替代,从而获得遮蔽角随距离的变化情况。

在上述分析中,还有一个问题需要解决。获取雷达站的站心坐标系中某一点的高程后,相应高程与其经纬度构成该点的大地坐标系的描述,需要将该点转换成站心坐标系的坐标,故需要利用博文覆盖范围分析中用到的坐标系统的方法进行处理。

3. 具体求解过程

1) 获取地物的大地坐标

假设要获取方位的遮蔽角变化。根据SRTM3数据是90米的距离间隔,不防选取距离步进为100米,即,以雷达站为原点,沿方位方向,每100米选取一点,来获取该点的大地坐标。

在站心坐标系中,在水平面上沿方位 每100米选取一点,假设当前计算的点为第点,,该点在站心坐标系中的坐标为,单位为米。

考虑平面显示的位置是根据兰勃特变换而来的,因此,若获得该点所对应的经纬度,需要利用利用逆兰勃特变换(详情参见博文兰勃特等角圆锥(Lambert Conformal Conic)投影正反变换)。

相应变换代码

/**
 * 单点Lambert反变换
 * @param Y 正东向坐标值,单位:米
 * @param X 正北向坐标值,单位:米
 * @param latitude0 参考点纬度,单位:度
 * @param longitude0 参考点经度,单位:度
 * @return 被变换点的经纬坐标,GroundPosition.x为经度,单位:度;
                             GroundPosition.y为纬度,单位:度。
 */
public static Point2D.Double arclambert(double Y, double X, 
        double latitude0, double longitude0) {
    //Y:纬度方向, X:经度方向
    double R_earth = 6378137.0;                 //地球半径
    double latitude1 = 25.0, latitude2 = 47.0;  //标准纬度(割纬度)
    Point2D.Double GroundPosition=new Point2D.Double() ;

    //圆锥常数: alpha = frac{ln(cos{latitude1})-ln(cos{latitude2})}
            ln(tan{45+latitude2/2})-ln(tan{45+latitude1/2})}
    double alpha = (Math.log(Math.cos(latitude1 * Math.PI / 180.0)) 
            - Math.log(Math.cos(latitude2 * Math.PI / 180.0)))
            / (Math.log(Math.tan((45.0 + latitude2 / 2.0) * Math.PI / 180.0)) 
            - Math.log(Math.tan((45.0 + latitude1 / 2.0) * Math.PI / 180.0)));
    //坐标原点对应纬圈投影圆弧半径
    double rho_original = (R_earth / alpha) * Math.cos(latitude1 * Math.PI / 180.0)
            * Math.pow(Math.tan((45.0 + latitude1 / 2.0) * Math.PI / 180.0) 
            / Math.tan((45.0 + latitude0 / 2.0) * Math.PI / 180.0), alpha);

    //需变换点对应纬圈投影圆弧半径
    double rho = Math.sqrt(Math.pow(X, 2) + Math.pow(Y, 2) 
            + Math.pow(rho_original, 2) - 2.0 * rho_original * Y);

    double latitude = (2 * Math.atan(Math.tan((45.0 + latitude1 / 2.0) * Math.PI / 180.0)
            / Math.pow(alpha * rho / (R_earth * Math.cos(latitude1 * Math.PI / 180.0)), 
            1 / alpha)) - Math.PI / 2.0) / Math.PI * 180.0;
    double longitude = longitude0 + Math.asin(X / rho) / alpha / Math.PI * 180.0;
    GroundPosition.x = longitude;
    GroundPosition.y = latitude;

    return GroundPosition;
}

根据点所对应的经纬度,采用最近邻法,可由SRTM3的文件中获得该点的高程数据(详情参见博文在Java GDAL环境中提取指定经纬度的高程数据),进而获得相应的大地坐标点。

2) 获取点的站心坐标

首先,由点的大地坐标计算出其在地心坐标系中的坐标,然后再变换为站心坐标.详情参见博文覆盖范围分析中用到的坐标系统。

这里需要注意的是,点的坐标与和点的坐标与之间存在差异,不完全相等,但差值并不会太大。

3) 计算遮蔽角

根据 点的站心坐标,计算各点相对雷达天线的俯仰角。

在获得所有点的俯仰角后,再由天线处出发,由近及远,两两比较俯仰角,并用较大的俯仰角更新较远点的俯仰角,作为该点的遮蔽角。文章来源地址https://www.toymoban.com/news/detail-558171.html

到了这里,关于如何计算遮蔽角的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 如何计算遮蔽角

    对于指定地点(E94.36488,N29.29815)安置的雷达系统,假设已知雷达天线的安装高度,如10米,根据DEM数据,计算某一方位的遮蔽角随距离的变化情况。 遮蔽角是指由于地物影响,导致相应地物远方特定俯仰角以下的区域无法被雷达探测,该角度被称为遮蔽角,如下图所示。

    2024年02月15日
    浏览(54)
  • 如何计算单元测试的覆盖率

    单元测试的覆盖率有:语句覆盖率(即行覆盖率)、分支覆盖率、条件覆盖率、分支条件覆盖率、路径覆盖率等。 语句覆盖率 所谓语句就是那些非分支、非判断的语句。 计算公式:程序执行到的语句总数 / 全部语句的总数 分支覆盖率 有判定语句的地方都会出现2个分支。

    2024年01月21日
    浏览(47)
  • 概念解析 | 电磁计算的新篇章:智能电磁计算

    注1:本文系“概念解析”系列之一,致力于简洁清晰地解释、辨析复杂而专业的概念。本次辨析的概念是:智能电磁计算。 随着人工智能的飞速发展,我们正处在一个信息爆炸的时代。各个领域,无论是自然科学还是社会科学,都在以前所未有的速度产生巨量的数据。这些

    2024年02月11日
    浏览(51)
  • 【电磁场计算】磁偶极子在外磁场中的受力

    最近在总结电磁场计算公式的时候发现,关于磁偶极子受力的计算公式,在不同文献中的说明有点出入,在此特地写一篇博客澄清一下。 通过磁偶极子在外磁场中的磁势能求梯度得出其受力: F = − ∇ ( U ) = − ∇ ( − m ⋅ B ) = m × ( ∇ × B ) + ( m ⋅ ∇ ) B begin{align*} F = -nab

    2024年02月08日
    浏览(47)
  • 基于HFSS的微带线特性阻抗仿真-与基于FDTD的计算电磁学方法对比(Matlab)

    工程下载: HFSS的微带线特性阻抗仿真工程文件(注意版本:HFSS 2023R2 ): https://download.csdn.net/download/weixin_44584198/88748285 基于FDTD的微带线特性阻抗仿真Matlab工程: https://download.csdn.net/download/weixin_44584198/88748215 微带线特性阻抗的仿真是电磁学里面一个比较基本的仿真了,很多

    2024年01月17日
    浏览(41)
  • 电磁兼容(EMC)设计如何融入产品研发流程~系统流程法

    系统流程法(System Flow Method) 产品工程师可以通过短期的培训以及通过积累经验基本掌握 EMC 设计的方法,但对于一个企业来讲,目前迫切的是建立一套规范的 EMC 设计流程,把电磁兼容要求融入产品设计中去,这样才能保证企业大多产品经过这样的流程顺利通过测试认证。如果

    2023年04月22日
    浏览(50)
  • sonar覆盖率、代码覆盖率、分支覆盖率的计算方式

    代码质量的覆盖率分为三种,覆盖率、代码覆盖率、分支覆盖率,那每一种的计算方式是怎么样的呢? 举例: 上面最有疑惑的是覆盖率,不知道怎么算出了来的,后面再说。 通过sonarqube可以分析出: 指标 值 可覆盖行(lines_to_cover) 13242 未覆盖的代码(uncovered_lines) 7943 可

    2024年02月06日
    浏览(40)
  • 利用Matab进行覆盖计算----战术计算

    每个点范围内的近似面积为resdistlat*resdistlon 三维视图: 数据的三维图: % 心型不完整,最上面是直线型的,不知道为什么 得出对应的XYZ坐标之后,放大,坐标转换到对应的地理坐标系中,然后加载到地图。 后记: Capacity Map Using SINR Data 注:对原数据进行修改后,再返回进行

    2024年02月01日
    浏览(40)
  • MindSpore AI科学计算系类 | VAE基于MindSpore Elec的特征域MT反演,提升大地电磁反演的精度和分辨率

    背景 昇思MindSpore联合清华大学、华为先进计算与存储实验室合作构建了基于变分自编码器的特征域大地电磁反演算法,通过灵活嵌入多物理先验知识有效提升了复杂介质地球物理反演的精度。 大地电磁(Magnetotelluric, MT)数据反演是通过地表测量的天然电磁场推断地下介质分

    2024年01月25日
    浏览(46)
  • ​5种常用于LLM的令牌遮蔽技术介绍以及Pytorch的实现

    本文将介绍大语言模型中使用的不同令牌遮蔽技术,并比较它们的优点,以及使用Pytorch实现以了解它们的底层工作原理。 令牌掩码Token Masking是一种广泛应用于语言模型分类变体和生成模型训练的策略。BERT语言模型首先使用,并被用于许多变体(RoBERTa, ALBERT, DeBERTa…)。 而Te

    2024年04月17日
    浏览(35)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包