在ArcGIS中查找栅格的中心点？

例如，如果我在ArcMap 10文档中加载了矩形且具有地理参考的栅格文件（.tiff w /关联的tfw），如何轻松找到其中心点并将其存储在点矢量层中？

另外，如果我的ArcMap文档中有多个这样的栅格，该如何将过程应用于所有栅格？

不幸的是，我有零Python经验。因此，以编程方式解决方案是可以的，但是我将需要有关如何将现有脚本加载到ArcGIS 10中以及如何在有问题的栅格上运行它的特定说明。（顺便说一句，栅格都在各自的图层中）。

虽然最初的问题是10.0，但是我将下面的代码更新为10.3.1。

将其复制粘贴到arcmap中的python窗口中以创建RasterCenter函数：

import arcpy, os
def RasterCenter(raster):
    #raster: string reference to raster
    raster = arcpy.Raster(raster)
    fcname = "{}_center".format(os.path.basename(str(raster)))
    x = raster.extent.XMin + (raster.extent.XMax - raster.extent.XMin)/2
    y = raster.extent.YMin + (raster.extent.YMax - raster.extent.YMin)/2
    featureclass = arcpy.CreateFeatureclass_management("in_memory", fcname, "POINT",spatial_reference = raster.spatialReference)
    with arcpy.da.InsertCursor(featureclass, ['SHAPE@XY']) as cursor:
        cursor.insertRow(((x, y),))
    mxd = arcpy.mapping.MapDocument("CURRENT")
    df = arcpy.mapping.ListDataFrames(mxd)[0]
    arcpy.MakeFeatureLayer_management(featureclass, fcname)
    layer = arcpy.mapping.Layer(fcname)
    arcpy.mapping.AddLayer(df, layer)

然后，您可以使用python窗口通过调用以下命令来创建要素类

RasterCenter("<reference to raster">)

因此，例如，如果您有一个名为DEM的栅格，则可以在python窗口中调用RasterCenter（“ dem”），它将添加一个名为“ dem_center”的图层，并且在栅格中心处有一个点。该图层存储在内存中，因此，如果要保留它，则将其导出。

要更进一步，您可以将脚本保存到.py文件，并将.py文件放置在python的搜索路径中。例如，将其另存为RasterCenter.py并将其放置在PYTHONPATH中（通常此位置为C：\ Python26 \ ArcGIS10.0 \ Lib）

然后，您可以执行以下操作：

import RasterCenter
RasterCenter.RasterCenter("<reference to raster">)

非常容易，简单地获取栅格属性并从最小，最大x和y得出中心点

MinX = arcpy.GetRasterProperties_management("raster", "LEFT")
MinY = arcpy.GetRasterProperties_management("raster", "BOTTOM")
MaxX = arcpy.GetRasterProperties_management("raster", "RIGHT")
MaxY = arcpy.GetRasterProperties_management("raster", "TOP")

centreX = (MaxX + MinX) / 2
centreY = (MaxY + MinY) / 2

和通常的错误检查等。

然后使用updateCursor添加到您的点表