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

jinlong1835 发表于 2017-9-18 18:39
我现在的问题有以下几个：第一，geodist计算的结果是什么呀？这个距离的单位是什么？为什么我自己用原始的 ...

请问您的问题解决了吗，我也不知道我的计算结果为什么一直出现偏差

要在Stata中使用`geodist`命令生成400多个地区间的所有两两组合的距离（形成一个距离矩阵），你需要按照以下步骤操作：

### 步骤 1: 确保安装 geodist

首先，确保你已经安装了`geodist`命令。如果还没有安装，可以在Stata中输入：
```stata
ssc install geodist
```

### 步骤 2: 准备数据

假设你的数据集有以下结构：

- `code`: 地区代码（唯一标识）
- `longitude`: 经度
- `latitude`: 纬度

### 步骤 3: 创建距离矩阵

使用`geodist`命令计算所有组合的距离。但是，`geodist`不能直接生成一个完整的两两地区间距离的矩阵，它主要用于比较两个特定点之间的距离。

为了解决这个问题，我们可以创建一个新的数据集来存储所有可能的地区对及其对应的经纬度值，然后使用`geodist`命令计算每一对地区的距离。这可以通过下面的步骤实现：

1. **复制并合并数据**：首先，你需要将你的数据集复制一份，以确保可以进行自我连接（即每个地区与自己以及其它所有地区比较）。
```stata
preserve // 保存原始数据状态

// 复制数据集两次以创建所有组合
tempfile temp1 temp2
save `temp1', replace
append using `temp1'
save `temp2', replace

use `temp1'
merge 1:1 code longitude latitude using `temp2'
drop _merge // 删除合并信息标记变量
```

2. **计算距离**：现在，你有了一张包含所有可能地区对的数据集。使用`geodist`命令来计算每一对地区的距离。
```stata
geodist dist, lon1(longitude) lat1(latitude) lon2(_merge_longitude) lat2(_merge_latitude)
```
请注意，这里假设在复制并合并数据后，经纬度字段被命名为`_merge_longitude`和`_merge_latitude`。如果Stata没有自动重命名这些变量，请使用适当的变量名。

3. **保存结果**：最后，将计算出的距离数据集保存或导出为需要的格式。
```stata
save "distance_matrix.dta", replace // 保存为dta文件
```

4. **恢复原始状态**：
```stata
restore // 恢复到原始数据集的状态
```

### 步骤 4: 验证结果

在完成上述步骤后，你将得到一个包含所有地区对及其距离的数据集。你可以通过查询该数据集中任意地区的组合来确认计算是否正确。

以上步骤可以生成一个包含所有两两地区间距离的矩阵形式的数据集，但请注意，由于数据集可能会变得非常大（尤其是对于400多个地区），处理和存储这个完整的距离矩阵可能需要一定的资源。

此文本由CAIE学术大模型生成，添加下方二维码，优先体验功能试用