import pyproj
import georaster
import pycrs
import toolbar
from mpl_toolkits.axes_grid1 import make_axes_locatable
pathNS = 'd27m12y2016S014815'
img = georaster.SingleBandRaster(pathNS + '\image.tif')
sigNS, LaNS, LoNS, thetaNS = toolbar.readFolder(pathNS)
extent = [135, 45, 165, 63]
# extent = [ 138, 55, 148, 60 ]
fig = plt.figure(figsize=(8, 6))
ax = fig.add_axes([0.1, 0.1, 0.8, 0.8])
# setup mercator map projection.
m = toolbar.makeMap(extent)
# draw parallels
m.drawparallels(np.arange(45, 65, 2), labels=[1, 0, 0, 1], color='grey')
# draw meridians
m.drawmeridians(np.arange(135, 165, 2), labels=[1, 0, 0, 1], color='grey')
plt.imshow(image.r, extent=image.extent)
# xm,ym = m(LoNS,LaNS)
# im = plt.scatter(xm,ym,30,sigNS, marker = '.',alpha = 1,cmap = 'jet',zorder=3)
# im = plt.scatter(xm,ym,5,sigNS, marker = '.',alpha = 1,cmap = 'jet',zorder=3)
# show colorbar
# divider = make_axes_locatable(ax)
# cax = divider.append_axes("right", size="5%", pad=0.05)
# S1_titles={'10.65 GHz V-Pol','10.65 GHz H-Pol','18.7 GHz V-Pol',
# '18.7 GHz H-Pol','23.8 GHz V-Pol','36.64 GHz V-Pol',
# '36.64 GHz H-Pol','89.0 GHz V-Pol ','89.0 GHz H-Pol'};
# S2_titles={'166.0 GHz V-Pol','166.0 GHz H-Pol','183.31 +/-3 GHz V-Pol','183.31 +/-7 GHz V-Pol'};
la1 = 40
la2 = 64
lo1 = 132
lo2 = 168
extent = [120, 130, 50, 60]
fig = plt.figure(figsize=(8, 6))
ax = fig.add_axes([0.1, 0.1, 0.8, 0.8])
# setup mercator map projection.
m = toolbar.makeMap([135, 45, 160, 63])
# draw parallels
m.drawparallels(np.arange(40, 80, 10), labels=[1, 1, 0, 1])
# draw meridians
m.drawmeridians(np.arange(130, 170, 10), labels=[1, 1, 0, 1])
# print(TcS1.shape)
# print(LaS1.shape)
xm, ym = m(LoS1, LaS1)
im = plt.scatter(xm, ym, 5, TcS1[:, 6], cmap='jet')
# plt.title(1)
# divider = make_axes_locatable(ax)
# cax = divider.append_axes("bottom", size="5%", pad=0.05)
plt.colorbar(im, ax=ax, orientation='vertical')
end = time.time()