怎么用stata geodist命令算出任意两地间的距离？

我有400多个地区的经纬度，数据变量分别是，地区，代码，经度，纬度。如何用stata geodist计算出任意两地之间的距离。就是需要创建出一个矩阵，里面有任意的两地之间的距离。谢谢大神！亲测有效追增论坛币！

1. * generate simulated data;
   
2. 
   
3. clear
   
4. set seed 1234
   
5. set obs 10
   
6. gen id1 = _n
   
7. gen double lat1 = 37 + (41 - 37) * uniform()
   
8. gen double lon1 = -109 + (109 - 102) * uniform()
   
9. 
   
10. * pairwise combination
    
11. preserve
    
12. rename (id1 lat1 lon1) (id2 lat2 lon2)
    
13. save geo_data, replace
    
14. restore
    
15. 
    
16. cross using geo_data
    
17. 
    
18. * 10*10 distance matrix
    
19. geodist lat1 lon1 lat2 lon2, generate(v)
    
20. keep id1 id2 v
    
21. reshape wide v, i(id1) j(id2)

复制代码

deem 发表于 2017-8-28 10:40

deem 发表于 2017-8-28 10:40

请问，能根据我的数据样子帮我编写一下吗？非常感谢！

雯雯517 发表于 2017-8-28 15:11
请问，能根据我的数据样子帮我编写一下吗？非常感谢！

你改下变量名就行了呀 这段代码挺清晰呀

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

deem 发表于 2017-8-28 16:17
你改下变量名就行了呀 这段代码挺清晰呀

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

雯雯517 发表于 2017-8-28 16:54
已用您的方法计算出来任意两地间的距离，谢谢。想追问一个，用这样的方法计算出来的数据，有重复值。比 ...

1. * generate data;
   
2. input str10 citycode longitude latitude
   
3. 110101 116 39
   
4. 110102 142 73
   
5. 110103 25 75
   
6. end
   
7. 
   
8. * make index var;
   
9. gen id2 = _n
   
10. save gist_data, replace
    
11. expand=_N-_n+1
    
12. sort citycode
    
13. bysort citycode: gen id1 = _n
    
14. rename (longitude latitude) (longitude2 latitude2)
    
15. 
    
16. merge m:1 id2 using gist_data
    
17. 
    
18. geodist longitude latitude longitude2 latitude2, generate(v)

复制代码

deem 发表于 2017-8-28 09:34

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

08zjliu 发表于 2017-9-17 12:26
您好，麻烦问一下，我想计算两组经纬度所代表地点之间的距离，stata应该怎么写程序啊？

1. geodist longitude latitude longitude2 latitude2, generate(v)

复制代码