私は同僚のIDLコードをPythonに書き直しているところですが、混乱しているいくつかの相違点があります。あなたがscipy.ndimage.interpolation.map_coordinates
を使用し、order=1
を指定すると、双線形補間を行うと思われる他のSOの質問とメーリングリストのスレッドによると、 IDLコード(GDLで実行)とpython(map_coordinates)の結果を比較すると、私は異なる結果を得ました。それから私はmpl_toolkits.basemap.interp
を使ってみました。私はIDLコードと同じ結果を得ました。以下は、何が間違っているかを示す簡単な例です。誰かが、私がmap_coordinates
で間違っていることを理解するのに役立つかもしれませんか、またはorder=1
は双線形ではありませんか?Scipy map_coordinates双線形補間とinterpおよびIDL補間との比較
from scipy.ndimage.interpolation import map_coordinates
from mpl_toolkits.basemap import interp
import numpy
in_data = numpy.array([[ 25.89125824, 25.88840675],[ 25.90930748, 25.90640068]], dtype=numpy.float32)
map_coordinates(in_data, [[0.0],[0.125]], order=1, mode='nearest')
# map_coordinates results in "25.89090157"
interp(in_data, numpy.array([0,1]), numpy.array([0,1]), numpy.array([0.0]), numpy.array([0.125]), order=1)
# interp results in "25.89351439", same as GDL's "25.8935" when printed
私は完全に罰金interp
を使用していますが、map_coordinates
は同じ結果を返しませんでした、なぜ私は興味がありました。私は、map_coordinates
のドキュメントにはbilinearについて言及していないことに気付きました。実際には双線形ですか?私は何が欠けていますか?
うわー、私はそのことを考えていなかったと信じてすることはできません。私はそれを試したと思った、ありがとうございました。 – daveydave400