Искажение высоты на сферическом изображении в Python

Я пытаюсь взять два прямоугольных изображения, одно из видимых элементов поверхности и одно, представляющее высоту, и сопоставить их с трехмерной сферой. Я знаю, как отображать объекты на сфере с помощью Cartopy, и я знаю, как делать карты поверхности рельефа, но я не могу найти простой способ объединить их, чтобы получить преувеличенную высоту на сферической проекции. Например, вот это сделано в MATLAB : Пример изображения

Кто-нибудь знает, есть ли простой способ сделать это в Python?


person vityav    schedule 01.09.2017    source источник


Ответы (1)


Мое решение не отвечает всем вашим требованиям. Но для начала это может быть хорошим стартером.

import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
import numpy as np
from matplotlib.cbook import get_sample_data
from matplotlib._png import read_png

# Use world image with shape (360 rows, 720 columns) 
pngfile = 'temperature_15-115.png'

fn = get_sample_data(pngfile, asfileobj=False)
img = read_png(fn)   # get array of color

# Some needed functions / constant
r = 5
pi = np.pi
cos = np.cos
sin = np.sin
sqrt = np.sqrt

# Prep values to match the image shape (360 rows, 720 columns)
phi, theta = np.mgrid[0:pi:360j, 0:2*pi:720j]

# Parametric eq for a distorted globe (for demo purposes)
x = r * sin(phi) * cos(theta)
y = r * sin(phi) * sin(theta)
z = r * cos(phi) + 0.5* sin(sqrt(x**2 + y**2)) * cos(2*theta)

fig = plt.figure()
fig.set_size_inches(9, 9)
ax = fig.add_subplot(111, projection='3d', label='axes1')

# Drape the image (img) on the globe's surface
sp = ax.plot_surface(x, y, z, \
                rstride=2, cstride=2, \
                facecolors=img)

ax.set_aspect(1)

plt.show()

Полученное изображение:

введите здесь описание изображения

person swatchai    schedule 23.09.2017
comment
У меня еще не было возможности попробовать, но, похоже, это именно то, что мне нужно. Спасибо! - person vityav; 03.10.2017