以下是我所做的:
- 我在ARCMAP中创建了一个形状文件(多边形几何体((在其属性表中有相关信息(
- 我将shapefile转换为raster.tif文件
- 我设法从ArcMap中将我需要的所有信息提取到spyder(python3(中。因此,我有一个光栅tif图像、一个dbf文件和一个csv,我用它来提取和存储一些坐标
- 我把tif图像转换成一个数组
im = Image.open(tif_filepath)
imarray = np.array(im)
这样我就可以处理像素值了。
我的问题是:如何将矩阵(imarray
(的每个条目与相应投影中的实际坐标值相关联?
例如,我想要一个函数,它接受条目I,j并返回x,y,其中x和y是Arcmap中定义良好的坐标。
GDAL是解决您的问题的最佳选择。不仅仅是光栅格式。您可以使用GDAL将多边形作为shp添加到python中。你可以找到许多带有GDAL的文档。
我希望这对你有帮助。