scanner.py 文件源码

def sphere2basemap(map, azimuthangle_polarangle_radialdistance):

    ## PT axis should not be givien to this, not same convention!

    azimuth = -1* (360 - (450. - np.rad2deg(azimuthangle_polarangle_radialdistance[0]))) # 450-
    takeoff =  90. - np.rad2deg(azimuthangle_polarangle_radialdistance[1]) #90. -
    #radius = azimuthangle_polarangle_radialdistance[2]

    while len(takeoff[takeoff>90.])>0 or len(takeoff[takeoff<-90.])>0 :

        azimuth[takeoff <-90.] += 180.
        takeoff[takeoff <-90.] = -180 - takeoff[takeoff <-90.]

        azimuth[takeoff  >90.] += 180.
        takeoff[takeoff  >90.] = 180 - takeoff[takeoff >90.]


    while len(azimuth[azimuth>360.])>0 or len(azimuth[azimuth<0.])>0 :

        azimuth[azimuth   <0.] += 360.
        azimuth[azimuth >360.] -= 360.


    return map( azimuth, takeoff )