CN116430346A - 一种sar图像像素点经纬度计算方法及装置 - Google Patents

一种sar图像像素点经纬度计算方法及装置 Download PDF

Info

Publication number
CN116430346A
CN116430346A CN202310692703.4A CN202310692703A CN116430346A CN 116430346 A CN116430346 A CN 116430346A CN 202310692703 A CN202310692703 A CN 202310692703A CN 116430346 A CN116430346 A CN 116430346A
Authority
CN
China
Prior art keywords
point
latitude
longitude
earth
distance
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202310692703.4A
Other languages
English (en)
Other versions
CN116430346B (zh
Inventor
张少雄
何宇
陈立坦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Lingkong Electronic Technology Co Ltd
Original Assignee
Xian Lingkong Electronic Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Lingkong Electronic Technology Co Ltd filed Critical Xian Lingkong Electronic Technology Co Ltd
Priority to CN202310692703.4A priority Critical patent/CN116430346B/zh
Publication of CN116430346A publication Critical patent/CN116430346A/zh
Application granted granted Critical
Publication of CN116430346B publication Critical patent/CN116430346B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/06—Systems determining position data of a target
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021—SAR image post-processing techniques
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/60—Analysis of geometric attributes
    • G06T7/66—Analysis of geometric attributes of image moments or centre of gravity
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/70—Determining position or orientation of objects or cameras
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00—Indexing scheme for image analysis or image enhancement
    • G06T2207/10—Image acquisition modality
    • G06T2207/10032—Satellite or aerial image; Remote sensing

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Electromagnetism (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)

Abstract

本申请公开了一种SAR图像像素点经纬度计算方法及装置，该方法包括：获取SAR图像参数；根据SAR图像参数确定待求点与已知点的横纵坐标之间的实际向量距离；根据椭球第一偏心率平方、地球赤道半径、已知点的经纬度和实际向量距离，确定待求点在地球坐标系中的坐标值；根据椭球第一偏心率平方、地球赤道半径和待求点在地球坐标系中的坐标值，确定待求点的经纬度。该方法在确定SAR图像中某个点的经纬度时，只需要一个已知点的像素坐标和经纬度，减少了对已知条件的要求，并且在计算某个点在地球坐标系中的坐标值和经纬度的具体步骤中，考虑了地球的实际形状，最终所确定的SAR图像中某个点的经纬度更加准确。

Description

一种SAR图像像素点经纬度计算方法及装置
技术领域
本申请涉及合成孔径技术领域，尤其涉及一种SAR图像像素点经纬度计算方法及装置。
背景技术
合成孔径雷达（英文：Synthetic Aperture Radar；简称：SAR）是一种主动式的对地观测系统，其在进行工作时，向地面发送连续的无线电脉冲，并接收和记录每个脉冲的回波，形成SAR图像。
目前，在计算SAR图像各像素经纬度时，将SAR图像按照平面图像处理方法，通过四个角的经纬度或四个角经纬度加中间点的经纬度去计算其他像素点的经纬度。目前所采用的方法至少需要SAR图像四个角的精准经纬度，对已知条件要求比较苛刻，并且由于没有考虑地球的实际形状，导致所得到的其他像素点的经纬度误差较大。
发明内容
本申请实施例通过提供一种SAR图像像素点经纬度计算方法及装置，解决了现有技术中计算SAR图像各像素经纬度时，对已知条件要求比较苛刻和所得到的其他像素点的经纬度误差较大的技术问题。
第一方面，本申请实施例提供了一种SAR图像像素点经纬度计算方法，所述方法包括：
获取SAR图像参数；所述SAR图像参数包括已知点的像素坐标和经纬度、分辨率、中心方位角、航迹角、侧视方向、以及待求点的像素坐标；
根据所述SAR图像参数确定所述待求点与所述已知点的横纵坐标之间的实际向量距离；
根据椭球第一偏心率平方、地球赤道半径、所述已知点的经纬度和所述实际向量距离，确定所述待求点在地球坐标系中的坐标值；
根据椭球第一偏心率平方和所述待求点在地球坐标系中的坐标值，确定所述待求点的经纬度。
结合第一方面，在一种可能的实现方式中，所述根据所述SAR图像参数确定所述待求点与所述已知点的横纵坐标之间的实际向量距离，包括：
根据所述已知点的像素坐标、所述待求点的像素坐标和所述分辨率，确定所述待求点与所述已知点的横纵坐标之间的理论距离；
根据所述理论距离、所述中心方位角和所述侧视方向，确定所述待求点与所述已知点的横纵坐标之间的实际距离；
根据所述航迹角和所述实际距离，确定所述待求点与所述已知点的横纵坐标之间的实际向量距离。
结合第一方面，在一种可能的实现方式中，所述确定所述待求点与所述已知点的横纵坐标之间的理论距离，包括：
根据所述已知点的像素坐标和所述待求点的像素坐标，计算所述待求点与所述已知点的横纵坐标之间像素点个数；
根据所述分辨率、以及所述待求点与所述已知点的横纵坐标之间像素点个数，计算所述待求点与所述已知点的横纵坐标之间的所述理论距离。
结合第一方面，在一种可能的实现方式中，所述确定所述待求点与所述已知点的横纵坐标之间的实际距离，包括：
根据所述理论距离、所述中心方位角和所述侧视方向，计算所述待求点与所述已知点的横纵坐标之间的误差距离；
根据所述理论距离和所述误差距离，计算所述待求点与所述已知点的横纵坐标之间的所述实际距离。
结合第一方面，在一种可能的实现方式中，所述确定所述待求点在地球坐标系中的坐标值，包括：
根据椭球第一偏心率平方、地球赤道半径和所述已知点的经纬度，计算所述已知点在地球坐标系中的坐标值；
根据所述实际向量距离和所述已知点的经纬度，计算所述已知点与所述待求点在地球坐标系中的差值；
根据所述已知点在地球坐标系中的坐标值、以及所述已知点与所述待求点在地球坐标系中的差值，计算所述待求点在地球坐标系中的坐标值。
结合第一方面，在一种可能的实现方式中，所述计算所述已知点在地球坐标系中的坐标值，包括：
根据椭球第一偏心率平方、地球赤道半径和所述已知点的纬度，计算所述已知点与地心的距离；
根据所述已知点与地心的距离、以及所述已知点的经纬度，计算所述已知点在地球坐标系中的X、Y轴方向的坐标值；
根据所述已知点与地心的距离、椭球第一偏心率平方和所述已知点的纬度，计算所述已知点在地球坐标系中的Z轴方向的坐标值。
结合第一方面，在一种可能的实现方式中，所述确定所述待求点的经纬度，包括：
根据所述待求点在地球坐标系中的坐标值，计算距离地球坐标系Z轴的垂直距离，以及计算所述待求点的经度：
根据椭球第一偏心率平方、所述待求点在地球坐标系中的坐标值、以及所述待求点距离地球坐标系Z轴的垂直距离，计算所述待求点的纬度。
第二方面，本申请实施例提供了一种SAR图像像素点经纬度计算装置，所述装置包括：
获取模块，用于获取SAR图像参数；所述SAR图像参数包括已知点的像素坐标和经纬度、分辨率、中心方位角、航迹角、侧视方向、以及待求点的像素坐标；
距离计算模块，用于根据所述SAR图像参数确定所述待求点与所述已知点的横纵坐标之间的实际向量距离；
坐标计算模块，用于根据椭球第一偏心率平方、地球赤道半径和所述实际向量距离，确定所述待求点在地球坐标系中的坐标值；
经纬计算模块，用于根据椭球第一偏心率平方、地球赤道半径和所述待求点在地球坐标系中的坐标值，确定所述待求点的经纬度。
第三方面，本申请实施例提供了一种SAR图像像素点经纬度计算设备，所述设备包括处理器，以及与所述处理器通信连接的存储器；
所述存储器存储有计算机可读指令；
所述处理器执行所述存储器存储的所述计算机可读指令时，实现如第一方面或第一方面任一种可能的实现方式所述的SAR图像像素点经纬度计算方法。
第四方面，本申请实施例提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机可读指令，所述计算机可读指令被处理器执行时实现如第一方面或第一方面任一种可能的实现方式所述的SAR图像像素点经纬度计算方法。
本申请实施例中提供的技术方案，至少具有如下技术效果：
本申请实施例提供了一种SAR图像像素点经纬度计算方法，该方法获取SAR图像参数后确定待求点与已知点的横纵坐标之间的实际向量距离，并根据椭球第一偏心率平方、地球赤道半径、已知点的经纬度和实际向量距离确定待求点在地球坐标系中的坐标值，而后确定待求点的经纬度。本申请实施例所提供的方法在确定SAR图像中某个点的经纬度时，只需要一个已知点的像素坐标和经纬度，而不需要SAR图像四个角的精准经纬度，减少了对已知条件的要求；并且在计算某个点在地球坐标系中的坐标值和经纬度的具体步骤中，考虑了地球的实际形状，最终所确定的SAR图像中某个点的经纬度更加准确。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对本申请实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本申请实施例提供的SAR图像像素点经纬度计算方法的流程图；
图2为本申请实施例提供的确定待求点与已知点的横纵坐标之间的实际向量距离的流程图；
图3为本申请实施例提供的计算待求点与已知点的横纵坐标之间的理论距离的流程图；
图4为本申请实施例提供的计算待求点与已知点的横纵坐标之间的实际距离的流程图；
图5为本申请实施例提供的计算待求点在地球坐标系中的坐标值的流程图；
图6为本申请实施例提供的计算已知点在地球坐标系中的坐标值；
图7为本申请实施例提供的计算待求点的经纬度的流程图；
图8为本申请实施例提供的SAR图像像素点经纬度计算装置的结构示意图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本申请实施例提供了一种SAR图像像素点经纬度计算方法，如图1所示，该方法包括步骤S101至S104。
S101：获取SAR图像参数。SAR图像参数包括已知点的像素坐标和经纬度、分辨率、中心方位角、航迹角、侧视方向、以及待求点的像素坐标。
分辨率包括方位像元尺寸和距离像元尺寸。方位像元尺寸是指SAR图像中一个像素在方位方向（即SAR系统的运动方向）上所代表的实际距离；距离像元尺寸是指SAR图像中一个像素在距离方向（即垂直于SAR系统的运动方向）上所代表的实际距离。
中心方位角是指SAR天线指向图像中央像素时的方位角，也即从正北方向开始逆时针旋转到SAR天线指向图像中央像素的方位角度数。
航迹角是指SAR系统的运动方向与正北方向之间的夹角，也即从正北方向开始逆时针旋转到SAR系统的运动方向的度数。
侧视方向是指雷达波束与地面法线之间的夹角，也即雷达朝向地面的方向。
S102：根据SAR图像参数确定待求点与已知点的横纵坐标之间的实际向量距离。
S103：根据椭球第一偏心率平方、地球赤道半径、已知点的经纬度和实际向量距离，确定待求点在地球坐标系中的坐标值。具体地，椭球第一偏心率平方为0.006694379995，地球赤道半径为6378137。其中，椭球第一偏心率平方体现了地球的实际形状，该步骤S103在确定待求点在地球坐标系中的坐标值时，考虑到了地球的实际形状，使得所求解的坐标值更加准确。
S104：根据椭球第一偏心率平方和待求点在地球坐标系中的坐标值，确定待求点的经纬度。该步骤S104同样考虑了地球的实际形状，使得所求解的经纬度更加准确。
如图2所示，S102在确定实际向量距离时具体包括步骤S201至S203，其描述了从SAR图像参数逐步推导实际向量距离的过程。
S201：根据已知点的像素坐标、待求点的像素坐标和分辨率，确定待求点与已知点的横纵坐标之间的理论距离。
S202：根据理论距离、中心方位角和侧视方向，确定待求点与已知点的横纵坐标之间的实际距离。
S203：根据航迹角和实际距离，确定待求点与已知点的横纵坐标之间的实际向量距离。
图3示出了S201在确定理论距离时的具体步骤S301和S302，详细描述及计算公式如下。
S301：根据已知点的像素坐标和待求点的像素坐标，计算待求点与已知点的横纵坐标之间像素点个数。
具体地，已知点的像素坐标值为
Figure SMS_1
，待求点的像素坐标值为
Figure SMS_2
。那 么，待求点与已知点的横坐标之间的像素点个数为：
Figure SMS_3
；待求点与已知点的纵坐标之 间的像素点个数为：
Figure SMS_4
。
S302：根据分辨率、以及待求点与已知点的横纵坐标之间像素点个数，计算待求点与已知点的横纵坐标之间的理论距离。具体地，理论距离通过以下公式计算：
Figure SMS_5
；
Figure SMS_6
；
Figure SMS_7
、
Figure SMS_8
分别表示待求点与已知点的横纵坐标之间的理论距离，
Figure SMS_9
表示距离像元 尺寸，
Figure SMS_10
表示方位像元尺寸。
图4示出S202在确定实际距离时的具体步骤S401和S402，详细描述及计算公式如下。
S401：根据理论距离、中心方位角和侧视方向，计算待求点与已知点的横纵坐标之间的误差距离。具体计算公式如下：
Figure SMS_11
；
Figure SMS_12
；
其中，
Figure SMS_13
、
Figure SMS_14
分别表示待求点与已知点的横纵坐标之间的误差距离，
Figure SMS_15
、
Figure SMS_16
分别 表示待求点与已知点的横纵坐标之间的理论距离，
Figure SMS_17
表示中心方位角，
Figure SMS_18
表示侧视方 向。
S402：根据理论距离和误差距离，计算待求点与已知点的横纵坐标之间的实际距离。具体计算公式如下：
Figure SMS_19
；
Figure SMS_20
；
其中，
Figure SMS_21
、
Figure SMS_22
表示待求点与已知点的横纵坐标之间的实际距离。
S203在确定实际向量距离时具体包括：根据航迹角、以及待求点与已知点的横纵坐标之间的实际距离，计算待求点与已知点的横纵坐标之间的实际向量距离。具体计算公式如下：
Figure SMS_23
；
Figure SMS_24
；
其中，
Figure SMS_25
、
Figure SMS_26
分别表示待求点与已知点的横纵坐标之间的实际向量距离，
Figure SMS_27
表示航迹角。
如图5所示，S103确定待求点在地球坐标系中的坐标值具体包括步骤S501至S503， 其中已知点的经纬度为
Figure SMS_28
。
S501：根据椭球第一偏心率平方、地球赤道半径和已知点的经纬度，计算已知点在地球坐标系中的坐标值。
S502：根据实际向量距离和所述已知点的经纬度，计算已知点与所述待求点在地球坐标系中的差值。
Figure SMS_29
Figure SMS_30
Figure SMS_31
；
其中，
Figure SMS_32
、
Figure SMS_33
、
Figure SMS_34
分别表示已知点与待求点在地球坐标系X、Y、Z轴方向的差值。
S503：根据已知点在地球坐标系中的坐标值、以及已知点与所述待求点在地球坐标系中的差值，计算待求点在地球坐标系中的坐标值。
Figure SMS_35
；
Figure SMS_36
；
Figure SMS_37
；
其中，
Figure SMS_38
、
Figure SMS_39
、
Figure SMS_40
分别表示已知点在地球坐标系中的X、Y、Z轴方向的坐标值，
Figure SMS_41
、
Figure SMS_42
、
Figure SMS_43
分别表示待求点在地球坐标系中X、Y、Z轴方向的坐标值。
如图6所示，S501计算已知点在地球坐标系中的坐标值具体包括S601至S603，详细描述及计算公式如下。
6501：根据椭球第一偏心率平方、地球赤道半径和已知点的纬度，计算已知点与地心的距离。具体计算公式如下：
Figure SMS_44
；
其中，r表示已知点与地心的距离。
S602：根据已知点与地心的距离、以及已知点的经纬度，计算已知点在地球坐标系中的X、Y轴方向的坐标值。具体计算公式如下：
Figure SMS_45
；
Figure SMS_46
；
其中，
Figure SMS_47
、
Figure SMS_48
分别表示已知点在地球坐标系中的X、Y轴方向的坐标值。
S603：根据已知点与地心的距离、椭球第一偏心率平方和已知点的纬度，计算已知点在地球坐标系中的Z轴方向的坐标值。具体计算公式如下：
Figure SMS_49
；
其中，
Figure SMS_50
表示已知点在地球坐标系中的Z轴方向的坐标值。
本申请实施例中的S601至S603在对已知点在地球坐标系中的坐标值进行求解时，考虑到了地球的实际形状，进而使得计算得出的坐标值更加准确。
如图7所示，S104在确定待求点的经纬度时具体包括步骤S701和S702。
S701：根据待求点在地球坐标系中的坐标值，计算距离地球坐标系Z轴的垂直距离，以及计算待求点的经度。具体计算公式如下：
Figure SMS_51
；其中，
Figure SMS_52
表示待求点距离地球坐标系Z轴的垂直距离；
Figure SMS_53
；其中，
Figure SMS_54
表示待求点的经度。
S702：根据椭球第一偏心率平方、待求点在地球坐标系中的坐标值、以及待求点距离地球坐标系Z轴的垂直距离，计算待求点的纬度。具体计算公式如下：
Figure SMS_55
；其中，
Figure SMS_56
表示待 求点的纬度。
图7中所示的S701和S702分别对待求点的经度和纬度进行了求解，并且在求解待求点的纬度时加入了椭球第一偏心率平方，考虑到了地球的实际形状，进而使得所确定的待求点的经纬度更加准确。
本申请实施例提供的方法在确定SAR图像中某个点的经纬度时，只需要一个已知点的像素坐标和经纬度，而不需要SAR图像四个角的精准经纬度，减少了对已知条件的要求；并且在计算某个点在地球坐标系中的坐标值和经纬度的具体步骤中，考虑了地球的实际形状，最终所确定的SAR图像中某个点的经纬度更加准确。
虽然本申请提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。本实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照本实施例或者附图所示的方法顺序执行或者并行执行（例如并行处理器或者多线程处理的环境）。
本申请实施例还提供了一种SAR图像像素点经纬度计算装置800，如图8所示，该SAR图像像素点经纬度计算装置800包括获取模块801、距离计算模块802、坐标计算模块803和经纬计算模块804。
获取模块801用于获取SAR图像参数。SAR图像参数包括已知点的像素坐标和经纬度、分辨率、中心方位角、航迹角、侧视方向、以及待求点的像素坐标。
距离计算模块802用于根据SAR图像参数确定待求点与已知点的横纵坐标之间的实际向量距离。坐标计算模块803用于根据椭球第一偏心率平方、地球赤道半径、已知点的经纬度和实际向量距离，确定待求点在地球坐标系中的坐标值。经纬计算模块804用于根据椭球第一偏心率平方和待求点在地球坐标系中的坐标值，确定待求点的经纬度。
其中，距离计算模块802包括理论距离计算模块、实际距离计算模块和实际向量距离计算模块。理论距离计算模块用于根据已知点的像素坐标、待求点的像素坐标和分辨率，确定待求点与已知点的横纵坐标之间的理论距离。实际距离计算模块用于理论距离、中心方位角和侧视方向，确定待求点与已知点的横纵坐标之间的实际距离。实际向量距离计算模块用于根据航迹角和实际距离，确定待求点与已知点的横纵坐标之间的实际向量距离。
理论距离计算模块具体用于：根据已知点的像素坐标和待求点的像素坐标，计算待求点与已知点的横纵坐标之间像素点个数；根据分辨率、以及待求点与已知点的横纵坐标之间像素点个数，计算待求点与已知点的横纵坐标之间的理论距离。具体计算公式参照本申请实施例以上所提供的方法所示。
实际距离计算模块具体用于：根据理论距离、中心方位角和侧视方向，计算待求点与已知点的横纵坐标之间的误差距离；根据理论距离和误差距离，计算待求点与已知点的横纵坐标之间的实际距离。具体计算公式参照本申请实施例以上所提供的方法所示。
坐标计算模块803具体用于：根据椭球第一偏心率平方、地球赤道半径和已知点的经纬度，计算已知点在地球坐标系中的坐标值；根据实际向量距离和已知点的经纬度，计算已知点与待求点在地球坐标系中的差值；根据已知点在地球坐标系中的坐标值、以及已知点与待求点在地球坐标系中的差值，计算待求点在地球坐标系中的坐标值。具体计算公式参照本申请实施例以上所提供的方法所示。
坐标计算模块803在计算已知点在地球坐标系中的坐标值时，具体包括：根据椭球第一偏心率平方、地球赤道半径和已知点的纬度，计算已知点与地心的距离；根据已知点与地心的距离、以及已知点的经纬度，计算已知点在地球坐标系中的X、Y轴方向的坐标值；根据已知点与地心的距离、椭球第一偏心率平方和已知点的纬度，计算已知点在地球坐标系中的Z轴方向的坐标值。
经纬计算模块804具体用于：根据待求点在地球坐标系中的坐标值，计算距离地球坐标系Z轴的垂直距离，以及计算待求点的经度；根据椭球第一偏心率平方、待求点在地球坐标系中的坐标值、以及待求点距离地球坐标系Z轴的垂直距离，计算待求点的纬度。具体计算公式参照本申请实施例以上所提供的方法所示。
上述实施例阐明的装置或模块，具体可以由计算机芯片或实体实现，或者由具有某种功能的产品来实现。为了描述的方便，描述以上装置时以功能分为各种模块分别描述。在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。当然，也可以将实现某功能的模块由多个子模块或子单元组合实现。
本申请所述装置中的部分模块可以在由计算机执行的计算机可执行指令的一般上下文中描述，例如程序模块。一般地，程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等。也可以在分布式计算环境中实施本申请，在这些分布式计算环境中，由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中，程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本申请实施例还提供了一种SAR图像像素点经纬度计算设备，设备包括处理器，以及与处理器通信连接的存储器。存储器存储有计算机可读指令；处理器执行存储器存储的计算机可读指令时，实现本申请实施例上述的SAR图像像素点经纬度计算方法。
本申请实施例还提供了一种计算机可读存储介质，计算机可读存储介质存储有计算机可读指令，计算机可读指令被处理器执行时实现上述的SAR图像像素点经纬度计算方法。
上述存储介质包括但不限于随机存取存储器（英文：Random Access Memory；简称：RAM）、只读存储器（英文：Read-Only Memory；简称：ROM）、缓存（英文：Cache）、硬盘（英文：Hard Disk Drive；简称：HDD）或者存储卡（英文：Memory Card）。
本说明书中的各个实施方式采用递进的方式描述，各个实施方式之间相同或相似的部分互相参见即可，每个实施方式重点说明的都是与其他实施方式的不同之处。本申请的全部或者部分可用于众多通用或专用的计算机系统环境或配置中。例如：个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、移动通信终端、多处理器系统、基于微处理器的系统、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
以上实施例仅用以说明本申请的技术方案，而非对本申请限制；尽管参照前述实施例对本申请进行了详细的说明，本领域普通技术人员应当理解：其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请技术方案的范围。

Claims (10)

1.一种SAR图像像素点经纬度计算方法，其特征在于，包括：
获取SAR图像参数；所述SAR图像参数包括已知点的像素坐标和经纬度、分辨率、中心方位角、航迹角、侧视方向、以及待求点的像素坐标；
根据所述SAR图像参数确定所述待求点与所述已知点的横纵坐标之间的实际向量距离；
根据椭球第一偏心率平方、地球赤道半径、所述已知点的经纬度和所述实际向量距离，确定所述待求点在地球坐标系中的坐标值；
根据椭球第一偏心率平方和所述待求点在地球坐标系中的坐标值，确定所述待求点的经纬度。
2.根据权利要求1所述的SAR图像像素点经纬度计算方法，其特征在于，所述根据所述SAR图像参数确定所述待求点与所述已知点的横纵坐标之间的实际向量距离，包括：
根据所述已知点的像素坐标、所述待求点的像素坐标和所述分辨率，确定所述待求点与所述已知点的横纵坐标之间的理论距离；
根据所述理论距离、所述中心方位角和所述侧视方向，确定所述待求点与所述已知点的横纵坐标之间的实际距离；
根据所述航迹角和所述实际距离，确定所述待求点与所述已知点的横纵坐标之间的实际向量距离。
3.根据权利要求2所述的SAR图像像素点经纬度计算方法，其特征在于，所述确定所述待求点与所述已知点的横纵坐标之间的理论距离，包括：
根据所述已知点的像素坐标和所述待求点的像素坐标，计算所述待求点与所述已知点的横纵坐标之间像素点个数；
根据所述分辨率、以及所述待求点与所述已知点的横纵坐标之间像素点个数，计算所述待求点与所述已知点的横纵坐标之间的所述理论距离。
4.根据权利要求3所述的SAR图像像素点经纬度计算方法，其特征在于，所述确定所述待求点与所述已知点的横纵坐标之间的实际距离，包括：
根据所述理论距离、所述中心方位角和所述侧视方向，计算所述待求点与所述已知点的横纵坐标之间的误差距离；
根据所述理论距离和所述误差距离，计算所述待求点与所述已知点的横纵坐标之间的所述实际距离。
5.根据权利要求1所述的SAR图像像素点经纬度计算方法，其特征在于，所述确定所述待求点在地球坐标系中的坐标值，包括：
根据椭球第一偏心率平方、地球赤道半径和所述已知点的经纬度，计算所述已知点在地球坐标系中的坐标值；
根据所述实际向量距离和所述已知点的经纬度，计算所述已知点与所述待求点在地球坐标系中的差值；
根据所述已知点在地球坐标系中的坐标值、以及所述已知点与所述待求点在地球坐标系中的差值，计算所述待求点在地球坐标系中的坐标值。
6.根据权利要求5所述的SAR图像像素点经纬度计算方法，其特征在于，所述计算所述已知点在地球坐标系中的坐标值，包括：
根据椭球第一偏心率平方、地球赤道半径和所述已知点的纬度，计算所述已知点与地心的距离；
根据所述已知点与地心的距离、以及所述已知点的经纬度，计算所述已知点在地球坐标系中的X、Y轴方向的坐标值；
根据所述已知点与地心的距离、椭球第一偏心率平方和所述已知点的纬度，计算所述已知点在地球坐标系中的Z轴方向的坐标值。
7.根据权利要求6所述的SAR图像像素点经纬度计算方法，其特征在于，所述确定所述待求点的经纬度，包括：
根据所述待求点在地球坐标系中的坐标值，计算距离地球坐标系Z轴的垂直距离，以及计算所述待求点的经度；
根据椭球第一偏心率平方、所述待求点在地球坐标系中的坐标值、以及所述待求点距离地球坐标系Z轴的垂直距离，计算所述待求点的纬度。
8.一种SAR图像像素点经纬度计算装置，其特征在于，包括：
获取模块，用于获取SAR图像参数；所述SAR图像参数包括已知点的像素坐标和经纬度、分辨率、中心方位角、航迹角、侧视方向、以及待求点的像素坐标；
距离计算模块，用于根据所述SAR图像参数确定所述待求点与所述已知点的横纵坐标之间的实际向量距离；
坐标计算模块，用于根据椭球第一偏心率平方、地球赤道半径、所述已知点的经纬度和所述实际向量距离，确定所述待求点在地球坐标系中的坐标值；
经纬计算模块，用于根据椭球第一偏心率平方和所述待求点在地球坐标系中的坐标值，确定所述待求点的经纬度。
9.一种SAR图像像素点经纬度计算设备，其特征在于，包括处理器，以及与所述处理器通信连接的存储器；
所述存储器存储有计算机可读指令；
所述处理器执行所述存储器存储的所述计算机可读指令时，实现如权利要求1-7中任一项所述的SAR图像像素点经纬度计算方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有计算机可读指令，所述计算机可读指令被处理器执行时实现如权利要求1-7中任一项所述的SAR图像像素点经纬度计算方法。
CN202310692703.4A 2023-06-13 2023-06-13 一种sar图像像素点经纬度计算方法及装置 Active CN116430346B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310692703.4A CN116430346B (zh) 2023-06-13 2023-06-13 一种sar图像像素点经纬度计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310692703.4A CN116430346B (zh) 2023-06-13 2023-06-13 一种sar图像像素点经纬度计算方法及装置

Publications (2)

Publication Number Publication Date
CN116430346A true CN116430346A (zh) 2023-07-14
CN116430346B CN116430346B (zh) 2023-10-17

Family

ID=87080095

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310692703.4A Active CN116430346B (zh) 2023-06-13 2023-06-13 一种sar图像像素点经纬度计算方法及装置

Country Status (1)

Country Link
CN (1) CN116430346B (zh)

Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6225942B1 (en) * 1999-07-30 2001-05-01 Litton Systems, Inc. Registration method for multiple sensor radar
EP2017647A1 (en) * 2007-07-19 2009-01-21 Consiglio Nazionale delle Ricerche Method for processing data sensed by a synthetic aperture radar (SAR) and related remote sensing system
CN102565797A (zh) * 2011-12-21 2012-07-11 北京航空航天大学 一种针对聚束模式星载sar图像的几何校正方法
CN102654576A (zh) * 2012-05-16 2012-09-05 西安电子科技大学 基于sar图像和dem数据的图像配准方法
CN107238824A (zh) * 2017-05-24 2017-10-10 西安电子科技大学 基于先验dem数据的星载sar图像几何精校正方法
CN107607947A (zh) * 2017-08-22 2018-01-19 西安电子科技大学 基于卡尔曼滤波的星载雷达成像参数在线估计方法
KR102028324B1 (ko) * 2019-02-26 2019-11-04 엘아이지넥스원 주식회사 영상 레이더의 영상 보정 방법 및 좌표 계산 방법
CN111652931A (zh) * 2020-04-09 2020-09-11 中国科学院空天信息创新研究院 一种地理定位方法、装置、设备及计算机可读存储介质
CN114325709A (zh) * 2022-03-14 2022-04-12 中国科学院空天信息创新研究院 多角度星载sar成像方法、装置、设备及介质
CN115685198A (zh) * 2022-10-16 2023-02-03 西安电子工程研究所 一种基于sar图像快速批量定位目标像素位置的方法
CN115932823A (zh) * 2023-01-09 2023-04-07 中国人民解放军国防科技大学 基于异源区域特征匹配的飞行器对地面目标定位方法
CN116224238A (zh) * 2023-03-08 2023-06-06 北京机电工程研究所 坐标系周期重建的sar多子区成像射频仿真方法及系统

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6225942B1 (en) * 1999-07-30 2001-05-01 Litton Systems, Inc. Registration method for multiple sensor radar
EP2017647A1 (en) * 2007-07-19 2009-01-21 Consiglio Nazionale delle Ricerche Method for processing data sensed by a synthetic aperture radar (SAR) and related remote sensing system
CN102565797A (zh) * 2011-12-21 2012-07-11 北京航空航天大学 一种针对聚束模式星载sar图像的几何校正方法
CN102654576A (zh) * 2012-05-16 2012-09-05 西安电子科技大学 基于sar图像和dem数据的图像配准方法
CN107238824A (zh) * 2017-05-24 2017-10-10 西安电子科技大学 基于先验dem数据的星载sar图像几何精校正方法
CN107607947A (zh) * 2017-08-22 2018-01-19 西安电子科技大学 基于卡尔曼滤波的星载雷达成像参数在线估计方法
KR102028324B1 (ko) * 2019-02-26 2019-11-04 엘아이지넥스원 주식회사 영상 레이더의 영상 보정 방법 및 좌표 계산 방법
CN111652931A (zh) * 2020-04-09 2020-09-11 中国科学院空天信息创新研究院 一种地理定位方法、装置、设备及计算机可读存储介质
CN114325709A (zh) * 2022-03-14 2022-04-12 中国科学院空天信息创新研究院 多角度星载sar成像方法、装置、设备及介质
CN115685198A (zh) * 2022-10-16 2023-02-03 西安电子工程研究所 一种基于sar图像快速批量定位目标像素位置的方法
CN115932823A (zh) * 2023-01-09 2023-04-07 中国人民解放军国防科技大学 基于异源区域特征匹配的飞行器对地面目标定位方法
CN116224238A (zh) * 2023-03-08 2023-06-06 北京机电工程研究所 坐标系周期重建的sar多子区成像射频仿真方法及系统

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CONG LI等: "Quantity Simulation Analysis Method based Novel RSA Timing Attack Algorithm for Single-Chip Microcomputer Platform", 2022 IEEE 22ND INTERNATIONAL CONFERENCE ON SOFTWARE QUALITY, RELIABILITY AND SECURITY *
WU ZHU等: "Performance Evaluation of Azimuth Offset Method for Mitigating the Ionospheric Effect on SAR Interferometry", HINDAWI JOURNAL OF SENSORS *
卫宇;: "考虑地球曲率情况下两点距离问题的求解", 航空兵器, no. 03 *
赵杰;林辉;: "网格坐标转换为雷达测量坐标的工程化算法", 火力与指挥控制, no. 07 *
郭庆;李明;: "机载条带SAR图像经纬度的计算方法", 现代雷达, no. 09 *

Also Published As

Publication number Publication date
CN116430346B (zh) 2023-10-17

Similar Documents

Publication Publication Date Title
CN110927708B (zh) 智能路侧单元的标定方法、装置及设备
US8264537B2 (en) Photogrammetric networks for positional accuracy
CN104902142A (zh) 一种移动终端视频的电子稳像方法
EP3136500A1 (en) Method, device and system for acquiring antenna engineering parameters
EP2759847A1 (en) Method and apparatus for determining equivalent velocity
US20230005109A1 (en) Remote Sensing Image Geometric Normalization Method and Apparatus
CN111460071A (zh) 高精地图的偏转方法、装置、设备及可读存储介质
CN114240805B (zh) 一种多角度sar动态成像检测方法及装置
CN111652931A (zh) 一种地理定位方法、装置、设备及计算机可读存储介质
CN116430346B (zh) 一种sar图像像素点经纬度计算方法及装置
Gu et al. Calibration of stereo rigs based on the backward projection process
CN115685202A (zh) 星载分布式干涉合成孔径雷达的基线估计方法及装置
CN114518585B (zh) 基于单颗卫星的目标定位方法、装置、终端及存储介质
CN112379344B (zh) 信号补偿方法及装置、设备、存储介质
CN114660927A (zh) 光电吊舱的控制方法和装置
CN115267711A (zh) 一种天气雷达天线指向性粗差校准方法及相关产品
CN112487532B (zh) 高精度自适应三维实时可视化打桩预警方法、装置及应用
CN114187332A (zh) 一种雷达图像配准方法和系统
CN113640799A (zh) 一种雷达波束中心照射点确定方法、装置及存储介质
CN113126058A (zh) 存储器、用于机载激光雷达系统的控制方法及装置
CN111681299B (zh) 基于InSAR解缠相位生成数字表面模型的方法及装置
CN113124816A (zh) 天线工参生成方法、装置、存储介质和计算机设备
CN113671501B (zh) 一种基于isar成像的方位向仿真校正方法和装置
CN114440885B (zh) 一种静止轨道遥感卫星定位方法及装置
CN117129956B (zh) 定位校正方法、装置、探测设备、计算机设备及存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant