rgdal-dev / rasterwise

Hard-won lessons! Don't lose 'em!
4 stars 0 forks source link

netcdf vicgrid coastline elevation #19

Open mdsumner opened 1 year ago

mdsumner commented 1 year ago

this one has huge lon/lat arrays in it but it's just VicGrid

http://imos-data.s3-website-ap-southeast-2.amazonaws.com/?prefix=Deakin_University/bathymetry/

the ncdump is


 ncdump -h Victorian-coast_Bathy_10m.nc
netcdf Victorian-coast_Bathy_10m {
dimensions:
        Y = 42454 ;
        X = 92606 ;
variables:
        double Y(Y) ;
                Y:long_name = "y coordinate of projection" ;
                Y:standard_name = "projection_y_coordinate" ;
                Y:units = "m" ;
                Y:axis = "Y" ;
                Y:valid_min = -40000000. ;
                Y:valid_max = 40000000. ;
                Y:reference_datum = "VICGRID94 projection (EPSG:3111)" ;
        double X(X) ;
                X:long_name = "x coordinate of projection" ;
                X:standard_name = "projection_x_coordinate" ;
                X:units = "m" ;
                X:axis = "X" ;
                X:valid_min = -40100000. ;
                X:valid_max = 40100000. ;
                X:reference_datum = "VICGRID94 projection (EPSG:3111)" ;
        double LATITUDE(Y, X) ;
                LATITUDE:long_name = "latitude" ;
                LATITUDE:standard_name = "latitude" ;
                LATITUDE:units = "degrees_north" ;
                LATITUDE:valid_min = -90. ;
                LATITUDE:valid_max = 90. ;
                LATITUDE:reference_datum = "WGS84 coordinate reference system (EPSG:4326)" ;
        double LONGITUDE(Y, X) ;
                LONGITUDE:long_name = "longitude" ;
                LONGITUDE:standard_name = "longitude" ;
                LONGITUDE:units = "degrees_east" ;
                LONGITUDE:valid_min = -180. ;
                LONGITUDE:valid_max = 180. ;
                LONGITUDE:reference_datum = "WGS84 coordinate reference system (EPSG:4326)" ;
        int CRS ;
                CRS:long_name = "WGS84 lat/lon coordinate reference system" ;
                CRS:grid_mapping_name = "latitude_longitude" ;
                CRS:epsg_code = "EPSG:4326" ;
                CRS:longitude_of_prime_meridian = 0. ;
                CRS:semi_major_axis = 6378137. ;
                CRS:inverse_flattening = 298.257223563 ;
        float HEIGHT(Y, X) ;
                HEIGHT:_FillValue = 99999.f ;
                HEIGHT:long_name = "height" ;
                HEIGHT:standard_name = "height" ;
                HEIGHT:units = "m" ;
                HEIGHT:positive = "up" ;
                HEIGHT:grid_mapping = "CRS" ;
                HEIGHT:coordinates = "LATITUDE LONGITUDE" ;
                HEIGHT:valid_min = -12000.f ;
                HEIGHT:valid_max = 9000.f ;
                HEIGHT:reference_datum = "AHD71 vertical datum (EPSG:5711)" ;

// global attributes:
                :project = "Cooperative Research Centre for Spatial Information (CRCSI)" ;
                :Conventions = "CF-1.6,IMOS-1.4" ;
                :standard_name_vocabulary = "NetCDF Climate and Forecast (CF) Metadata Convention Standard Name Table Version 29" ;
                :title = "Victorian Coastal Digital Elevation Model (VCDEM 2017) 10m resolution." ;
                :institution = "Cooperative Research Centre for Spatial Information (CRCSI)" ;
                :institution_references = "http://www.crcsi.com.au/" ;
                :date_created = "2018-03-08T02:31:20Z" ;
                :abstract = "A gap free Digital Elevation Model (DEM) for the Victorian Coastal region created from the Updated High Resolution 2.5m VCDEM as well as parts of the Shuttle Radar Topography Mission (SRTM) 1 Second DEM, the Australian Bathymetry and Topography grid and a 10m DEM of Port Phillip Bay (PPB) that were provided for this project. This product uses a combination of low resolution data and interpolation techniques to provide a DEM that appears representative of the bathymetry in the Victorian Coastal Region." ;
                :lineage = "The Continuous Seamless 10m DEM is derived from the Updated High Resolution 2.5m VCDEM as well as parts of the SRTM 1 second DEM data and the Australian Bathymetry and Topography grid with all gaps interpolated smoothly to create a comprehensive gap free DEM for the Victorian Coastal Region.\nPositional Accuracy: The final product is derived from various data sources - from DELWP, Port of Melbourne, Deakin University, CRCSI, Royal Australian Navy, and Geoscience Australia.\nThe 2010 VCDEM product consisted of a 1m Topography DEM and 2.5m Bathymetry product that covered 90% of the Updated High Resolution DEM. The 1m Topography DEM is stated to have a horizontal accuracy of ± 35cm and a vertical accuracy of ± 10cm @ 68% confidence interval. The 2.5m Bathymetry DEM accuracy meets the International Hydrographic Organization (IHO) Order 1 specifications of ± 50cm vertical and ± 3.17m horizontal accuracy @ 2σ.\nThe multibeam data obtained from Deakin University was acquired using a Kongsberg EM 2040 MBES with a capability of achieving IHO Order of Special for hydrographic standard. The final accuracy of the data acquired exceeded that required for Order 1b of the IHO S-44 standard for hydrographic surveys.\nMBES and LiDAR data obtained from Port of Melbourne was acquired using systems designed to meet IHO Order 1 for horizontal and vertical accuracy.\nReprojection, aggregation, resampling and slight re-alignment of datasets will have had a negative impact on the final accuracy.\nOriginal Datums were as follows:\nVertical Datum: AHD71 (EPSG:5711)\nHorizontal Datum GDA94 - Projection: VicGrid 94 (EPSG:3111).\n" ;
                :naming_authority = "IMOS" ;
                :geospatial_lat_min = -40.9126005119486 ;
                :geospatial_lat_max = -36.9737584013878 ;
                :geospatial_lat_units = "degrees_north" ;
                :geospatial_lon_min = 139.444242120472 ;
                :geospatial_lon_max = 150.400025934403 ;
                :geospatial_lon_units = "degrees_east" ;
                :geospatial_vertical_min = -118.f ;
                :geospatial_vertical_max = 753.9601f ;
                :geospatial_vertical_positive = "up" ;
                :geospatial_vertical_units = "metres" ;
                :data_centre = "Australian Ocean Data Network (AODN)" ;
                :data_centre_email = "info@aodn.org.au" ;
                :author = "Quadros, Nathan" ;
                :author_email = "nquadros@crcsi.com.au" ;
                :point_of_contact = "Ierodiaconou, Daniel" ;
                :point_of_contact_email = "iero@deakin.edu.au" ;
                :citation = "The citation in a list of references is: \"Cooperative Research Centre for Spatial Information (CRCSI), [year-of-data-download], Victorian Coastal Digital Elevation Model (VCDEM 2017), [data-access-url], accessed [date-of-access].\"" ;
                :acknowledgement = "Any users of CRCSI data are required to clearly acknowledge the source of the material in the format: \"Data was sourced from the Cooperative Research Centre for Spatial Information (CRCSI). The following organisations were involved in the collection of the data: Department of Environment, Land, Water and Planning (DELWP), Victorian Government; Port of Melbourne; Deakin University; Royal Australian Navy (RAN); and Commonwealth of Australia (Geoscience Australia).\"" ;
                :disclaimer = "This data is not suitable for navigational purposes." ;
                :credit = "Department of Environment, Land, Water and Planning (DELWP), Victorian Government; Port of Melbourne; Deakin University; Royal Australian Navy (RAN); Commonwealth of Australia (Geoscience Australia)" ;
                :license = "https://creativecommons.org/licenses/by-nc/4.0/" ;
                :cdm_data_type = "Grid" ;
}