GeoGridUnitArea[proj,loc]
位置 loc の周囲の小さい地理領域の境界で評価された,投影法 proj で入手された地理格子上の単位面積に対応する実際の地理面積を与える.
GeoGridUnitArea
GeoGridUnitArea[proj,loc]
位置 loc の周囲の小さい地理領域の境界で評価された,投影法 proj で入手された地理格子上の単位面積に対応する実際の地理面積を与える.
詳細とオプション
- GeoGridUnitAreaは,与えられた場所の周りの地理投影による局所的な面積の歪みを説明する.
- GeoGridUnitAreaは,大域的名義スケール(地理モデルと地図の基準モデルの縮小係数,伝統的に1:125000およびこれに類似したもの)と地理投影から導かれる面積の局所的歪みを組み合せる.
- GeoGridUnitArea[…]の結果は,(地球あるいはその他の天体の)地理モデル上の地理領域のQuantityと投影された地理格子の無次元領域の比に対応する.
- 地理投影は,デフォルトのパラメータを持った名前付きの投影法"proj"として,あるいは{"proj",params}として与えることができる.ただし,"proj"はGeoProjectionDataの任意の実体,params は"StandardParallels"->{33,60}のようなパラメータの規則である.GeoProjectionData["proj"]は,投影法"proj"のパラメータのデフォルト値を与える.
- 位置 loc は,度を単位とした座標のペア{lat,lon},GeoPosition[…]あるいはGeoGridPosition[…]のような地理位置オブジェクト,あるいは地理実体Entity[…]として与えることができる.
- GeoGridUnitAreaはその位置と方向の引数に縫い込まれる.
- 次は,GeoGridUnitAreaの使用可能なオプションである.
-
GeoModel Automatic 地球あるいは天体のモデル UnitSystem $UnitSystem 結果に使う単位系
例題
すべて開く すべて閉じる例 (1)
エッフェル塔の位置におけるメルカトル図法による地理格子の単位面積を計算する:
scale = GeoGridUnitArea["Mercator", Entity["Building", "EiffelTower::5h9w8"]]投影座標における辺の長さが0.01の小さい正方形の面積は以下のようになる:
0.01 ^ 2scalexy = GeoGridPosition[Entity["Building", "EiffelTower::5h9w8"], "Mercator"]["GridXY"]GeoGraphics[{Opacity[0.2], Rectangle[xy - {0.005, 0.005}, xy + {0.005, 0.005}]}, GeoProjection -> "Mercator", Frame -> True, PlotRangePadding -> 0.003]スコープ (5)
"VanDerGrinten"投影法によるユーザの地理位置の地理格子の単位面積を計算する:
GeoGridUnitArea["VanDerGrinten", Here]次は"VanDerGrinten"投影法のパラメータのデフォルト値である:
GeoProjectionData["VanDerGrinten"]GeoGridUnitArea[{"VanDerGrinten", "ReferenceModel" -> 1000}, Here]GeoGridUnitArea[{"VanDerGrinten", "Centering" -> {-40, 150}}, Here]{lat,lon}の度を単位とするペアを使って位置を指定する:
GeoGridUnitArea["Orthographic", {-80.5, 30.32}]GeoGridUnitArea["Orthographic", GeoPosition[{-80.5, 30.32}]]GeoGridUnitArea["Orthographic", GeoPositionXYZ[{911681, 533170, -6268998}]]GeoGridUnitArea["Orthographic", GeoGridPosition[{0.714285, -0.753047}, "Albers"]]地理Entityオブジェクトを使って位置を指定する:
GeoGridUnitArea["Albers", Entity["City", {"Toronto", "Ontario", "Canada"}]]GeoGridUnitArea["Mollweide", {GeoPosition[{10, 40}], GeoGridPosition[{-30.3, 20}, "Mercator"], {85, 178}}]QuantityArray出力を正規形に変換する:
%//NormalGeoGridUnitAreaは,多数の場所についての値を有効に処理する:
locs = RandomGeoPosition["World", 10 ^ 6];GeoGridUnitArea["Mollweide", locs]//AbsoluteTimingオプション (2)
GeoModel (1)
デフォルトで,GeoGridUnitAreaは地球についての値を返す:
GeoGridUnitArea["Bonne", {50, -100}]月の対応する点について同じ計算を行うと,より小さいスケールが返される:
GeoGridUnitArea["Bonne", {50, -100}, GeoModel -> "Moon"]GeoGridUnitArea["Bonne", {50, -100}, GeoModel -> Quantity[2000, "Kilometers"]]アプリケーション (3)
Plot3D[GeoGridUnitArea["Mercator", {lat, lon}], {lat, -89, 89}, {lon, -180, 180}]GeoGridUnitAreaを使って,地図上に投影された狭い領域の地理領域を近似する:
geopol = EntityValue[Entity["City", {"SanFrancisco", "California", "UnitedStates"}], "Polygon"]pol = GeoGridPosition[geopol, "Albers"]Area[pol /. GeoGridPosition[expr_, proj_] :> expr]領域の中心的な位置における地理格子の単位面積を結果に掛ける:
% GeoGridUnitArea["Albers", Entity["City", {"SanFrancisco", "California", "UnitedStates"}]]geopol//GeoAreaGeoGraphics[{Opacity[0.5], Disk[{0, 0}, 1]}, GeoRange -> "World", GeoProjection -> "ConicEquidistant", Frame -> True, Axes -> True]area[x_Real, y_Real] := QuantityMagnitude[GeoGridUnitArea["ConicEquidistant", GeoGridPosition[{x, y}, "ConicEquidistant"], GeoModel -> "ITRF00"], "Kilometers" ^ 2];Quantity[NIntegrate[area[x, y], {x, y}∈Disk[{0, 0}, 1]], "Kilometers" ^ 2]GeoGridPosition[N@CirclePoints[10000], "ConicEquidistant"]GeoArea[Polygon[GeoPosition[%]], UnitSystem -> "Metric"]特性と関係 (3)
loc = Entity["Lake", "LakeVictoria::b8795"];proj = "WagnerI";ある位置に収束する,投影された徐々に小さくなる正方形族を構築する:
{x, y} = GeoGridPosition[loc, proj]["GridXY"]square[s_] := Polygon[{{x - s, y + s}, {x + s, y + s}, {x + s, y - s}, {x - s, y - s}}]squares = square /@ PowerRange[1, 1 / 2000, 1 / 2];GeoGraphics[{EdgeForm[Black], FaceForm[], squares}, GeoRange -> "World", GeoProjection -> proj, Frame -> True]これらの正方形について,投影されていない面積と投影された面積の比を計算する:
ratio[pol_] := GeoArea[GeoPosition[GeoGridPosition[#, proj]& /@ pol]] / Area[pol]ratio /@ squaresこれらの値は,その位置における投影についてのGeoGridUnitAreaで与えられた値に収束する:
GeoGridUnitArea[proj, loc, GeoModel -> "ITRF00"]投影された単位面積は地理モデルサイズの正方形によって異なる:
proj = GeoProjectionData["Mollweide"]GeoGridUnitArea["Mollweide", GeoPosition[{10, 20}], GeoModel -> #]& /@ Quantity[{10, 100, 1000}, "Kilometers"]投影された単位面積は基準モデルの逆二乗または中心スケール因子につれて変化する:
GeoGridUnitArea[{"Mollweide", "ReferenceModel" -> #}, GeoPosition[{10, 20}]]& /@ {10, 100, 1000}GeoGridUnitArea[{"Mollweide", "CentralScaleFactor" -> #}, GeoPosition[{10, 20}]]& /@ {10, 100, 1000}equalarea = GeoProjectionData["EqualArea"]世界中のランダムな位置の集合について,投影された単位面積区分の幅が0であることを確かめる:
points = RandomGeoPosition["World", 1000];MinMax /@ QuantityMagnitude[GeoGridUnitArea[{#, "ReferenceModel" -> 1}, points, GeoModel -> Quantity[1, "Meters"]]& /@ equalarea, "Meters" ^ 2]考えられる問題 (1)
地理位置が投影できなければ,その位置についての投影された単位面積も計算できない:
GeoGridPosition[GeoPosition[{20, 120}], "Orthographic"]GeoGridUnitArea["Orthographic", GeoPosition[{20, 120}]]この位置は,デフォルトの中心を使った"Orthographic"図法による地球の半分には含まれない:
GeoGraphics[GeoRange -> "World", GeoProjection -> "Orthographic", GeoGridLines -> Automatic]テクニカルノート
テキスト
Wolfram Research (2019), GeoGridUnitArea, Wolfram言語関数, https://reference.wolfram.com/language/ref/GeoGridUnitArea.html.
CMS
Wolfram Language. 2019. "GeoGridUnitArea." Wolfram Language & System Documentation Center. Wolfram Research. https://reference.wolfram.com/language/ref/GeoGridUnitArea.html.
APA
Wolfram Language. (2019). GeoGridUnitArea. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/GeoGridUnitArea.html
BibTeX
@misc{reference.wolfram_2026_geogridunitarea, author="Wolfram Research", title="{GeoGridUnitArea}", year="2019", howpublished="\url{https://reference.wolfram.com/language/ref/GeoGridUnitArea.html}", note=[Accessed: 09-July-2026]}
BibLaTeX
@online{reference.wolfram_2026_geogridunitarea, organization={Wolfram Research}, title={GeoGridUnitArea}, year={2019}, url={https://reference.wolfram.com/language/ref/GeoGridUnitArea.html}, note=[Accessed: 09-July-2026]}