python gdal 矢量裁剪柵格

下面的代碼是python gdal 利用矢量文件在柵格圖上裁剪出對應區域的實現，首先說明裁剪出來的大小肯定是完全一致的，但是我仔細看了下，貌似有點偏移，這點偏移應該是無法避免的。
代碼有兩版，第一版是完全把圖讀進去了，問題就是圖像太大就沒辦法，後來想gdal本來就可以不讀出圖只獲取部分圖像，然後改了下，最後得到第二版，完美。

# -*- coding: utf-8 -*-
import os
import numpy as np
import gdal
from osgeo import gdal, gdalnumeric, ogr, osr, gdal_array
gdal.UseExceptions()

def world2Pixel(