factory_25831 = RGeo::Cartesian.factory(:srid => 25831)
point_25831 = factory_25831.point(505500.421875, 4699669.4375)
factory_4326 = RGeo::Geographic.spherical_factory(:srid => 4326)
point_4326 = RGeo::Feature.cast(point_25831, :factory => factory_4326, :project => true)
=> #<RGeo::Geographic::SphericalPointImpl:0xe35f6c "POINT (60.421875 90.0)">