Scikitimage从tif读取卫星图像并删除具有nodata值的所有行/列

Question

我试图在一个储存为tif的卫星衍生光栅上运行一个阈值算法。在算法中，场景边界上出现的无数据值（-3.40282306e+38）将被考虑在内，并导致不希望的行为。我在寻找一种方法，在运行阈值算法之前删除所有没有数据值的行/列。在

我试图将所有无数据值的单元格设置为np.nan公司但是我没有得到我想要的结果。在

我的代码：

import numpy as np
from skimage import filters
from skimage import exposure
from skimage.io import imread, show
from skimage.filters import try_all_threshold

no_data_value = -3.40282306e+38

ndwi = imread(<'my.tif'>)
ndwi[ndwi == no_data_value] = np.nan
ndwi = ndwi[~np.isnan(ndwi).any(axis=1)]
val = filters.threshold_otsu(ndwi)

Answer 1

filters.threshold_otsu不关心传入的数组的形状，因此您也可以执行以下操作：

threshold_otsu(z[z > -1])