Changeset 14481

Show
Ignore:
Timestamp:
05/16/08 17:24:40 (2 months ago)
Author:
rouault
Message:

Add test for GetGeoTransform? and GetMetadata? on dataset level

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved
  • trunk/autotest/gdrivers/lcp.py

    r14480 r14481  
    4848        gdaltest.post_reason('wrong number of bands') 
    4949        return 'fail' 
     50 
     51    metadata = [ ('LATITUDE', '49'), 
     52                 ('LINEAR_UNIT', 'Meters'), 
     53                 ('DESCRIPTION', 'This is a test LCP file created with FARSITE 4.1.054, using data downloaded from the USGS \r\nNational Map for LANDFIRE (2008-05-06). Data were reprojected to UTM zone 12 on NAD83 \r\nusing gdalwarp (GDAL 1.4.2).\r\n') ] 
     54    md = ds.GetMetadata() 
     55    for item in metadata: 
     56        if md[item[0]] != item[1]: 
     57            gdaltest.post_reason('wrong metadataitem for dataset. md[\'%s\']=\'%s\', expected \'%s\'' % (item[0], md[item[0]], item[1])) 
     58            return 'fail' 
     59 
     60    check_gt = (285807.932887174887583,30,0,5379230.386217921040952,0,-30) 
     61    new_gt = ds.GetGeoTransform() 
     62    for i in range(6): 
     63        if abs(new_gt[i]-check_gt[i]) > 1e-5: 
     64            print 
     65            print 'old = ', check_gt 
     66            print 'new = ', new_gt 
     67            gdaltest.post_reason( 'Geotransform differs.' ) 
     68            return 'fail' 
    5069 
    5170    dataPerBand = [ ( 18645, [  ('ELEVATION_UNIT', '0' ), 
     
    125144        return 'fail' 
    126145 
     146    metadata = [ ('LATITUDE', '48'), 
     147                 ('LINEAR_UNIT', 'Meters'), 
     148                 ('DESCRIPTION', '') ] 
     149    md = ds.GetMetadata() 
     150    for item in metadata: 
     151        if md[item[0]] != item[1]: 
     152            gdaltest.post_reason('wrong metadataitem for dataset. md[\'%s\']=\'%s\', expected \'%s\'' % (item[0], md[item[0]], item[1])) 
     153            return 'fail' 
     154 
     155    check_gt = (-1328145,30,0,2961735,0,-30) 
     156    new_gt = ds.GetGeoTransform() 
     157    for i in range(6): 
     158        if abs(new_gt[i]-check_gt[i]) > 1e-5: 
     159            print 
     160            print 'old = ', check_gt 
     161            print 'new = ', new_gt 
     162            gdaltest.post_reason( 'Geotransform differs.' ) 
     163            return 'fail' 
     164 
    127165    dataPerBand = [ ( 28381, [  ('ELEVATION_UNIT', '0' ), 
    128166                                ('ELEVATION_UNIT_NAME', 'Meters' ),