forked from xiaoxuzhao/FSRS_old
/
fsrs_map_all_1year_1000days.py
66 lines (65 loc) · 2.54 KB
/
fsrs_map_all_1year_1000days.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
# -*- coding: utf-8 -*-
"""
Created on Thu Jan 18 10:12:29 2018
@author: xiaoxu
map with different sites
"""
from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
import numpy as np
from pandas import read_csv
#####################
#HARDCODES
input_dir='/home/zdong/xiaoxu/FSRS/all_files/'
save_dir='/home/zdong/xiaoxu/FSRS/figure/'
#####################
df_fsrs=read_csv(input_dir+'fsrs_sites(>1 km,year).csv') #fsrs_sites(>1 km,year).py
df_getfsrs=read_csv(input_dir+'getfsrs.csv') #fsrs_sites(>1 km,year).py
df_fsrs_1000=read_csv(input_dir+'fsrs_sites(>1000days).csv') #fsrs_sites(>some days).py
fig = plt.figure()
a=fig.add_subplot(1,3,1)
my_map = Basemap(projection='merc', lat_0 = 45, lon_0 = -63,
resolution = 'h', area_thresh = 0.3,
llcrnrlon=-68, llcrnrlat=42.0,
urcrnrlon=-59, urcrnrlat=48)
my_map.drawcoastlines()
my_map.drawcountries()
my_map.fillcontinents(color = 'coral')
my_map.drawmapboundary()
x,y=my_map(df_getfsrs['Longitude'].values,df_getfsrs['Latitude'].values)
my_map.plot(x, y, 'bo', markersize=4)
a.set_title('All sites with raw data',fontsize=10)
my_map.drawparallels(np.arange(40,80,3),labels=[1,0,0,0])
my_map.drawmeridians(np.arange(-180,180,3),labels=[1,1,0,1])
a=fig.add_subplot(1,3,2)
my_map = Basemap(projection='merc', lat_0 = 45, lon_0 = -63,
resolution = 'h', area_thresh = 0.3,
llcrnrlon=-68, llcrnrlat=42.0,
urcrnrlon=-59, urcrnrlat=48)
my_map.drawcoastlines()
my_map.drawcountries()
my_map.fillcontinents(color = 'coral')
my_map.drawmapboundary()
x,y=my_map(df_fsrs['Longitude'].values,df_fsrs['Latitude'].values)
my_map.plot(x, y, 'bo', markersize=4)
a.set_title('Longterm sites (>1 year)',fontsize=10)
my_map.drawparallels(np.arange(40,80,3),labels=[0,0,0,1])
my_map.drawmeridians(np.arange(-180,180,3),labels=[1,1,0,1])
plt.savefig(save_dir+'compare_getfsrs_map(all sites,>1 year)',dpi=200)
plt.show()
a=fig.add_subplot(1,3,3)
my_map = Basemap(projection='merc', lat_0 = 45, lon_0 = -63,
resolution = 'h', area_thresh = 0.3,
llcrnrlon=-68, llcrnrlat=42.0,
urcrnrlon=-59, urcrnrlat=48)
my_map.drawcoastlines()
my_map.drawcountries()
my_map.fillcontinents(color = 'coral')
my_map.drawmapboundary()
x,y=my_map(df_fsrs_1000['Longitude'].values,df_fsrs_1000['Latitude'].values)
my_map.plot(x, y, 'bo', markersize=4)
a.set_title('Longterm sites (>1000 days)',fontsize=10)
my_map.drawparallels(np.arange(40,80,3),labels=[0,0,0,1])
my_map.drawmeridians(np.arange(-180,180,3),labels=[1,1,0,1])
plt.savefig(save_dir+'compare_getfsrs_map(all sites,>1 year,> 1000 days)',dpi=200)
plt.show()