базовая карта matplotlib и сечение гипоцентров - PullRequest
Новая
       5

базовая карта matplotlib и сечение гипоцентров

0 голосов
/ 15 мая 2018

Мне удалось нанести эпицентры (из файла широты, долготы, глубины, величины) на двухмерную карту XY, но я не могу найти способ их сечения. (Имеется в виду выбор азимута и диапазона ширины и нанесение гипоцентров на график xz).

GMT имеет возможность делать такие сечения, используя проекцию, как описано здесь: http://mygeologypage.ucdavis.edu/gps/GMT/LONG_VALLEY/hypocenter.html

Есть идеи, как это можно сделать?

Пример, аналогичный тому, что я делал до сих пор: 

import requests
from csv import DictReader 
DATA_URL='http://earthquake.usgs.gov/earthquakes/feed/v1.0/summary/4.5_month.csv'
print("Downloading", DATA_URL)
resp = requests.get(DATA_URL)
quakes = list(DictReader(resp.text.splitlines()))

lngs = [float(q['longitude']) for q in quakes]
lats = [float(q['latitude']) for q in quakes]
depths=[float(q['depth']) for q in quakes]
mags = [2 ** float(q['mag']) for q in quakes]

import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
plt.figure(figsize=(14, 8))
earth = Basemap()
earth.bluemarble(alpha=0.42)
earth.drawcoastlines(color='#555566', linewidth=1)  
plt.scatter(lngs, lats, mags, 
            c='red',alpha=0.5, zorder=10)
...