已用您的方法计算出来任意两地间的距离，谢谢。想追问一个，用这样的方法计算出来的数据，有重复值。比如说：citycode1=110101 citycode2=110102 计算出来一个距离，citycode1=110102 citycode2=110101,又会计算出来一个相同的距离。我试过用duplicates drop dist,force。但是这样删掉的重复值不对。怎么仅保留一个距离呢？非常感谢！

您好，麻烦问一下，我想计算两组经纬度所代表地点之间的距离，stata应该怎么写程序啊？

非常感谢您的回答，不过我的问题更复杂些。id1和id2分别代表两组数据，怎么写程序求id2每个地方与id1中所有地方之间的距离？非常感激，愿意付出论坛币求教。

id1	lat1	lon1	id2	lat2	lon2
1	36.4345005273	115.9884387433	1	22.5408438779	114.1073437386
2	38.3044167994	116.8386942891	2	22.5739411063	114.1289712414
3	39.7268362964	116.3415800533	3	22.5388238085	114.1150023445
4	38.1639583261	114.6507245991	4	22.494806421	113.9253
5	34.8254568173	112.4451136144	5	22.5384661673	114.1197
6	39.1258031075	116.3916925454	6	22.5443198085	114.1113
7	37.1914757238	119.9634814343	7	22.5425299264	114.0873
8	38.0002309195	115.5600181231	8	22.5426343514	114.1229
9	34.8310757279	110.3291709124	9	22.5402387121	114.0911
10	36.2913151777	120.0203177897	10	22.5413547121	114.0942
11	29.2468388525	119.9110785558	11	22.5408483514	114.1196
12	24.9496872654	118.5135601716	12	22.4927325903	113.9183
13	35.8009112101	112.9183777177	13	22.5762051063	114.1295
14	36.2913151777	120.0203177897	14	22.5728342526	114.1441
15	34.8255486701	115.5421477633	15	22.5658913926	114.102
16	34.8406677240	119.1773057054	16	22.5384361647	113.9962
17	38.1461892667	114.5709685257	17	22.5762925302	114.1225
18	36.6845776316	118.4796830112	18	22.5419188085	114.1184

非常感谢您的回答，我的问题更复杂些。假设id1和id2代表两组地方经纬度，如何写程序求id2中每个地方与id2每个地方的距离。非常感谢，愿意付100币求教！
id1              lat1                              lon1              id2                lat2                        lon2
1        36.4345005273        115.9884387433        1        22.5408438779        114.1073437386
2        38.3044167994        116.8386942891        2        22.5739411063        114.1289712414
3        39.7268362964        116.3415800533        3        22.5388238085        114.1150023445
4        38.1639583261        114.6507245991        4        22.494806421        113.9252521
5        34.8254568173        112.4451136144        5        22.5384661673        114.1196634
6        39.1258031075        116.3916925454        6        22.5443198085        114.1112563
7        37.1914757238        119.9634814343        7        22.5425299264        114.0872721
8        38.0002309195        115.5600181231        8        22.5426343514        114.1228738
9        34.8310757279        110.3291709124        9        22.5402387121        114.0910752
10        36.2913151777        120.0203177897        10        22.5413547121        114.0942172
11        29.2468388525        119.9110785558        11        22.5408483514        114.1195538
12        24.9496872654        118.5135601716        12        22.4927325903        113.9183497
13        35.8009112101        112.9183777177        13        22.5762051063        114.1294792
14        36.2913151777        120.0203177897        14        22.5728342526        114.1440686
15        34.8255486701        115.5421477633        15        22.5658913926        114.102039
16        34.8406677240        119.1773057054        16        22.5384361647        113.996173
17        38.1461892667        114.5709685257        17        22.5762925302        114.1225164
18        36.6845776316        118.4796830112        18        22.5419188085        114.1183853

非常感谢您的回答，不过我的问题比较复杂，就是两组经纬度分别有很多数据，假设第一组有10组经纬度，第二组有10组，如果写程序求第一组每组经纬度与第二组所有经纬度之间的距离？非常感谢，愿意付100币求教！麻烦您了

不太明白 我在8楼的回答应该可以吧

id1 lat1 lon1 id2 lat2 lon2
  1   a1  b1    4   c1   d1
  2   a2  b2    5  c2    d2
  3   a3 b3     6  c3    d3
假设id1、id2分别代表两组经纬度，a、c分别纬度，b、d代表纬度，如何求id1中1、2、3分别与id2中所有地点4、5、6的距离，也就是求1、4之间距离，1、5之间距离，1、6之间距离；2、4之间距离,2、5之间距离,2、6之间距离；3、4之间距离，3、5之间距离，3、6之间距离？ 这样能明白吗？

非常感谢，不用您之前讲的公式了吗？就是这些公式gen double lat1 = 37 + (41 - 37) * uniform()
gen double lon1 = -109 + (109 - 102) * uniform()

那个是生成模拟数据的

明白了，非常感谢啊。不过可能是我的数据太大了，用stata的时候，运行cross using dat2，会出现no room to add more observations，不知道怎么办

最后的结果一般不需要算出来任意两点的距离，你要想清楚自己最终要的是什么

您在啊，还想麻烦您一个问题，就是按照下面的程序，我需要不停的按下去，您看能帮忙写个程序吗。上次您帮忙，欠您100论坛币，怎么给您，或者您随便发给东西标价100币，我付给您。非常感谢。
. geodist lat1 lon1  var105 var106, generate(h)

. geodist lat1 lon1  var107 var108, generate(j)

. geodist lat1 lon1  var109 var110, generate(k)

. geodist lat1 lon1  var111 var112, generate(l)

. geodist lat1 lon1  var113 var114, generate(z)

. geodist lat1 lon1  var115 var116, generate(x)

. geodist lat1 lon1  var117 var118, generate(c)

. geodist lat1 lon1  var119 var120, generate(v)

. geodist lat1 lon1  var121 var122, generate(b)

. geodist lat1 lon1  var123 var124, generate(n)

请问你的问题解决了吗？就是一个坐标对应同组内其他坐标的距离？到底用哪个代码实现呀？

这个楼主早解决了，不过我的问题不一样

您好，有个问题向你请教？我想实现一个地区内任意两两坐标组合的距离之和，请问该怎么实现呀？比如 A B C D E是同一个地区内的五个不同的坐标，我想得出A相对其他四个位置，B相对其他四个位置（AB AC AD AE），C以此类推，距离之和咋么实现呀？

同组内任意两个坐标之和吗？

按照deemg给写的就可以搞定的

geodist计算出来的结果是什么呀？千米吗？

你有什么问题啊？

我现在的问题有以下几个：第一，geodist计算的结果是什么呀？这个距离的单位是什么？为什么我自己用原始的三角函数计算公式编写的代码算出的结果跟这个程序对不上那？第二，把同样的数据分别列相同两行，然后执行geodist代码就可以算观测值跟任意一个坐标的距离之和了吗？

结果计算的是两地之间的距离，单位应该是公里。跟三角函数计算肯定有差距，估计这个是曲线的距离，因为地球是圆形的，但是差距应该不会太大。