跟踪太阳的经度和纬度

Question

我有一张4096 * 4096的太阳图像。我想用MATLAB跟踪它的纬度和经度。我试过下面的代码。虽然它在轴上工作正常（纬度= 0/90经度= 0/90），但它在边缘处失败（例如纬度= 30经度= 90）

x=round((b*s*cos((l*pi)/180))+xC)
y=round((-l*s)+yC)

where x=horizontal displacement
y=vertical displacement
l=latitude
b=longitude
(xC,yC)=center coordinates of the Sun
s=scaling factor.

提前致谢 AU

我正在绘制给定纬度的所有经度点。图1是纬度= 0而图2和3分别是纬度= 45和-45。随着我增加纬度，线变得更偏斜（纬度= 70时见图4）。

请参阅以下链接以获取图片：https://drive.google.com/drive/folders/0B2X4xzsGcarCbGM5aVM1UXBRX2M?usp=sharing

我还尝试绘制纬度以获得恒定的经度和类似的趋势（经度增加时的偏度）。