o tBhb6@sdZddlZddlZddlmZddlmZddlmZddl m Z m Z m Z ddl mZmZdd l mZdd lmZdd lmZdd lmZdd lmZGddde e e ZdS)zIsomap for manifold learningN)issparse) shortest_path)connected_components) BaseEstimatorTransformerMixin _ClassNamePrefixFeaturesOutMixin)NearestNeighborskneighbors_graph)radius_neighbors_graph)check_is_fitted) KernelPCA)KernelCenterer)_fix_connected_componentsc @s`eZdZdZddddddddddddd d d Zd d Zd dZdddZdddZddZ dS)IsomapaIsomap Embedding. Non-linear dimensionality reduction through Isometric Mapping Read more in the :ref:`User Guide `. Parameters ---------- n_neighbors : int or None, default=5 Number of neighbors to consider for each point. If `n_neighbors` is an int, then `radius` must be `None`. radius : float or None, default=None Limiting distance of neighbors to return. If `radius` is a float, then `n_neighbors` must be set to `None`. .. versionadded:: 1.1 n_components : int, default=2 Number of coordinates for the manifold. eigen_solver : {'auto', 'arpack', 'dense'}, default='auto' 'auto' : Attempt to choose the most efficient solver for the given problem. 'arpack' : Use Arnoldi decomposition to find the eigenvalues and eigenvectors. 'dense' : Use a direct solver (i.e. LAPACK) for the eigenvalue decomposition. tol : float, default=0 Convergence tolerance passed to arpack or lobpcg. not used if eigen_solver == 'dense'. max_iter : int, default=None Maximum number of iterations for the arpack solver. not used if eigen_solver == 'dense'. path_method : {'auto', 'FW', 'D'}, default='auto' Method to use in finding shortest path. 'auto' : attempt to choose the best algorithm automatically. 'FW' : Floyd-Warshall algorithm. 'D' : Dijkstra's algorithm. neighbors_algorithm : {'auto', 'brute', 'kd_tree', 'ball_tree'}, default='auto' Algorithm to use for nearest neighbors search, passed to neighbors.NearestNeighbors instance. n_jobs : int or None, default=None The number of parallel jobs to run. ``None`` means 1 unless in a :obj:`joblib.parallel_backend` context. ``-1`` means using all processors. See :term:`Glossary ` for more details. metric : str, or callable, default="minkowski" The metric to use when calculating distance between instances in a feature array. If metric is a string or callable, it must be one of the options allowed by :func:`sklearn.metrics.pairwise_distances` for its metric parameter. If metric is "precomputed", X is assumed to be a distance matrix and must be square. X may be a :term:`Glossary `. .. versionadded:: 0.22 p : int, default=2 Parameter for the Minkowski metric from sklearn.metrics.pairwise.pairwise_distances. When p = 1, this is equivalent to using manhattan_distance (l1), and euclidean_distance (l2) for p = 2. For arbitrary p, minkowski_distance (l_p) is used. .. versionadded:: 0.22 metric_params : dict, default=None Additional keyword arguments for the metric function. .. versionadded:: 0.22 Attributes ---------- embedding_ : array-like, shape (n_samples, n_components) Stores the embedding vectors. kernel_pca_ : object :class:`~sklearn.decomposition.KernelPCA` object used to implement the embedding. nbrs_ : sklearn.neighbors.NearestNeighbors instance Stores nearest neighbors instance, including BallTree or KDtree if applicable. dist_matrix_ : array-like, shape (n_samples, n_samples) Stores the geodesic distance matrix of training data. n_features_in_ : int Number of features seen during :term:`fit`. .. versionadded:: 0.24 feature_names_in_ : ndarray of shape (`n_features_in_`,) Names of features seen during :term:`fit`. Defined only when `X` has feature names that are all strings. .. versionadded:: 1.0 See Also -------- sklearn.decomposition.PCA : Principal component analysis that is a linear dimensionality reduction method. sklearn.decomposition.KernelPCA : Non-linear dimensionality reduction using kernels and PCA. MDS : Manifold learning using multidimensional scaling. TSNE : T-distributed Stochastic Neighbor Embedding. LocallyLinearEmbedding : Manifold learning using Locally Linear Embedding. SpectralEmbedding : Spectral embedding for non-linear dimensionality. References ---------- .. [1] Tenenbaum, J.B.; De Silva, V.; & Langford, J.C. A global geometric framework for nonlinear dimensionality reduction. Science 290 (5500) Examples -------- >>> from sklearn.datasets import load_digits >>> from sklearn.manifold import Isomap >>> X, _ = load_digits(return_X_y=True) >>> X.shape (1797, 64) >>> embedding = Isomap(n_components=2) >>> X_transformed = embedding.fit_transform(X[:100]) >>> X_transformed.shape (100, 2) Nrautor minkowski n_neighborsradius n_components eigen_solvertolmax_iter path_methodneighbors_algorithmn_jobsmetricp metric_paramsc CsL||_||_||_||_||_||_||_||_| |_| |_ | |_ | |_ dSNr) selfrrrrrrrrrrrr r#o/var/www/html/riverr-enterprise-integrations-main/venv/lib/python3.10/site-packages/sklearn/manifold/_isomap.py__init__s zIsomap.__init__c Cs|jdur|jdurtd|jdt|j|j|j|j|j|j|jd|_ |j ||j j |_ t |j dr<|j j |_ t|jd|j|j|j|jd|_|jdurct|j |j|j|j|jd|jd}nt|j |j|j|j|jd|jd }t|\}}|d kr|jdkrt|rtd |d tjd |d ddtd|j j|||d|j jd|j j}t||j dd|_!|j!d}|d9}|j"||_#|j#j$d |_%dS)Nz 1. The graph cannot be completed with metric='precomputed', and Isomap cannot befitted. Increase the number of neighbors to avoid this issue, or precompute the full distance matrix instead of passing a sparse neighbors graph.zm > 1. Completing the graph to fit Isomap might be slow. Increase the number of neighbors to avoid this issue.r) stacklevel)Xgraphn_connected_componentscomponent_labelsr+rF)methoddirectedr#)&rr ValueErrorr rrrr rnbrs_fitn_features_in_hasattrr'r rrrr kernel_pca_r r rr RuntimeErrorwarningswarnr_fit_Xeffective_metric_effective_metric_params_rr dist_matrix_ fit_transform embedding_shape_n_features_out)r"r.nbgr0labelsGr#r#r$_fit_transforms        zIsomap._fit_transformcCsNd|jd}t|}|jj}tt|dt|d|jdS)a(Compute the reconstruction error for the embedding. Returns ------- reconstruction_error : float Reconstruction error. Notes ----- The cost function of an isomap embedding is ``E = frobenius_norm[K(D) - K(D_fit)] / n_samples`` Where D is the matrix of distances for the input data X, D_fit is the matrix of distances for the output embedding X_fit, and K is the isomap kernel: ``K(D) = -0.5 * (I - 1/n_samples) * D^2 * (I - 1/n_samples)`` r4rr) rArrBr: eigenvalues_npsqrtsumrD)r"rHG_centerevalsr#r#r$reconstruction_errors ,zIsomap.reconstruction_errorcCs|||S)a Compute the embedding vectors for data X. Parameters ---------- X : {array-like, sparse graph, BallTree, KDTree, NearestNeighbors} Sample data, shape = (n_samples, n_features), in the form of a numpy array, sparse graph, precomputed tree, or NearestNeighbors object. y : Ignored Not used, present for API consistency by convention. Returns ------- self : object Returns a fitted instance of self. )rIr"r.yr#r#r$r76s z Isomap.fitcCs|||jS)aFit the model from data in X and transform X. Parameters ---------- X : {array-like, sparse graph, BallTree, KDTree} Training vector, where `n_samples` is the number of samples and `n_features` is the number of features. y : Ignored Not used, present for API consistency by convention. Returns ------- X_new : array-like, shape (n_samples, n_components) X transformed in the new space. )rIrCrQr#r#r$rBKs zIsomap.fit_transformcCst||jdur|jj|dd\}}n |jj|dd\}}|jj}|jd}t||f}t |D]}t |j ||||dddfd||<q2|dC}|d9}|j |S)aTransform X. This is implemented by linking the points X into the graph of geodesic distances of the training data. First the `n_neighbors` nearest neighbors of X are found in the training data, and from these the shortest geodesic distances from each point in X to each point in the training data are computed in order to construct the kernel. The embedding of X is the projection of this kernel onto the embedding vectors of the training set. Parameters ---------- X : array-like, shape (n_queries, n_features) If neighbors_algorithm='precomputed', X is assumed to be a distance matrix or a sparse graph of shape (n_queries, n_samples_fit). Returns ------- X_new : array-like, shape (n_queries, n_components) X transformed in the new space. NT)return_distancerrr4)r rr6 kneighborsradius_neighborsn_samples_fit_rDrKzerosrangeminrAr: transform)r"r. distancesindices n_samples_fit n_queriesG_Xir#r#r$rZ_s   0 zIsomap.transformr!) __name__ __module__ __qualname____doc__r%rIrPr7rBrZr#r#r#r$rs* ^   r)rdr<numpyrK scipy.sparserscipy.sparse.csgraphrrbaserrr neighborsr r r utils.validationr decompositionr preprocessingr utils.graphrrr#r#r#r$s