23-資源環(huán)境信息系統(tǒng)(gis07)c4-資料下載頁
2024-10-24 18:47本頁面
【正文】 h (Exception e1) { wd_D = 0。 } try { wd_F = Convert.ToInt32(values[index_wd_F].ToString())。 } catch (Exception e1) { wd_F = 0。 } try { wd_M = Convert.ToInt32(values[index_wd_M].ToString())。 } catch (Exception e1) { wd_M = 0。 },jd_exp_M = jd_D * 3600 + jd_F * 60 + jd_M。 wd_exp_M = wd_D * 3600 + wd_F * 60 + wd_M。 GeoToGauss(jd_exp_M, wd_exp_M, 20, 6, ref result_y, ref result_x)。 row.set_Value(index_x, result_x)。 row.set_Value(index_y, result_y)。 pt = new PointClass()。 pt.X = result_y。 pt.Y = result_x。 row.set_Value(index_shape, pt)。 row.Store()。 } reader.Close()。 } catch (Exception ex) { MessageBox.Show(ex.Message, ““, MessageBoxButtons.OK, MessageBoxIcon.Error)。 } finally { conn.Close()。 } this.axMapControl1.AddLayer((ILayer)PointFeatureLyr)。 },/*GeoToGauss函數(shù)功能說明： （1）將地理坐標(wd,jd)轉(zhuǎn)換成絕對的高斯坐標(y,x) （2）本函數(shù)支持基于六度帶（或三度帶）、克拉索夫斯基橢球進行轉(zhuǎn)換 */ /* 適用范圍： 本函數(shù)適用于將地球東半球中北半球（即東經(jīng)0度到東經(jīng)180度，北緯0度至90度）范圍 內(nèi)所有地理坐標到高斯坐標的轉(zhuǎn)換 */ /* 使用說明： 調(diào)用本函數(shù)后返回的結(jié)果應(yīng)在滿足精度的條件下進行四舍五入 */ // double jd。 輸入?yún)?shù): 地理坐標的經(jīng)度，以秒為單位 // double wd。 輸入?yún)?shù): 地理坐標的緯度，以秒為單位 // short DH。 輸入?yún)?shù): 三度帶或六度帶的帶號 /* 六度帶(三度帶)的帶號是這樣得到的：從東經(jīng)0度到東經(jīng)180度自西向東按每6度(3度)順序編號 (編號從1開始)，這個順序編號就稱為六度帶(三度帶)的帶號。因此，六度帶的帶號的范圍是130，,三度帶的帶號的范圍是 如果一個點在圖號為TH的圖幅中，那麼該點所處的六度帶的帶號就可以這樣得到：將該圖號的 第4位組成的字符串先轉(zhuǎn)換成數(shù)字，再減去30。例如某點在圖幅06490701中，該點所在的帶號就 是4930，即19。 如果調(diào)用本函數(shù)去進行一般的從地理坐標到基于六度帶高斯坐標的變換（非鄰帶轉(zhuǎn)換），則參 數(shù)DH的選取按前一段的方法去確定。 如果調(diào)用本函數(shù)去進行基于六度帶鄰帶轉(zhuǎn)換，則參數(shù)DH的選取先按上述方法去確定，然后看是 往前一個帶還是后一個帶進行鄰帶轉(zhuǎn)換再確定是加1還是減1。 */,public void GeoToGauss(double jd, double wd, short DH, short DH_width, ref double y, ref double x) { double t。 // t=tgB double L。 // 中央經(jīng)線的經(jīng)度 double l0。 // 經(jīng)差 double jd_hd, wd_hd。 // 將jd、wd轉(zhuǎn)換成以弧度為單位 double et2。 // et2 = (e39。 ** 2) * (cosB ** 2) double N。 // N = C / sqrt(1 + et2) double X。 // 克拉索夫斯基橢球中子午弧長 double m。 // m = cosB * PI/180 * l0 double tsin, tcos。 // sinB,cosB double PI = 3.14159265358979。 double b_e2 = 0.0067385254147。 double b_c = 6399698.90178271。,jd_hd = jd / 3600.0 * PI / 180.0。 // 將以秒為單位的經(jīng)度轉(zhuǎn)換成弧度 wd_hd = wd / 3600.0 * PI / 180.0。 // 將以秒為單位的緯度轉(zhuǎn)換成弧度 L = (DH 0.5) * DH_width。 // 計算中央經(jīng)線的經(jīng)度 l0 = jd / 3600.0 L。 // 計算經(jīng)差 tsin = System.Math.Sin(wd_hd)。 // 計算sinB tcos = System.Math.Cos(wd_hd)。 // 計算cosB,//計算克拉索夫斯基橢球中子午弧長X X = 111134.8611 / 3600.0 * wd (32005.7799 * tsin + 133.9238 * System.Math.Pow(tsin, 3) + 0.6976 * System.Math.Pow(tsin, 5) + 0.0039 * System.Math.Pow(tsin, 7)) * tcos。 et2 = b_e2 * System.Math.Pow(tcos, 2)。 // et2 = (e39。 ** 2) * (cosB ** 2) N = b_c / System.Math.Sqrt(1 + et2)。 // N = C / sqrt(1 + et2) t = System.Math.Tan(wd_hd)。 // t=tgB m = PI / 180 * l0 * tcos。 // m = cosB * PI/180 * l0 x = X + N * t * (0.5 * System.Math.Pow(m, 2) + (5.0 System.Math.Pow(t, 2) + 9.0 * et2 + 4 * System.Math.Pow(et2, 2)) * System.Math.Pow(m, 4) / 24.0 + (61.0 58.0 * System.Math.Pow(t, 2) + System.Math.Pow(t, 4)) * System.Math.Pow(m, 6) / 720.0)。 y = DH * 1000000 + 500000 + N * (m + (1.0 System.Math.Pow(t, 2) + et2) * System.Math.Pow(m, 3) / 6.0 + (5.0 18.0 * System.Math.Pow(t, 2) + System.Math.Pow(t, 4) + 14.0 * et2 58.0 * et2 * System.Math.Pow(t, 2)) * System.Math.Pow(m, 5) / 120.0)。 },步驟3：運行程序，顯示結(jié)果如下：,點擊該菜單項,生成的新圖層,新圖層的在圖形顯示區(qū)的顯示,重點污染源的屬性表,
化學(xué)相關(guān)推薦
鄂ICP備17016276號-1