[netcdf-java] no gridded fields in netcdf file which is generated by GDAL

• To: netcdf-java@xxxxxxxxxxxxxxxx
• Subject: [netcdf-java] no gridded fields in netcdf file which is generated by GDAL
• From: Hendy Sumanto <hsumanto@xxxxxxxx>
• Date: Fri, 26 Aug 2011 11:42:14 +1000 (EST)

Hi Everyone,

I am using gdal_translate to generate a netcdf from a USGS LANDSAT geotiff 
data, but I'm not seeing why the grid 
described below is not recognized as a grid by NetCDF Tools (4.2).  Could 
someone please point out the problem?

Thanks,

Hendy

netcdf L71092086_08620110129_B50-variant3 {
dimensions:
        x = 8121 ;
        y = 7121 ;
variables:
        char transverse_mercator ;
                transverse_mercator:Northernmost_Northing = -4040685. ;
                transverse_mercator:Southernmost_Northing = -4254315. ;
                transverse_mercator:Easternmost_Easting = 525015. ;
                transverse_mercator:Westernmost_Easting = 281385. ;
                transverse_mercator:spatial_ref = "PROJCS[\"WGS 84 / UTM zone 
55N\",GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 
84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433],AUTHORITY[\"EPSG\",\"4326\"]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",0],PARAMETER[\"central_meridian\",147],PARAMETER[\"scale_factor\",0.9996],PARAMETER[\"false_easting\",500000],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AUTHORITY[\"EPSG\",\"32655\"]]"
 ;
                transverse_mercator:GeoTransform = "281384.9999999999 30 0 
-4040685 0 -30 " ;
                transverse_mercator:grid_mapping_name = "transverse_mercator" ;
                transverse_mercator:latitude_of_projection_origin = 0.f ;
                transverse_mercator:longitude_of_central_meridian = 147.f ;
                transverse_mercator:scale_factor_at_central_meridian = 0.9996f ;
                transverse_mercator:false_easting = 500000.f ;
                transverse_mercator:false_northing = 0.f ;
        byte Band1(y, x) ;
                Band1:grid_mapping = "transverse_mercator" ;
                Band1:long_name = "GDAL Band Number 1" ;

// global attributes:
                :Conventions = "CF-1.0" ;
}

• Follow-Ups:
  • Re: [netcdf-java] no gridded fields in netcdf file which is generated by GDAL
    • From: John Caron