Finding SRID from WKT out of GRIB file?

Question

I have a grib2 file which I'm trying to import into PostGIS 2.0 Raster and process. Trouble is, I can't figure out what the SRID is supposed to be, or one close enough for government work. Running gdalinfo on the file gives me the output below. I've tried 3785 just as a shot in the dark, but when I transform the coordinates back to 4326 they don't match up. I really would like to work with 4326 in my code, especially since that's how people will be signing up for subscriptions and the like.

If I do something like

select st_value(rast, st_makepoint(-277757.789, 4547961.333)) from test_grib_table;

I get a result. (Note that that's the center point specified below.) But I can't figure out how to transform 4326 lat/long degree coords into this (unknown to me) SRS.

Gdalinfo output:

Driver: GRIB/GRIdded Binary (.grb)
Files: 130-2013.050614_000_TMPK.grib2
Size is 451, 337
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["Coordinate System imported from GRIB file",
        DATUM["unknown",
            SPHEROID["Sphere",6371229,0]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Conformal_Conic_2SP"],
    PARAMETER["standard_parallel_1",25],
    PARAMETER["standard_parallel_2",25],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",265],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0]]
Origin = (-3332155.288903323933482,6830293.833488883450627)
Pixel Size = (13545.000000000000000,-13545.000000000000000)
Corner Coordinates:
Upper Left  (-3332155.289, 6830293.833) (139d51'22.04"W, 54d10'20.71"N)
Lower Left  (-3332155.289, 2265628.833) (126d 6'34.06"W, 16d 9'49.48"N)
Upper Right ( 2776639.711, 6830293.833) ( 57d12'21.76"W, 55d27'10.73"N)
Lower Right ( 2776639.711, 2265628.833) ( 68d56'16.73"W, 17d11'55.33"N)
Center      ( -277757.789, 4547961.333) ( 98d 8'30.73"W, 39d54'5.40"N)
Band 1 Block=451x1 Type=Float64, ColorInterp=Undefined
  Description = 2[m] HTGL="Specified height level above ground"
  Metadata:
    GRIB_UNIT=[K]
    GRIB_COMMENT=Temperature [K]
    GRIB_ELEMENT=TMP
    GRIB_SHORT_NAME=2-HTGL
    GRIB_REF_TIME=  1367848800 sec UTC
    GRIB_VALID_TIME=  1367848800 sec UTC
    GRIB_FORECAST_SECONDS=0 sec

Answer 1

If there is an EPSG code for this projection, you should find it on spatialreference.org. Search for Lambert_Conformal_Conic_2SP.

Otherwise you'll have to define a custom Proj.4 projection string and add it to the spatial_ref_sys table in your PostGIS database. The Proj.4 string would probably look like the following (not tested!):

+proj=lcc +lat_1=25 +lat_2=25 +lat_0=0 +lon_0=-95 +x_0=0 +y_0=0 +ellps=sphere +a=6371229 +b=6371229

Still the reprojection won't work unless the sphere's origin matches the origin of the WGS84 ellipsoid (used in EPSG:4326) by chance. The translation vector can be specified in the Proj.4 string with the +towgs parameter.

Answer 2

If it says it's Lambert Conic Conformal 2SP then its EPSG Code should be 9802. You can check for remotesensing.org for more details.

---

Update

The documents from remotesensing.org have been moved to http://geotiff.maptools.org/proj_list/lambert_conic_conformal_2sp.html