def map_pixel_inv(row, col, cellx, celly, xmin, ymax):
col = np.asarray(col)<EOL>row = np.asarray(row)<EOL>point_x = xmin+col*cellx<EOL>point_y = ymax+row*celly<EOL>return point_x, point_y<EOL>
Usage:
    map_pixel(xcoord, ycoord, x_cell_size, y_cell_size, xmin, ymax)


where:
        xmin is leftmost X coordinate in system

        ymax is topmost Y coordinate in system

Example:
        raster = HMISea.tif

        ndv, xsize, ysize, geot, projection, datatype = get_geo_info(raster)

        row, col = map_pixel(x,y,geot[1],geot[-1], geot[0],geot[3])