import os
import cv2
import sys
import skimage.io as io
import shapefile
import numpy as np
from osgeo import gdal
from PIL import Image
Image.MAX_IMAGE_PIXELS=100000000000000
import xlsxwriter
def coord_trans(adfGeoTransform,coord):
px = adfGeoTransform[0] + coord[0] * adfGeoTransform[1] + coord[1] * adfGeoTransform[2]
py = adfGeoTransform