Add query option for ?f=coveragejson

[C-Loftus]

Describe the solution you'd like

In EDR some endpoints return either geojson or coveragejson.

• For instance, /locations returns geojson but /locations/{id} returns coveragejson.

A user may want to get a different formatting of json output other than the default. However, it is my understanding that pygeoapi processes ?f= before it reaches the provider so it is not possible to switch on the value of the format and return a user-specified format.

• We can get both json and jsonld just since it is hard coded as another option, but this is not possible for coveragejson

Describe alternatives you've considered

• pass the format query param (f=) as an argument for the provider that the provider can optionally use to change the type of json
• hard code covjson similar to json vs jsonld as another json option

Additional context

Similar to how having the option to output raw json is useful to debug jsonld, having the output to output geojson alongside covjson as another format is also useful for debugging.

SInce coveragejson is less common and more complicated, it is useful for a reference to be able to query both the geojson and coveragejson representation of the same resource and the differences in expressiveness.

cc @webb-ben

[C-Loftus]

Believe this is related to #1821

[github-actions]

This Issue has been inactive for 90 days. As per RFC4, in order to manage maintenance burden, it will be automatically closed in 7 days.