实验三：重力模型用于出行分布预测
一、实验的输入与输出文件
1、输入文件
? 地理文件：小区地理文件（njZone.dbd）、线类型地理文件
（njroad.dbd）；
? 矩阵文件：现状出行分布矩阵（BaseOD.mtx）； 2、输出文件
? 地理文件：质心地理文件（Center.dbd）； ? 网络文件：道路网络（njroad.net）；
? 矩阵文件：小区质心间最短出行时间阻抗矩阵（SPMAT.mtx）、出行分
布预测矩阵（CGRAV.MTX）；
? 数据表文件：重力模型参数标定输出数据表文件（Summary.bin）；
二、实验过程
1、创建质心连杆和交通网络
（1）创建质心地理文件
在TransCAD中打开njZone.dbd小区地理文件，选择Tools/Export菜单项，弹出Export Zone Geography对话框。在该对话框中，To后选择Standard Geographic File，ID Field后选择ZoneID，同时勾选Exports as Centroid points，点击OK进行确定，并将其保存为Center.dbd地理文件。
（2）质心连杆
点击顶部工具栏
按钮，弹出Layers对话框，点击Add layer按钮，分别
打开njroad.dbd和Center.dbd地理文件，选中njroad point图层，点击Show Layer按钮，完成后点击Close按钮，关闭对话框。
将njroad point 置为当前层，选择Dataview/Modify Table菜单项，在弹出的对话框中为节点图层地理文件添加一个名为Index（Integer型）的属性字段。
将njZone置为当前层，选择Tools/Map Editing/Connect菜单项，弹出Connect对话框。在该对话框中，选择Fill选项卡，在Node Field后的下拉列表框中选择Index，同时勾选IDs from Zone layer单选框，点击OK，完成小区质心与现有路网的连接。
（3）设置质心连接线属性
将njroad置为当前图层，点击
按钮，弹出路段图层的数据表文件，此
时选择Selection/Selection by Condition 菜单项，在弹出的对话框中键入
Capacity=null，然后点击OK，筛选出6条质心连接线。鼠标点击Capacity字段，选中Capacity一列，右键选择Fill命令，在弹出的对话框Single Value后输入100000。类似的方法为Time一列输入0.01，Speed一列输入40，关闭数据表视图，完成质心连接线的属性设置。
（4）创建道路网路
将njroad 置为当前层，选择Networks/Paths菜单下的Create命令，弹出Create Network对话框。在该对话框的左列表框中选中Capacity和Time两项，然后点击OK按钮，并将创建的网络文件保存为njroad.net。
2、生成小区间阻抗矩阵及矩阵索引转换
（1）小区间出行阻抗矩阵的生成
将njroad point 置为当前层，选择Selection/Select by Condition菜单项，在弹出的对话框中Enter a Condition后键入Index>0，点击OK按钮，筛选出质心节点。选择Networks/Paths菜单下的Multiple Paths 菜单项，弹出Multiple Shortest Paths对话框。在该对话框中，Minimize后下拉列表框中选择Time字段，From和To后均选择Selection字段，点击OK，并将该矩阵文件保存为SPMAT.mtx，即为小区间的最短出行时间矩阵。
（2）阻抗矩阵对角线0值处理
由上面步骤生成的阻抗矩阵对角线上的阻抗值（区内出行阻抗）为0，不能用于包含区内出行的出行分布过程，因而需要对其进行相应的处理。处理方法为将小区拟化为一个圆，圆的半径即为平均区内出行距离，进而求出小区内部的出行时间阻抗，将其填充到SPMAT.mtx矩阵的对角线中。（该值计算过程和方法参考小区内部出行阻抗计算.xlsx）
（3）阻抗矩阵索引转换
生成的阻抗矩阵的索引号与交通小区的编号并不一致，不能直接用于交通分布，因而需要对其进行转换。方法为在阻抗矩阵视图的某一位置单击鼠标右
键，选择Indices命令，弹出Matrix Indices对话框。点击Add Index按钮，弹出Add Add Index对话框，在Original Ro-Column IDs from框架下的Dataview后下拉列表框中选择njroad point 层，Field字段选择ID，在Row-Column IDs from框架下的Field字段选择Index，Selection后下拉列表框选择Selection。点击OK完成矩阵索引的添加并返回到Matrix Indices对话框。在Current Indices框架下Rows和Columns后的下拉列表框均选择New，完成后点击Close按钮，完成矩阵索引的转换。
