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

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

新規登録して質問してみよう
ただいま回答率
85.50%
Python 3.x

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

Q&A

0回答

1826閲覧

basemapにおける『cyl projection cannot cross pole』のエラーの対処

R.Shigemori

総合スコア3376

Python 3.x

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

0グッド

0クリップ

投稿2017/12/13 02:42

###前提・実現したいこと
取得したデータを地図上にプロットしたくbasemapのインストールを行いました。
インストールできていることを確認するため、単純に日本地図を描画するコードを実行したところ、『cyl projection cannot cross pole』というエラーが出ました。

###発生している問題・エラーメッセージ

ValueError Traceback (most recent call last) <ipython-input-6-d7b0840e2c1b> in <module>() 6 7 # Basemapの定義 ----> 8 m = Basemap(llcrnrlat=south,urcrnrlat=north, llcrnrlon=west,urcrnrlon=east,resolution='h') 9 # 海岸線を引く 10 m.drawcoastlines() C:\Users\9047247\AppData\Local\Continuum\Anaconda3\lib\site-packages\mpl_toolkits\basemap\__init__.py in __init__(self, llcrnrlon, llcrnrlat, urcrnrlon, urcrnrlat, llcrnrx, llcrnry, urcrnrx, urcrnry, width, height, projection, resolution, area_thresh, rsphere, ellps, lat_ts, lat_1, lat_2, lat_0, lon_0, lon_1, lon_2, o_lon_p, o_lat_p, k_0, no_rot, suppress_ticks, satellite_height, boundinglat, fix_aspect, anchor, celestial, round, epsg, ax) 1107 if self.resolution is not None: 1108 self.coastsegs, self.coastpolygontypes =\ -> 1109 self._readboundarydata('gshhs',as_polygons=True) 1110 # reformat for use in matplotlib.patches.Polygon. 1111 self.coastpolygons = [] C:\Users\9047247\AppData\Local\Continuum\Anaconda3\lib\site-packages\mpl_toolkits\basemap\__init__.py in _readboundarydata(self, name, as_polygons) 1239 if containsPole and\ 1240 self.projection in _cylproj + _pseudocyl + ['geos']: -> 1241 raise ValueError('%s projection cannot cross pole'%(self.projection)) 1242 # make sure some projections have has containsPole=True 1243 # we will compute the intersections in stereographic ValueError: cyl projection cannot cross pole

うまく描画用データが読み込めていないような感じがしますが、わかりません。

###該当のソースコード

python

1import numpy as np 2import pandas as pd 3import matplotlib.pyplot as plt 4import seaborn as sns 5import _geoslib 6from mpl_toolkits.basemap import Basemap 7 8# 緯度経度で範囲を指定する 9north = 46. 10south = 30. 11east = 147. 12west = 128. 13 14# Basemapの定義 15m = Basemap(llcrnrlat=south,urcrnrlat=north, llcrnrlon=west,urcrnrlon=east,resolution='h') 16# 海岸線を引く 17m.drawcoastlines()

###補足情報(言語/FW/ツール等のバージョンなど)
os : windows(32bit)
python : 3.6.1
basemap: 1.1.0

basemapのインストールはhttp://www.lfd.uci.edu/~gohlke/pythonlibs/より該当するwhlファイルをダウンロードしてpipで行いました。
また、basemapの前提となるpyprojのインストールは、Anaconda Navigatorから行いました。
(versionは1.9.5.1です)

###確認したことなど
stackoverflowを検索しましたが、検索の仕方が悪いのか解決策を見つけることが出来ませんでした。
別のPCのwindows(32bit)+python3.5+basemap1.1.0は、うまく動いています。また、windows(64bit)+python3.6+basemap1.1.0も同様にうまく動いています。python3.5にすれば対応できることはわかっているのですがinstall資源がなく不可能な状態です。
以前に描画解像度(resolution='h')を指定したところ『データをアップデートしろ』という主旨のエラーメッセージが出たことがあるので、同じことをすればいいかもと思いつつ、当時に行ったことを記録しておらず試すことができていません。

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

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

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

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

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

guest

あなたの回答

tips

太字

斜体

打ち消し線

見出し

引用テキストの挿入

コードの挿入

リンクの挿入

リストの挿入

番号リストの挿入

表の挿入

水平線の挿入

プレビュー

まだ回答がついていません

会員登録して回答してみよう

アカウントをお持ちの方は

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

ただいまの回答率
85.50%

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

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

質問する

関連した質問