問題タブ [cartopy]

カートピーで作成された地図に等高線をプロットしたいと思います。次の場所で例を実行すると:

http://scitools.org.uk/cartopy/docs/latest/matplotlib/introductory_examples/03.contours.html

github からの最新コードのチェックアウトを使用すると、次のようになります。

これをどのように行うべきですか？

(元の発表http://matplotlib.1069221.n5.nabble.com/ANN-Cartopy-A-new-mapping-library-using-matplotlib-td39886.htmlへの応答として質問)

cartopy の最新バージョンを使用して極平射図上に等高線をプロットすることは可能ですか? 私は自分でそれを解決するのに苦労しているので、これがどのように行われるかの例を見たいです!

カートピーがさまざまな地図投影法で自動的にプロットできるというアイデアが本当に気に入っています。しかし、アイリスキューブの扱い方がわかりませんでした。その姉妹プロジェクトとして、私はできるかもしれないと思っていました。このようなことは可能ですか？

ありがとう

私は次の電子メールを受け取り、この質問への回答がすべての人に利用可能であることを確認したいと思いました。

やあ、

日付線を横切り、左側に東アジア、右側に北アメリカの西を示すカートピーを使用して、単純な緯度経度マップを設定したいと思います。次のグーグルマップは大まかに私が求めているものです：

https://maps.google.co.uk/?ll=56.559482,-175.253906&spn=47.333523,133.066406&t=m&z=4

これはCartopyで実行できますか？

太平洋地域を強調できるように、プロットの下部に180を使用して、北半球の極心平射投影プロットをプロットしたいと思います。gitの最新のカートピーを使用しており、極心平射投影法のプロットを問題なく作成できますが、プロットの下部にある経度を変更する方法がわかりません。経度の範囲を[-180、180]に設定しようとしましたが、これは役に立ちません。また、NorthPolarStereo（）はcentral_longitudeなどのキーワード引数を受け入れません。これは現在可能ですか？

私は cartopy を使用して平射図法で北極の地図を作成し、その上に (断面の位置を示すために) 線をプロットしています。次のコードを使用すると、線は極を通る直線ではなく、緯度の線に沿って進みます。

これを回避するには、x と y を次のように変更する必要があります。

北極に 2 つのデータ ポイントがあるようにします。これに対するより良い解決策はありますか？

編集：画像添付

ありがとう、

ティム

高さと幅を指定して (地理) 軸のセットを定義した場合、プロットがこれらの軸を埋めるようにするにはどうすればよいですか?

これにより、マップの適切な縦横比を持つプロットが生成されますが、代わりに軸を埋めるためにプロットが幅よりも高くなっています。これは単なる例ですが、これを行うことが重要な現実のシナリオがあります。

編集

念のために言うと、実際には結果を歪ませたいので、この例では、幅よりも高さが大きいグローバルな範囲を持つマップが本当に必要です。の使用はとの射影でax.set_aspect('auto')機能するように見えますが、おそらくすべてで機能するとは限りません (たとえば)。PlateCarreeNorthPolarStereographicOSGB

Cartopy を使用して表示したい衛星画像データがあります。ここで説明されている画像の例にうまく従いました。このコードの結果：

このコードは、次の画像を生成します

私の問題は、衛星画像データがプレートカリー図法ではなく、メルカトル図法であることです。

しかし、軸オブジェクトを取得すると

海岸線を失います。

ここで報告された問題を見ました。しかし

この画像の結果：

データは存在せず、サンディエゴは間違った場所にあります。また、緯度/経度の範囲が変更されました。私は何を間違っていますか？

ポスト ディスカッションの更新

主な問題は、メソッドを使用してターゲット プロジェクションのイメージ範囲を適切に指定していなかったことtransform_pointsです。imshowPhil が提案するように、メソッドの座標参照系についても具体的に説明する必要がありました。正しいコードは次のとおりです。

この正しく地理投影された衛星画像が得られます。

海氷データのプロットを作成しようとしています。データは EASE-North グリッドで配信されます。サンプル ファイル (HDF4) は次の場所からダウンロードできます。

EASE-Grid 用のカスタム プロジェクション クラスを作成しましたが、機能しているようです (海岸線はデータとうまく一致しています)。

Natural Earth 機能を追加しようとすると、空の Matplotlib フィギュアが返されます。

上記のスクリプトは正常に動作し、次の結果が生成されます。

しかし、コメントを外すと、ax.add_feature(land)エラーなしで失敗し、空の図のみが返されます。明らかな何かが欠けていますか？

IPython ノートブックは次のとおりです: http://nbviewer.ipython.org/6779935

私の Cartopy ビルドは、Christoph Gohlke の Web サイト (ありがとう!) のバージョン 0.9 です。

編集：

Figure を保存しようとすると、例外がスローされます。

「土地」を調べてcartopy.featureも問題はありません。すべてのポリゴンが通過し.isvalid()、すべてのリング (ext en int) が 4 つ以上のタプルで構成されています。したがって、入力形状は問題ではないようです (そして、PlateCaree() で正常に動作します)。

EASE_North に変換した後、一部のリング (南半球など) が「破損」する可能性がありますか?

編集2:

組み込みの NE 機能を削除して同じシェープファイルをロードすると (ただし、40N 未満のものはすべてクリップされます)、機能します。したがって、ある種の再投影の問題のようです。