k2q.m

来自「StaMps最新测试版」· M 代码 · 共 15 行

function [q]=K2q(K,range_pixel)
%[q]=K2q(K)
%convert from K to topo error
if nargin<2
    range_pixel=2300;
end
lambda =0.0565647;
rho = 831265+range_pixel*7.9048902811596;
h=794141;
re=6345245;
%load parms

alpha=pi-acos((rho.^2+re^2-(re+h)^2)/2./rho/re);

q=K*lambda.*rho.*sin(alpha)/4/pi;