python中如何输入坐标_如何使用Python将数组从球形坐标插入到笛卡尔坐标？

这样的事情应该有效：

import scipy.interpolate

rflat=scipy.array( r_coord.flat )

tflat=scipy.array( theta_coord.flat )

pflat=scipy.array( phi_coord.flat )

coordpoints=scipy.concatenate( [ rflat[:, scipy.newaxis], tflat[:,scipy.newaxis], pflat[:,scipy.newaxis] ] , axis=1 )

rtpinterpolator=scipy.interpolate.linearNDInterpolate( coordppoints, density.flat )

def xyz2rtp( x,y,z):

r=scipy.sqrt( x**2+y**2+z**2)

t=scipy.acos( z/r )

p=scipy.atan2( y, x )

return (r,t,p)

# now you can get the interpolated value for any (x,y,z) coordinate you want.

val=rtpinterpolator( xyz2rtp( x,y,z) )

关键点：

>使用现有的scipy多维插值,>将xyz坐标传递给插值器时将其转换为rtp.