Help with mapfile for NetCDF-lambert conic conformal

Jeff Bullard jeffbullard at VERIZON.NET
Fri Dec 1 20:32:27 EST 2006


I am able to get Geotiff rasters to diplay with Mapserver but I'm
struggling with netCDF CF-1 compliant files.

If you have time to look at my mapfile (astro5.map), you can see my html
index at 
http://bigbend.homeip.net/

My netCDF files have the Lambert Conic Conformal projection.
gdal_translate won't handle them and gdalinfo gives a warning that the
lat/long was not found.  This is strange because Ferret will handle them
and ncdump shows lat/long variables with the correct units
"degrees_east" and "degrees_north".  Here is the output:

$gdalinfo 200612010300_sky.nc

Warning 1: xdim/x/lon/longitude or ydim/y/lat/latitude variable(s) not
found!
Driver: netCDF/network Common Data Format
Size is 689, 1
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            TOWGS84[0,0,0,0,0,0,0],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9108"]],
        AXIS["Lat",NORTH],
        AXIS["Long",EAST],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Lambert_Conformal_Conic_1SP"],
    PARAMETER["latitude_of_origin",25],
    PARAMETER["central_meridian",-95],
    PARAMETER["scale_factor",25],
    PARAMETER["false_easting",0],
    PARAMETER["false_northing",0]]
Origin = (-263789.468708,3.000000)
Pixel Size = (5072.03385776,0.00000000)
Metadata:
  NC_GLOBAL#title=GRIB Data translated to NetCDF
  NC_GLOBAL#Conventions=CF-1.0
  NC_GLOBAL#CreatedBy=degrib
  NC_GLOBAL#comment=Used degrib with: -NetCDF 3
  NC_GLOBAL#references=http://www.nws.noaa.gov/mdl/NDFD_GRIB2Decoder/
  Sky_SFC#long_name=Total cloud cover
  Sky_SFC#units=%
  Sky_SFC#coordinates=longitude latitude
  Sky_SFC#ReferenceTime=1164931200
  Sky_SFC#ReferenceTimeString=2006-12-01 00:00:00 00:00
  Sky_SFC#gridType=SCALAR
  Sky_SFC#level=SFC
  Sky_SFC#DescriptiveLevel=0[-] SFC="Ground or water surface"
  Sky_SFC#OriginatingCenter=865535
  Sky_SFC#precision=0
  Sky_SFC#valid_min=0.000000e+00
  Sky_SFC#valid_max=1.000000e+02
  Sky_SFC#NumMissing=375560
  Sky_SFC#_FillValue=9.999000e+03
  Sky_SFC#missing_value=9.999000e+03
  Sky_SFC#grid_mapping=MapProjection
  Sky_SFC#GRIBMessageVersion=2
  Sky_SFC#GRIBMessageUnits=%
  Sky_SFC#GRIB2_Discipline=0
  Sky_SFC#GRIB2_Sect4TemplateNumber=0
  Sky_SFC#GRIB2_Sect4ParameterCategory=6
  Sky_SFC#GRIB2_Sect4ParameterNumber=1
  Sky_SFC#GRIB2_GenerationProcess=2
  Sky_SFC#GRIB2_FstFixedSurfaceType=1
  Sky_SFC#GRIB2_FstFixedSurfaceValue=0
  Sky_SFC#GRIB2_SndFixedSurfaceType=255
  Sky_SFC#GRIB2_SndFixedSurfaceValue=0
  ProjectionHr#long_name=Projection Hour After ReferenceTime
  ProjectionHr#units=hours since 2006-12-01 00:00:00 00:00
  ProjectionHr#_CoordinateAxisType=Time
  ProjectionHr#FirstValidTime=1164942000
  ProjectionHr#FirstValidTimeString=2006-12-01 03:00:00 00:00
  ProjectionHr#FirstProjectionHR=3
  YCells#units=m
  YCells#long_name=y coordinate of projection
  YCells#standard_name=projection_y_coordinate
  YCells#grid_spacing=5079.406000 m
  YCells#TrueGridLength_Latitude=25
  YCells#_CoordinateAxisType=GeoY
  XCells#units=m
  XCells#long_name=x coordinate of projection
  XCells#standard_name=projection_x_coordinate
  XCells#grid_spacing=5079.406000 m
  XCells#TrueGridLength_Latitude=25
  XCells#_CoordinateAxisType=GeoX
  MapProjection#long_name=MapProjection
  MapProjection#grid_mapping_name=lambert_conformal_conic
  MapProjection#EarthShape=sphere
  MapProjection#EarthRadius_KM=6371.2
  MapProjection#EarthRadius=6371.200000 km
  MapProjection#standard_parallel=25
  MapProjection#longitude_of_central_meridian=-95
  MapProjection#latitude_of_projection_origin=25
  MapProjection#LowerLeftLatitude=20.192
  MapProjection#LowerLeftLongitude=-121.554
Corner Coordinates:
Upper Left  ( -263789.469,       3.000) ( 95d 6'16.28"W, 24d59'59.87"N)
Lower Left  ( -263789.469,       3.000) ( 95d 6'16.28"W, 24d59'59.87"N)
Upper Right ( 3230841.859,       3.000) ( 93d43'11.51"W, 24d59'40.17"N)
Lower Right ( 3230841.859,       3.000) ( 93d43'11.51"W, 24d59'40.17"N)
Center      ( 1483526.195,       3.000) ( 94d24'43.84"W, 24d59'55.82"N)
Band 1 Block=689x1 Type=Float32, ColorInterp=Undefined
  NoData Value=9999
Metadata:
  NETCDF_VARNAME=Sky_SFC
  NETCDF_DIMENSION_xcells=1
Band 2 Block=689x1 Type=Float32, ColorInterp=Undefined
  NoData Value=9999
Metadata:
  NETCDF_VARNAME=Sky_SFC
  NETCDF_DIMENSION_xcells=2
Band 3 Block=689x1 Type=Float32, ColorInterp=Undefined
  NoData Value=9999
Metadata:
  NETCDF_VARNAME=Sky_SFC
  NETCDF_DIMENSION_xcells=3
Band 4 Block=689x1 Type=Float32, ColorInterp=Undefined
  NoData Value=9999
Metadata:
  NETCDF_VARNAME=Sky_SFC



More information about the mapserver-users mailing list