In [1]:
import numpy as np
In [2]:
arr=np.array([[2,45,6,34], [7,4,12,76], [3,5,0,13]])
arr, arr.shape
Out[2]:
Finn dette i arr:
- Første element?
- Første rad?
- Tredje kolonne?
- Første tre kolonner i to siste rader?
In [3]:
arr[0,0], arr[0,:], arr[:,2],
arr[1:, :3]
Out[3]:
Bruke funksjoner:¶
Husk at brakkeparantes brukes for indeksering, mens vanlig parantes brukes til funksjoner.
In [4]:
# Gjennomsnitt av hele "arr"
np.mean(arr)
Out[4]:
In [5]:
# Gjennomsnitt av første rad
np.mean(arr[0,:])
Out[5]:
In [6]:
# Gjennomsnitt av rad nr 2 og 3, kolonne nr 1.
np.mean(arr[1:,:2])
Out[6]:
Laste inn data og lagre dataene til variabler¶
In [7]:
import sys
# Bytt ut stien under slik at den peker på hvor på din PC du har
# lagret mappen "Funksjoner".
sys.path.append(r"W:\Work\Documents\EkteData\EkteData\Funksjoner")
from EkteDataFunk import loadData # for filinnlasting
In [8]:
file='TempBergenYearlyNonan.txt'
data = loadData(file, rydd='N')
data.shape, data
Out[8]:
In [9]:
tid=data[:,0]
temp=data[:,1]
Velge ut en spesifikk tidsperiode¶
fra en tidsserie med årlig oppløsning.
In [10]:
np.where(tid==1862)
Out[10]:
In [11]:
tid[1]
Out[11]:
In [12]:
np.where(tid==1873)
Out[12]:
In [13]:
tid[12]
Out[13]:
In [14]:
np.where((tid>=1870) & (tid<=1893))
Out[14]:
In [15]:
tid[9:32+1]
Out[15]:
Lage figurer¶
In [16]:
import matplotlib.pyplot as plt
In [17]:
plt.plot(tid, temp)
plt.title('Temperatur Bergen')
plt.xlabel('År')
plt.ylabel('Temperatur [\u2103]')
plt.show()
Plot et spesifikk tidsperiode
In [18]:
plt.plot(tid[1:12], temp[1:12])
plt.title('Temperatur Bergen')
plt.xlabel('År')
plt.ylabel('Temperatur [\u2103]')
plt.show()
Plot flere tidsserier:
- hele tidsserien,
- hver femte verdi,
- verdiene fra indeks 50 til og med 149.
In [19]:
plt.plot(tid, temp)
plt.plot(tid[::5], temp[::5])
plt.plot(tid[50:150], temp[50:150])
plt.title('Temperatur Bergen')
plt.xlabel('År')
plt.ylabel('Temperatur [\u2103]')
plt.show()
Feilmeldinger¶
In [20]:
tid[2.0:4]
In [21]:
tid(2:4)
In [22]:
tid[2:4
In [23]:
plt.plot(tid[3:],temp[5:])
In [24]:
plt.plot(tid,temp)
plt.xlabel('tid)
plt.ylabel('temp')
plt.title('Årlig temperatur fra 1861 til 2020')
Dele en tidserie inn i like lange bolker¶
In [25]:
ind=[
0,
len(tid)*1/4,
len(tid)*2/4,
len(tid)*3/4,
len(tid)
]
In [26]:
ind
Out[26]:
In [27]:
ind=[
0,
int(len(tid)*1/4),
int(len(tid)*2/4),
int(len(tid)*3/4),
len(tid)
]
ind
Out[27]:
In [28]:
plt.plot(tid[ind[0]:ind[1]], temp[ind[0]:ind[1]])
plt.plot(tid[ind[1]:ind[2]], temp[ind[1]:ind[2]])
plt.plot(tid[ind[2]:ind[3]], temp[ind[2]:ind[3]])
plt.plot(tid[ind[3]:ind[4]], temp[ind[3]:ind[4]])
plt.title('Temperatur Bergen')
plt.xlabel('År')
plt.ylabel('Temperatur [\u2103]')
plt.show()
In [ ]:
Last ned og prøv selv:
Forfatter
Vår Iren Hjorth Dundas