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(point_x, point_y, 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])
    '''
    point_x = np.asarray(point_x)
    point_y = np.asarray(point_y)
    col = np.floor((point_x - xmin) / cellx).astype(int)
    row = np.floor((point_y - ymax) / celly).astype(int)
    return row, col