[发明专利]一种采用北斗网格位置编码的地理坐标解算方法

说明书

本发明公开了一种采用北斗网格位置编码的地理坐标解算方法，根据国家标准《北斗网格位置码》(GB/T 394092020)对狭义北斗网格位置码的描述，确定北斗网格位置码和经纬度之间的数学关系；利用计算机数学运算速度快、可靠性高的特点，利用python语言描述北斗网格位置码和经纬度之间的数学关系，进而计算出狭义北斗网格位置码。

技术领域

本发明属于地理坐标信息计算领域，具体涉及一种采用北斗网格位置编码的地理坐标解算方法。

背景技术

北斗卫星导航系统是我国着眼于国家安全和经济社会发展需要，自主建设运行的全球卫星导航系统，是为全球用户提供全天候、全天时、高精度的定位导航授时服务的国家重要时空基础设施。北斗系统的建成运用为相关技术与产业的发展提供了基础，北斗网格码就是其中之一。简化北斗网格位置编码的地理坐标解算方法和优化表现形式，以便于更好的利用和推广北斗网格码。

发明内容

发明目的：本发明提供一种简化北斗网格位置编码的地理坐标解算方法和优化北斗网格码的表现形式。

技术方案：一种卫星飞临时刻与拍摄范围的计算方法，包括以下步骤：

步骤1，计算L1级北斗网格码，按照1：100万图幅划分，经度方向用1～60编码，纬度方向分南北半球(N,S)，按照A～V编码，大小6°×4°。

步骤2，将1:100万6°×4°网格，从左下角划分(东北半球)分成6x4个网格，一个方向用A～X表示，相当于GeoSOT第9层1°×1°网格，相当于128km×128km网格。

步骤3，将GeoSOT 1°×1°网格，分成4×6个网格，一个方向分别用A～P表示，相当于1:5万地图图幅15’×10’网格，相当于30km×20km网格。

步骤4，将1:5万地图图幅15’×10’网格，从左下角划分(东北半球)分成15×10个网格，两个方向分别用0-E，0～9表示，相当于GeoSOT第15层1'×1'网格，相当于2km×2km网格。

步骤5，将GeoSOT第15层1×1’网格，从左下角划分(东北半球)分成15x15个网格，一个方向分别用0～E表示，相当于GeoSOT第19层4”×4”网格，相当于128m×128m网格。

步骤6，将GeoSOT第19层4'×4’网格，从左下角划分(东北半球)分成2×2个网格，一个方向分别用0～3表示，相当于GeoSOT第20层2×2网格，相当于64m×64m网格。

步骤7，将GeoSOT第20层2×2网格，从左下角划分(东北半球)分成8x8个网格，两个方向分别用0～7表示,相当于GeoSOT第23层1/4×1/4，相当于8m×8m网格。

步骤8，将GeoSOT第23层1/4×1/4网格，从左下角划分(东北半球)分成8x8个网格，两个方向分别用0～7表示，相当于GeoSOT第26层1/32×1/32，相当于1m×1m网格。

附图说明

图1是L3-L7级二维北斗位置码示意图；

图2是L3-L7级三维北斗位置码示意图；

图3为经纬度转北斗网格码程序截图；

具体实施方式

下面结合附图和具体实施方式，进一步阐明本发明。

如图1-3所示，一种卫星飞临时刻与拍摄范围的计算方法，包括以下步骤：

步骤1，计算L1级北斗网格码，按照1：100万图幅划分，经度方向用1～60编码，纬度方向分南北半球(N,S)，按照A～V编码，大小6°×4°。通过图三的经纬度转北斗网格码程序进行转换，之后网格建模，分别在二三维中展示。