/
spatial_transformer_grid.py
169 lines (135 loc) · 5.72 KB
/
spatial_transformer_grid.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
import numpy
import chainer
from chainer.backends import cuda
from chainer import function
from chainer.utils import argument
from chainer.utils import type_check
if cuda.cudnn_enabled:
cudnn = cuda.cudnn
libcudnn = cuda.cuda.cudnn
_sampler_type = libcudnn.CUDNN_SAMPLER_BILINEAR
class SpatialTransformerGrid(function.Function):
def __init__(self, output_shape):
self.output_shape = output_shape
def check_type_forward(self, in_types):
n_in = in_types.size()
type_check.expect(n_in == 1)
theta_type = in_types[0]
type_check.expect(
theta_type.dtype.char == 'f',
theta_type.ndim == 3,
theta_type.shape[1] == 2,
theta_type.shape[2] == 3,
)
def forward_cpu(self, inputs):
return self._forward(inputs)
def forward_gpu(self, inputs):
if not chainer.should_use_cudnn('>=auto', 5000):
return self._forward(inputs)
theta, = inputs
B, _, _ = theta.shape
H, W = self.output_shape
grid_t = cuda.cupy.empty((B, H, W, 2), dtype=theta.dtype)
# Unlike spatial_transformer_sampler,
# channel size can be anything in this case.
shape = numpy.array((B, 1, H, W), dtype=numpy.int32)
theta = cuda.cupy.ascontiguousarray(theta)
handle = cudnn.get_handle()
self.st_desc =\
cuda.cupy.cudnn.create_spatial_transformer_descriptor(
_sampler_type, grid_t.dtype, len(shape), shape.ctypes.data)
libcudnn.spatialTfGridGeneratorForward(
handle, self.st_desc.value, theta.data.ptr, grid_t.data.ptr)
grid = cuda.cupy.transpose(grid_t, (0, 3, 1, 2))
return grid,
def _forward(self, inputs):
theta, = inputs
H, W = self.output_shape
B, _, _ = theta.shape
xp = cuda.get_array_module(theta)
ys, xs = xp.meshgrid(
xp.linspace(-1, 1, H, dtype=numpy.float32),
xp.linspace(-1, 1, W, dtype=numpy.float32), indexing='ij',
copy=False
)
coords = xp.concatenate(
[xs[None], ys[None], xp.ones((1, H, W), dtype=numpy.float32)],
axis=0)
grid = theta.dot(coords.reshape(3, H * W)).reshape(B, 2, H, W)
return grid,
def backward_cpu(self, inputs, grad_outputs):
return self._backward(inputs, grad_outputs)
def backward_gpu(self, inputs, grad_outputs):
if not chainer.should_use_cudnn('>=auto', 5000):
return self._backward(inputs, grad_outputs)
theta, = inputs
ggrid, = grad_outputs
ggrid_t = cuda.cupy.transpose(ggrid, (0, 2, 3, 1))
gtheta = cuda.cupy.empty_like(theta)
handle = cudnn.get_handle()
ggrid_t = cuda.cupy.ascontiguousarray(ggrid_t)
libcudnn.spatialTfGridGeneratorBackward(
handle, self.st_desc.value, ggrid_t.data.ptr, gtheta.data.ptr)
return gtheta,
def _backward(self, inputs, grad_outputs):
theta, = inputs
ggrid, = grad_outputs
H, W = self.output_shape
B, _, _ = theta.shape
xp = cuda.get_array_module(theta)
ys, xs = xp.meshgrid(
xp.linspace(-1, 1, H, dtype=numpy.float32),
xp.linspace(-1, 1, W, dtype=numpy.float32), indexing='ij',
copy=False
)
coords = xp.concatenate(
[xs[None], ys[None], xp.ones((1, H, W), dtype=numpy.float32)],
axis=0)
coords_T = coords.reshape(3, H * W).transpose(1, 0)
ggrid = ggrid.reshape(B, 2, H * W)
gtheta = ggrid.dot(coords_T).reshape(B, 2, 3)
return gtheta,
def spatial_transformer_grid(theta, output_shape, **kwargs):
"""2D Spatial Transformer grid.
This function generates coordinates of the points sampled from an image
to perform warping described in `Spatial Transformer Networks \
<https://arxiv.org/abs/1506.02025>`_.
Given a coordinate in the warped image :math:`(x_i^t, y_i^t)`, the point
sampled from the source image :math:`(x_i^s, y_i^s)` are calculated
by the following equation.
.. note::
cuDNN supports SpatialTransformerGrid from version 5.0.0.
.. math::
\\left(\\begin{matrix} x_i^s \\\\
y_i^s \\end{matrix}\\right)
=
\\left(\\begin{matrix} \\theta_{11} & \\theta_{12} & \\theta_{13} \\\\
\\theta_{21} & \\theta_{22} & \\theta_{23} \\end{matrix}\\right)
\\left(\\begin{matrix} x_i^t \\\\
y_i^t \\\\
1 \\end{matrix}\\right)
Notatition: here is a notation for dimensionalities.
- :math:`n` is the batch size.
- :math:`h_O` and :math:`w_O` are the height and the width of the output
image.
Args:
theta (~chainer.Variable): An array of shape :math:`(n, 2, 3)`.
This is a batch of :math:`2 \\times 3` matrix used for
the warping described above.
output_shape (tuple): A tuple of 2 elements: :math:`h_O, w_O`.
Returns:
~chainer.Variable: A variable of shape :math:`(n, 2, h_O, w_O)`.
In the 2nd dimension, the first element is the coordinate along the
x axis, and the second element is the coordinate along the y axis.
All the coordinates in the image are scaled to fit range
:math:`[-1, 1]`.
This means that the coordinate :math:`(-1, -1)` corresponds to
the upper-left corner of the input image.
"""
argument.check_unexpected_kwargs(
kwargs, use_cudnn="The argument \"use_cudnn\" is not "
"supported anymore. "
"Use chainer.using_config('use_cudnn', value) "
"context where value can be `always`, `never`, or `auto`.")
argument.assert_kwargs_empty(kwargs)
return SpatialTransformerGrid(output_shape)(theta)