5公里范围内 经纬度判断_「技能技巧」Python根据经纬度批量提取遥感图像光谱信息...

以下代码，可以根据经纬度批量提取遥感影像光谱反射率信息，其实也就是提取特定经纬度的像元值。我找了一幅遥感影像，实验了一下，效果还是可以的，分享一下，如需自取。

实验材料：

1.一幅(也可多幅，不需拼接)经过辐射定标和QUAC大气校正的遥感影像(Landsat8，ENVI处理)，包含DAT文件和HDR头文件；

2.一个包含经纬度信息的csv文件(虚拟的点)；

如下图所示：

接下来，运行文章末尾的代码，运行过程和结果如下图：

将record文件内的提取内容经过整理后，与ENVI提取的光谱曲线对比，结果一致，如下所示，其中只提取了4个点，因为第5个点是在图像外，被过滤掉了：

---

注意事项：

1.代码运行过程中对存在于图像边缘的点可能会存在范围上的误判，不过一般不影响结果。对于出现范围误判的点，在结果文件中会有“AN ERROR”字样出现，对个别误判的，手动查一下即可，一般没什么意义。本例中没有误判点，没法演示；