使用Bresenham算法做视域分析时报Series对象无法转整数的TypeError

直接报错原因 #

你触发类型错误的核心原因是x1、y1的数据类型不对：

• 你通过x1 = obs['latitude']、y1 = obs['longitude']拿到的是Pandas的Series对象（哪怕观测点CSV里只有一行数据，取值默认也是序列而非单个数值），而range()函数要求输入整数，所以直接抛出类型错误。

其他需要修正的问题 #

• 你当前用的Bresenham算法仅支持整数栅格坐标输入，你现在传入的是地理坐标（浮点型的经纬度/投影坐标），哪怕解决了Series的问题，也会因为浮点值无法传入range()报错。你需要先把所有地理坐标映射为对应栅格的整数行列号之后，再调用算法。
• 你的bresenham(x1, y1, x2, y2)调用写在缓冲区点遍历循环的外部，实际只会处理最后一个缓冲区边界点，需要缩进放到for循环内部。
• 坐标系设置存在逻辑矛盾：EPSG:32620是UTM投影坐标系，单位为米，你设置的0.08缓冲区仅为8厘米，不符合常规分析需求；如果你的CSV存储的是WGS84经纬度坐标，初始CRS应该设为EPSG:4326，再重投影到EPSG:32620做缓冲区分析。
• 你当前实现的Bresenham算法仅支持x递增、斜率在0到1之间的线段，没有处理所有方向的线段情况，实际视域分析会出现漏线。

修正后的核心代码片段 #

if __name__=='__main__':
    # 取单个观测点的数值，注意先把地理坐标转成栅格行列号（整数）再传入算法
    x1 = int(obs['latitude'].iloc[0]) # 实际使用这里替换为对应栅格的行号
    y1 = int(obs['longitude'].iloc[0]) # 实际使用这里替换为对应栅格的列号

    for i in range(len(lat)):
        x2 = int(lat[i]) # 实际使用替换为缓冲区点对应的栅格行号
        y2 = int(long[i]) # 实际使用替换为缓冲区点对应的栅格列号
        # 缩进放到循环内部，处理每个缓冲区点
        bresenham(x1, y1, x2, y2)

额外优化建议 #

如果要实现完整的栅格视域分析，建议直接调用GDAL的视域分析接口gdal.ViewshedGenerate()，不需要自己从零实现Bresenham算法，稳定性和效率更高。

内容的提问来源于stack exchange，提问作者cer.as