← Back to branch summary

~ubuntu-branches/ubuntu/trusty/liblas/trusty-proposed

« back to all changes in this revision

Viewing changes to python/tests/SRS-GDAL.txt

  • Committer: Package Import Robot
  • Author(s): Francesco Paolo Lovergine
  • Date: 2014-01-05 17:00:29 UTC
  • mfrom: (7.1.2 sid)
  • Revision ID: package-import@ubuntu.com-20140105170029-ddtp0j63x5jvck2u
Tags: 1.7.0+dfsg-2
http://bugs.debian.org/733282
Fixed missing linking of system boost component.
(closes: #733282)

Show diffs side-by-side

added added

removed removed

Lines of Context:
python/tests/SRS-GDAL.txt
 
1
  >>> from liblas import srs
 
2
  >>> from liblas import point
 
3
  >>> from liblas import header
 
4
  
 
5
  >>> import liblas
 
6
  >>> s = srs.SRS()
 
7
  >>> s.proj4
 
8
  ''
 
9
  >>> s.proj4 = '+proj=utm +zone=15 +datum=WGS84 +units=m +no_defs'
 
10
  >>> s.proj4 == '+proj=utm +zone=15 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs '
 
11
  True
 
12
  
 
13
  >>> s = srs.SRS()
 
14
  >>> s.set_userinput('EPSG:4326')
 
15
  True
 
16
  >>> s.proj4 == '+proj=longlat +datum=WGS84 +no_defs '
 
17
  True
 
18
  
 
19
  >>> from liblas import file
 
20
  >>> f = file.File('../test/data/1.2_3.las',mode='r')
 
21
  >>> s = f.header.srs
 
22
  >>> s.wkt == """PROJCS["NAD83 / UTM zone 15N",GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS 1980",6378137,298.2572221010002,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4269"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-93],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","26915"]]"""
 
23
  True
 
24
  
 
25
  >>> s2 = srs.SRS()
 
26
  >>> s2.wkt = """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"]]"""
 
27
  >>> p = f.read(0)
 
28
  >>> p.x
 
29
  470692.44
 
30
  >>> f.set_srs(s2)
 
31
  True
 
32
  >>> p = f.read(0)
 
33
  >>> s2.vlr_count()
 
34
  4
 
35
  >>> s2.GetVLR(0).recordlength
 
36
  64
 
37
  >>> int(round(p.x))
 
38
  -93
 
39
  >>> int(round(p.y))
 
40
  42
 
41
 
 
42
  >>> del f
 
43
 
 
44
 
 
45
# -93.3515625902 41.5771483954
 
46
 
 
47
  >>> def new_offset(old_scale, new_scale, old_offset, x):
 
48
  ...       return (new_scale*(x - old_offset) - old_scale*x)/(-1.0*old_scale)
 
49
  
 
50
  >>> utm_wkt = """PROJCS["NAD83 / UTM zone 15N",GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS 1980",6378137,298.2572221010002,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4269"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-93],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","26915"]]"""
 
51
 
 
52
  >>> dd_wkt = """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"]]"""
 
53
  >>> s_dd = srs.SRS()
 
54
  >>> s_dd.wkt = dd_wkt
 
55
  >>> s_utm = srs.SRS()
 
56
  >>> s_utm.wkt = utm_wkt
 
57
 
 
58
  >>> f = file.File('../test/data/1.2_3.las',mode='r')
 
59
  >>> dd_header = liblas.header.Header(handle=f.header.handle,copy=True)
 
60
  >>> utm_header = liblas.header.Header(handle=f.header.handle,copy=True)
 
61
  >>> del f
 
62
 
 
63
  >>> utm_header.offset
 
64
  [0.0, 0.0, 0.0]
 
65
  >>> utm_header.scale
 
66
  [0.01, 0.01, 0.01]
 
67
  
 
68
#  >>> new_offset(0.01, 0.0001, 0.0, 470692.44)
 
69
#  >>> utm_header.offset = [offset+1.0/0.000001 for offset in utm_header.offset]
 
70
#  >>> utm_header.offset
 
71
  [1000000.0, 1000000.0, 1000000.0]
 
72
#  >>> utm_header.scale = [0.000001,0.000001,0.000001]
 
73
  >>> utm_header.srs = s_utm
 
74
  
 
75
#  >>> dd_header.scale = [0.000001,0.000001,0.000001]
 
76
  >>> dd_header.srs = s_dd
 
77
  
 
78
  >>> f = file.File('../test/data/1.2_3.las',mode='r', header = utm_header)
 
79
  >>> f.header.data_offset
 
80
  438L
 
81
  >>> f.header.scale
 
82
  [0.01, 0.01, 0.01]
 
83
  >>> p = f.read(0)
 
84
  >>> origx, origy = p.x, p.y
 
85
  >>> origx, origy
 
86
  (470692.44, 4602888.9000000004)
 
87
  >>> f.set_srs(s_dd)
 
88
  True
 
89
  >>> p = f.read(0)
 
90
 
 
91
We only get truncated values because our header scale 
 
92
values are 0.01
 
93
 
 
94
  >>> p.x, p.y
 
95
  (-93.350000000000009, 41.579999999999998)
 
96
 
 
97
#real values
 
98
#  (-93.351562590199833, 41.577148395415108)
 
99
  
 
100
  
 
101
 
 
102
  >>> f_project = file.File('junk_srs_project.las',mode='w',header=dd_header)
 
103
 
 
104
  >>> p.header = dd_header
 
105
  >>> p.x, p.y
 
106
  (-93.350000000000009, 41.579999999999998)
 
107
 
 
108
  >>> dd_header.srs.proj4
 
109
  '+proj=longlat +datum=WGS84 +no_defs '
 
110
  >>> f_project.write(p);f_project.write(p);f_project.write(p)
 
111
  >>> f_project.close()
 
112
  >>> del f_project
 
113
  >>> f3 = file.File('junk_srs_project.las')
 
114
  >>> f3.header.data_offset
 
115
  789L
 
116
 
 
117
  >>> s_utm = srs.SRS()
 
118
  >>> s_utm.wkt = utm_wkt
 
119
  >>> p3 = f3.read(1)
 
120
  >>> int(round(p3.x)), int(round(p3.y))
 
121
  (-93, 42)
 
122
  >>> p3 = f3.read(0)
 
123
  >>> int(round(p3.x)), int(round(p3.y))
 
124
  (-93, 42)
 
125
 
 
126
  >>> import os
 
127
  >>> os.remove('junk_srs_project.las')
 
128
 
 
129
 
 
130
  >>> f = file.File('../test/data/srs_vertcs.las',mode='r')
 
131
  >>> s = f.header.srs
 
132
  >>> s.get_wkt_compoundok() == """COMPD_CS["unknown",PROJCS["WGS 84 / UTM zone 17N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-81],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","32617"]],VERT_CS["NAVD88 height",VERT_DATUM["North American Vertical Datum 1988",2005,AUTHORITY["EPSG","5103"],EXTENSION["PROJ4_GRIDS","g2003conus.gtx,g2003alaska.gtx,g2003h01.gtx,g2003p01.gtx"]],AXIS["Up",UP],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","5703"]]]"""
 
133
  True
 
134
 
 
135
  >>> s2 = srs.SRS()
 
136
  >>> s2.wkt = """PROJCS["WGS 84 / UTM zone 17N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-81],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","32617"]]"""
 
137
  >>> s2.set_verticalcs( 5703, 'abc', 5103, 9001 )
 
138
  True
 
139
  >>> s2.get_wkt_compoundok()
 
140
  'COMPD_CS["unknown",PROJCS["WGS 84 / UTM zone 17N",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",-81],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","32617"]],VERT_CS["NAVD88 height",VERT_DATUM["North American Vertical Datum 1988",2005,AUTHORITY["EPSG","5103"],EXTENSION["PROJ4_GRIDS","g2003conus.gtx,g2003alaska.gtx,g2003h01.gtx,g2003p01.gtx"]],AXIS["Up",UP],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","5703"]]]'