基于Python在Spyder中获取栅格文件内商店与房屋的像素坐标

嘿，我来帮你搞定这个问题！在Spyder里用Python处理这个300×300的栅格文件其实很简单，咱们用numpy库就能轻松定位商店和房屋的坐标，具体步骤如下：

1. 导入必要的库 #

首先，我们需要导入numpy——它处理数值数组和栅格数据特别顺手。在Spyder的控制台或脚本里输入：

import numpy as np

2. 读取栅格文件 #

接下来把你的文本栅格文件加载成numpy数组。假设你的文件每行是一行栅格数据，数值用空格分隔（如果是逗号分隔，记得加delimiter=','参数）：

# 替换成你的栅格文件实际路径，比如'C:/data/grid.txt'或者直接文件名（如果在同目录）
grid = np.loadtxt('your_grid_file.txt')

3. 查找商店的坐标 #

商店对应数值是1，我们用np.where()函数找到所有值为1的位置。注意numpy数组的索引是行在前，列在后，如果你的栅格左上角是(0,0)原点，列对应x坐标，行对应y坐标，所以需要转换一下：

# 获取商店的行、列索引
shop_rows, shop_cols = np.where(grid == 1)
# 转换成(x, y)坐标列表
shop_coords = list(zip(shop_cols, shop_rows))

# 打印所有商店坐标
print("所有商店的坐标：")
for coord in shop_coords:
    print(coord)

4. 查找房屋的坐标 #

房屋对应数值是10到250之间，我们用范围条件来筛选：

# 获取房屋的行、列索引
house_rows, house_cols = np.where((grid >= 10) & (grid <= 250))
# 转换成(x, y)坐标列表
house_coords = list(zip(house_cols, house_rows))

# 打印所有房屋坐标
print("\n所有房屋的坐标：")
for coord in house_coords:
    print(coord)

一些注意事项 #

• 坐标原点问题：如果你的栅格系统原点在右下角，那么y坐标需要反转，比如300 - 1 - shop_rows（因为300行的索引是0到299）。
• 文件格式问题：如果你的文本文件有注释行，比如开头几行是说明，可以用np.loadtxt(..., skiprows=N)跳过前N行。
• 路径问题：如果文件不在Spyder当前工作目录，一定要写完整的绝对路径，避免找不到文件。

内容的提问来源于stack exchange，提问作者A. Kor