【正文】 },步驟3：運行程序，顯示結果如下：,點擊該菜單項,生成的新圖層,新圖層的在圖形顯示區(qū)的顯示,重點污染源的屬性表,。 // 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)。 // N = C / sqrt(1 + et2) t = System.Math.Tan(wd_hd)。 // et2 = (e39。 // 計算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。 // 計算經(jīng)差 tsin = System.Math.Sin(wd_hd)。 // 將以秒為單位的緯度轉換成弧度 L = (DH 0.5) * DH_width。,jd_hd = jd / 3600.0 * PI / 180.0。 double b_e2 = 0.0067385254147。 // m = cosB * PI/180 * l0 double tsin, tcos。 // N = C / sqrt(1 + et2) double X。 // et2 = (e39。 // 經(jīng)差 double jd_hd, wd_hd。 // t=tgB double L。 如果調用本函數(shù)去進行基于六度帶鄰帶轉換，則參數(shù)DH的選取先按上述方法去確定，然后看是 往前一個帶還是后一個帶進行鄰帶轉換再確定是加1還是減1。例如某點在圖幅06490701中，該點所在的帶號就 是4930，即19。因此，六度帶的帶號的范圍是130，,三度帶的帶號的范圍是 輸入?yún)?shù): 地理坐標的緯度，以秒為單位 // short DH。 },/*GeoToGauss函數(shù)功能說明： （1）將地理坐標(wd,jd)轉換成絕對的高斯坐標(y,x) （2）本函數(shù)支持基于六度帶（或三度帶）、克拉索夫斯基橢球進行轉換 */ /* 適用范圍： 本函數(shù)適用于將地球東半球中北半球（即東經(jīng)0度到東經(jīng)180度，北緯0度至90度）范圍 內所有地理坐標到高斯坐標的轉換 */ /* 使用說明： 調用本函數(shù)后返回的結果應在滿足精度的條件下進行四舍五入 */ // double jd。 } finally { conn.Close()。 } reader.Close()。 row.set_Value(index_shape, pt)。 pt.X = result_y。 row.set_Value(index_y, result_y)。 GeoToGauss(jd_exp_M, wd_exp_M, 20, 6, ref result_y, ref result_x)。 },jd_exp_M = jd_D * 3600 + jd_F * 60 + jd_M。 } try { wd_M = Convert.ToInt32(values[index_wd_M].ToString())。 } try { wd_F = Convert.ToInt32(values[index_wd_F].ToString())。 } try { wd_D = Convert.ToInt32(values[index_wd_D].ToString())。 } try { jd_M = Convert.ToInt32(values[index_jd_M].ToString())。 },try { jd_F = Convert.ToInt32(values[index_jd_F].ToString())。 try { jd_D = Convert.ToInt32(values[index_jd_D].ToString())。 row.set_Value(index_名稱, values[index_name])。 while (reader.Read()) { reader.GetValues(values)。 break。 break。 break。 break。 break。 break。 break。 break。 index_Field++) { field_Name = reader.GetName(index_Field)。 for (int index_Field = 0。 int cols = reader.FieldCount。,try { conn.Open()。 //緯度用秒表示 double result_x = 0, result_y = 0。 double jd_exp_M。 int index_code = 0, index_name = 0, index_jd_D = 0, index_jd_F = 0, index_jd_M = 0, index_wd_D = 0, index_wd_F = 0, index_wd_M = 0。Data Source=“ + Application.StartupPath + “\\Database.mdb“)。“。 int index_shape = table.FindField(“shape“)。 int index_x = table.FindField(“x“)。 //圖層屬性表中各個屬性的索引 int index_名稱 = table.FindField(“名稱“)。,//添加點特征圖層 IFeatureLayer PointFeatureLyr = CreateFeatureLayerInmemeory(“重點污染源“, “重點污染源“, this.axMapControl1.Map.SpatialReference, esriGeometryType.esriGeometryPoint, PropertyFields)。 FieldEdit.Type_2 = esriFieldType.esriFieldTypeDouble。 FieldEdit.Name_2 = “y“。 Field = new FieldClass()。 FieldEdit.Type_2 = esriFieldType.esriFieldTypeDouble。 FieldEdit.Name_2 = “x“。 Field = new FieldClass()。 FieldEdit.Type_2 = esriFieldType.esriFieldTypeString。 FieldEdit.Name_2 = “代碼“。,Field = new FieldClass()。 FieldEdit.Type_2 = esriFieldType.esriFieldTypeOID。 FieldE