このページは機械翻訳を使用して翻訳されました。元の英語を参照するには、ここをクリックします。

projfwd

緯度経度座標をx-yの地図座標に投影します

構文

[x,y] = projfwd(proj,lat,lon)

説明

[x,y] = projfwd(proj,lat,lon) は、lat と lon で指定された緯度経度座標を、proj で指定された投影座標参照系の x と y マップ座標に変換。proj は、projcrs オブジェクト (R2020b 以降)、地図投影構造、またはGeoTIFF情報構造を使用して指定します。

例

緯度経度座標を x-y に投影します

地図投影を指定して、緯度経度座標を x-y 座標に投影します。次に、投影された座標を地図上に表示します。

これを行うには、まずボストンのランドマークの緯度と経度の座標を指定します。NAD83 地理 CRS で座標を指定します。

lat = [42.3604 42.3691 42.3469 42.3480 42.3612];
lon = [-71.0580 -71.0710 -71.0623 -71.0968 -71.0941];

次に、ボストンのGeoTIFFイメージを配列とマップ セル参照オブジェクトとしてインポートします。参照オブジェクトの ProjectedCRS プロパティを照会して、地図投影に関する情報を取得します。投影された CRS の基礎となる地理 CRS が NAD83 であることを確認します。

[A,R] = readgeoraster("boston.tif");
proj = R.ProjectedCRS;
proj.GeographicCRS.Name

ans = 
"NAD83"

GeoTIFF画像と同じ投影された CRS を使用して、緯度経度座標を x-y 座標に投影します。

[x,y] = projfwd(proj,lat,lon);

GeoTIFF画像と投影座標を同じ地図上に表示します。座標のマーカー シンボルと色を変更して、より目立つようにします。次に、軸ラベルを追加します。

mapshow(A,R)
mapshow(x,y,DisplayType="point",Marker="o", ...
    MarkerFaceColor="y",MarkerEdgeColor="none")
xlabel("x (survey feet)")
ylabel("y (survey feet)")

入力引数

proj — 地図投影

projcrs オブジェクト | スカラー地図投影構造 | スカラーGeoTIFF情報構造

地図投影。projcrs オブジェクト (R2020b 以降)、スカラー地図投影構造(mstruct)、またはGeoTIFF情報構造として指定されます。地図投影構造の詳細については、defaultm を参照してください。GeoTIFF情報構造の詳細については、geotiffinfo を参照してください。

データ型: struct

lat — 測地緯度

スカラー値 | ベクトル | 行列 | N 次元配列

測地緯度。度単位でスカラー値、ベクトル、行列、または ND 配列として指定されます。lat と lon のサイズは一致する必要があります。

データ型: single | double

lon — 測地経度

スカラー値 | ベクトル | 行列 | N 次元配列

測地経度。度単位で、スカラー値、ベクトル、行列、または ND 配列として指定されます。lat と lon のサイズは一致する必要があります。

データ型: single | double

出力引数

x — 投影されたx座標

スカラー値 | ベクトル | 行列 | N 次元配列

投影された x 座標。スカラー値、ベクトル、行列、または ND 配列として返されます。

y — 投影されたy座標

スカラー値 | ベクトル | 行列 | N 次元配列

投影された y 座標。スカラー値、ベクトル、行列、または ND 配列として返されます。

ヒント

lat と lon の地理 CRS が proj の地理 CRS と一致しない場合、x と y の値は不正確になる可能性があります。proj が projcrs オブジェクトの場合、GeographicCRS プロパティを照会することで、その地理 CRS を見つけることができます。たとえば、このコードは、EPSG コード 32610 から projcrs オブジェクトを作成し、関連する地理 CRS を見つける方法を示しています。

proj = projcrs(32610);
proj.GeographicCRS.Name

ans = 

    "WGS 84"

バージョン履歴

R2006a より前に導入

R2022b: projfwdは、投影を繰り返し使用することでWindowsのパフォーマンスが向上しました。

projfwd 関数と projinv 関数は、同じ projcrs オブジェクトまたは地図投影構造を入力として同じ MATLAB^® セッションで複数回使用すると、Windows^® でパフォーマンスが向上します。

たとえば、このコードは以前のリリースよりも約 1.67 倍高速です。コード:

• DTEDファイルのグリッド座標を取得します

• projcrs オブジェクトを使用して座標を投影および逆投影します

• 別のDTEDファイルのグリッド座標を取得します

• 同じprojcrsオブジェクトを使用して座標を投影および逆投影するのに必要な時間を測定します

function timingTest
    info1 = georasterinfo("n39_w106_3arc_v2.dt1");
    R1 = info1.RasterReference;
    [latGrid1,lonGrid1] = geographicGrid(R1);

    p = projcrs(4088);
    [xfwd1,yfwd1] = projfwd(p,latGrid1,lonGrid1);
    [latinv1,loninv1] = projinv(p,xfwd1,yfwd1);

    info2 = georasterinfo("n40_w106_3arc_v2.dt1");
    R2 = info2.RasterReference;
    [latGrid2,lonGrid2] = geographicGrid(R2);

    tic
    [xfwd2,yfwd2] = projfwd(p,latGrid2,lonGrid2);
    [latinv2,loninv2] = projinv(p,xfwd2,yfwd2);
    toc
end

おおよその実行時間は次のとおりです。

• R2022a: 0.70秒

• R2022b: 0.42 秒

このコードは、関数 timingTest を呼び出すことによって、Windows 10 Intel^® Xeon^® CPU W-2133 @ 3.6 GHzテストシステムで計測されました。

R2020b: projcrs オブジェクトと追加の地図投影構造を使用して座標を投影する

projcrs オブジェクトを指定して座標を投影します。projcrs オブジェクトを使用すると、ホーチン斜メルカトル図法など、以前のリリースでは利用できなかった投影法を使用して座標を変換できます。

さらに、有効な mapprojection フィールドを使用して地図投影構造を指定して座標を投影することもできます。mapprojection フィールド オプションのリストについては、maplist または maps 関数を使用してください。

参考

関数

• geotiffinfo | projinv | projlist

オブジェクト

• projcrs

トピック

• Project and Display Raster Data
• Transform Coordinates to a Different Projected CRS