GDAL GeoTiff -> NetCDF Testing results
======================================

Albers Equal Area: OK
------------------------

*CF-1 compliance?:* YES
*Loads in IDV?*: YES

(Works when we set 'longitude_of_center' -> 'longitude_of_central_meridian')

Test EPSG: 3577 ("GDA94 / Australian Albers")

======================================  =============================================
GDAL:                                   CF-1 Equivalents:
======================================  =============================================
PROJECTION["Albers_Conic_Equal_Area"],  grid_mapping_name = albers_conical_equal_area
PARAMETER["standard_parallel_1",-18],   standard_parallel=-18, -36  
PARAMETER["standard_parallel_2",-36],   
PARAMETER["latitude_of_center",0],      latitude_of_projection_origin=0
PARAMETER["longitude_of_center",132],   longitude_of_central_meridian=132
======================================  =============================================

Azimuthal equidistant: Partial
------------------------------

*CF-1 compliance?:* YES
*Loads in IDV?*: NO (says it doesn't contain gridded data)

Notes:
 - (Problem loading into IDV? At least with params below)
 - (Problem loading into IDV? - complains its 'not gridded data'?)

TODO:- maybe test further in NetCDF Java.

=========================================== ===============================================
GDAL:                                       CF-1 Equivalents
=========================================== ===============================================
PROJECTION["Azimuthal_Equidistant"],        grid_mapping_name=azimuthal_equidistant
PARAMETER["latitude_of_center",0],          latitude_of_projection_origin=0
PARAMETER["longitude_of_center",134],       longitude_of_projection_origin=134
=========================================== ===============================================


Lambert azimuthal equal area: OK
--------------------------------

*CF-1 compliance?:* YES
*Loads in IDV?*: YES (Reprojection very sensitive to input params)

(Works when we set 'longitude_of_center' -> 'longitude_of_projection_origin')

=========================================== ===============================================
GDAL:                                       CF-1 Equivalents
=========================================== ===============================================
PROJECTION["Lambert_Azimuthal_Equal_Area"], grid_mapping_name=lambert_azimuthal_equal_area
PARAMETER["latitude_of_center",-37],        latitude_of_projection_origin=-37
PARAMETER["longitude_of_center",145],       longitude_of_projection_origin=145
=========================================== ===============================================

Lambert conformal: OK
-----------------------

*CF-1 compliance?:* YES
*Loads in IDV?*: YES

Test EPSG: 3112 ("GDA94 / Geoscience Australia Lambert")

==========================================  =============================================
GDAL:                                       CF-1 Equivalents
==========================================  =============================================
PROJECTION["Lambert_Conformal_Conic_2SP"],  grid_mapping_name=lambert_conformal_conic
PARAMETER["standard_parallel_1",-18],       standard_parallel=-18, -36
PARAMETER["standard_parallel_2",-36],
PARAMETER["latitude_of_origin",0],          latitude_of_projection_origin=0
PARAMETER["central_meridian",134],          longitude_of_central_meridian=134
==========================================  =============================================

Lambert Cylindrical Equal Area: Partial
---------------------------------------

*CF-1 compliance?:* YES
*Loads in IDV?*: NO (Claims its not gridded data)

(works - note also has a 'scale_factor_at_projection_origin' alternative to 'standard_parallel')

=========================================== ===============================================
GDAL:                                       CF-1 Equivalents
=========================================== ===============================================
PROJECTION["Cylindrical_Equal_Area"],       
PARAMETER["standard_parallel_1",-37],       standard_parallel=-37
PARAMETER["central_meridian",145],          longitude_of_central_meridian=145
=========================================== ===============================================

Mercator-1SP: No
----------------

*CF-1 compliance?:* YES
*Loads in IDV?*: NO (Claims its not gridded data)

Note:
 - Needs 'central_meridian' -> 'longitude_of_projection_origin'
 - Needs 'scale_factor' -> 'scale_factor_at_projection_origin'

=========================================== ===============================================
GDAL:                                       CF-1 Equivalents
=========================================== ===============================================
PROJECTION["Mercator_1SP"],                 grid_mapping_name = mercator
PARAMETER["central_meridian",150],          longitude_of_projection_origin
PARAMETER["scale_factor",1],                scale_factor_at_projection_origin
=========================================== ===============================================

Mercator-2SP: GDAL Prob
-----------------------

Not added as yet :- seems to be problems as GDAL refuses to create Mercator-2SP files ??

TODO: perhaps file a bug

=========================================== ===============================================
GDAL:                                       CF-1 Equivalents
=========================================== ===============================================
PROJECTION["Mercator_2SP"],                 grid_mapping_name = mercator
PARAMETER["standard_parallel_1",-45],       standard_parallel (? 2 values, to cover lat of origin?)  
PARAMETER["latitude_of_origin",0],          
PARAMETER["central_meridian",150],          longitude_of_projection_origin
=========================================== ===============================================

Orthographic: No
----------------

Need 'central_meridian' -> 'longitude_of_projection_origin'?

=========================================== ===============================================
GDAL:                                       CF-1 Equivalents
=========================================== ===============================================
PROJECTION["Orthographic"],                 grid_mapping_name=orthographic
PARAMETER["latitude_of_origin",-37],        latitude_of_projection_origin=-37
PARAMETER["central_meridian",145],          longitude_of_projection_origin=145
=========================================== ===============================================

Polar stereographic: No
-----------------------

Issues:
 - need "scale_factor" -> 'scale_factor_at_projection_origin'
 - think need "central_meridian" -> 'straight_vertical_longitude_from_pole'
 - not sure if need "latitude_of_origin" -> latitude_of_projection_origin? (see below)

=========================================== ===============================================
GDAL:                                       CF-1 Equivalents
=========================================== ===============================================
PROJECTION["Polar_Stereographic"],          
PARAMETER["latitude_of_origin",0],          (?) - in proj4/WKT, this is extra to +90/-90
PARAMETER["central_meridian",134],          straight_vertical_longitude_from_pole (?)
PARAMETER["scale_factor",1],                scale_factor_at_projection_origin
                                            latitude_of_projection_origin (?)
=========================================== ===============================================

Rotated pole: No GDAL Equivalent?
---------------------------------

Not sure if there's a GDAL Equivalent of this?

Stereographic: No
-----------------

Notes:
 - need "central_meridian" -> 'longitude_of_projection_origin'
 - need 'scale_factor' -> 'scale_factor_at_projection_origin'

=========================================== ===============================================
GDAL:                                       CF-1 Equivalents
=========================================== ===============================================
PROJECTION["Oblique_Stereographic"],        grid_mapping_name = stereographic
PARAMETER["latitude_of_origin",-37],        latitude_of_projection_origin
PARAMETER["central_meridian",134],          longitude_of_projection_origin
PARAMETER["scale_factor",1],                scale_factor_at_projection_origin
=========================================== ===============================================

Transverse Mercator: OK
-----------------------

=========================================== ===============================================
GDAL:                                       CF-1 Equivalents
=========================================== ===============================================
PROJECTION["Transverse_Mercator"],          grid_mapping_name=transverse_mercator
PARAMETER["latitude_of_origin",0],          latitude_of_projection_origin=0
PARAMETER["central_meridian",147],          longitude_of_central_meridian=147
=========================================== ===============================================

