版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (3972)
>虫友互识 (905)
>导师招生 (170)
>文献求助 (154)
>论文投稿 (146)
>硕博家园 (84)
>休闲灌水 (82)
>考博 (77)
>博后之家 (68)
>招聘信息布告栏 (62)
>论文道贺祈福 (56)
>基金申请 (48)
>公派出国 (48)
>教师之家 (46)
>考研 (35)
>找工作 (25)

返回列表

【奖励】本帖被评价5次，作者3s_studio增加金币 4.5 个

当前主题已经存档。

3s_studio

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 980.2
帖子: 266
在线: 40.1小时
虫号: 491779

[资源] 【分享】由经纬度坐标求距离的VB程序（附源代码）

经纬度坐标的地图中求距离的VB6.0程序，距离单位为M，利用大地正反算白塞尔公式编制而成。

下载：

程序及源代码下载

[ Last edited by 3s_studio on 2008-12-13 at 09:38 ]

回复此楼

» 猜你喜欢

祈福，祝自己今年好运已经有100人回复
关于会评已经有1人回复
地球化学论文润色/翻译怎么收费? 已经有256人回复
据说！出现这些情况，国基金会评本子会被拿下！已经有19人回复
无人打捞，直接陪跑已经有16人回复
定了定了，明天出结果！2025 国自然基金结果即将揭晓，查询攻略看这里已经有107人回复
2025国家自然科学基金放榜时间，根据官方回应，8月15日基本没戏了已经有5人回复
各位未来院士帮忙看看评审意见已经有9人回复
【招聘】科研助理+五险一金已经有18人回复
《风吹过日常的缝隙》已经有0人回复

1楼2008-06-23 23:33:47

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

cherishing

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 2977.3
帖子: 1340
在线: 270.3小时
虫号: 315141

★★★★★ 五星级,优秀推荐

谢谢共享，收藏！

回复此楼

2楼2008-06-24 09:53:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

cherishing

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 2977.3
帖子: 1340
在线: 270.3小时
虫号: 315141

但是楼主是否测试程序？你的内容是否不全？里面怎么是空的？？请解释一下好不？

赞一下

回复此楼

3楼2008-06-24 09:59:28

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

3s_studio

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 980.2
帖子: 266
在线: 40.1小时
虫号: 491779

cherishing，不好意思，我看了，这个文件确实不正确。我已经修改完毕，测试可以用，请到http://www.sciencenet.cn/blog/user_content.aspx?id=30185免费下载，谢谢。

[ Last edited by 3s_studio on 2008-6-24 at 12:31 ]

赞一下

回复此楼

4楼2008-06-24 10:50:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

cherishing

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 2977.3
帖子: 1340
在线: 270.3小时
虫号: 315141

好的，已经下载，测试可用！谢谢楼主分享！

赞一下

回复此楼

5楼2008-06-24 13:47:42

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

3s_studio

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 980.2
帖子: 266
在线: 40.1小时
虫号: 491779

呵呵，不用客气。谢谢大家关注！

赞一下

回复此楼

6楼2008-06-24 20:57:13

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

★★★★★ 五星级,优秀推荐

谢谢

回复此楼

7楼2008-06-25 09:00:44

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

serena84

应助: (幼儿园)
在线:
虫号: 575585

谢谢楼主哈

回复此楼

8楼2008-06-25 17:05:14

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

3s_studio

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 980.2
帖子: 266
在线: 40.1小时
虫号: 491779

光说不练，呵呵。

回复此楼

9楼2008-06-25 18:10:01

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

moonlight3988

木虫 (正式写手)

应助: 10 (幼儿园)
金币: 4120
帖子: 959
在线: 116.2小时
虫号: 566974

★★★ 三星级,支持鼓励

我有C#版的源码，谁需要？

赞一下

回复此楼

10楼2008-06-25 22:59:00

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

3s_studio

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 980.2
帖子: 266
在线: 40.1小时
虫号: 491779

哦，发我一份吧，我也看看你的算法。
邮箱 zjwang2008@gmail.com
谢谢。

赞一下

回复此楼

11楼2008-06-26 10:34:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

moonlight3988

木虫 (正式写手)

应助: 10 (幼儿园)
金币: 4120
帖子: 959
在线: 116.2小时
虫号: 566974

更详细源码请看http://bbs.i-location.org/thread-294-1-1.html

///

/// 贝赛尔大地问题正解
///

      /// 已知点纬度
      /// 已知点经度
      /// 大地线长
      /// 大地方位角
      /// 待求点纬度
      /// 待求点经度
      /// 大地反方位角
      /// 参考椭球长半轴
      /// 参考椭球扁率倒数
      public static void Bessel_PntSA_Pnt(double B1, double L1, double S, double A1, out double B2, out double L2, out double A2, double a, double f)
      {
         double ee = (2 * f - 1) / f / f;    //第一偏心率的平方
         double ee2 = ee / (1 - ee);          //第二偏心率的平方

         double u1 = Math.Atan(Math.Sqrt(1 - ee) * Math.Tan(B1 * Math.PI / 180));
         double m = Math.Cos(u1) * Math.Sin(A1);
         m = Math.Atan(m / Math.Sqrt(1 - m * m));

         double M_ = Math.Atan(Math.Tan(u1) / Math.Cos(A1));
         if (M_ < 0) M_ += Math.PI;

         double KK = ee2 * Math.Pow(Math.Cos(m), 2);
         double alpha = Math.Sqrt(1 + ee2) * (1 - KK / 4 + 7 * KK * KK / 64 - 15 * Math.Pow(KK, 3) / 256) / a;
         double beta = KK / 4 - KK * KK / 8 + 37 * Math.Pow(KK, 3) / 512;
         double gama = KK * KK * (1 - KK) / 128;

         double sigma, temp;
         sigma = alpha * S;
         temp = 0;
         while (Math.Abs(temp - sigma) > 0.0000000001)
         {
            temp = sigma;
            sigma = alpha * S + beta * Math.Sin(sigma) * Math.Cos(2 * M_ + sigma) + gama * Math.Sin(2 * sigma) * Math.Cos(4 * M_ + 2 * sigma);
         }

         A2 = Math.Atan(Math.Tan(m) / Math.Cos(M_ + sigma));
         if (Math.Cos(M_) > 0) A2 += Math.PI;
         if (m < 0 && Math.Cos(M_) < 0) A2 += 2 * Math.PI;

         double u2 = Math.Atan(-Math.Cos(A2) * Math.Tan(M_ + sigma));

         double lamda1;
         lamda1 = Math.Atan(Math.Sin(u2) * Math.Tan(A2));
         if (Math.Cos(M_) < 0) lamda1 += Math.PI;
         if (m < 0 && Math.Cos(M_) > 0) lamda1 += 2 * Math.PI;

         double lamda2;
         lamda2 = Math.Atan(Math.Sin(u1) * Math.Tan(A1));
         if (m > 0)
         {
            if (Math.Cos(M_ + sigma) < 0)
                  lamda2 += Math.PI;
            else if (Math.Sin(M_ + sigma) < 0)
                  lamda2 += 2 * Math.PI;
         }
         else
         {
            if (Math.Cos(2 * Math.PI - M_ - sigma) < 0)
                  lamda2 += Math.PI;
            else if (Math.Sin(M_ + sigma) < 0)
                  lamda2 += 2 * Math.PI;
         }

         B2 = Math.Atan(Math.Sqrt(1 + ee2) * Math.Tan(u2)) * 180 / Math.PI;

         KK = ee * Math.Pow(Math.Cos(m), 2);
         alpha = ee / 2 + ee * ee / 8 + Math.Pow(ee, 3) / 16 - ee * (1 + ee) * KK / 16 + 3 * ee * KK * KK / 128;
         gama = ee * KK * KK / 256;

         L2 = lamda1 + lamda2 - Math.Sin(m) * (alpha * sigma + beta * Math.Sin(sigma) * Math.Cos(2 * M_ + sigma) + gama * Math.Sin(2 * sigma) * Math.Cos(4 * M_ + 2 * sigma));
         L2 = L1 + L2 * 180 / Math.PI;
      }

      ///

/// 由两点大地坐标求解大地方位角
///

      /// 起点纬度
      /// 起点经度
      /// 末点纬度
      /// 末点经度
      /// 参考椭球长半轴
      /// 参考椭球扁率倒数
      /// 两点间的大地方位角
      public static double Bessel_BL_A(double B1, double L1, double B2, double L2, double a, double f)
      {
         if (L1 == L2) return 0;

         double ee = (2 * f - 1) / f / f;    //第一偏心率的平方

         double u1 = Math.Atan(Math.Sqrt(1 - ee) * Math.Tan(B1 * Math.PI / 180));
         double u2 = Math.Atan(Math.Sqrt(1 - ee) * Math.Tan(B2 * Math.PI / 180));
         double dL = (L2 - L1) * Math.PI / 180;

         double sigma = Math.Sin(u1) * Math.Sin(u2) + Math.Cos(u1) * Math.Cos(u2) * Math.Cos(dL);
         sigma = Math.Atan(Math.Sqrt(1 - sigma * sigma) / sigma);
         if (sigma <= 0) sigma += Math.PI;

         double m = Math.Cos(u1) * Math.Cos(u2) * Math.Sin(dL) / Math.Sin(sigma);
         m = Math.Atan(m / Math.Sqrt(1 - m * m));

         double KK = ee * Math.Pow(Math.Cos(m), 2);
         double alpha = ee / 2 + ee * ee / 8 + Math.Pow(ee, 3) / 16 - ee * (1 + ee) * KK / 16 + 3 * ee * KK * KK / 128;
         double lamda = dL + alpha * sigma * Math.Sin(m);

         sigma += Math.Sin(m) * (alpha * sigma * Math.Sin(m));

         m = Math.Cos(u1) * Math.Cos(u2) * Math.Sin(lamda) / Math.Sin(sigma);
         m = Math.Atan(m / Math.Sqrt(1 - m * m));

         double A = Math.Atan(Math.Sin(lamda) / (Math.Cos(u1) * Math.Tan(u2) - Math.Sin(u1) * Math.Cos(lamda)));
         if (A <= 0) A += Math.PI;
         if (m <= 0) A += Math.PI;

         double M_ = Math.Atan(Math.Sin(u1) * Math.Tan(A) / Math.Sin(m));
         if (M_ <= 0) M_ += Math.PI;

         KK = ee * Math.Pow(Math.Cos(m), 2);
         alpha = ee / 2 + ee * ee / 8 + Math.Pow(ee, 3) / 16 - ee * (1 + ee) * KK / 16 + 3 * ee * KK * KK / 128;
         double beta = ee * (1 + ee) * KK / 16 - ee * KK * KK / 32;
         lamda = dL + Math.Sin(m) * (alpha * sigma + beta * Math.Sin(sigma) * Math.Cos(2 * M_ + sigma));

         A = Math.Atan(Math.Sin(lamda) / (Math.Cos(u1) * Math.Tan(u2) - Math.Sin(u1) * Math.Cos(lamda)));
         if (A <= 0) A += Math.PI;
         if (m <= 0) A += Math.PI;

         return A * 180 / Math.PI;
      }

      ///

/// 由大地坐标计算大地线长
///

      /// 起点纬度
      /// 起点经度
      /// 末点纬度
      /// 末点经度
      /// 参考椭球长半轴
      /// 参考椭球扁率倒数
      /// 两点间的大地线长
      public static double Bessel_BL_S(double B1, double L1, double B2, double L2, double a, double f)
      {
         if (L1 == L2) return Math.Abs(MeridianLength(B1, a, f) - MeridianLength(B2, a, f));

         double ee = (2 * f - 1) / f / f;    //第一偏心率的平方
         double ee2 = ee / (1 - ee);          //第二偏心率的平方

         double u1 = Math.Atan(Math.Sqrt(1 - ee) * Math.Tan(B1 * Math.PI / 180));
         double u2 = Math.Atan(Math.Sqrt(1 - ee) * Math.Tan(B2 * Math.PI / 180));
         double dL = (L2 - L1) * Math.PI / 180;

         double sigma = Math.Sin(u1) * Math.Sin(u2) + Math.Cos(u1) * Math.Cos(u2) * Math.Cos(dL);
         sigma = Math.Atan(Math.Sqrt(1 - sigma * sigma) / sigma);
         if (sigma <= 0) sigma += Math.PI;

         double m = Math.Cos(u1) * Math.Cos(u2) * Math.Sin(dL) / Math.Sin(sigma);
         m = Math.Atan(m / Math.Sqrt(1 - m * m));

         double KK = ee * Math.Pow(Math.Cos(m), 2);
         double alpha = ee / 2 + ee * ee / 8 + Math.Pow(ee, 3) / 16 - ee * (1 + ee) * KK / 16 + 3 * ee * KK * KK / 128;
         double lamda = dL + alpha * sigma * Math.Sin(m);

         sigma += Math.Sin(m) * (alpha * sigma * Math.Sin(m));

         m = Math.Cos(u1) * Math.Cos(u2) * Math.Sin(lamda) / Math.Sin(sigma);
         m = Math.Atan(m / Math.Sqrt(1 - m * m));

         double A = Math.Atan(Math.Sin(lamda) / (Math.Cos(u1) * Math.Tan(u2) - Math.Sin(u1) * Math.Cos(lamda)));
         if (A <= 0) A += Math.PI;
         if (m <= 0) A += Math.PI;

         double M_ = Math.Atan(Math.Sin(u1) * Math.Tan(A) / Math.Sin(m));
         if (M_ <= 0) M_ += Math.PI;

         KK = ee * Math.Pow(Math.Cos(m), 2);
         alpha = ee / 2 + ee * ee / 8 + Math.Pow(ee, 3) / 16 - ee * (1 + ee) * KK / 16 + 3 * ee * KK * KK / 128;
         double beta = ee * (1 + ee) * KK / 16 - ee * KK * KK / 32;

         lamda = dL + Math.Sin(m) * (alpha * sigma + beta * Math.Sin(sigma) * Math.Cos(2 * M_ + sigma));
         sigma = Math.Sin(u1) * Math.Sin(u2) + Math.Cos(u1) * Math.Cos(u2) * Math.Cos(lamda);
         sigma = Math.Atan(Math.Sqrt(1 - sigma * sigma) / sigma);
         if (sigma < 0) sigma += Math.PI;

         KK = ee2 * Math.Pow(Math.Cos(m), 2);
         alpha = Math.Sqrt(1 + ee2) / a * (1 - KK / 4 + 7 * KK * KK / 64 - 15 * Math.Pow(KK, 3) / 256);
         beta = KK / 4 - KK * KK / 8 + 37 * Math.Pow(KK, 3) / 512;
         double gama = KK * KK * (1 - KK) / 128;

         double S = (sigma - beta * Math.Sin(sigma) * Math.Cos(2 * M_ + sigma) - gama * Math.Sin(2 * gama) * Math.Cos(4 * M_ + 2 * sigma)) / alpha;

         return S;
      }

[ Last edited by moonlight3988 on 2008-7-29 at 21:06 ]

赞一下

回复此楼

12楼2008-06-26 12:49:33

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖