[ncl-talk] hdf to netcdf coords
Dennis Shea
shea at ucar.edu
Thu Feb 14 08:02:14 MST 2019
These numeric types are supported by netCDF-4
%> *ncl_convert2nc* MOD11_L2.A2008009.0435.006.2015337174400.hdf *-nc4c*
On Thu, Feb 14, 2019 at 5:45 AM Vanúcia Schumacher <
vanucia-schumacher at hotmail.com> wrote:
> Hi NCL team,
>
> I would like to convert hdf data to netcdf.
> I tried to use ncl_convert2nc, but it did not work:
>
> Classic mode NetCDf does not support unsigned integer types: LST will be
> written as type short
> Classic mode NetCDf does not support unsigned integer types: QC will be
> written as type short
> Classic mode NetCDf does not support unsigned integer types: Error_LST
> will be written as type byte
> Classic mode NetCDf does not support unsigned integer types: Emis_31 will
> be written as type byte
> Classic mode NetCDf does not support unsigned integer types: Emis_32 will
> be written as type byte
> Classic mode NetCDf does not support unsigned integer types: View_angle
> will be written as type byte
> Classic mode NetCDf does not support unsigned integer types: View_time
> will be written as type byte
>
> I also tried this script and it did not work
>
> fils = ("MOD11_L2.A2008009.0435.006.2015337174400.hdf")
> f = addfile(fils, "r")
> var = short2flt_hdf(f->LST)
> lat2d = f->Latitude
> lon2d = f->Longitude
>
> var!0 = "latitude"
> var!1 = "longitude"
> var&latitude = lat2d(:,0)
> var&longitude = lon2d(0,:)
> var at lat2d = lat2d
> var at lon2d = lon2d
>
> newfile = "of.nc"
> new_var = var
> file_create = addfile(newfile, "c")
> file_create->LST= new_var
>
> fatal:Coordinate variables must have one dimension only
> fatal:No coordinate variable exists for dimension (latitude) in variable
> (var)
>
> print(f):
> identifier_product_doi : 10.5067/MODIS/MOD11_L2.006
> identifier_product_doi_authority : http://dx.doi.org
> dimensions:
> Coarse_swath_lines_5km_MOD_Swath_LST = 410
> Coarse_swath_pixels_5km_MOD_Swath_LST = 271
> Along_swath_lines_1km_MOD_Swath_LST = 2050
> Cross_swath_pixels_1km_MOD_Swath_LST = 1354
> variables:
> float Latitude ( Coarse_swath_lines_5km_MOD_Swath_LST,
> Coarse_swath_pixels_5km_MOD_Swath_LST )
> long_name : Latitude of every 5 scan lines and 5 pixels
> units : degree
> Number_Type : float32
> valid_range : ( -90, 90 )
> _FillValue : -999
> Cell_Along_Swath_Sampling : 3, 2048, 5
> Cell_Across_Swath_Sampling : 3, 1353, 5
> hdf_name : Latitude
>
> float Longitude ( Coarse_swath_lines_5km_MOD_Swath_LST,
> Coarse_swath_pixels_5km_MOD_Swath_LST )
> long_name : Longitude of every 5 scan lines and 5 pixels
> units : degree
> Number_Type : float32
> valid_range : ( -180, 180 )
> _FillValue : -999
> Cell_Along_Swath_Sampling : 3, 2048, 5
> Cell_Across_Swath_Sampling : 3, 1353, 5
> hdf_name : Longitude
>
> ushort LST ( Along_swath_lines_1km_MOD_Swath_LST,
> Cross_swath_pixels_1km_MOD_Swath_LST )
> long_name : Land-surface Temperature
> units : K
> Number_Type : uint16
> valid_range : ( 7500, 65535 )
> _FillValue : 0
> LST : LST data * scale_factor
> scale_factor : 0.02
> add_offset : 0
> scale_factor_err : 0
> add_offset_err : 0
> calibrated_nt : 5
> hdf_name : LST
>
> ushort QC ( Along_swath_lines_1km_MOD_Swath_LST,
> Cross_swath_pixels_1km_MOD_Swath_LST )
> long_name : Quality control for daytime LST and emissivity
> Number_Type : uint16
> valid_range : ( 0, 65535 )
> hdf_name : QC
>
> ubyte Error_LST ( Along_swath_lines_1km_MOD_Swath_LST,
> Cross_swath_pixels_1km_MOD_Swath_LST )
> long_name : Land-surface Temperature Error
> units : K
> Number_Type : uint8
> valid_range : ( 1, 255 )
> _FillValue : 0
> Error_LST : Error_LST data * scale_factor
> scale_factor : 0.04
> add_offset : 0
> scale_factor_err : 0
> add_offset_err : 0
> calibrated_nt : 5
> hdf_name : Error_LST
>
> ubyte Emis_31 ( Along_swath_lines_1km_MOD_Swath_LST,
> Cross_swath_pixels_1km_MOD_Swath_LST )
> long_name : Band 31 emissivity
> Number_Type : uint8
> valid_range : ( 1, 255 )
> _FillValue : 0
> Emis_31 : Emis_31 data * scale_factor + add_offset
> scale_factor : 0.002
> add_offset : 0.49
> scale_factor_err : 0
> add_offset_err : 0
> calibrated_nt : 5
> hdf_name : Emis_31
>
> ubyte Emis_32 ( Along_swath_lines_1km_MOD_Swath_LST,
> Cross_swath_pixels_1km_MOD_Swath_LST )
> long_name : Band 32 emissivity
> Number_Type : uint8
> valid_range : ( 1, 255 )
> _FillValue : 0
> Emis_32 : Emis_32 data * scale_factor + add_offset
> scale_factor : 0.002
> add_offset : 0.49
> scale_factor_err : 0
> add_offset_err : 0
> calibrated_nt : 5
> hdf_name : Emis_32
>
> ubyte View_angle ( Along_swath_lines_1km_MOD_Swath_LST,
> Cross_swath_pixels_1km_MOD_Swath_LST )
> long_name : zenith angle of MODIS viewing at the pixel
> units : deg
> Number_Type : uint8
> valid_range : ( 0, 180 )
> _FillValue : 255
> View_angle : View_angle data * scale_factor
> scale_factor : 0.5
> add_offset : 0
> scale_factor_err : 0
> add_offset_err : 0
> calibrated_nt : 5
> hdf_name : View_angle
>
> ubyte View_time ( Along_swath_lines_1km_MOD_Swath_LST,
> Cross_swath_pixels_1km_MOD_Swath_LST )
> long_name : Time of Land-surface Temperature observation
> units : hrs
> Number_Type : uint8
> valid_range : ( 0, 240 )
> _FillValue : 255
> View_time : View_time data * scale_factor
> scale_factor : 0.1
> add_offset : 0
> scale_factor_err : 0
> add_offset_err : 0
> calibrated_nt : 5
> hdf_name : View_time
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190214/c3563a10/attachment.html>
More information about the ncl-talk
mailing list