Python DEM高程精度计算__新建空的元组，存入

import gdal
import sys
import os
import math

        
        controlFile=open("D:\\DEM2\\11.txt","r")
        lstcon=controlFile.readlines()
        pNum=0
        sum=0

        id1= []

        for point in lstcon:
            lstvalue=point.split(",")
            ID=lstvalue[0]
            x=float(lstvalue[1])
            y=float(lstvalue[2])
            z=float(lstvalue[3])

            if x>originX and x<finalX and y>finalY and y<originY:
                xOffset=(x-originX)/pixelWidth
                yOffset=(y-originY)/pixelHeight
                data=dataset.ReadAsArray(xOffset,yOffset,1,1)
                deviationZ=data-z
                dt=deviationZ*deviationZ
                sum += dt
                pNum+=1
                deZ =float(deviationZ)
                id1.append(ID+str(x).rjust(20)+str(y).rjust(20)+str(z).rjust(20)+str(deZ).rjust(25)+"\n")#右对齐
        id_str=''.join(id1).strip("\n")
               
        m=math.sqrt(sum/pNum)
        m_file.write("点数："+str(pNum)+"\n"+"中误差："+str(m)+"\n"+"ID              X                   Y                   Z                       dz   "+"\n"+id_str)


        m_file.close()