Intentar una traducción directa del código Matlab a python no es lo mejor en este caso.
Ya que tienes los datos en una excel, puedes beneficiarte de la biblioteca pandas
para leerlos:
import pandas as pd
xe = pd.read_excel("anexo.xlsx", header=None, sheet_name="xe")
nee = pd.read_excel("anexo.xlsx", header=None, sheet_name="nee")
Una vez los tenemos leídos, hay que procesarlos como deseas antes de pintarlos. Aquí es donde la cosa se vuelve un poco confusa, pues el código Matlab que pones no lo entiendo (y tampoco su traducción a Python), pues pareces estar queriendo pintar una superficie 3D por un lado, pero en realidad un scatter plot por otro.
Voy a darte la implementación de lo que he entendido que quieres hacer.
Entiendo que la hoja nee
contiene en realidad índices a coordenadas que están en la hoja xe
, pero que deben ser transformados, con la fórmula que has indicado: multiplicar por 3 el valor encontrado en nee
y restar 3 para obtener el índice de la coordenada x
deseada, o bien restar 2 para obtener el de la coordenada y
deseada. Estos índices se usan sobre el vector guardado en xe
para extraer los valores x,y
a plotear.
Pandas permite hacer operaciones "vectoriales", esto es, que operan con todos los valores de una columna "a la vez" (en realidad internamente pandas itera por esos valores, pero tu código python no necesita bucles, sino que usa una sintaxis "vectorial").
Es decir, podríamos agregar al dataframe nee
un par de columnas adicionales que serían los índices a usar en xe
para sacar x
e y
, con solo dos líneas de código:
nee["ix"] = nee.iloc[:,0]*3-3
nee["iy"] = nee.iloc[:,0]*3-2
Aquí, nee.loc[]
permite acceder a elementos concretos de la tabla ne
según sus coordenadas. Por ejemplo nee.iloc[i,j]
iría al elemento de la fila i
, columna j
. Pero también permite acceder a "columnas completas" de una sola vez, si en el índice de la fila pones :
. Así pues nee.iloc[:,0]
extrae la columna 0 completa.
Como ves esa columna se usa en la operación *3-3
, y pandas multiplicará por 3 y restará 3 a todos los elementos de esa columna. El resultado lo guardamos en una columna nueva que hemos llamado ix
(por "índices de x"). Análogamente computamos la columna nueva iy
("índices de y").
Una vez tenemos éstas, podemos extraer las coordenadas x
necesarias de la tabla xe
con esta sintaxis:
x = xe.iloc[nee.ix]
y = xe.iloc[nee.iy]
De nuevo esto es una operación "vectorial", que actúa sobre toda la tabla "a la vez". El resultado es que tanto x
como y
son vectores columna, que se pueden usar ya directamente en un plot: plt.scatter(x, y)
.
Todo junto
Ahora bien, la tabla nee
tiene tres columnas, y ahí es donde me pierdo. Aparentemente en tu código no distingues de ninguna forma de qué columna vienen los datos, y te limitas a iterar por ellas.
Para mayor claridad yo voy a iterar también sobre las tres columnas (0, 1 y 2) de la tabla nee
, pero voy a usar un color diferente (y una leyenda diferente) para los puntos extraídos de cada columna.
El código completo sería entonces:
import matplotlib.pyplot as plt
import pandas as pd
xe = pd.read_excel("anexo.xlsx", header=None, sheet_name="xe")
nee = pd.read_excel("anexo.xlsx", header=None, sheet_name="nee")
for i in range(3):
nee["ix"] = nee.loc[:,i]*3-3
nee["iy"] = nee.loc[:,i]*3-2
x = xe.iloc[nee.ix]
y = xe.iloc[nee.iy]
plt.scatter(x, y, label=f"Column {i}", marker="s")
plt.gcf().set_size_inches(10,4)
plt.legend()
Y este es el resultado:

Como ves, los puntos azules (de la columna 0) no se ven, entiendo que porque quedan tapados por los verdes.
Bonus
Otra versión aún más breve si realmente no quieres separar por colores las columnas.
En ese caso no necesitamos ni siquiera iterar por esas columnas, y podemos aplicar la transformación *3-3
a la tabla completa, o mejor a sus valores. El resultado será otra tabla (de 112x3 como la tabla original) en la que los valores ya están multiplicados por 3 y restados. Mediante la función .reshape(-1)
podemos convertir esa matriz de 112x3 a un solo vector de 336 elementos, que podemos usar como índices sobre xe
.
Podemos usar pandas para leer la excel, pero seguidamente extraer solo la matriz de valores (con .values
). El resultado es un array numpy, que nos permite una sintaxis más corta al evitar tener que usar iloc
.
Todo el código se reduce ahora a tres líneas:
xe = pd.read_excel("anexo.xlsx", header=None, sheet_name="xe").values
nee = pd.read_excel("anexo.xlsx", header=None, sheet_name="nee").values
plt.scatter(xe[(nee*3-3).reshape(-1)], xe[(nee*3-2).reshape(-1)])
Lo que produce la gráfica:
