我現在有分佈在經緯度上的數據資料,一個經緯度可能會對應到很多組資料或都沒有資料
目前有使用MatrixReal模組,部份程式碼如下:
for ($i=1; $i<=$line; $i++) {
$count_matrix = new Math::MatrixReal(180,360);
$vtec_matrix = new Math::MatrixReal(180,360);
$grid_x[$i]=1+int($lon[$i]+180);
$grid_y[$i]=1+int(90-$lat[$i]);
$count_matrix->assign($grid_y[$i],$grid_x[$i],1);
$vtec_matrix->assign($grid_y[$i],$grid_x[$i],$vtec[$i]);
}
其中$vtec[$i]是我的資料,$count_matrix是用來計算在每個網格上的資料筆數,
$vtec_matrix是要將所有資料加起來,最後再將$vtec_matrix資料總和除以$count_matrix
相對應的資料筆數,得到全球網格的平均數據。目前的問題是無法對$count和$vtec這兩個
零矩陣內元素做運算,print出來的結果依然是原始給定的零矩陣,請問該如何解決?