| | 117 | ############################################################################### |
|---|
| | 118 | # Verify that prime meridians are preserved when round tripping. (#1940) |
|---|
| | 119 | # |
|---|
| | 120 | |
|---|
| | 121 | def 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 | |
|---|