REFAC: update multi-baseline for large input tests (same projection, but different wkt strings) and removed out-of-date baselines

The application Superimpose crashed with a segmentation fault during
resampling. Investigation revealed that the interpolated index was out
of image bounds. The application uses the PHR mode and resample the XS
image with a otb::GridResampleImageFilter. This filter doesn't check
explicitly that input indexes are in the input buffered region. It
uses a m_ReachableOutputRegion to crop the output region
processed. The problem appears because an output pixel is right on the
border of the input image extent, and when someone uses a nearest neighbor
interpolator, this one converts the continuous index to an
out-of-bound index.

ENH: Add new Extendded Filename 'epsg' to set the projection of the output image with an EPSG code + tests + doc

BUG: updated multi baselines for ioTvMultiDatasetReading 1 and 2, 'Not specified based on custom spheroid' is not formatted is the same way between different gdal versions

TEST: added multi baseline for ioTvMultiDatasetReading tests when using gdal3 (projection layout is different)

BUG: check geometry validity before computing intersection to avoid errors due to self intersecting polygons in geos 3.8