Transformation error using EPSG:23031

[fitz]

I got an exception trying to transform the following bounding box to EPSG:23031: 3.8142776, 51.285914,5.0541355,51.7975699,EPSG:4326. The exception message was "Infinite longitude".

I traced the problem into org.osgeo.proj4j.datum.Datum.transformToGeocentricFromWgs84. The code checks the length transform array to see if it is 3 or 7 parameter transformation.

Regardless of the type datum transformation, however, the array always has 7 elements. The org.osgeo.proj4j.parser.Proj4Parser.parseDatumTransform() does a split(",") on the string, resulting in an array with 7 values in it based on the format of the epsg file (i.e. +towgs84=-87,-98,-121,0,0,0,0 produces an array with 7 elements not 3).

[mdavis]

Fixed by detecting zero transform entries and creating 3-param transformation where appropriate.