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])

def map_pixel_inv(row, col, cellx, celly, xmin, ymax):
    '''
    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])
    '''
    col = np.asarray(col)
    row = np.asarray(row)
    point_x = xmin+col*cellx
    point_y = ymax+row*celly
    return point_x, point_y