def METHOD_NAME(self):
    import lenstronomy.Util.util as util
    numPix = 10
    theta = 0.5
    deltaPix = 0.05
    subgrid_res = 3
    Mpix2a = (
        np.array([[np.cos(theta), -np.sin(theta)], [np.sin(theta), np.cos(theta)]])
        * deltaPix
    )
    x_grid, y_grid = util.make_grid_transformed(numPix, Mpix2a)
    coords = Coordinates(Mpix2a, ra_at_xy_0=x_grid[0], dec_at_xy_0=y_grid[0])
    x_grid_high_res, y_grid_high_res = util.make_subgrid(
        x_grid, y_grid, subgrid_res=subgrid_res
    )
    coords_sub = Coordinates(
        Mpix2a / subgrid_res,
        ra_at_xy_0=x_grid_high_res[0],
        dec_at_xy_0=y_grid_high_res[0],
    )
    x, y = coords_sub.map_coord2pix(x_grid[1], y_grid[1])
    npt.assert_almost_equal(x, 4, decimal=10)
    npt.assert_almost_equal(y, 1, decimal=10)
    x, y = coords_sub.map_coord2pix(x_grid[0], y_grid[0])
    npt.assert_almost_equal(x, 1, decimal=10)
    npt.assert_almost_equal(y, 1, decimal=10)
    ra, dec = coords_sub.map_pix2coord(1, 1)
    npt.assert_almost_equal(ra, x_grid[0], decimal=10)
    npt.assert_almost_equal(dec, y_grid[0], decimal=10)
    ra, dec = coords_sub.map_pix2coord(1 + 2 * subgrid_res, 1)
    npt.assert_almost_equal(ra, x_grid[2], decimal=10)
    npt.assert_almost_equal(dec, y_grid[2], decimal=10)
    x_2d = util.array2image(x_grid)
    y_2d = util.array2image(y_grid)
    ra, dec = coords_sub.map_pix2coord(1 + 2 * subgrid_res, 1 + 3 * subgrid_res)
    npt.assert_almost_equal(ra, x_2d[3, 2], decimal=10)
    npt.assert_almost_equal(dec, y_2d[3, 2], decimal=10)
    ra, dec = coords.map_pix2coord(2, 3)
    npt.assert_almost_equal(ra, x_2d[3, 2], decimal=10)
    npt.assert_almost_equal(dec, y_2d[3, 2], decimal=10)