r - RでマップとXY ggplotチャートを組み合わせる

Question

遺跡を表すマップと、さまざまな考古学的オブジェクトの XY ggplot チャートを組み合わせる必要があります。マップは tiff ファイルにあり、その比率を尊重する必要があります。

まず、これはマップで、基準縮尺が赤で強調表示されています (たとえば、X 軸では、-6000 から -4000 までの距離は 20 メートル、Y 軸では、900 から 2100 までは 12 メートルです)。

私のggplotチャートは、次のコードを実行することで取得されます:

archaeo <- ggplot() + 
    geom_ellipsis(data=Unit_H, 
        aes(x0 = X, y0 = Y, a = Diameter_E.W/2+250, b = Diameter_N.S/2+250, angle = 0), 
        lwd=0, col="darkgray", fill="gray", alpha=0.15) +     
    geom_ellipsis(data=Unit_H, 
        aes(x0 = X, y0 = Y, a = Diameter_E.W/2+120, b = Diameter_N.S/2+120, angle = 0), 
        lwd=0, col="darkgray", fill="gray", alpha=0.25) + 
    geom_ellipsis(data=Unit_H, 
        aes(x0 = X, y0 = Y, a = Diameter_E.W/2, b = Diameter_N.S/2, angle = 0), 
        lwd=0.5, col="darkgray", fill="gray", alpha=0.75) + 
    geom_point(data=Unit_H, aes(X, Y), size = 0.5) + 
    geom_point(data=Refits_H_trans, aes(x,y,group=sample, colour=factor(sample))) + 
    geom_line(data=Refits_H_trans, lwd=0.2, lty=1, aes(x,y, group=sample, colour=factor(sample))) + 
    coord_fixed() + 
    theme_bw() + 
    theme(legend.position="none") + 
    ggtitle("Unit H") + 
    xlim(-6600,-3800) + 
    ylim(400,2400)

結果のチャートは次のとおりです。

さて、私の問題は、ggplot の背景としてマップを含めることです。ggpubr の background_image() を使用すると、次の結果が得られました。

map_levelH <- readPNG("Planta H-I.png")

Map.archaeo <- ggplot() + 
    background_image(map_levelH) + 
    geom_ellipsis(data=Unit_H, 
        aes(x0 = X, y0 = Y, a = Diameter_E.W/2+250, b = Diameter_N.S/2+250, angle = 0), 
        lwd=0, col="darkgray", fill="gray", alpha=0.15) +     
    geom_ellipsis(data=Unit_H, 
        aes(x0 = X, y0 = Y, a = Diameter_E.W/2+120, b = Diameter_N.S/2+120, angle = 0), 
        lwd=0, col="darkgray", fill="gray", alpha=0.25) + 
    geom_ellipsis(data=Unit_H, 
        aes(x0 = X, y0 = Y, a = Diameter_E.W/2, b = Diameter_N.S/2, angle = 0), 
        lwd=0.5, col="darkgray", fill="gray", alpha=0.75) + 
    geom_point(data=Unit_H, aes(X, Y), size = 0.5) + 
    geom_point(data=Refits_H_trans, aes(x,y,group=sample, colour=factor(sample))) + 
    geom_line(data=Refits_H_trans, lwd=0.2, lty=1, aes(x,y, group=sample, colour=factor(sample))) + 
    coord_fixed() + 
    theme_bw() + 
    theme(legend.position="none") + 
    ggtitle("Unit H") + 
    xlim(-6600,-3800) + 
    ylim(400,2400)

ご覧のとおり、ggplot とマップの縮尺が一致していません。だから、私の質問は次のとおりです。

• ggplot の X 軸と Y 軸の値を使用してマップをジオリファレンスするにはどうすればよいですか?
• 画像が歪まないように、画像の比率を維持する必要があります。どうすればいいですか？xlim 値を変更すると、画像も変更され、その比率が変更されるため、これを求めています。