土地荒漠化敏感度监测

[「已注销」]

用机器学习的方法进行逻辑回归分析，如何编写程序实现基于地理位置的栅格数据的条件读写呢？
也就是说整幅图上有不规则的分布的点，这些点有位置信息，我们如何根据这一位置信息来读相应位置的栅格图像像元的数值并写入一个新的文本文档？

[cobinzhang]

对于一个有投影的地理栅格数据，如果有栅格点的定位方法，比如经纬度，必然可以找到该栅格点的像素值。

gt=ras_ds.GetGeoTransform()
mx,my 为经纬度
px = int((mx - gt[0]) / gt[1]) #x pixel 像素点行列号
py = int((my - gt[3]) / gt[5]) #y pixel

[「已注销」]

引用 1 楼 cobinzhang 的回复:

对于一个有投影的地理栅格数据，如果有栅格点的定位方法，比如经纬度，必然可以找到该栅格点的像素值。

gt=ras_ds.GetGeoTransform()
mx,my 为经纬度
px = int((mx - gt[0]) / gt[1]) #x pixel 像素点行列号
py = int((my - gt[3]) / gt[5]) #y pixel

谢谢啦，我之后试一试，请问有没有这方面的具体教程可以推荐的呢？