使用循環量的加權插值法

Question

我正在做一個加權插值項目。 每個站點在地圖上都有一個坐標點，如下所示。

var stationCoor = [[408,352],[525,348],[535,495],[420,400],[272,145],[175,195],[197,335]];

我正在取位於湖中的點，並使用這些點為來自這些站的輸入創建加權平均值。 這是我確定加權數字的函數。

function findWeightSpeed(xPos, yPos){
   var totalHypt = 0;
   var arrHpyt = [];
   var arrWeight = [];
   for(var l=0;l<7;l++){
       var xDis = Math.abs(xPos-stationCoor[l][0]);
       var yDis = Math.abs(yPos-stationCoor[l][1]);
       var hptSq = Math.pow(xDis,2)+Math.pow(yDis,2);
       var hypt = Math.sqrt(hptSq);
       totalHypt = totalHypt+hypt;
       arrHpyt.push(hypt);
   }
   for(var j=0;j<7;j++){
       arrWeight.push(arrHpyt[j]/totalHypt)
   }
   return arrWeight;
}

這將找到點 (xPos,yPos) 和站之間的斜邊。 然后將數據相加並將每個站點除以產生加權數字的總數。

我需要使用這些點來衡量這些站的風向。 我正在使用下面的函數來計算點數的平均值。

function averageAngles(){
    var x = 0;
    var y = 0;
    var pi = 22/7;
    var angle = [2.7925,2.8797,2.9670,3.0543, 0.0872]; // 310,320,330,340,10
    for(var i = 0; i < angle.length; i++) {
        x += Math.cos(angle[i]);
        y += Math.sin(angle[i]);
    }
    var average_angle = Math.atan2(y, x);
    console.log((average_angle/pi)*360);
}

這為我提供了所有點的加權平均值為 0.20 的准確信息。 然而，對於7個站（如圖下面地圖上）的加權平均點類似於[0.1076839005418769, 0.08051796093187284, 0.003987308213631277, 0.08458358029618485, 0.2463427297217639, 0.26463834002675196, 0.21224618026791833]

我會如何做一個函數，從取加權平均數字findWeightSpeed()並用它來重量的圓弧數量在averageAngles()

我用了這個如何計算一組循環數據的平均值？ 制作角度平均函數。

非常感謝您提供的任何建議。

Answer 1

這是我在網上找到的一個鏈接，它解釋了整個過程。

計算風速和風向的加權平均值

代碼與此類似。

function weightAllData(xPos,yPos,windData){
                var uVecSum = 0;
                var vVecSum = 0;
                var arrayWeightSpeed = findWeightSpeed(xPos, yPos); //using weighted interpolation based on distance
                var arrayWindSpeed = [WSData];
                var arrayWindDirection = [WDData];
                for(var m=0;m<7;m++){
                    uVecSum = uVecSum + (arrayWeightSpeed[m] * getUVector(arrayWindSpeed[m],(arrayWindDirection[m]/180)*Math.PI));
                    vVecSum = vVecSum + (arrayWeightSpeed[m] * getVVector(arrayWindSpeed[m],(arrayWindDirection[m]/180)*Math.PI));
                }

                var weightWS = Math.sqrt(Math.pow(uVecSum,2)+Math.pow(vVecSum,2));
                if(vVecSum!=0){
                    weightWDRad = Math.atan(uVecSum/vVecSum);
                }
                if(vVecSum==0){
                    weightWDRad = Math.atan(uVecSum/(0.0001+vVecSum));
                }
                if(weightWDRad<0){
                    weightWDRad = weightWDRad + Math.PI
                }
                weightWD = (weightWDRad * (180/Math.PI));
}

如果您需要解釋，請告訴我