使用ArcGIS for Desktop从30 x30 km间隔的点集中制作县级shapefile？

我有一个CSV格式的数据集，例如：

latitude, longitude, value
-45, 45, 10
....

我想使用ArcGIS或R + grass从这些数据计算县级均值。

我具有ArcGIS USA Counties图层，但是无法成功完成空间连接。

我怎样才能做到这一点？

— 安倍
source

Answers:

用于从ArcGIS中的xyz数据表中生成每个县中均值（z）的县级形状文件的工作流

请注意，这是基于@MLowry的建议，并添加了几乎逐步的说明。

打开ArcCatalog

创建新的个人地理数据库（例如foo.mdb）
文件->新->个人地理数据库
添加数据
- 右键单击map.csv->导出到地理数据库（单个）-> foo.mdb
- 右键单击counties.lyr->导出到地理数据库（单个）-> foo.mdb
从xy表创建要素类（.shp文件）
- 打开foo.mdb
- 右键单击map.csv表->从xytable创建要素类
- 输入字段x = lon，y = lat，z = yield
- 输入的坐标系的坐标导入->从counties.lyr（或等效地，->地理..->世界-> wgs_1984->打开->确定）
将要素类导入地理数据库（多个）
- 选择上一个创建的shp文件。步骤->添加->好的

打开ArcMap

新增资料
- 选择foo.mdb-> CTL +选择表->打开
空间连接
- ArcToolbox->分析工具->叠加->空间连接->
- 目标功能：counties.lyr
- 加盟特色：xyz表图
- 输出要素类：filename_spatialjoin
- 加入操作：JOIN_ONE_TO_ONE
- 匹配选项：最近
- 连接特征的字段图：删除不需要的字段（仅STATE_FIPS，CNTY_FIPS，必需的z值
- 搜索半径：30km
- 点击 OK
- Z值＝产率；选择“平均值”（或替代统计量）

— 大卫
source

@David感谢您的详细答案...将不得不检查出来。

— 日安倍晋三

@Abe随时进行编辑。我做到了CW

— David LeBauer 2011年

@David在这种情况下，IDW将为区域统计生成不正确的结果：请参阅@scw回复后的讨论。

— ub

@whuber。感谢您指出了这一点。当我结束时，我已经忘记了最初的要点是与县之间进行空间连接。固定。

— David LeBauer 2011年

@大卫谢谢。但是现在我很困惑。首先，在标题中，您实际上在工作流中的哪个位置创建了“栅格图层”？其次，它是如何，3层产生的？CSV文件仅表示一组点（x，y），每个点都有一个数字属性（z）。

— ub

听起来您可以将.csv转换为事件主题，然后导出为.shp，然后使用新的.shp进行空间连接到County层，而BOOM则很适合进行计算。

— 低级
source

一种通用方法是将规则间隔的点数据转换为栅格（在Arc * 中将 XYZ转换为栅格；在GRASS 中将v.in.xyz转换为栅格），然后执行区域统计运算符以汇总每个县内的像元值并计算统计信息（平均值是计算得出的标准统计数据之一）。要执行统计，请使用Arc *中的区域统计，或GRASS中的v.rast.stats。

如果您需要进行更高级的统计，则可以使用raster和来在R中进行分析，sp但难度更高。

— w
source

好主意，但是除非点数组的坐标正确定位，否则在转换为栅格时将发生重采样，并且您将受到重采样算法的支配：无论如何，平均值都可能会偏离一点点。

— whuber

正确，如果像元大小设置为粗略，可能会发生重采样。作为近似，您应该能够将像元大小设置为任意两点之间最小距离的一半，并避免出现此问题。这种方法确实具有必须处理重采样的缺点，但是我发现它比点数大（数万个或更多）要好得多。

— scw

无论单元格大小如何，都会进行@scw重采样（除非数据已经与网格完美对齐）。当使用更好的像元大小时，您可能会得到更差的答案，因为大多数网格值是在原始图像之间进行插值的。除非使用最近邻插值，否则这将导致极值显示不足。均值只会受到轻微影响；其他统计信息（例如标准差）可能会受到更大的影响。

— ub

@whuber我想我知道我们的误解是从哪里来的-我不主张将数据转换为连续的表面，而只是将点值（如果存在）放置到常规的栅格网格中。因此，绝大多数网格将保持NULL，只有在单元格位置与点位置一致的偶然值。

— scw

@scw啊，很好！谢谢你的澄清。实际上，您正在执行@MLowry建议的基于网格的空间连接。

— ub