Erreur d’exécution 7 mémoire insuffisante fichiers lourds

Question

Pour une application de métrologie, j'ai une série de cordonnées X et Y de 3000 points environ et je souhaite calculer le plus petit cercle circonscrit qui englobe tout les points. Je sais qu'un cercle est entièrement défini par 3 points donc j'ai réalisé une macro qui calcule tout les cercles possibles, qui sélectionne ceux qui englobent tout les points puis qui affiche le plus petit de tous. La macro marche nickel pour un nombre de points à 200-300 mais pour plus de points cela ne fonctionne plus j'ai un problème de "mémoire insuffisante". Bon je comprends bien que le calcul soit pas évident car le système doit calculer les diamètres de tous les trios de points donc une matrice 3000x3000x3000 soit 27 milliards de données. Avant de revoir entièrement mon code pour qu'il soit (beaucoup) plus rapide et plus léger, je souhaite savoir si quelqu'un parmi vous aurait une solution plus simple permettant de "débloquer" la mémoire ou je ne sais pas quoi.. Merci d'avance. Lucas...

AlgoPlus · Answer

Sans voir ni données ni code, pas facile de proposer quelque chose . Puisque l'on a les coordonnées de tous les points: le centre du cercle doit se trouver au milieu de la différence des ordonnées extrêmes et au milieu de la différence des abscisses extrêmes ? et son rayon la différence la plus grande ...

AlgoPlus · Answer

Après une lecture rapide :

Peut-être modifier la structure des test avec des goto . mais ce n'est pas ce qui va accélérer le traitement.

Pas compris l'utilité de tableaux à 3 dimensions.

Pour rendre le traitement plus rapide il faudrait passer par un tableau plutôt que de faire de multiples lectures/écritures dans la feuille, et copier le tableau final des résultats dans la feuille en fin de traitement

on peut appliquer ce que j'ai proposé sur un échantillon (de 1350 valeurs). Cet échantillon se déplaçant de 0.8 ( une cellule vers le bas).

A chaque itération les coordonnées du centre et le diamètre sont stockés dans un tableau.

En fin on choisit le plus grand ???

Pour comprendre ce que tu fais, il faut que je me familiarise avec tes index paramétrés.

A+

AlgoPlus · Answer

Un premier jet, qui ne prend pas tout en compte:

Sub Cercles()
Dim DL As Long, Plage As Range, T, DerIt As Long, V As Integer, TF()
Dim OK As Boolean, x As Long, D1 As Integer, D2 As Integer, D As Integer
Dim MinX As Integer, MinY As Integer, MaxX As Integer, MaxY As Integer
V = 1350 'nombre de valeurs dans l'échantillon
T = Range("C5:E" & Range("C" & Rows.Count).End(xlUp).Row) ' tableau des données
DerIt = UBound(T, 1) - 1350 'détermine la dernière itération

For i = 1 To DerIt
    MinX = 10000: MinY = 10000: MaxX = -10000: MaxY = -10000
    For j = 0 To V - 1
        If Not IsError(T(i + j, 2)) And Not IsError(T(i + j, 3)) Then
            If T(i + j, 2) > MaxX Then MaxX = T(i + j, 2)
            If T(i + j, 3) > MaxY Then MaxY = T(i + j, 3)
            If T(i + j, 2) < MinX Then MinX = T(i + j, 2)
            If T(i + j, 3) < MinY Then MinY = T(i + j, 3)

        End If
    Next
    x = x + 1 'augmentation taille du tableau final
    ReDim Preserve TF(1 To 5, 1 To x)
    D1 = MaxX - MinX
    D2 = MaxY - MinY
    If D1 > D2 Then
        D = D1
    Else
        D = D2
    End If
    TF(1, x) = i
    TF(2, x) = (MaxX + MinX) / 2
    TF(3, x) = (MaxY + MinY) / 2
    TF(4, x) = D
Next

Range("G5").Resize(UBound(TF, 2), UBound(TF, 1)) = Application.Transpose(TF)
End Sub

Edit :pas vu le message précédent avant d'envoyer celui ci

le cercle n'est pas toujours formé par les points extrêmes par rapport à l'origine

et mon code va donc à la poubelle !

galopin01 · Answer

Pour mois un tableau à 2 dimensions semble suffisant : (3000 lignes et 9 colonnes) pour appliquer les calculs du tableau de ton pdf. Si on pouvait disposer de toute la série de points... Après quand on dispose de ce tableau complet pour trouver le plus petit cercle, je ne suis sur de rien, mais ça ne me parait pas infaisable... En tout cas 3 dimensions me parait injustifié. A+...

eriiic · Answer

Oulalaaa, ça me parait bien compliqué tout ça. Pas besoin de calculer tous les cercles. Calcul le barycentre de tous tes points, ça sera le centre de ton cercle circonscrit. Sur le plan c'est simple. X= (moyenne des x), Y = (moyenne des y) Et pour avoir son rayon tu calcules sa distance à tous les points et tu prends le maxi. Distance de 2 points : =RACINE((x2-X)^2+(y2-Y)^2) eric Edit : enfin je pense... tu confirmeras...

eriiic · Answer

Effectivement, c'était un mauvais souvenir. Valable uniquement pour 3 points (triangle) Il te faut la valeur exacte ou une approchée suffit ...

eriiic · Answer

La valeur exacte peut faire plusieurs milliers de décimales

Si une valeur approchée à 5-6 décimales te suffit tu peux voir avec le solveur.

Tu peux changer la convergence dans les options

Tu peux aussi tester le modèle évolutionnaire mais c'est plus long. Le GRC a l'air de bien convenir.

Test avec 129 points

eric

11classeur1.xlsx (50.75 Ko)

AlgoPlus · Answer

On pourrait diminuer le nombre d'itérations (et ainsi la mémoire utilisée ?) dans la détermination des points d'un cercle.

avec

    For A = 4 + T To 4 + T + V
        For B = 4 + T To 4 + T + V
            For C = 4 + T To 4 + T + V
                If A = B Or B = C Or C = A Then

on va déterminer toutes les permutations possibles, donc on traitera, par exemple, Point 1-point 2-Point3, mais aussi Point 1-point 3-Point2, et Point 3-point 2-Point1, ....

Il suffit d'avoir 1 fois ces 3 point qu'importe l'ordre où on les trouve pour déterminer le cercle.

Avec

    For A = 4 + T To 4 + T + V
        For B =A + 1 To 4 + T + V
            For C =B + 1 To 4 + T + V

on élimine ces redondances.

Est ce que ça ça sera suffisant ?