10/11/24, 18:89 PC3_SPATIAL_JATENG_LISMASARI.
ipynb - Colab
PM
!pip install xarray
!pip install matplotlib
!pip install eofs
!pip install cartopy
!pip install Basemap
!pip install netcdf4
[Link] 1/9
10/11/24, 18:89 PC3_SPATIAL_JATENG_LISMASARI.ipynb - Colab
PM
Requirement already satisfied: xarray in /usr/local/lib/python3.10/dist-packages (2024.10.0)
Requirement already satisfied: numpy>=1.24 in /usr/local/lib/python3.10/dist-packages (from xarray) (1.26.4)
Requirement already satisfied: packaging>=23.1 in /usr/local/lib/python3.10/dist-packages (from xarray) (24.1)
Requirement already satisfied: pandas>=2.1 in /usr/local/lib/python3.10/dist-packages (from xarray) (2.2.2)
Requirement already satisfied: python-dateutil>=2.8.2 in /usr/local/lib/python3.10/dist-packages (from pandas>=2.1->xarray) (2.8.2)
Requirement already satisfied: pytz>=2020.1 in /usr/local/lib/python3.10/dist-packages (from pandas>=2.1->xarray) (2024.2)
Requirement already satisfied: tzdata>=2022.7 in /usr/local/lib/python3.10/dist-packages (from pandas>=2.1->xarray) (2024.2)
Requirement already satisfied: six>=1.5 in /usr/local/lib/python3.10/dist-packages (from python-dateutil>=2.8.2->pandas>=2.1->xarray)
Requirement already satisfied: matplotlib in /usr/local/lib/python3.10/dist-packages (3.8.0)
Requirement already satisfied: contourpy>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (1.3.0)
Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (0.12.1)
Requirement already satisfied: fonttools>=4.22.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (4.54.1)
Requirement already satisfied: kiwisolver>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (1.4.7)
Requirement already satisfied: numpy<2,>=1.21 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (1.26.4)
Requirement already satisfied: packaging>=20.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (24.1)
Requirement already satisfied: pillow>=6.2.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (10.4.0)
Requirement already satisfied: pyparsing>=2.3.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (3.2.0)
Requirement already satisfied: python-dateutil>=2.7 in /usr/local/lib/python3.10/dist-packages (from matplotlib) (2.8.2)
Requirement already satisfied: six>=1.5 in /usr/local/lib/python3.10/dist-packages (from python-dateutil>=2.7->matplotlib) (1.16.0)
Collecting eofs
Downloading [Link] (3.6 kB)
Requirement already satisfied: numpy in /usr/local/lib/python3.10/dist-packages (from eofs) (1.26.4)
Downloading [Link] (1.1 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 1.1/1.1 MB 6.3 MB/s eta [Link]
Installing collected packages: eofs
Successfully installed eofs-1.4.1
Collecting cartopy
Downloading Cartopy-0.24.1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.[Link] (7.9 kB)
Requirement already satisfied: numpy>=1.23 in /usr/local/lib/python3.10/dist-packages (from cartopy) (1.26.4)
Requirement already satisfied: matplotlib>=3.6 in /usr/local/lib/python3.10/dist-packages (from cartopy) (3.8.0)
Requirement already satisfied: shapely>=1.8 in /usr/local/lib/python3.10/dist-packages (from cartopy) (2.0.6)
Requirement already satisfied: packaging>=21 in /usr/local/lib/python3.10/dist-packages (from cartopy) (24.1)
Requirement already satisfied: pyshp>=2.3 in /usr/local/lib/python3.10/dist-packages (from cartopy) (2.3.1)
Requirement already satisfied: pyproj>=3.3.1 in /usr/local/lib/python3.10/dist-packages (from cartopy) (3.7.0)
Requirement already satisfied: contourpy>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (1.3.0)
Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (0.12.1)
Requirement already satisfied: fonttools>=4.22.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (4.54.1)
Requirement already satisfied: kiwisolver>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (1.4.7)
Requirement already satisfied: pillow>=6.2.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (10.4.0)
Requirement already satisfied: pyparsing>=2.3.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (3.2.0)
Requirement already satisfied: python-dateutil>=2.7 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=3.6->cartopy) (2.8.2
Requirement already satisfied: certifi in /usr/local/lib/python3.10/dist-packages (from pyproj>=3.3.1->cartopy) (2024.8.30)
Requirement already satisfied: six>=1.5 in /usr/local/lib/python3.10/dist-packages (from python-dateutil>=2.7->matplotlib>=3.6->carto
Downloading Cartopy-0.24.1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (11.7 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 11.7/11.7 MB 59.2 MB/s eta [Link]
Installing collected packages: cartopy
Successfully installed cartopy-0.24.1
Collecting Basemap
Downloading basemap-1.4.1-cp310-cp310-manylinux1_x86_64.[Link] (9.1 kB)
Collecting basemap-data<1.4,>=1.3.2 (from Basemap)
Downloading basemap_data-[Link] (2.7 kB)
Requirement already satisfied: pyshp<2.4,>=1.2 in /usr/local/lib/python3.10/dist-packages (from Basemap) (2.3.1)
Requirement already satisfied: matplotlib<3.9,>=1.5 in /usr/local/lib/python3.10/dist-packages (from Basemap) (3.8.0)
Collecting pyproj<3.7.0,>=1.9.3 (from Basemap)
Downloading pyproj-3.6.1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.[Link] (31 kB)
Collecting packaging<24.0,>=16.0 (from Basemap)
Downloading [Link] (3.2 kB)
Requirement already satisfied: numpy<1.27,>=1.21 in /usr/local/lib/python3.10/dist-packages (from Basemap) (1.26.4)
Requirement already satisfied: contourpy>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (1.3.
Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (0.12.1)
Requirement already satisfied: fonttools>=4.22.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (4.5
Requirement already satisfied: kiwisolver>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (1.4
Requirement already satisfied: pillow>=6.2.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (10.4.0)
Requirement already satisfied: pyparsing>=2.3.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (3.2.
Requirement already satisfied: python-dateutil>=2.7 in /usr/local/lib/python3.10/dist-packages (from matplotlib<3.9,>=1.5->Basemap) (
Requirement already satisfied: certifi in /usr/local/lib/python3.10/dist-packages (from pyproj<3.7.0,>=1.9.3->Basemap) (2024.8.30)
Requirement already satisfied: six>=1.5 in /usr/local/lib/python3.10/dist-packages (from python-dateutil>=2.7->matplotlib<3.9,>=1.5->
Downloading basemap-1.4.1-cp310-cp310-manylinux1_x86_64.whl (935 kB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 936.0/936.0 kB 16.7 MB/s eta [Link]
Downloading basemap_data-[Link] (30.5 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 30.5/30.5 MB 31.3 MB/s eta [Link]
Downloading [Link] (53 kB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 53.0/53.0 kB 3.4 MB/s eta [Link]
Downloading pyproj-3.6.1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (8.3 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 8.3/8.3 MB 62.4 MB/s eta [Link]
Installing collected packages: pyproj, packaging, basemap-data, Basemap
Attempting uninstall: pyproj
Found existing installation: pyproj 3.7.0
Uninstalling pyproj-3.7.0:
Successfully uninstalled pyproj-3.7.0
Attempting uninstall: packaging
Found existing installation: packaging 24.1
Uninstalling packaging-24.1:
Successfully uninstalled packaging-24.1
[Link] 2/9
10/11/24, 18:89 PC3_SPATIAL_JATENG_LISMASARI.ipynb - Colab
PM Successfully installed Basemap-1.4.1 basemap-data-1.3.2 packaging-23.2 pyproj-3.6.1
Collecting netcdf4
Downloading netCDF4-1.7.2-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.[Link] (1.8 kB)
Collecting cftime (from netcdf4)
Downloading cftime-1.6.4.post1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.[Link] (8.7 kB)
Requirement already satisfied: certifi in /usr/local/lib/python3.10/dist-packages (from netcdf4) (2024.8.30)
Requirement already satisfied: numpy in /usr/local/lib/python3.10/dist-packages (from netcdf4) (1.26.4)
Downloading netCDF4-1.7.2-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (9.1 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 9.1/9.1 MB 49.3 MB/s eta [Link]
Downloading cftime-1.6.4.post1-cp310-cp310-manylinux_2_17_x86_64.manylinux2014_x86_64.whl (1.3 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 1.3/1.3 MB 57.9 MB/s eta [Link]
Installing collected packages: cftime, netcdf4
import numpy as np
import [Link] as ma
import xarray as xr
import [Link] as plt
import matplotlib as mpl
import [Link] as colors
import [Link] as mpatches
from mpl_toolkits.basemap import Basemap
from [Link] import fft, fftfreq, ifft
import pandas as pd
from [Link] import *
#import [Link] as ccrs
#import [Link] as cfeature
#from scipy import signal
#from [Link] import Eof
#[Link]['[Link]'] = 'deJavu Sans'
#[Link]['[Link]'] = 35 #22
#[Link]['[Link]'] = 'bold' #22
#[Link]['[Link]'] = [35,10]
import [Link] as ticker
from [Link] import drive
[Link]('/content/drive')
Mounted at /content/drive
pcs=xr.open_dataset("/content/drive/MyDrive/Fisika Komputasi/Mutiara_Jateng_PC3.nc")
pcs
[Link]
► Dimensions: (time: 371, latitude: 281, longitude: 841)
▼ Coordinates:
time (time) datetime64[ns] 1993-01-01 ... 2023-01-12
latitude (latitude) float32 35.0 34.75 34.5 ... -34.75 -35.0
longitude (longitude) float32 40.0 40.25 40.5 ... 249.8 250.0
▼ Data variables:
rr_spasial (time, latitude, longitude) float64 ...
► Indexes: (3)
▼ Attributes:
RR : Data Rata Rata Stasiun
ds_sst_slp_wind_era5=xr.open_dataset("/content/drive/MyDrive/Fisika Komputasi/FULL_LONG-[Link]")
t_start='1993-01-01'
t_stop='2023-12-01'
min_lon = 40
max_lon = 250
min_lat = -35
max_lat = 35
[Link] 3/9
10/11/24, 18:89 PC3_SPATIAL_JATENG_LISMASARI.ipynb - Colab
PM #ds_uwind_sel = ds_uwind.sel(time=slice(t_start, t_stop), lat=slice(max_lat,min_lat), lon=slice(min_lon,max_lon))
#ds_vwind_sel = ds_vwind.sel(time=slice(t_start, t_stop), lat=slice(max_lat,min_lat), lon=slice(min_lon,max_lon))
#ds_slp_sel = ds_slp.sel(time=slice(t_start, t_stop), lat=slice(max_lat,min_lat), lon=slice(min_lon,max_lon))
ds_era5_sel = ds_sst_slp_wind_era5.sel(time=slice(t_start, t_stop), latitude=slice(max_lat,min_lat), longitude=slice(min_lon,max_lon))
#ds_era5_sel
#ds_era5_sel_climo = ds_era5_sel.groupby('[Link]').mean('time')
#ds_era5_sel_anoms = (ds_era5_sel.groupby('[Link]') - ds_era5_sel_climo)
ds_era5_sel
[Link]
► Dimensions: (time: 360, latitude: 281, longitude: 841)
▼ Coordinates:
longitude (longitude) float32 40.0 40.25 40.5 ... 249.8 250.0
latitude (latitude) float32 35.0 34.75 34.5 ... -34.75 -35.0
time (time) datetime64[ns] 1993-01-01 ... 2022-12-01
▼ Data variables:
u10 (time, latitude, longitude) float32 ...
v10 (time, latitude, longitude) float32 ...
msl (time, latitude, longitude) float32 ...
sst (time, latitude, longitude) float32 ...
► Indexes: (3)
▼ Attributes:
Conventions : CF-1.6
history : 2023-03-10 [Link] GMT by grib_to_netcdf-2.25.1: /opt/ecmwf/mars-client/bin/grib_to
_netcdf.bin -S param -o /cache/data6/[Link]-1678431328.9136436-787
[Link] /cache/tmp/7e6d2ff1-c136-4379-91b1
-[Link]
BUAT SEASON (ASO dan NDJ) UNTUK MASING-MASING DATA
s_y=1993
e_y=2023
s_ASO='-08'
e_ASO='-10'
list_sASO=[]#str(s_y)+'-01']
list_eASO=[]#str(s_y)+'-02']
list_sesASO=[]
for yr in range(s_y,e_y):
list_sASO.append(str(yr)+s_ASO)
list_eASO.append(str(yr)+e_ASO)
list_sesASO.append('ASO'+str(yr))
#list_sJJA.append([Link])
#list_eJJA.append([Link])
#list_sesJJA.append([Link])
s_NDJ='-11'
e_NDJ='-01'
list_sNDJ=[]#str(s_y)+'-01']
list_eNDJ=[]#str(s_y)+'-02']
list_sesNDJ=[]
for yr in range(s_y,e_y):
list_sNDJ.append(str(yr)+s_NDJ)
if (yr==(e_y-1)):
list_eNDJ.append(str(yr)+'-12')
else:
list_eNDJ.append(str(yr+1)+e_NDJ)
list_sesNDJ.append('NDJ'+str(yr))
[Link] 4/9
10/11/24, 18:89 PC3_SPATIAL_JATENG_LISMASARI.ipynb - Colab
PM #list_sNDJ.append([Link])
#list_eNDJ.append([Link])
#list_sesNDJ.append([Link])
df_list=[Link](list(zip(list_sASO,list_eASO,list_sesASO,
list_sNDJ,list_eNDJ,list_sesNDJ)),columns=['sASO','eASO','sesASO',
'sNDJ','eNDJ','sesNDJ'])
df_list
sASO eASO sesASO sNDJ eNDJ sesNDJ
0 1993-08 1993-10 ASO1993 1993-11 1994-01 NDJ1993
1 1994-08 1994-10 ASO1994 1994-11 1995-01 NDJ1994
2 1995-08 1995-10 ASO1995 1995-11 1996-01 NDJ1995
3 1996-08 1996-10 ASO1996 1996-11 1997-01 NDJ1996
4 1997-08 1997-10 ASO1997 1997-11 1998-01 NDJ1997
5 1998-08 1998-10 ASO1998 1998-11 1999-01 NDJ1998
6 1999-08 1999-10 ASO1999 1999-11 2000-01 NDJ1999
7 2000-08 2000-10 ASO2000 2000-11 2001-01 NDJ2000
8 2001-08 2001-10 ASO2001 2001-11 2002-01 NDJ2001
9 2002-08 2002-10 ASO2002 2002-11 2003-01 NDJ2002
10 2003-08 2003-10 ASO2003 2003-11 2004-01 NDJ2003
11 2004-08 2004-10 ASO2004 2004-11 2005-01 NDJ2004
12 2005-08 2005-10 ASO2005 2005-11 2006-01 NDJ2005
13 2006-08 2006-10 ASO2006 2006-11 2007-01 NDJ2006
14 2007-08 2007-10 ASO2007 2007-11 2008-01 NDJ2007
15 2008-08 2008-10 ASO2008 2008-11 2009-01 NDJ2008
16 2009-08 2009-10 ASO2009 2009-11 2010-01 NDJ2009
17 2010-08 2010-10 ASO2010 2010-11 2011-01 NDJ2010
18 2011-08 2011-10 ASO2011 2011-11 2012-01 NDJ2011
19 2012-08 2012-10 ASO2012 2012-11 2013-01 NDJ2012
20 2013-08 2013-10 ASO2013 2013-11 2014-01 NDJ2013
21 2014-08 2014-10 ASO2014 2014-11 2015-01 NDJ2014
22 2015-08 2015-10 ASO2015 2015-11 2016-01 NDJ2015
23 2016-08 2016-10 ASO2016 2016-11 2017-01 NDJ2016
24 2017-08 2017-10 ASO2017 2017-11 2018-01 NDJ2017
25 2018-08 2018-10 ASO2018 2018-11 2019-01 NDJ2018
26 2019-08 2019-10 ASO2019 2019-11 2020-01 NDJ2019
27 2020-08 2020-10 ASO2020 2020-11 2021-01 NDJ2020
28 2021-08 2021-10 ASO2021 2021-11 2022-01 NDJ2021
29 2022-08 2022-10 ASO2022 2022-11 2022-12 NDJ2022
SEASON data SST
d1=ds_era5_sel.sel(time=slice(df_list.iloc[0]['sASO'], df_list.iloc[0]['eASO'])).mean(dim='time')
d11 = d1.assign_coords({'season':df_list.iloc[0]['sesASO']})
seas_SST_ASO_TS = d11.expand_dims({"season": 1})
for i in range(len(df_list)):
d1=ds_era5_sel.sel(time=slice(df_list.iloc[i]['sASO'], df_list.iloc[i]['eASO'])).mean(dim='time')
d11 = d1.assign_coords({'season':df_list.iloc[i]['sesASO']})
d111 = d11.expand_dims({"season": 1})
seas_SST_ASO_TS=[Link]([seas_SST_ASO_TS, d111], compat="no_conflicts", join="outer")
[Link] 5/9
10/11/24, 18:89 PC3_SPATIAL_JATENG_LISMASARI.ipynb - Colab
PM d2=ds_era5_sel.sel(time=slice(df_list.iloc[0]['sNDJ'], df_list.iloc[0]['eNDJ'])).mean(dim='time')
d22 = d2.assign_coords({'season':df_list.iloc[0]['sesNDJ']})
seas_SST_NDJ_TS = d22.expand_dims({"season": 1})
for i in range(1,len(df_list)-1):
d2=ds_era5_sel.sel(time=slice(df_list.iloc[i]['sNDJ'], df_list.iloc[i]['eNDJ'])).mean(dim='time')
d22 = d2.assign_coords({'season':df_list.iloc[i]['sesNDJ']})
d222 = d22.expand_dims({"season": 1})
seas_SST_NDJ_TS=[Link]([seas_SST_NDJ_TS, d222], compat="no_conflicts", join="outer")
SEASON PC
d5=[Link](time=slice(df_list.iloc[0]['sASO'], df_list.iloc[0]['eASO'])).mean(dim='time')
d55 = d5.assign_coords({'season':df_list.iloc[0]['sesASO']})
seas_PC_ASO_TS = d55.expand_dims({"season": 1})
for i in range(len(df_list)):
d5=[Link](time=slice(df_list.iloc[i]['sASO'], df_list.iloc[i]['eASO'])).mean(dim='time')
d55 = d5.assign_coords({'season':df_list.iloc[i]['sesASO']})
d555 = d55.expand_dims({"season": 1})
seas_PC_ASO_TS=[Link]([seas_PC_ASO_TS, d555], compat="no_conflicts", join="outer")
d6=[Link](time=slice(df_list.iloc[0]['sNDJ'], df_list.iloc[0]['eNDJ'])).mean(dim='time')
d66 = d6.assign_coords({'season':df_list.iloc[0]['sesNDJ']})
seas_PC_NDJ_TS = d66.expand_dims({"season": 1})
for i in range(len(df_list)):
d6=[Link](time=slice(df_list.iloc[i]['sNDJ'], df_list.iloc[i]['eNDJ'])).mean(dim='time')
d66 = d6.assign_coords({'season':df_list.iloc[i]['sesNDJ']})
d666 = d66.expand_dims({"season": 1})
seas_PC_NDJ_TS=[Link]([seas_PC_NDJ_TS, d666], compat="no_conflicts", join="outer")
HITUNG KORELASI
from [Link] import t
import math as mt
import numpy as np
def val_p_val(r,n,sig_lev):
if([Link](r)):
return [Link]
else:
p_tabel=1-sig_lev
t_hit = (r*[Link]((n-2))) / [Link]((1-r**2))
dff=n-2
pval=(1 - [Link](x=t_hit, df=dff))*2
if (pval>p_tabel):
return r
else:
return [Link]
cor_pc3_sst=[Link](pcs.rr_spasial, ds_era5_sel.sst-273, dim='time')
cor_pc3_sst
[Link] 6/9
10/11/24, 18:89 PC3_SPATIAL_JATENG_LISMASARI.ipynb - Colab
PM
[Link] (latitude: 281, longitude: 841)
array([[ nan, nan, nan, ..., nan,
nan, nan],
[ nan, nan, nan, ..., nan,
nan, nan],
[ nan, nan, nan, ..., nan,
nan, nan],
...,
[ 0.2097926 , 0.18411746, 0.18185603, ..., -0.12998584,
-0.13428626, -0.12888835],
[ 0.19159282, 0.17804667, 0.17573168, ..., -0.10697331,
-0.11233453, -0.11247877],
[ 0.20367306, 0.19419666, 0.17535955, ..., -0.08243089,
-0.08828885, -0.09047806]])
▼ Coordinates:
latitude (latitude) float32 35.0 34.75 34.5 ... -34.75 -35.0
longitude (longitude) float32 40.0 40.25 40.5 ... 249.8 250.0
► Indexes: (2)
► Attributes: (0)
[Link]['[Link]'] = 'deJavu Sans'
[Link]['[Link]'] = 10 #22
[Link]['[Link]'] = 'bold' #22
# create new figure, axes instances.
fig=[Link](figsize=(7,4))
ax=fig.add_axes([0.1,0.1,0.8,0.8])
# setup mercator map projection.
m = Basemap(llcrnrlon=min_lon,llcrnrlat=min_lat,urcrnrlon=max_lon,urcrnrlat=max_lat,\
resolution='i', projection='cyl')
[Link]()
[Link]()
# draw parallels
[Link]([Link](min_lat,max_lat,15),labels=[1,0,0,0])
# draw meridians
[Link]([Link](min_lon,max_lon,20),labels=[0,0,0,1])
###################################################################
###################################################################
sekmen=7
#max_sst=cor_pc1_sst.max().values*1
#min_sst=cor_pc1_sst.min().values*1
max_sst=1
min_sst=-1
px,py = [Link](cor_pc3_sst['longitude'], cor_pc3_sst['latitude'])
levels1 = [Link](min_sst, max_sst, sekmen)
x, y = m(px, py)
cmap = plt.get_cmap('RdYlBu_r', sekmen)
[Link](x, y, cor_pc3_sst, cmap = cmap, levels=levels1, vmin = min_sst, vmax = max_sst)
bounds = [-4,-2,-1.0, -0.5, 0.5, 1.0,2,4]
norm = [Link](levels1, cmap.N)
[Link](
[Link](cmap=cmap, norm=norm),
#cax=cax,
#extend='both',
extendfrac='auto',
ticks=levels1,
spacing='uniform',
shrink=0.2, pad=0.1
)
[Link]("Correlation of SM Maritim Tanjung Emas and SST")
[Link] 7/9
10/11/24, 18:89 PC3_SPATIAL_JATENG_LISMASARI.ipynb - Colab
PM
Text(0.5, 1.0, 'Correlation of SM Maritim Tanjung Emas and SST')
[Link]['[Link]'] = 'deJavu Sans'
[Link]['[Link]'] = 10 #22
[Link]['[Link]'] = 'bold' #22
# create new figure, axes instances.
fig=[Link](figsize=(7,4))
ax=fig.add_axes([0.1,0.1,0.8,0.8])
# setup mercator map projection.
m = Basemap(llcrnrlon=min_lon,llcrnrlat=min_lat,urcrnrlon=max_lon,urcrnrlat=max_lat,\
resolution='i', projection='cyl')
[Link]()
[Link]()
# draw parallels
[Link]([Link](min_lat,max_lat,15),labels=[1,0,0,0])
# draw meridians
[Link]([Link](min_lon,max_lon,20),labels=[0,0,0,1])
###################################################################
###################################################################
sekmen=7
max_sst=cor_pc3_sst.max().values*1
min_sst=cor_pc3_sst.min().values*1
px,py = [Link](cor_pc3_sst['longitude'], cor_pc3_sst['latitude'])
levels1 = [Link](min_sst, max_sst, sekmen)
x, y = m(px, py)
cmap = plt.get_cmap('RdYlBu_r', sekmen)
[Link](x, y, cor_pc3_sst, cmap = cmap, levels=levels1, vmin = min_sst, vmax = max_sst)
bounds = [-4,-2,-1.0, -0.5, 0.5, 1.0,2,4]
norm = [Link](levels1, cmap.N)
[Link](
[Link](cmap=cmap, norm=norm),
#cax=cax,
#extend='both',
extendfrac='auto',
ticks=levels1,
spacing='uniform',
shrink=0.2, pad=0.1
)
[Link]("Correlation of SM Maritim Tanjung Emas and SST")
Text(0.5, 1.0, 'Correlation of SM Maritim Tanjung Emas and SST')
n_all=len(ds_era5_sel.time)
sig_lev=0.95
for i in range(len(cor_pc3_sst.latitude)):
for j in range(len(cor_pc3_sst.longitude)):
[Link] 8/9
10/11/24, 18:89 PC3_SPATIAL_JATENG_LISMASARI.ipynb - Colab
PM r=cor_pc3_sst[i][j]
[Link]['[Link]'] = 'deJavu Sans'
[Link]['[Link]'] = 10 #22
[Link]['[Link]'] = 'bold' #22
# create new figure, axes instances.
fig=[Link](figsize=(7,4))
ax=fig.add_axes([0.1,0.1,0.8,0.8])
# setup mercator map projection.
m = Basemap(llcrnrlon=min_lon,llcrnrlat=min_lat,urcrnrlon=max_lon,urcrnrlat=max_lat,\
resolution='i', projection='cyl')
[Link]()
[Link]()
# draw parallels
[Link]([Link](min_lat,max_lat,15),labels=[1,0,0,0])
# draw meridians
[Link]([Link](min_lon,max_lon,20),labels=[0,0,0,1])
###################################################################
###################################################################
sekmen=7
max_sst=cor_pc3_sst.max().values*1
min_sst=cor_pc3_sst.min().values*1
px,py = [Link](cor_pc3_sst['longitude'], cor_pc3_sst['latitude'])
levels1 = [Link](min_sst, max_sst, sekmen)
x, y = m(px, py)
cmap = plt.get_cmap('RdYlBu_r', sekmen)
[Link](x, y, cor_pc3_sst, cmap = cmap, levels=levels1, vmin = min_sst, vmax = max_sst)
bounds = [-4,-2,-1.0, -0.5, 0.5, 1.0,2,4]
norm = [Link](levels1, cmap.N)
[Link](
[Link](cmap=cmap, norm=norm),
#cax=cax,
[Link] 9/9
