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

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

ただいまの
回答率

90.33%

  • Python

    9232questions

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

basemapを使って地図に駅と利用者を描画したい

解決済

回答 1

投稿 編集

  • 評価
  • クリップ 0
  • VIEW 65

kenken0023

score 2

前提・実現したいこと

https://qiita.com/duonys/items/c941bc2818abe5cc1da7に記載されている、basemapを用いて山手線の駅の座標に利用者数を表した円を表示するというもので、実際に自分で行った際にエラーが出て実現できていない。

ちなみに私はwindows10を使っているのですが、記載されているプログラムの環境はMacです。

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

記載されているプログラムは2つに分かれており、それらを一つにまとめた際、よくわからないエラーが出てしまう。

エラーメッセージ


Traceback (most recent call last):
File "mapping3.py", line 92, in <module>
bg_img = Image.open(file_name)
File "C:\Users\mikihome2015\Anaconda3\lib\site-packages\PIL\Image.py", line 2477, in open
fp = builtins.open(filename, "rb")
FileNotFoundError: [Errno 2] No such file or directory: 'osm.png'

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File "mapping3.py", line 94, in <module>
bg_img = get_osm_img(minlat=minlat, minlon=minlon, maxlat=maxlat, maxlon=maxlon, scale=60000)
File "mapping3.py", line 80, in get_osm_img
return Image.open(BytesIO(response.content))
File "C:\Users\mikihome2015\Anaconda3\lib\site-packages\PIL\Image.py", line 2519, in open
% (filename if filename else fp))
OSError: cannot identify image file <_io.BytesIO object at 0x0000019AF7C1CF68>

該当のソースコード

import pandas as pd
from http.client import IncompleteRead
import time
#-----------------------------------------------
import requests
from io import BytesIO
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
from PIL import Image
#%matplotlib inline
#もともとこの上の%matplotlib inlineの行はコメントアウトしていなかったがsegmentErrorのためコメントアウト
#-----------------------------------------------


# Retrieve station's lat/lon
def fetch_latlon():
    url = 'http://qiita.com/butchi_y/items/3a6b70b38e13dc56ef13'
    data = pd.read_html(url)[0]
    data.columns = ['name', 'lat', 'lon']
    return data

# Retrieve the number of passenger
def fetch_passenger():
    # fetch data
    urls = ['http://www.jreast.co.jp/passenger/index.html', 'http://www.jreast.co.jp/passenger/2014_01.html']
    data = pd.DataFrame()
    for url in urls:
        df = None
        for i in range(5):
            try:
                df = pd.read_html(url)
                break
            except IncompleteRead as ir:
                time.sleep(0.1)
        # Sometimes fetching passenger fails
        if df is None:
            print('Failed to fetch passenger.')
            sys.exit(1)

        for i in range(2):
            # Delete rank column in the first table
            if url == urls[0]:
                df[i] = df[i].iloc[:, 1:]
                df[i].columns = [col for col in range(len(df[i].columns))]
            data = data.append(df[i])
    # arrange data
    data = data.iloc[2:, [0, 3]]
    data.columns = ['name', 'passenger']
    data.reset_index(inplace=True, drop=True)
    return data

# Retrieve lat/lon and passenger info and merge them into one data
def fetch_data():
    latlon = fetch_latlon()
    passenger = fetch_passenger()
    data = latlon.merge(passenger, how='inner', on='name')
    return data

data = fetch_data()

#-------------------------------------------------------------------

# min/max circle size of plot
min_size = 200
max_size = 2000

# Retrive static OpenStreetMap
def get_osm_img(minlat, minlon, maxlat, maxlon, scale=60000, img_format='png'):
    url = 'http://www.openstreetmap.org/export/finish'
    payload = {
        'mapnik_format': img_format, 
        'mapnik_scale': scale, 
        'minlon': minlon, 
        'minlat': minlat, 
        'maxlon': maxlon, 
        'maxlat': maxlat, 
        'format': 'mapnik'
    }
    response = requests.post(url, payload)
    return Image.open(BytesIO(response.content))

fig = plt.figure(figsize=(15, 15))

minlat, minlon, maxlat, maxlon = 35.61, 139.67, 35.75, 139.80
bmap = Basemap(projection='merc', llcrnrlat=minlat, urcrnrlat=maxlat, llcrnrlon=minlon, urcrnrlon=maxlon, lat_ts=0, resolution='l')

x, y = bmap(data['lon'].values, data['lat'].values)

file_name = 'osm.png'
bg_img = None
try:
    bg_img = Image.open(file_name)
except FileNotFoundError as fnfe:
    bg_img = get_osm_img(minlat=minlat, minlon=minlon, maxlat=maxlat, maxlon=maxlon, scale=60000)
    bg_img.save(file_name)
bmap.imshow(bg_img, origin='upper')
bmap.scatter(x, y, c=data['passenger'], cmap=plt.cm.get_cmap('seismic'), alpha=0.5, s=data['passenger'].map(lambda x: (x - data['passenger'].min()) / (data['passenger'].max() - data['passenger'].min()) * (max_size - min_size) + min_size))

# plt.colorbar()
plt.show()
# plt.savefig('out.png')

試したこと

まず2つあったプログラムの内の1つ(山手線の駅の座標データをスクレイピングするもの)だけでプログラムを動かしてみて、エラーが出ないことをか確認した。

その後、プログラムをつなげた際エラーが出たのでpythonに必要なライブラリとbasemapがきちんとインストールされているか確認した。

補足情報(FW/ツールのバージョンなど)

python3.7.0
プログラム上から11行目にある%matplotlib inlineの%の意味も分かりませんでした。

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

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

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

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

質問への追記・修正、ベストアンサー選択の依頼

  • t_obara

    2018/12/07 11:26

    FileNotFoundError: [Errno 2] No such file or directory: 'osm.png' というエラーについてはいかがですか?

    キャンセル

  • kenken0023

    2018/12/07 11:43

    返信ありがとうございます。このエラーについてはfile_name='osm.png'という時点ではただの何もない画像ファイルですがtry、except文で恐らくfile_nameを開くことが出来なければ、exceptに入ってbs_omg_img関数に具体的なステータスを入れて画像データとして入手するという構造になっているんだと思っています。try,exceptについては私自身使ったことがないのでどのような挙動になるか調べてみようと思います。

    キャンセル

  • kenken0023

    2018/12/07 12:24

    try,except文をなくして、直接もともとexcept文に書いてる内容を書くとNo such fileのエラーは消えました。

    キャンセル

回答 1

checkベストアンサー

+1

Open Street MapはWebAPIを介したアクセスが許容されていないそうです。そのため、get_osm_imgの結果である地図画像がないため、エラーになっています。どうしても地図画像が必要であればマニュアル操作で取得することを考えるべきでしょう。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2018/12/07 11:29

    回答ありがとうございます。
    マニュアル操作というのはいわゆる手動でということでしょうか?

    キャンセル

  • 2018/12/07 12:36

    Open Street Map のサーバーに多数のアクセスをかけるとIPアドレスが、アクセスブロックされたりするので、作成した Web サービスで使うのであれば、OpenStreetMap のサーバーを自前で構築する必要があると思います。

    https://dev.classmethod.jp/cloud/aws/install-openstreetmap-on-ec2-ubuntu14_04/ や https://blog.mori-soft.com/entry/2017/12/28/104556 等を参考に...

    キャンセル

  • 2018/12/10 20:24

    回答ありがとうございます。少し遅くなりましたが、今からやってみます。

    キャンセル

同じタグがついた質問を見る

  • Python

    9232questions

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