基于OSGEO平台的遥感图像动态区域切割技术

 def rasterCrop(inputRaster='', cropOutputAbspath='', cropSize=300, xPiexlStart=50, yPiexlStart=50):

    
     """
    
     inputRaster需要裁剪的遥感影像栅格数据的绝对路径
    
     cropSize裁剪的尺寸，默认是裁剪正方形的
    
     xPiexlStart, yPiexlStart为裁剪起点的行列号坐标，起点也相当于从（0,0）点的偏移距离
    
     """
    
     ds = gdal.Open(inputRaster)
    
     geoTranform = ds.GetGeoTransform()
    
     projectionInform = ds.GetProjection()
    
     # 获取栅格数据行列数
    
     numX, numY = ds.RasterXSize, ds.RasterYSize
    
     if xPiexlStart + cropSize > numX:
    
     xPiexlStart = int(numX - cropSize)
    
     if yPiexlStart + crop