pythonで作成した二次元強度マップのz（強度）をxで微分したいです。

Question

前提・実現したいこと

ある実験を行い、xとyの二つのパラメータを変化させて、zを測定しました。
その後、pythonを用いて、画像のような図を作成しました。
xが横軸に、yが縦軸に対応しています。色の変化がzに対応しています。

ここから、zをxで微分して、そのグラフを作成したいのですが、うまくいきません。

発生している問題・エラーメッセージ

エラーメッセージとして、次の文章が表示されます。

ValueError                                Traceback (most recent call last)
<ipython-input-4-67aedc53eb19> in <module>()
     13 y = numerical_diff(f,x)
     14 
---> 15 plt.plot(x, y)

~\Anaconda3\envs\Test1\lib\site-packages\matplotlib\pyplot.py in plot(*args, **kwargs)
   3356                       mplDeprecation)
   3357     try:
-> 3358         ret = ax.plot(*args, **kwargs)
   3359     finally:
   3360         ax._hold = washold

~\Anaconda3\envs\Test1\lib\site-packages\matplotlib\__init__.py in inner(ax, *args, **kwargs)
   1853                         "the Matplotlib list!)" % (label_namer, func.__name__),
   1854                         RuntimeWarning, stacklevel=2)
-> 1855             return func(ax, *args, **kwargs)
   1856 
   1857         inner.__doc__ = _add_data_doc(inner.__doc__,

~\Anaconda3\envs\Test1\lib\site-packages\matplotlib\axes\_axes.py in plot(self, *args, **kwargs)
   1525         kwargs = cbook.normalize_kwargs(kwargs, _alias_map)
   1526 
-> 1527         for line in self._get_lines(*args, **kwargs):
   1528             self.add_line(line)
   1529             lines.append(line)

~\Anaconda3\envs\Test1\lib\site-packages\matplotlib\axes\_base.py in _grab_next_args(self, *args, **kwargs)
    404                 this += args[0],
    405                 args = args[1:]
--> 406             for seg in self._plot_args(this, kwargs):
    407                 yield seg
    408 

~\Anaconda3\envs\Test1\lib\site-packages\matplotlib\axes\_base.py in _plot_args(self, tup, kwargs)
    381             x, y = index_of(tup[-1])
    382 
--> 383         x, y = self._xy_from_xy(x, y)
    384 
    385         if self.command == 'plot':

~\Anaconda3\envs\Test1\lib\site-packages\matplotlib\axes\_base.py in _xy_from_xy(self, x, y)
    240         if x.shape[0] != y.shape[0]:
    241             raise ValueError("x and y must have same first dimension, but "
--> 242                              "have shapes {} and {}".format(x.shape, y.shape))
    243         if x.ndim > 2 or y.ndim > 2:
    244             raise ValueError("x and y can be no greater than 2-D, but have "

ValueError: x and y must have same first dimension, but have shapes (1001,) and (17, 1001

該当のソースコード

``python

以下を行ったところ、エラーが出ました。

import numpy as np
import matplotlib.pyplot as plt

def numerical_diff(f,x):
    h = 1e-4
    return (f(x+h) - f(x-h)) / (2*h)

def f(x):
    return zdata

x = xdata

y = numerical_diff(f,x)

plt.plot(x, y)

### 試したこと

ここに問題に対して試したことを記載してください。

### 補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

Answer 1

過去に同じエラー文で質問をしている人がいました。
https://teratail.com/questions/94289

要するにxとyの次元の数がちがうために起こるエラーだそうです。
コードがかなりはしょっているようで確認できませんが心当たりはありますか？

Comments

[kid39]

ご回答ありがとうございます。
調べたところ、xの次元が1、yの次元が2と、両者の次元が異なっていました。

二つのパラメータを変化させて測定したデータを、パラメータの一つ(x)で微分したものをyと置きました。そのためにyの次元が2になったのだと思います。

xとyの次元を、同じ次元にする方法はありますでしょうか？

[Ryupe]

とりあえず質問文にコード全体を載せてもらえますか？
でないと、こちらは答えようがありません。

補足としてプログラムにおける次元というのは、
1次元…[1,2,3,4,5,6...]
2次元…[[1,2,3,4,5],[6,7,8,9,10]]
このような配列のことを指しています。

https://note.nkmk.me/python-numpy-ndarray-ndim-shape-size/

数学的な y^2 のようなことは言っていません。

[kid39]

承知しました。ご回答、ありがとうございました。