python处理DEM坡度坡向

Question

您好，可否咨询一下问题，您的这个程序运行的结果是正确的，但我跟教程上对照公式，发现dy,dx好像反了，可否指导一下，谢谢！https://blog.csdn.net/weixin_40501429/article/details/115302922

s_dx[i,j] = ((npgrid[i-1,j+1]+2npgrid[i,j+1]+npgrid[i+1,j+1])-(npgrid[i-1,j-1]+2npgrid[i,j-1]+npgrid[i+1,j-1])) / (8 * sizex)
s_dy[i, j] = ((npgrid[i+1, j-1] + 2 * npgrid[i+1, j] + npgrid[i+1,j+1])-(npgrid[i-1,j-1]+2 * npgrid[i-1,j] + npgrid[i-1,j+1])) / (8 * sizey)
s_dx的计算脚标为 j+1，j-1之差，即y方向一致的数据，对应就是ghi与abc了，相当于dy
s_dy的计算脚标为 i+1, i-1之差，即x方向一致的数据，对应就是cfi与adg了，相当于dx吧
难道DEM中x，y方向与编号的IJ方向对应是反的？

Answer 1

大地坐标系中，xy方向一般x指向北，y指向东。与一般坐标系中的xy轴是相反的。