質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

新規登録して質問してみよう
ただいま回答率
85.35%
cartopy

cartopyは、イギリス気象局(MetOffice)が開発している、地図の描画や地理空間データ解析を行うためのPythonライブラリです。点や線、画像などを地図座標系に変換できるなどの特徴を持ちます。

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

Q&A

解決済

1回答

1580閲覧

python:cartopyを用いて指定した領域の地図を出力したい

napo5221

総合スコア7

cartopy

cartopyは、イギリス気象局(MetOffice)が開発している、地図の描画や地理空間データ解析を行うためのPythonライブラリです。点や線、画像などを地図座標系に変換できるなどの特徴を持ちます。

Python 3.x

Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Python

Pythonは、コードの読みやすさが特徴的なプログラミング言語の1つです。 強い型付け、動的型付けに対応しており、後方互換性がないバージョン2系とバージョン3系が使用されています。 商用製品の開発にも無料で使用でき、OSだけでなく仮想環境にも対応。Unicodeによる文字列操作をサポートしているため、日本語処理も標準で可能です。

0グッド

0クリップ

投稿2023/02/02 09:50

編集2023/02/02 09:55

実現したいこと

cartopyを用いて指定した領域の地図を出力したいです。
また、適当な緯度経度ごとに罫線を入れて、その緯度経度のラベルも表示させたいです。

前提

pythonのcartopyを用いていろいろ試していますが、なかなか想定通りのものが表示されません。
■■な機能を実装中に以下のエラーメッセージが発生しました。

該当のソースコード

Python

1import matplotlib.pyplot as plt 2import cartopy.crs as ccrs 3import matplotlib.ticker as mticker 4 5fig, ax = plt.subplots(figsize=(10, 8), subplot_kw={'projection': ccrs.PlateCarree()}) 6 7ax.set_extent([85, 205, -60, 60]) 8 9gl = ax.gridlines(crs=ccrs.PlateCarree(), draw_labels=True) 10gl.xlocator = mticker.FixedLocator(range(85, 205, 30)) 11gl.ylocator = mticker.FixedLocator(range(-60, 60, 30)) 12 13ax.coastlines() 14plt.show()

試したこと

今回は東経85度-西経155度、北緯60度-南緯60度の120度×120度の範囲で30度ずつの罫線を表示させようとしてますが、子午線をまたいでいるからか、下の図のようになってしまっています。なかなか解決方法も見つからず、困っております。
ご存知の方いらっしゃいましたら、教えていただきたいです。
イメージ説明
Jupyterlabを用いています。

気になる質問をクリップする

クリップした質問は、後からいつでもMYページで確認できます。

またクリップした質問に回答があった際、通知やメールを受け取ることができます。

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

guest

回答1

0

ベストアンサー

python

1import matplotlib.pyplot as plt 2import cartopy.crs as ccrs 3import matplotlib.ticker as mticker 4 5fig, ax = plt.subplots(figsize=(10, 8), subplot_kw={'projection': ccrs.PlateCarree(central_longitude=145)}) 6 7ax.set_extent([85, 206, -61, 61], crs=ccrs.PlateCarree()) 8 9gl = ax.gridlines(crs=ccrs.PlateCarree(), draw_labels=True) 10gl.xlocator = mticker.FixedLocator([*range(85, 180, 30), -155]) 11gl.ylocator = mticker.FixedLocator(range(-60, 61, 30)) 12 13ax.coastlines() 14plt.tight_layout() 15plt.show()

イメージ説明

投稿2023/02/02 17:33

編集2023/02/02 18:24
melian

総合スコア20655

バッドをするには、ログインかつ

こちらの条件を満たす必要があります。

napo5221

2023/02/02 18:10

ありがとうございます。私が実行したところ、一番右の155度のラベルだけ表示されません。ax.set_extent([85, 206, -60, 60])のように少し範囲をずらすと表示されます。 また追加の質問で申し訳なのですが、緯度の範囲を-60,60に限定して上下の無駄な範囲をカットすることはできるのでしょうか。
melian

2023/02/02 18:28

明示的に crs=ccrs.PlateCarree() を追加することで -60, 60 に限定することができました。 > ax.set_extent([85, 206, -60, 60])のように少し範囲をずらすと表示されます。 こちらの環境では 205 でも155度のラベルが表示されるのですが、緯度の範囲指定では -60, 60 のラベルが表示されません。
napo5221

2023/02/05 09:49

環境によってかなり変わるんでしょうか、、いろいろ試してみたいと思います。ありがとうございました。
guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

15分調べてもわからないことは
teratailで質問しよう!

ただいまの回答率
85.35%

質問をまとめることで
思考を整理して素早く解決

テンプレート機能で
簡単に質問をまとめる

質問する

関連した質問