[问题] 矩阵运算问题

楼主: ncuatm (VC #15)   2012-04-09 17:16:55
我现在有分布在经纬度上的数据资料,一个经纬度可能会对应到很多组资料或都没有资料
目前有使用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出来的结果依然是原始给定的零矩阵,请问该如何解决?

Links booklink

Contact Us: admin [ a t ] ucptt.com