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