Gdal随笔

开始研究遥感影像的分割，尝试去通过gdal对影像进行读取并对图像进行分割。过程中遇到了许多问题，现总结如下。

1.error1:range must be finite

　　报这个错误是因为数组中有Nan值。一开始很纳闷，因为readasarray的图像在envi里打开并没有问题，黑边也作了处理。查询satckoverflow后发现是因为影像太大，无法读取，崩溃导致nan。通过采用分块读取的方法解决了问题。

但是仍有不足，这里的分块读取分割成四部分，不够intelligence。后续将作改进。

2.idea:采用polygonize方法可以用mask来提取栅格数据并且转换成shp文件。gdal:polygonize。

3.GetRasterBand(n)方法读取波段，readasarray转换成数组阵列。

2019/11/30

今天发现一些问题，以前对range must be finite的报错给出的答案是错误的，所以来进行修改。其实在读取遥感影像时，一般小图幅是不需要分块处理的，只有大到了一定程度才会需要去分块读取。今天用data = gdal_array.LoadFile(raster)读取了8000*8000的数据，并且是配置较差的电脑完成的，所以猜测ReadAsArray应该也是可以直接读取的（的确可以），这种图幅完全不需要分块读取。那么经过排查发现，报错问题在于ndwi值范围异常，-inf~1341的范围毫无疑问是错的，所以这里会导致pillow图像处理报错。最终发现错误原因在于辐射定标做错了，这里贴个图记录一下。下次不要再错了。

总结：如果遇到range must be finite，那么十之八九是你的辐射定标、大气校正存在问题，导致计算的ndwi值有巨大偏差。