Changeset 12536

Show
Ignore:
Timestamp:
10/25/07 23:13:27 (8 months ago)
Author:
warmerdam
Message:

added prime meridian round trip test (#1940)

Files:

Legend:

Unmodified
Added
Removed
Modified
Copied
Moved
  • trunk/autotest/osr/osr_proj4.py

    r11819 r12536  
    115115    return 'fail' 
    116116 
     117############################################################################### 
     118# Verify that prime meridians are preserved when round tripping. (#1940) 
     119# 
     120 
     121def osr_proj4_5(): 
     122     
     123    srs = osr.SpatialReference() 
     124 
     125    srs.ImportFromProj4( '+proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs' ) 
     126 
     127    if abs(float(srs.GetAttrValue('PRIMEM',1)) - 2.3372291667) > 0.00000001: 
     128        gdaltest.post_reason('prime meridian lost?') 
     129        return 'fail' 
     130 
     131    if abs(srs.GetProjParm('central_meridian')) != 0.0: 
     132        gdaltest.post_reason( 'central meridian altered?' ) 
     133        return 'fail' 
     134 
     135    p4 = srs.ExportToProj4() 
     136    srs2 = osr.SpatialReference() 
     137    srs2.ImportFromProj4( p4 ) 
     138 
     139    if not srs.IsSame(srs2): 
     140        gdaltest.post_reason( 'round trip via PROJ.4 damaged srs?' ) 
     141        print srs.ExportToPrettyWkt() 
     142        print srs2.ExportToPrettyWkt() 
     143     
     144    return 'success' 
     145 
    117146gdaltest_list = [  
    118147    osr_proj4_1, 
     
    120149    osr_proj4_3, 
    121150    osr_proj4_4, 
     151    osr_proj4_5, 
    122152    None ] 
    123153