C 高斯正算 高斯反算 高斯換帶等

2021-08-26 18:33:05 字數 977 閱讀 9295

首先,你要確定橢球引數:

a = 6378140;  //西安80橢球 iga75

e2 = 0.006694384999588;

m0 = a * (1 - e2);

m2 = 3.0 / 2 * e2 * m0;

m4 = 5.0 / 4 * e2 * m2;

m6 = 7.0 / 6 * e2 * m4;

m8 = 9.0 / 8 * e2 * m6;

a0 = m0 + m2 / 2 + (3.0 / 8.0) * m4 + (5.0 / 16.0) * m6 + (35.0 / 128.0) * m8;

a2 = m2 / 2 + m4 / 2 + 15.0 / 32 * m6 + 7.0 / 16 * m8;

a4 = m4 / 8 + 3.0 / 16 * m6 + 7.0 / 32 * m8;

a6 = m6 / 32 + m8 / 16;

a8 = m8 / 128;

xx = 0;

yy = 0;

_x = 0;

_y = 0;

bb = 0;

ll = 0;

下面才開始正題:

高斯正算:把經緯度座標轉換為平面座標

void gausspositive(double b, double l, double l0)

高斯反算:把平面座標轉換為經緯度座標:

void gaussnegative(double x, double y, double l0)

裡面涉及到的弧長反算:

double hcfansuan(double px)

return bf1;

}

高斯換帶就比較簡單了:

void gausszone(double x, double y, double l0, double l0new)

高斯座標自然值計算 高斯模糊

之前看馮樂樂女神的 unity shader入門精要 高斯模糊這一節只是大概略過,昨晚重新看了一下,感受頗深,隨即搞了乙個效果試一下。老規矩先上圖 原圖 高斯模糊 毛玻璃的既視感 清晰的,各個畫素之間會有明顯的過渡,而如果各個畫素之間的差距不是很大,那麼影象就會模糊了 鑑於這個定義,我們就可以通過 ...

計算高斯積分與插值

print np.sqrt np.pi integrate.quad lambda x np.exp x 2 np.inf,np.inf 差值即在資料集已知資料點之間填補空白 scipy.interpolate函式可以根據實驗資料進行插值。interpld類可以建立線性插值或者三次差值 import...

用矩陣解方程式的運算 高斯消元法

x 2y 4z 7 3x 7y 2z 11 2x 3y 3z 1 我們暫且將 矩陣看為多元方程組,將新知識與已有知識建立鏈結是記憶的好方法 x 2y 4z7 3x 7y 2z 11 2x 3y 3z1 上面的方程式寫成矩陣,左邊是方程式,第一列代表xi的係數,第二列代表y,第三列代表z,矩陣中最右邊...