id,summary,reporter,owner,description,type,status,priority,milestone,component,version,severity,resolution,keywords,cc 5191,PROJ.4 parameters for LCC (2SP) misinterpreted,tomyun,Kyle Shannon,"A PROJ.4 syntax for '''Lambert Conic Conformal with 2 standard parallels''' seems to be incorrectly parsed as '''1 standard parallel''' parameters. {{{ srs = osr.SpatialReference() srs.ImportFromProj4('+proj=lcc +lat_1=30 +lat_2=60 +lat_0=30 +lon_0=126 +x_0=0 +y_0=0 +a=6378160 +b=6356775 +units=m +no_defs') srs.ExportToProj4() }}} At this point, `srs` shows different parameters from what it has originally imported. {{{ '+proj=lcc +lat_1=30 +lat_0=30 +lon_0=126 +k_0=1 +x_0=0 +y_0=0 +a=6378160 +b=6356775 +units=m +no_defs ' }}} However, if we use a bit verbose WKT syntax for the same projection, {{{ srs.ImportFromWkt(''' PROJCS[""unnamed"", GEOGCS[""Coordinate System imported from GRIB file"", DATUM[""unknown"", SPHEROID[""Spheroid imported from GRIB file"",6378160,298.2539162964695]], PRIMEM[""Greenwich"",0], UNIT[""degree"",0.0174532925199433]], PROJECTION[""Lambert_Conformal_Conic_2SP""], PARAMETER[""standard_parallel_1"",30], PARAMETER[""standard_parallel_2"",60], PARAMETER[""latitude_of_origin"",30], PARAMETER[""central_meridian"",126], PARAMETER[""false_easting"",0], PARAMETER[""false_northing"",0]] ''') srs.ExportToProj4() }}} it now matches with the intended parameters of LCC (2SP). {{{ '+proj=lcc +lat_1=30 +lat_2=60 +lat_0=30 +lon_0=126 +x_0=0 +y_0=0 +a=6378160 +b=6356775 +units=m +no_defs ' }}} ",defect,closed,normal,1.11.0,OGR_SRS,1.10.0,normal,fixed,,