实验设计指南 > 空间填充设计 > 快速灵活填充设计 > 为地图形状创建空间填充设计
发布日期: 08/07/2020

为地图形状创建空间填充设计

可以使用快速灵活填充 (FFF) 设计来设计要为跨地图形状之类的区域执行的实验。例如,若您有两种样本收集方法(可能用于空气或土壤样本),希望在美国某州分布的各个位置进行比较,则考虑 FFF 设计。

假定相关的美国州是佐治亚州。您需要 40 个样本,每种抽样方法 20 个样本,遍布全州。若要为地图形状创建 FFF 设计,您需要一个具有地图边界纬度和经度的文件(请参见背景地图)。

在下例中,您将确定美国佐治亚州的经度和纬度,定义该州的形状,然后创建设计。该设计将为 40 个样本中的每一个样本提供位置和抽样方法。

查找佐治亚州的经度和纬度

1. 导航到 Maps 目录,通常会在此找到:

Windows: C:/Program Files/SAS/JMP/<版本>/Maps

Macintosh: /Library/Application Support/JMP/<版本>/Maps

2. 打开 US-State-XY.jmpUS-State-Name.jmp

请注意,在 US-State-Name.jmp 中,Shape ID 11 与美国佐治亚州对应。在本例中,这就是相关州。

3. US-State-XY.jmp 中,选择行 > 行选择 > 选择符合条件的行

4. 选择形状 ID 并在文本框中输入 11,然后点击确定

这将选定形状 ID 为 11(代表佐治亚州)的所有行。

5. 选择表 > 子集并点击确定以获取 US-State-XY.jmp 表中包含佐治亚州纬度和经度数据的子集(选定行)。

6. 关闭 US-State-XY.jmpUS-State-Name.jmp

获取值范围

要定义空间填充设计的空间,需要知道佐治亚州纬度和经度的值范围。然后,您可以使用这些值创建设计因子。

1. 选择文件 > 新建 > 脚本。在脚本窗口中,右击并选择显示嵌入的日志

2. 复制以下内容并将其粘贴到新脚本中:

dt = current data table();
mymap = dt << get as matrix({X,Y});
xx = mymap[0,1];
yy = mymap[0,2];
show(min(xx), max(xx));
show(min(yy), max(yy));
show(xx,yy);

3. US-State-XY.jmp 表的子集处于活动状态时,运行该脚本。

佐治亚州纬度和经度的值范围以及 xx yy 矩阵显示在日志中。

“范围查找”脚本和日志中的结果 

添加因子

1. 选择实验设计 > 特殊目的 > 空间填充设计

2. 向设计中添加两个连续因子。在添加因子数文本框中输入 2,然后点击连续型

3. 重命名因子经度和纬度。按如下方式输入最小值和最大值:

对于经度,输入 -85.6 和 -80.9 作为经度值。这些是日志中 x 的最小和最大经度值。

对于纬度,输入 30.4 和 35.0 作为纬度值。这些是日志中 y 的最小和最大纬度值。

4. 添加一个分类因子。点击分类 > 2 水平。重命名因子方法。

因子 

5. 点击继续

定义形状

要进一步定义空间填充设计的空间,需要使用多边形定义州的形状。

1. 选择使用不被允许的组合脚本打开“不被允许的组合表达式”文本框。

经度和纬度范围定义一个正方形。佐治亚州是一个多边形,而不是正方形。您想排除落在佐治亚州多边形之外的点。

2. 转至脚本窗口并从脚本日志中复制 xx yy 矩阵。将其粘贴到“不被允许的组合表达式”文本框。

3. 最后,添加以下行:

!In Polygon(longitude, latitude , xx, yy);

这将排除 xxyy 坐标所定义的多边形中未包含的点。

不被允许的组合表达式 

制作设计

1. 试验次数改为 40 并点击快速灵活填充

2. 点击制表

3. 从您的设计表中,选择图形 > 图形生成器

4. 将以下列拖入以下区域:

经度拖入 X 区域

纬度拖入 Y 区域

方法拖入颜色区域

5. 在图形之上,点击“平滑线”元素图标将其关闭。

6. 要添加州形状,右击该图形并选择图形 > 背景地图

7. 边界下,选择美国各州并点击确定

美国佐治亚州的 FFF 设计示例 

注意:由于 FFF 算法的随机起点的原因,地图上点的准确位置可能与所示位置不同。

地图上标识出您应该收集样本的 40 个地点。抽样方法(L1 或 L2)由点的颜色指示。

“图形”>“背景地图”窗口包含以下附加选项:

要确定用于测试的郡县,请在边界下选择美国各郡

要在街道一级绘制测试位置,请在图像下方选择街道图服务

需要更多信息?有问题?从 JMP 用户社区得到解答 (community.jmp.com).