Foros Club Delphi

Foros Club Delphi (https://www.clubdelphi.com/foros/index.php)
-   Python (https://www.clubdelphi.com/foros/forumdisplay.php?f=46)
-   -   Lista bidimensional de python a txt (https://www.clubdelphi.com/foros/showthread.php?t=96640)

gusspagano 11-03-2024 19:19:39

Lista bidimensional de python a txt
 
Antes que nada, gracias por su tiempo y su generosa ayuda.

Tengo un archivo txt conformado por 344,160 lineas. Cada linea tiene 3 valores: coordenada x, coordenada y, dato. Estos 3 valores estan separados por espacio. Ej: 338619.4247572054 2304584.9999724245 -10.0.

Lo que intento es meter estos 344,160 datos en una lista bidimensional de tamaño 100 x 100. Utilizo las coordenadas x, y de cada linea para calcular la columna y renglon en la cual debo colocar cada dato. Dado que la lista unicame tiente 10,000 posiciones (100x100 = 10,000) y el txt tiene 344,160 datos, muchos datos van a caer dentro de una misma posición dentro de la lista, por lo cual estoy haciendo un promedio. Por otra parte, habran posiciones dentro de la lista en los cuales no caerá ningún dato. A estas posiciones de la lista les agrego el valor -9999 para indicar que es un dato nulo.

El proceso que sigo es el siguiente:
  • creo una lista vacia (lista = []) a la cual posteriormente agrego 100 veces dato nulo (-9999)
  • creo otra lista vacia (listadatos = []) a la cual agrego 100 veces la lista anterior
  • viciamos la lista "lista":lista = [] y la llenamos ahora con 100 datos = 0
  • creamos otra lista vacia llamada "listacant" la cual contiene la cantidad de datos que han caido en cada posición de la lista "listadatos"
  • llenamos listacant con 100 veces la lista "lista"


De este modo tengo dos lista. listadatos que contiene hasta ahora 10,000 valores iguales a -9999 y listacant que cotiene 10,000 valores iguales a cero.
para promediar los datos y colocarlos en listadatos, verifico en primer lugar si la posición ya contiene valores distintos a -9999 o no, y deacuerdo con esto hago el promedio.



Código:

for i in range (344160):
            linea = temp_arc.readline()                                            # lee una linea del txt 
            xdato, ydato, zdato = linea.split(' ')                                # se obtienen coordenadas y dato por separado
            ncol = math.trunc ( ( float ( xdato ) - xmin ) / rsl )        # calculo de columna. rsl = resolución
            nfil = math.trunc ( ( ymax - float ( ydato ) ) / rsl )        # calculo de fila
            if listadatos[ncol][nfil] == -9999:                                #
                listadatos[ncol][nfil] = float ( zdato )
            else:
                listadatos[ncol][nfil] = ( ( listadatos[ncol][nfil]
                                        * listacant[ncol][nfil]
                                        + float (zdato) )
                                          / ( listacant[ncol][nfil] + 1 ) )
            listacant[ncol][nfil] += 1

El problema es que al verificar el resutado, el codigo me pone 100 veces el mismo dato para una misma fila. la siguiente fila tiene un dato distinto que igualmente se repite 100 veces y así suscecivamente.
antes habia hecho el codigo con delphi y corria sin problemas, pero aquí no entiendo que pasa.

Me disculpo por lo extenso de mi consulta y ojala alguien pueda darme una pista de lo que estoy haciendo mal.

Saludos y gracias nuevamente.

duilioisola 12-03-2024 14:16:04

No es necesario calcular el promedio cada vez que modificas una celda.
Creo que lo mejor es dividirlo en varios pasos.
Además no creo que debas utilizar el valor -9999 como indicador de nulo. Yo simplemente dejaría la celda con un 0.

Te dejo un pseudocódigo Delphi:
Código Delphi [-]
lista_val = array[344160];
lista_Can = array[100, 100];
lista_Sum = array[100, 100];


// Inicializamos listas ----------
RellenarListaValores(lista_val);
RellenarConCeros(lista_Can);
RellenarConCeros(lista_Sum);


// Acumulamos valores ------------
for i = 0 to lista_val.Count do
begin
  // Dividimos la linea leida en 3 valores
  ObtenerDatos(lista_val[i], x, y, z);
  
  // Averiguamos en qué posición debe guardar el valor
  CalcularPosicion(x, y, columna, fila);
  
  // incrementamos la cantidad de 
  lista_Can[columna, fila] := lista_Can[columna, fila] + 1;
  
  // Sumamos el valor leido
  lista_Sum[columna, fila] := lista_Sum[columna, fila] + z;
end;

  
// Calculamos promedios ----------
for Columna = 1 to 100 do
begin
  for Fila = 1 to 100 do
    if lista_Can[columna, fila] <> 0 then
       // Calculo Promedio
       lista_Sum[columna, fila] := lista_Sum[columna, fila] / lista_Can[columna, fila]
    else     
       // No hay datos y se considera nulo
       lista_Sum[columna, fila] := -9999;
end;

gusspagano 16-03-2024 01:48:09

Cierto, no es necesario promediar cada vez.

Tengo un código en delhi que corre sin problemas. La cosa es que quiero hacerlo con python porque, posterirmente voy a utlizar este código en QGis. En python no hace lo que debería y no entiendo por qué.

Gracias por tu ayuda!


La franja horaria es GMT +2. Ahora son las 02:37:14.

Powered by vBulletin® Version 3.6.8
Copyright ©2000 - 2024, Jelsoft Enterprises Ltd.
Traducción al castellano por el equipo de moderadores del Club Delphi