In [1]:
from charistools.hypsometry import Hypsometry
In [2]:
%cd /vagrant/source/charistools/test/compare
%ls
/vagrant/source/charistools/test/compare
00readme.txt
compare.SOI_by_elev.2001.0100m.txt
compare.test_basin.500.correct_sDDF.ALBEDO_MCD.0035.ice_melt.dat
compare.test_basin.500.correct_sDDF.ALBEDO_MCD.0035.snow_off_ice_melt.dat
compare.test_basin.500.correct_sDDF.ALBEDO_MCD.0035.snow_off_ice_melt.partial.dat
compare.test_basin.500.correct_sDDF.ALBEDO_MCD.0035.snow_on_ice_melt.dat
fromIDL/
IN_basinmask_h23v05.tif
IN_basinmask_h23v05.tif.aux.xml
IN_basinmask_h23v06.tif
IN_basinmask_h23v06.tif.aux.xml
IN_basinmask_h24v05.tif
IN_basinmask_h24v05.tif.aux.xml
IN_basinmask_h24v06.tif
IN_basinmask_h24v06.tif.aux.xml
IN_basinmask_h25v05.tif
IN_basinmask_h25v05.tif.aux.xml
IN_basinmask_h25v06.tif
IN_basinmask_h25v06.tif.aux.xml
IN_basin_outline.cpg
IN_basin_outline.dbf
IN_basin_outline.prj
IN_basin_outline.sbn
IN_basin_outline.sbx
IN_basin_outline.shp
IN_basin_outline.shx
IN_basin_outline_sinu.dbf
IN_basin_outline_sinu.prj
IN_basin_outline_sinu.shp
IN_basin_outline_sinu.shx
IN_Hunza_at_Danyour.0100m.area_by_elev.txt
IN_Hunza_at_Danyour.2001.0100m.ALBEDO_MOD10A1.0046.ablation_area.dat
IN_Hunza_at_Danyour.2001.0100m.ALBEDO_MOD10A1.0046.ice_melt.dat
IN_Hunza_at_Danyour.2001.0100m.ALBEDO_MOD10A1.0046.snow_off_ice_area.dat
IN_Hunza_at_Danyour.2001.0100m.ALBEDO_MOD10A1.0046.snow_off_ice_melt.dat
IN_Hunza_at_Danyour.2001.0100m.ALBEDO_MOD10A1.0046.snow_on_ice_area.dat
IN_Hunza_at_Danyour.2001.0100m.ALBEDO_MOD10A1.0046.snow_on_ice_melt.dat
IN_Hunza_at_Danyour.2001.0100m.GRSIZE_SCAG.0200.ablation_area.txt
IN_Hunza_at_Danyour.2001.0100m.GRSIZE_SCAG.0200.snow_off_ice_area.txt
IN_Hunza_at_Danyour.2001.0100m.GRSIZE_SCAG.0200.snow_on_ice_area.txt
In [3]:
newsoi = Hypsometry(filename="compare.SOI_by_elev.2001.0100m.txt")
In [14]:
newsoi.data.columns
Out[14]:
Index([u'2900.0', u'3000.0', u'3100.0', u'3200.0', u'3300.0', u'3400.0',
u'3500.0', u'3600.0', u'3700.0', u'3800.0', u'3900.0', u'4000.0',
u'4100.0', u'4200.0', u'4300.0', u'4400.0', u'4500.0', u'4600.0',
u'4700.0', u'4800.0', u'4900.0', u'5000.0', u'5100.0', u'5200.0',
u'5300.0', u'5400.0', u'5500.0', u'5600.0', u'5700.0', u'5800.0',
u'5900.0', u'6000.0', u'6100.0', u'6200.0', u'6300.0', u'6400.0',
u'6500.0', u'6600.0', u'6700.0', u'6800.0', u'6900.0', u'7000.0',
u'7100.0', u'7200.0', u'7300.0', u'7400.0', u'7500.0', u'7600.0',
u'7700.0'],
dtype='object')
In [5]:
newsoi.data.transpose()
Out[5]:
Date
2001-01-01 00:00:00
2001-01-02 00:00:00
2001-01-03 00:00:00
2900.0
0.18
0.17
0.15
3000.0
0.38
0.38
0.36
3100.0
0.94
0.70
0.67
3200.0
2.26
1.87
1.52
3300.0
2.68
2.24
2.18
3400.0
2.87
2.86
2.46
3500.0
3.16
2.96
2.34
3600.0
5.14
4.65
3.97
3700.0
2.78
2.70
2.64
3800.0
4.76
4.81
4.07
3900.0
5.52
4.66
3.48
4000.0
8.08
7.63
6.75
4100.0
7.62
7.49
7.20
4200.0
10.53
9.32
8.81
4300.0
16.93
15.41
13.90
4400.0
26.06
23.06
20.64
4500.0
41.56
40.83
37.92
4600.0
55.93
52.48
51.02
4700.0
62.25
60.79
59.20
4800.0
85.20
81.01
76.42
4900.0
105.36
98.21
95.69
5000.0
143.73
132.41
129.27
5100.0
195.48
183.05
175.84
5200.0
204.26
193.41
185.92
5300.0
197.90
182.02
173.97
5400.0
164.18
152.03
148.29
5500.0
120.23
115.63
109.18
5600.0
97.39
89.90
87.86
5700.0
61.68
57.69
56.29
5800.0
42.07
41.05
39.63
5900.0
36.15
32.79
31.11
6000.0
25.07
23.92
22.95
6100.0
14.09
13.89
14.14
6200.0
18.05
17.64
17.02
6300.0
12.57
12.38
11.92
6400.0
9.93
9.90
9.99
6500.0
9.55
9.12
8.56
6600.0
6.83
6.94
6.96
6700.0
7.19
7.29
7.84
6800.0
5.26
5.13
5.30
6900.0
4.80
4.75
4.73
7000.0
2.26
2.27
2.30
7100.0
3.83
3.62
3.63
7200.0
2.08
2.07
2.07
7300.0
1.61
1.62
1.63
7400.0
0.16
0.16
0.17
7500.0
1.32
1.27
1.23
7600.0
0.50
0.49
0.49
7700.0
0.00
0.00
0.00
In [6]:
%cd fromIDL
%ls
/vagrant/source/charistools/test/compare/fromIDL
compare.test_basin.500.sDDF.ALBEDO_MCD.0035.ice_melt.dat
compare.test_basin.500.sDDF.ALBEDO_MCD.0035.snow_off_ice_melt.dat
compare.test_basin.500.sDDF.ALBEDO_MCD.0035.snow_on_ice_melt.dat
IN_Hunza_at_Danyour.2001.0100m.ALBEDO_MOD10A1.0046.ablation_area.dat
IN_Hunza_at_Danyour.2001.0100m.ALBEDO_MOD10A1.0046.snow_on_ice_area.1-3.dat
IN_Hunza_at_Danyour.2001.0100m.ALBEDO_MOD10A1.0046.snow_on_ice_area.1-3.dat~
IN_Hunza_at_Danyour.2001.0100m.ALBEDO_MOD10A1.0046.snow_on_ice_area.dat
IN_Hunza_at_Danyour.2001.0100m.GRSIZE_SCAG.0200.ablation_area.dat
IN_Hunza_at_Danyour.2001.0100m.GRSIZE_SCAG.0200.snow_off_ice_area.dat
IN_Hunza_at_Danyour.2001.0100m.GRSIZE_SCAG.0200.snow_on_ice_area.dat
In [7]:
origsoi = Hypsometry(filename="IN_Hunza_at_Danyour.2001.0100m.ALBEDO_MOD10A1.0046.snow_on_ice_area.1-3.dat")
In [15]:
origsoi.data.columns
origsoi.data.columns = ["%s0" % (i) for i in origsoi.data.columns]
origsoi.data.columns
Out[15]:
Index([u'2900.0', u'3000.0', u'3100.0', u'3200.0', u'3300.0', u'3400.0',
u'3500.0', u'3600.0', u'3700.0', u'3800.0', u'3900.0', u'4000.0',
u'4100.0', u'4200.0', u'4300.0', u'4400.0', u'4500.0', u'4600.0',
u'4700.0', u'4800.0', u'4900.0', u'5000.0', u'5100.0', u'5200.0',
u'5300.0', u'5400.0', u'5500.0', u'5600.0', u'5700.0', u'5800.0',
u'5900.0', u'6000.0', u'6100.0', u'6200.0', u'6300.0', u'6400.0',
u'6500.0', u'6600.0', u'6700.0', u'6800.0', u'6900.0', u'7000.0',
u'7100.0', u'7200.0', u'7300.0', u'7400.0', u'7500.0', u'7600.0',
u'7700.0'],
dtype='object')
In [16]:
origsoi.data.transpose()
Out[16]:
Date
2001-01-01 00:00:00
2001-01-02 00:00:00
2001-01-03 00:00:00
2900.0
0.18
0.17
0.15
3000.0
0.38
0.38
0.36
3100.0
0.94
0.88
0.67
3200.0
2.26
1.87
1.52
3300.0
2.68
2.24
2.18
3400.0
2.87
2.86
2.66
3500.0
3.34
3.08
2.34
3600.0
5.14
4.65
3.97
3700.0
2.89
2.70
2.64
3800.0
5.09
4.95
4.07
3900.0
5.52
4.66
3.96
4000.0
8.65
8.10
6.88
4100.0
8.15
7.64
7.52
4200.0
10.53
9.65
8.81
4300.0
17.67
15.82
14.06
4400.0
26.26
23.22
20.97
4500.0
41.56
41.17
38.10
4600.0
56.27
52.74
51.20
4700.0
63.36
62.27
60.00
4800.0
86.58
81.88
77.61
4900.0
105.91
99.34
96.48
5000.0
147.27
134.59
131.48
5100.0
199.12
187.44
178.67
5200.0
207.44
196.26
190.63
5300.0
201.51
189.24
178.94
5400.0
168.63
156.56
153.41
5500.0
123.26
117.68
112.73
5600.0
100.00
92.14
89.51
5700.0
62.83
58.69
57.92
5800.0
43.57
41.82
40.95
5900.0
36.69
33.46
32.09
6000.0
25.25
24.40
23.11
6100.0
14.59
14.27
14.43
6200.0
18.21
17.79
17.02
6300.0
12.90
12.69
12.10
6400.0
10.11
9.90
10.16
6500.0
9.55
9.12
8.71
6600.0
6.83
6.94
7.29
6700.0
7.35
7.29
7.84
6800.0
5.26
5.13
5.30
6900.0
4.80
4.75
4.73
7000.0
2.26
2.27
2.30
7100.0
3.83
3.82
3.63
7200.0
2.08
2.07
2.07
7300.0
1.76
1.62
1.63
7400.0
0.16
0.16
0.17
7500.0
1.32
1.27
1.23
7600.0
0.50
0.49
0.49
7700.0
0.00
0.00
0.00
In [17]:
diff = newsoi.data - origsoi.data
In [18]:
diff.transpose()
Out[18]:
Date
2001-01-01 00:00:00
2001-01-02 00:00:00
2001-01-03 00:00:00
2900.0
0.00
0.00
0.00
3000.0
0.00
0.00
0.00
3100.0
0.00
-0.18
0.00
3200.0
0.00
0.00
0.00
3300.0
0.00
0.00
0.00
3400.0
0.00
0.00
-0.20
3500.0
-0.18
-0.12
0.00
3600.0
0.00
0.00
0.00
3700.0
-0.11
0.00
0.00
3800.0
-0.33
-0.14
0.00
3900.0
0.00
0.00
-0.48
4000.0
-0.57
-0.47
-0.13
4100.0
-0.53
-0.15
-0.32
4200.0
0.00
-0.33
0.00
4300.0
-0.74
-0.41
-0.16
4400.0
-0.20
-0.16
-0.33
4500.0
0.00
-0.34
-0.18
4600.0
-0.34
-0.26
-0.18
4700.0
-1.11
-1.48
-0.80
4800.0
-1.38
-0.87
-1.19
4900.0
-0.55
-1.13
-0.79
5000.0
-3.54
-2.18
-2.21
5100.0
-3.64
-4.39
-2.83
5200.0
-3.18
-2.85
-4.71
5300.0
-3.61
-7.22
-4.97
5400.0
-4.45
-4.53
-5.12
5500.0
-3.03
-2.05
-3.55
5600.0
-2.61
-2.24
-1.65
5700.0
-1.15
-1.00
-1.63
5800.0
-1.50
-0.77
-1.32
5900.0
-0.54
-0.67
-0.98
6000.0
-0.18
-0.48
-0.16
6100.0
-0.50
-0.38
-0.29
6200.0
-0.16
-0.15
0.00
6300.0
-0.33
-0.31
-0.18
6400.0
-0.18
0.00
-0.17
6500.0
0.00
0.00
-0.15
6600.0
0.00
0.00
-0.33
6700.0
-0.16
0.00
0.00
6800.0
0.00
0.00
0.00
6900.0
0.00
0.00
0.00
7000.0
0.00
0.00
0.00
7100.0
0.00
-0.20
0.00
7200.0
0.00
0.00
0.00
7300.0
-0.15
0.00
0.00
7400.0
0.00
0.00
0.00
7500.0
0.00
0.00
0.00
7600.0
0.00
0.00
0.00
7700.0
0.00
0.00
0.00
In [ ]:
Content source: mjbrodzik/ipython_notebooks
Similar notebooks: