| | 366 | ############################################################################### |
|---|
| | 367 | # Test automatic morphing of ESRI-style LCC WKT prefixed with 'ESRI::' |
|---|
| | 368 | |
|---|
| | 369 | def osr_esri_12(): |
|---|
| | 370 | |
|---|
| | 371 | srs = osr.SpatialReference() |
|---|
| | 372 | srs.SetFromUserInput( 'ESRI::PROJCS["Lambert Conformal Conic",GEOGCS["grs80",DATUM["D_North_American_1983",SPHEROID["Geodetic_Reference_System_1980",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["standard_parallel_1",34.33333333333334],PARAMETER["standard_parallel_2",36.16666666666666],PARAMETER["latitude_of_origin",33.75],PARAMETER["central_meridian",-79],PARAMETER["false_easting",609601.22],PARAMETER["false_northing",0],UNIT["Meter",1]]' ) |
|---|
| | 373 | |
|---|
| | 374 | # No MorphFromESRI() is required |
|---|
| | 375 | |
|---|
| | 376 | if srs.GetAttrValue( 'PROJECTION' ) != 'Lambert_Conformal_Conic_2SP': |
|---|
| | 377 | gdaltest.post_reason( 'Got wrong PROJECTION name (%s) after ESRI morph.' % \ |
|---|
| | 378 | srs.GetAttrValue( 'PROJECTION' ) ) |
|---|
| | 379 | return 'fail' |
|---|
| | 380 | |
|---|
| | 381 | if abs( srs.GetProjParm('standard_parallel_1') - 34.333333333 ) > 0.00001: |
|---|
| | 382 | gdaltest.post_reason( 'Got wrong parameter value (%g) after ESRI morph.' % \ |
|---|
| | 383 | srs.GetProjParm('standard_parallel_1') ) |
|---|
| | 384 | return 'fail' |
|---|
| | 385 | |
|---|
| | 386 | if srs.GetAttrValue( 'DATUM' ) != 'North_American_Datum_1983': |
|---|
| | 387 | gdaltest.post_reason( 'Got wrong DATUM name (%s) after ESRI morph.' % \ |
|---|
| | 388 | srs.GetAttrValue( 'DATUM' ) ) |
|---|
| | 389 | |
|---|
| | 390 | if srs.GetAttrValue( 'UNIT' ) != 'Meter': |
|---|
| | 391 | gdaltest.post_reason( 'Got wrong UNIT name (%s) after ESRI morph.' % \ |
|---|
| | 392 | srs.GetAttrValue( 'UNIT' ) ) |
|---|
| | 393 | return 'fail' |
|---|
| | 394 | |
|---|
| | 395 | return 'success' |
|---|
| | 396 | |
|---|
| | 397 | ############################################################################### |
|---|
| | 398 | # Test automatic morphing of ESRI-style LCC WKT prefixed with 'ESRI::' |
|---|
| | 399 | # but read directly from file. |
|---|
| | 400 | |
|---|
| | 401 | def osr_esri_13(): |
|---|
| | 402 | |
|---|
| | 403 | srs = osr.SpatialReference() |
|---|
| | 404 | srs.SetFromUserInput( 'data/lcc_esri.prj' ) |
|---|
| | 405 | |
|---|
| | 406 | # No MorphFromESRI() is required |
|---|
| | 407 | |
|---|
| | 408 | if srs.GetAttrValue( 'PROJECTION' ) != 'Lambert_Conformal_Conic_2SP': |
|---|
| | 409 | gdaltest.post_reason( 'Got wrong PROJECTION name (%s) after ESRI morph.' % \ |
|---|
| | 410 | srs.GetAttrValue( 'PROJECTION' ) ) |
|---|
| | 411 | return 'fail' |
|---|
| | 412 | |
|---|
| | 413 | if abs( srs.GetProjParm('standard_parallel_1') - 34.333333333 ) > 0.00001: |
|---|
| | 414 | gdaltest.post_reason( 'Got wrong parameter value (%g) after ESRI morph.' % \ |
|---|
| | 415 | srs.GetProjParm('standard_parallel_1') ) |
|---|
| | 416 | return 'fail' |
|---|
| | 417 | |
|---|
| | 418 | if srs.GetAttrValue( 'DATUM' ) != 'North_American_Datum_1983': |
|---|
| | 419 | gdaltest.post_reason( 'Got wrong DATUM name (%s) after ESRI morph.' % \ |
|---|
| | 420 | srs.GetAttrValue( 'DATUM' ) ) |
|---|
| | 421 | |
|---|
| | 422 | if srs.GetAttrValue( 'UNIT' ) != 'Meter': |
|---|
| | 423 | gdaltest.post_reason( 'Got wrong UNIT name (%s) after ESRI morph.' % \ |
|---|
| | 424 | srs.GetAttrValue( 'UNIT' ) ) |
|---|
| | 425 | return 'fail' |
|---|
| | 426 | |
|---|
| | 427 | return 'success' |
|---|
| | 428 | |
|---|
| | 429 | |
|---|
| | 430 | ############################################################################### |
|---|
| | 431 | |
|---|