def _ortho_pixelize(self, data_source, field, bounds, size, antialias,
dim, periodic):
buff = pixelize_aitoff(data_source["py"], data_source["pdy"],
data_source["px"], data_source["pdx"],
size, data_source[field], None,
None, theta_offset = 0,
phi_offset = 0).transpose()
return buff
def _ortho_pixelize(self, data_source, field, bounds, size, antialias, dim,
periodic):
# For axis=2, x axis will be latitude, y axis will be longitude
px = (data_source["px"].d + 90) * np.pi / 180
pdx = data_source["pdx"].d * np.pi / 180
py = (data_source["py"].d + 180) * np.pi / 180
pdy = data_source["pdy"].d * np.pi / 180
# First one in needs to be the equivalent of "theta", which is
# longitude
buff = pixelize_aitoff(py, pdy, px, pdx, size, data_source[field],
None, None).transpose()
return buff
def _ortho_pixelize(self, data_source, field, bounds, size, antialias, dim,
periodic):
# For a radial axis, px will correspond to longitude and py will
# correspond to latitude.
px = (data_source["px"].d + 180) * np.pi / 180
pdx = data_source["pdx"].d * np.pi / 180
py = (data_source["py"].d + 90) * np.pi / 180
pdy = data_source["pdy"].d * np.pi / 180
# First one in needs to be the equivalent of "theta", which is
# longitude
buff = pixelize_aitoff(px, pdx, py, pdy, size, data_source[field],
None, None).transpose()
return buff
def _ortho_pixelize(self, data_source, field, bounds, size, antialias,
dim, periodic):
# For axis=2, x axis will be latitude, y axis will be longitude
px = (data_source["px"].d + 90) * np.pi/180
pdx = data_source["pdx"].d * np.pi/180
py = (data_source["py"].d + 180) * np.pi/180
pdy = data_source["pdy"].d * np.pi/180
# First one in needs to be the equivalent of "theta", which is
# longitude
buff = pixelize_aitoff(py, pdy, px, pdx,
size, data_source[field], None,
None).transpose()
return buff
def _ortho_pixelize(self, data_source, field, bounds, size, antialias, dim, periodic):
buff = pixelize_aitoff(
data_source["py"],
data_source["pdy"],
data_source["px"],
data_source["pdx"],
size,
data_source[field],
None,
None,
theta_offset=0,
phi_offset=0,
).transpose()
return buff
def _ortho_pixelize(self, data_source, field, bounds, size, antialias, dim,
periodic):
# use Aitoff projection
# http://paulbourke.net/geometry/transformationprojection/
bounds = tuple(_.ndview for _ in self._aitoff_bounds)
buff = pixelize_aitoff(
azimuth=data_source["py"],
dazimuth=data_source["pdy"],
colatitude=data_source["px"],
dcolatitude=data_source["pdx"],
buff_size=size,
field=data_source[field],
bounds=bounds,
input_img=None,
azimuth_offset=0,
colatitude_offset=0,
).transpose()
return buff
