已经实现将所有大于阈值的像素赋值为1，其他为0并输出tif的功能，但是观察结果发现貌似单纯凭借这个阈值得到的并非是城市范围，某些年份很大，某些年份很小，城市面积与年份对应关系不大，甚至相反加两行代码，试试提取在最佳阈值情况下聚类之后的最大斑块