地理的な場所の格子
新関数のCoordinateBoundsArrayとCoordinateBoundingBoxArrayを使って,指定された刻み幅の多次元数値配列を作ることができる.
アメリカ合衆国の経緯度の境界を計算する.
In[1]:=
![Click for copyable input](assets.ja/grid-of-geographical-locations/In_13.png)
{lat, lon} = GeoBounds[Entity["Country", "UnitedStates"]]
Out[1]=
![](assets.ja/grid-of-geographical-locations/O_12.png)
そこから,刻み幅が2度の座標の配列を生成する.
In[2]:=
![Click for copyable input](assets.ja/grid-of-geographical-locations/In_14.png)
arr = CoordinateBoundsArray[QuantityMagnitude@{lon, lat}, 2];
配列から地理格子を作成する.
In[3]:=
![Click for copyable input](assets.ja/grid-of-geographical-locations/In_15.png)
grid = GeoGridPosition[arr, "Equirectangular"];
In[4]:=
![Click for copyable input](assets.ja/grid-of-geographical-locations/In_16.png)
pos = GeoPosition[grid]
Out[4]=
![](assets.ja/grid-of-geographical-locations/O_13.png)
地図上にその座標をプロットする.
In[5]:=
![Click for copyable input](assets.ja/grid-of-geographical-locations/In_17.png)
graphics = GeoGraphics[Point[pos], GeoGridLines -> Automatic]
Out[5]=
![](assets.ja/grid-of-geographical-locations/O_14.png)
他の投影法で点をプロットする.
In[6]:=
![Click for copyable input](assets.ja/grid-of-geographical-locations/In_18.png)
Show[graphics, GeoProjection -> "Mollweide"]
Out[6]=
![](assets.ja/grid-of-geographical-locations/O_15.png)
In[7]:=
![Click for copyable input](assets.ja/grid-of-geographical-locations/In_19.png)
Show[graphics, GeoProjection -> "Mercator"]
Out[7]=
![](assets.ja/grid-of-geographical-locations/O_16.png)