淺析太陽影子定位的數學模型

劉瀟（山東科技大學，山東青島 266590）

淺析太陽影子定位的數學模型

劉瀟
（山東科技大學，山東青島 266590）

目前，相機的應用非常普遍，從圖像或視頻中估算經緯度是目前計算機視覺領域的研究熱點問題之一。本文通過對問題的層層深入分析，引入太陽方位角、時角的概念，建立了物體與影子長度及太陽高度角之間的三角函數關系，而太陽高度角與緯度，太陽赤緯角和時角有關，所以確定了影子長度關于各個參數變化的數學模型，并分析了影子長度關于各個參數的變化規律。通過對某地直桿影子的頂點坐標進行分析，建立了數學模型來確定直桿所處的地理位置。

影子定位 數學模型 matlab軟件 拋物線擬合

1 概述

如何確定視頻的拍攝地點是視頻數據分析的重要方面，太陽影子定位技術就是通過分析視頻中物體的太陽影子變化，確定視頻拍攝的地點的一種方法。本文主要主要解決了兩個問題：

（1）本文通過引入各種與影子長度有關的參數，建立了一天中影子長度關于各個參數變化的數學模型，并分析了影子長度關于各個參數的變化規律。

（2）在太陽下豎立一根直桿，用相機錄下直桿的影子變化，通過分析影子的變化規律來確定直桿所在的地點。

本文中所用到的數學符號列表如表1。

2 建立影子長度變化規律的數學模型

表1

在探討影子長度與各因素之間的關系時，忽略陽光在大氣中的折射而引起的角度變化。由于光在同種均勻物質中延直線傳播，所以當光遇到不透明的物體時，就會在物體后面形成影子。物體形成影子的示意圖如圖1。

由三角形關系有：

H為太陽高度角，對于地球上的某個地點，太陽高度角是指太陽光的入射方向和地平面之間的夾角，專業上講太陽高度角是指某地太陽光線與該地作垂直于地心的地表切線的夾角。由于地球的自轉，造成太陽東升西落，太陽高度角在一日內不斷發生變化。太陽高度角隨著地方時和太陽的赤緯的變化而變化。太陽赤緯（與太陽直射點緯度相等）以δ表示，觀測地地理緯度用φ表示（太陽赤緯與地理緯度都是北緯為正，南緯為負），地方時（時角）以 Ω表示，有太陽高度角的計算公式：

式中，N表示積日，即日期在一年內的順序號；INT表示取整，即最終結果只保留整數部分即可，根據（3）、（4）、（5）、（6）式即可求得太陽赤緯。經過以上分析，影子長度關于各個參數變化的數學模型為：

由以上數學模型可知影響影子長度變化的參數為物體的高度L、物體所處的緯度 φ、太陽赤緯 δ以及物體所處地方的時角 Ω。對以上各參數分別求偏導得：

由以上對于各個參數的偏導函數，運用matlab對各偏導進行分析，我們可以得出影子長度關于各個參數的變化規律如下：

（1）當物體所處地點，日期以及所處地方的地方時都確定的條件下，物體的高度越高，在太陽照射下產生的影子越長。

（2）對于一個確定的物體，其高度是定值，當它所處的日期，以及所處地方的地方時都確定的條件下，它所處地方的緯度距離當時太陽直射點的距離越遠時，物體在太陽下形成的影子越長。

（3）對于一個高度確定的物體，當它所處的緯度以及地方時確定時，當太陽赤緯最小時即太陽赤緯為-23°26＇，南回歸線以南地區的地區影子長度達到最小，北回歸線以北地區影子長度達到最大；隨著太陽赤緯角的增大，南回歸線以南地區影子長度增長，北回歸線以北地區影子長度縮短。

（4）對于一個確定高度的物體，當它所處的緯度以及太陽赤緯確定時，即緯度與日期確定時，物體在一天中不同時刻形成影子的長度關于地方時的正午十二點呈對稱關系，在早上天將亮和傍晚天將黑時影子最長，從早上隨著時間的推移，影子不斷變短，并且縮短速度越來越慢，在正午十二點時影子長度最短，由中午到傍晚，影子長度又不斷拉長，拉長速度越來越快。

3 由影子頂點坐標變化求直桿所在地點

設固定直桿所在地點的經緯度為（φ，λ），以直桿底端為原點，水平地面 xoy為平面，直桿影長端點坐標記為 （x， y），記錄在一段時間內影子端點的多個坐標 （x1， y1）， （x2，y2）...（xn，yn），則影子長為

影子長是隨著時間變化的，利用在一段時間內不同時刻記錄的影子長度，對影子長度與時間進行拋物線擬合，得到影子長度與時間的函數關系如下：

接下來引入太陽方位角［1］的概念確定固定直桿所在地的緯度，太陽方位角即太陽所在的方位，指太陽光線在地平面上的投影與當地子午線的夾角，可近似地看作是豎立在地面上的直線在陽光下的陰影與正南方的夾角。方位角以目標物正北方向為零，順時針方向逐漸變大，其取值范圍是0°～360°。因此太陽方位角一般是以目標物的正北方向為起始方向，以太陽光的入射方向為終止方向，按順時針方向所測量的角度［2］。設太陽方位角為A則有

TT為真太陽時，TC為北京時間，CL為經度修正，如果地方子午圈在北京子午圈的東邊，則為正，反之為負，QE 為時差。

綜上，我們可以得到固定直桿所在地點經緯度的數學模型為：

4 結語

在太陽下豎立一根直桿，用相機拍攝下直桿影子變化的視頻。Matlab軟件工具箱中自帶的函數可以方便讀取AVI 文件的相關屬性和數據，利用這些屬性和數據可以提取出視頻文件中的任何一幀。然后針對每一幅圖像，采用二值化處理圖像求出分割閾值，進而求出最優閾值，進而求出影子的長度，利用所建立的數學模型對影子的頂點坐標進行處理，即可確定出直桿所處的地理位置。

［1］王國安.太陽高度角和日出日落時刻太陽方位角一年變化范圍的計算.氣象與環境科學，第30卷增刊:第163頁，2007年9月.

［2］汪和平.太陽方位角的公式求解及其應用.中國數學，2008年第9 版.高中版.

劉瀟（1994—），女，本科在讀，研究方向：機械設計制造及其自動化，就職于山東科技大學。