Python中利用numpy库创建矩阵，并根据DS,的分辨率和范围计算出矩阵的格网大小和数量来建立一个与DSM大小相同的矩阵，以便储存DSM中各个格网点的可见性

代码示例：

import numpy as np

# 假设DSM的分辨率为1米/像素，范围为1000*1000米
dsm_resolution = 1  # 分辨率，单位：米/像素
dsm_range = [1000, 1000]  # 范围，单位：米

# 计算矩阵的格网大小和数量
grid_size = dsm_resolution  # 格网大小，与分辨率相同
grid_num = [int(np.ceil(dsm_range[0] / grid_size)), int(np.ceil(dsm_range[1] / grid_size))]  # 格网数量

# 创建矩阵
dsm_matrix = np.zeros(grid_num)

# 示例：在矩阵中设置第(100, 200)个格网点的值为1
dsm_matrix[100, 200] = 1