source: cows/trunk/cows/service/imps/csmlbackend/csmltests/test_csmllayer.py @ 4270

Subversion URL: http://proj.badc.rl.ac.uk/svn/ndg/cows/trunk/cows/service/imps/csmlbackend/csmltests/test_csmllayer.py@4270
Revision 4270, 2.4 KB checked in by domlowe, 11 years ago (diff)

moving wms_csmllayer code

Line 
1import wms_csmllayer
2import csml
3
4
5def test_layerMapper():
6    lm=wms_csmllayer.CSMLLayerMapper()
7    print 'building layermap'
8    file='hadgem1'
9    layers=lm.map(fileoruri=file)
10   
11    #for key, layer in layers.items():
12        #print '*********************'
13        #print 'name: %s'%key  #name is currently gml:id..
14        #print 'title: %s'%layer.title
15        #print 'abstract: %s'%layer.abstract
16        #print 'dimensions: %s'%layer.dimensions
17        #for dim in layer.dimensions:
18            #print layer.dimensions[dim].units
19            #print layer.dimensions[dim].extent
20        #print 'units: %s'%layer.units
21        #print 'crss: %s'%layer.crss
22        #print '*********************'
23       
24    #should be 35 features in this file
25    layerid='UkndO02d'
26    assert len(layers) == 35
27    assert layers[layerid].title == 'm1s3i223'
28    assert layers[layerid].abstract=='SURF & BL TOTL MOISTURE FLUX KG/M2/S'
29    assert layers[layerid].units == ['degrees_east', 'degrees_north', 'none', 'days_since_1970-01-01_00:00:0']
30    assert layers[layerid].crss == ['EPSG:4326']
31    assert layers[layerid].dimensions.keys() == ['latitude', 'z0_surface', 'longitude', 'time']
32    #layerid='EAY38o0Z'   #total precip
33    layerid='DiXD05pN' #solar_1
34    print 'getting slab...'
35    slab= layers[layerid].getSlab('EPSG:4324', dimValues={'time':'1999-12-06T12:00:00.0', 'z0_surface':'0'})
36    img=slab.getImage((-80, -160, 80, 160) ,320, 320)
37    img.save("testimage.jpeg", "JPEG")
38    print 'getting slab2...'
39    slab2= layers[layerid].getSlab('EPSG:4324', dimValues={'time':'1998-12-06T12:00:00.0', 'z0_surface':'0'})
40    img2=slab2.getImage((-80, -160, 80, 160) ,320, 320)
41    img2.save("testimage2.jpeg", "JPEG")
42   
43    #Temperature:
44    print 'getting slab3...'
45    layerid = 'RJ1QcBwz'
46    slab3= layers[layerid].getSlab('EPSG:4324', dimValues={'time':'1995-12-06T12:00:00.0', 'height':'0'})
47    img3=slab3.getImage((-80, -160, 80, 160) ,320, 320)
48    img3.save("testimage3.jpeg", "JPEG")
49   
50    print 'getting lots of images'
51    miny=-90
52    minx=-180
53    maxy=90
54    maxx=180
55    for i in range(8):
56        miny=miny+10
57        minx=minx+10
58        maxy=maxy -10
59        maxx=maxx-10
60        img=slab3.getImage((miny, minx, maxy, maxy) ,320, 320)
61        img.save("sequence%s.jpeg"%str(i), "JPEG")
62   
63   
64   
65   
66   
67def main():
68    test_layerMapper()
69   
70if __name__=="__main__":
71    main()
Note: See TracBrowser for help on using the repository browser.