-
Notifications
You must be signed in to change notification settings - Fork 1
/
gdal_getstacksample.py
110 lines (107 loc) · 8.97 KB
/
gdal_getstacksample.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
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
# -*- coding: utf-8 -*-
import mygdal_old2
from mydas import *
import numpy
if __name__ == '__main__':
samples = from_recarray(numpy.recfromcsv('../data/samples_new2.csv',
dtype=[('id', 'i4'),
('longitude', 'f8'),
('latitude', 'f8'),
('from', 'datetime64[D]'),
('to', 'datetime64[D]'),
('label', 'U14')]))
samples['geolocs'] = numpy.array([samples['longitude'], samples['latitude']]).T
from_srs = 'PROJCS["WGS 84 / UTM zone 21N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-57],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],AUTHORITY["EPSG","32621"],AXIS["Easting",EAST],AXIS["Northing",NORTH]]'
to_srs = 'GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]'
samples['geolocs_wgs84'] = mygdal_old2.__reproject__(samples['geolocs'], from_srs, to_srs)
samples['longitude'] = samples['geolocs_wgs84'][:, 0]
samples['latitude'] = samples['geolocs_wgs84'][:, 1]
timeline = numpy.array(["2000-02-18", "2000-03-05", "2000-03-21", "2000-04-06",
"2000-04-22", "2000-05-08", "2000-05-24", "2000-06-09",
"2000-06-25", "2000-07-11", "2000-07-27", "2000-08-12",
"2000-08-28", "2000-09-13", "2000-09-29", "2000-10-15",
"2000-10-31", "2000-11-16", "2000-12-02", "2000-12-18",
"2001-01-01", "2001-01-17", "2001-02-02", "2001-02-18",
"2001-03-06", "2001-03-22", "2001-04-07", "2001-04-23",
"2001-05-09", "2001-05-25", "2001-06-10", "2001-06-26",
"2001-07-12", "2001-07-28", "2001-08-13", "2001-08-29",
"2001-09-14", "2001-09-30", "2001-10-16", "2001-11-01",
"2001-11-17", "2001-12-03", "2001-12-19", "2002-01-01",
"2002-01-17", "2002-02-02", "2002-02-18", "2002-03-06",
"2002-03-22", "2002-04-07", "2002-04-23", "2002-05-09",
"2002-05-25", "2002-06-10", "2002-06-26", "2002-07-12",
"2002-07-28", "2002-08-13", "2002-08-29", "2002-09-14",
"2002-09-30", "2002-10-16", "2002-11-01", "2002-11-17",
"2002-12-03", "2002-12-19", "2003-01-01", "2003-01-17",
"2003-02-02", "2003-02-18", "2003-03-06", "2003-03-22",
"2003-04-07", "2003-04-23", "2003-05-09", "2003-05-25",
"2003-06-10", "2003-06-26", "2003-07-12", "2003-07-28",
"2003-08-13", "2003-08-29", "2003-09-14", "2003-09-30",
"2003-10-16", "2003-11-01", "2003-11-17", "2003-12-03",
"2003-12-19", "2004-01-01", "2004-01-17", "2004-02-02",
"2004-02-18", "2004-03-05", "2004-03-21", "2004-04-06",
"2004-04-22", "2004-05-08", "2004-05-24", "2004-06-09",
"2004-06-25", "2004-07-11", "2004-07-27", "2004-08-12",
"2004-08-28", "2004-09-13", "2004-09-29", "2004-10-15",
"2004-10-31", "2004-11-16", "2004-12-02", "2004-12-18",
"2005-01-01", "2005-01-17", "2005-02-02", "2005-02-18",
"2005-03-06", "2005-03-22", "2005-04-07", "2005-04-23",
"2005-05-09", "2005-05-25", "2005-06-10", "2005-06-26",
"2005-07-12", "2005-07-28", "2005-08-13", "2005-08-29",
"2005-09-14", "2005-09-30", "2005-10-16", "2005-11-01",
"2005-11-17", "2005-12-03", "2005-12-19", "2006-01-01",
"2006-01-17", "2006-02-02", "2006-02-18", "2006-03-06",
"2006-03-22", "2006-04-07", "2006-04-23", "2006-05-09",
"2006-05-25", "2006-06-10", "2006-06-26", "2006-07-12",
"2006-07-28", "2006-08-13", "2006-08-29", "2006-09-14",
"2006-09-30", "2006-10-16", "2006-11-01", "2006-11-17",
"2006-12-03", "2006-12-19", "2007-01-01", "2007-01-17",
"2007-02-02", "2007-02-18", "2007-03-06", "2007-03-22",
"2007-04-07", "2007-04-23", "2007-05-09", "2007-05-25",
"2007-06-10", "2007-06-26", "2007-07-12", "2007-07-28",
"2007-08-13", "2007-08-29", "2007-09-14", "2007-09-30",
"2007-10-16", "2007-11-01", "2007-11-17", "2007-12-03",
"2007-12-19", "2008-01-01", "2008-01-17", "2008-02-02",
"2008-02-18", "2008-03-05", "2008-03-21", "2008-04-06",
"2008-04-22", "2008-05-08", "2008-05-24", "2008-06-09",
"2008-06-25", "2008-07-11", "2008-07-27", "2008-08-12",
"2008-08-28", "2008-09-13", "2008-09-29", "2008-10-15",
"2008-10-31", "2008-11-16", "2008-12-02", "2008-12-18",
"2009-01-01", "2009-01-17", "2009-02-02", "2009-02-18",
"2009-03-06", "2009-03-22", "2009-04-07", "2009-04-23",
"2009-05-09", "2009-05-25", "2009-06-10", "2009-06-26",
"2009-07-12", "2009-07-28", "2009-08-13", "2009-08-29",
"2009-09-14", "2009-09-30", "2009-10-16", "2009-11-01",
"2009-11-17", "2009-12-03", "2009-12-19", "2010-01-01",
"2010-01-17", "2010-02-02", "2010-02-18", "2010-03-06",
"2010-03-22", "2010-04-07", "2010-04-23", "2010-05-09",
"2010-05-25", "2010-06-10", "2010-06-26", "2010-07-12",
"2010-07-28", "2010-08-13", "2010-08-29", "2010-09-14",
"2010-09-30", "2010-10-16", "2010-11-01", "2010-11-17",
"2010-12-03", "2010-12-19", "2011-01-01", "2011-01-17",
"2011-02-02", "2011-02-18", "2011-03-06", "2011-03-22",
"2011-04-07", "2011-04-23", "2011-05-09", "2011-05-25",
"2011-06-10", "2011-06-26", "2011-07-12", "2011-07-28",
"2011-08-13", "2011-08-29", "2011-09-14", "2011-09-30",
"2011-10-16", "2011-11-01", "2011-11-17", "2011-12-03",
"2011-12-19", "2012-01-01", "2012-01-17", "2012-02-02",
"2012-02-18", "2012-03-05", "2012-03-21", "2012-04-06",
"2012-04-22", "2012-05-08", "2012-05-24", "2012-06-09",
"2012-06-25", "2012-07-11", "2012-07-27", "2012-08-12",
"2012-08-28", "2012-09-13", "2012-09-29", "2012-10-15",
"2012-10-31", "2012-11-16", "2012-12-02", "2012-12-18",
"2013-01-01", "2013-01-17", "2013-02-02", "2013-02-18",
"2013-03-06", "2013-03-22", "2013-04-07", "2013-04-23",
"2013-05-09", "2013-05-25", "2013-06-10", "2013-06-26",
"2013-07-12", "2013-07-28", "2013-08-13", "2013-08-29",
"2013-09-14", "2013-09-30", "2013-10-16", "2013-11-01",
"2013-11-17", "2013-12-03", "2013-12-19", "2014-01-01",
"2014-01-17", "2014-02-02", "2014-02-18", "2014-03-06",
"2014-03-22", "2014-04-07", "2014-04-23", "2014-05-09"], dtype='M8[D]')
for i in range(rows(samples)):
samples['from'][i] = timeline[timeline <= samples['from'][i]][-1]
for i in range(rows(samples)):
samples['to'][i] = timeline[timeline >= samples['to'][i]][0]
r = to_recarray(samples, ['id', 'longitude', 'latitude', 'from', 'to', 'label'])
numpy.savetxt('../data/samples_wgs84.csv', r,
header='id,longitude,latitude,from,to,label', fmt='%d,%.18f,%.18f,%s,%s,%s')