NetCDF in Mapserver...

[Norman Barker]

-----Original Message-----
From: UMN MapServer Users List [mailto:MAPSERVER-USERS at LISTS.UMN.EDU]On
Behalf Of Andrea Marelli
Sent: Wednesday, November 16, 2005 12:02 PM
To: MAPSERVER-USERS at LISTS.UMN.EDU
Subject: [UMN_MAPSERVER-USERS] NetCDF in Mapserver...


Hi all,

I can't see anything (only a white image) when I try to use mapserver with
NetCDF files.
I need some suggestions...

The server has installed Mapservser 4.6.1 (with support for
Proj/Geos/Gdal/Postgis/Wfs/Wcs), Gdal 1.3.1 (support for
Postgres/Geos/Proj4/Netcd) and NetCDF 3.6.0-p1, ...

The NetCDF file's Header is:


netcdf  N20001832000213_MO {
dimensions:
        lon = 1480 ;
        lat = 1207 ;
variables:
        float lon(lon) ;
                lon:units = "Degrees" ;
                lon:long_name = "Longitude" ;
        float lat(lat) ;
                lat:units = "Degrees" ;
                lat:long_name = "Latitude" ;
        float chlor_a(lat, lon) ;
                chlor_a:units = "mg.m-3" ;
                chlor_a:long_name = "Chlorophyll a concentration, OC4v4 Algorith
m" ;
                chlor_a:missing_value = -9999.f ;
        float K_490(lat, lon) ;
                K_490:units = "m-1" ;
                K_490:long_name = "Diffuse attenuation coefficient at 490 nm" ;
                K_490:missing_value = -9999.f ;
        float nLw_412(lat, lon) ;
                nLw_412:units = "mW.cm-2.sr-1.um-1" ;
                nLw_412:long_name = "Normalized water-leaving radiance at 412 nm
" ;
                nLw_412:missing_value = -9999.f ;
        float nLw_443(lat, lon) ;
                nLw_443:units = "mW.cm-2.sr-1.um-1" ;
                nLw_443:long_name = "Normalized water-leaving radiance at 443 nm
" ;
                nLw_443:missing_value = -9999.f ;
        float nLw_490(lat, lon) ;
                nLw_490:units = "mW.cm-2.sr-1.um-1" ;
                nLw_490:long_name = "Normalized water-leaving radiance at 490 nm
" ;
                nLw_490:missing_value = -9999.f ;
        float nLw_510(lat, lon) ;
                nLw_510:units = "mW.cm-2.sr-1.um-1" ;
                nLw_510:long_name = "Normalized water-leaving radiance at 510 nm
" ;
                nLw_510:missing_value = -9999.f ;
        float nLw_555(lat, lon) ;
                nLw_555:units = "mW.cm-2.sr-1.um-1" ;
                nLw_555:long_name = "Normalized water-leaving radiance at 555 nm
" ;
                nLw_555:missing_value = -9999.f ;
        float nLw_670(lat, lon) ;
                nLw_670:units = "mW.cm-2.sr-1.um-1" ;
                nLw_670:long_name = "Normalized water-leaving radiance at 670 nm
" ;
                nLw_670:missing_value = -9999.f ;
        float l3m_data(lat, lon) ;
                l3m_data:units = "Einstein m^-2 Day" ;
                l3m_data:long_name = "Photosynthetically Available Radiation" ;
                l3m_data:missing_value = -9999.f ;
        float sst(lat, lon) ;
                sst:units = "Degrees" ;
                sst:long_name = "Temperature" ;
                sst:missing_value = -9999.f ;

// global attributes:
                :Year = 2000s ;
                :Month = 7s ;
                :Southernmost_Latitude = 22.f ;
                :Easternmost_Latitude = -20.f ;
                :Northernmost_Latitude = 75.f ;
                :Westernmost_Latitude = 45.f ;
                :Number_of_Columns = 1480s ;
                :Number_of_Rows = 1207s ;
}
#########  end of NetCDF file's Header ########




and this is the mapfile:


MAP

NAME "Test"
EXTENT 22 -20 75 45
SIZE 600 600
SHAPEPATH "/raid/emis/satellite_test"

LAYER
   NAME "chlor_a"
   TYPE raster
   STATUS DEFAULT
   PROCESSING "BANDS=3"
   DATA "N20001832000213_MO.ncdf"
   CLASS
   NAME "chlor_a"
     STYLE
       COLOR 255 0 0
     END
   END
END

END
#########  end of Map file ########


So, why I can't see anything?
Is there anything wrong?

Please Help me..
Thanks!!!

Andrea

Hi Andrea,

are you trying to serve netCDF using a WMS, or a WCS?  If it is a WCS then there is a tutorial here http://ms.gis.umn.edu:8081/ms_plone/docs/howto/WCSServerFormatHowTo and I would be grateful for any feedback on how to improve this documentation.

Norman