首先,你要確定橢球引數:
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,矩陣中最右邊...