RSer工作時不可避免會用到大型的遙感影像,由于分辨率過高、區域過大、波段信息過多等原因,都會導致數據非常的大。這個時候我們在進行一些簡單的操作,如計算NDVI、二值化、分類等時,計算機的內存都會溢出。因此今天跟大家分享一下我平時分塊的方法,中間如何計算就按照自己的需求來即可。
原創作者:RS迷途小書童
博客地址:https://blog.csdn.net/m0_56729804?type=blog
1 忽略警告信息
由于GDAL版本的更新,現在運行代碼都會有一條警告Warning,但不影響程序,我有強迫癥所以必須沒有警告。另外就是指定Proj庫的路徑,以免坐標系進行警告。
gdal.DontUseExceptions()
os.environ['PROJ_LIB'] = 'G:/Anaconda/envs/pyDL/Lib/site-packages/osgeo/data/proj'
os.environ['GDAL_DATA'] = 'G:/Anaconda/envs/pyDL/Lib/site-packages/osgeo/data