卫片坐标系统检查

"""""""""""检查坐标不正确卫片"""""""""""""""""""""
import gdal
import numpy as np
import os
import shutil
path="I:/building/zhejiang_anji/anji_tif/"
os.chdir(path)
ps=[i for i in os.listdir(path) if i.endswith(".tif")]
a=[]
for p in ps:
    dataset = gdal.Open(p)  # 打开文件
    im_proj = dataset.GetProjection()  # 地图投影信息
    if im_proj.find("WGS_1984_Web_Mercator_Auxiliary_Sphere")!=-1:
        print(p)
        #os.remove(p)
        #shutil.copy(p,os.path.join("F:\\4-imp_shidian\\yangshan\\",p))
        a.append(p)