train_classifier_img_large.py 4.7 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221
  1. #!/usr/bin/env python3
  2. import time
  3. import random
  4. import numpy as np
  5. import gym
  6. from gym_minigrid.register import env_list
  7. from gym_minigrid.minigrid import Grid, OBJECT_TO_IDX
  8. import babyai
  9. import matplotlib.pyplot as plt
  10. import torch
  11. import torch.nn as nn
  12. import torch.optim as optim
  13. import torch.nn.functional as F
  14. from torch.autograd import Variable
  15. import torchvision
  16. import numpy as np
  17. import cv2
  18. import PIL
  19. ##############################################################################
  20. def make_var(arr):
  21. arr = np.ascontiguousarray(arr)
  22. #arr = torch.from_numpy(arr).float()
  23. arr = torch.from_numpy(arr)
  24. arr = Variable(arr)
  25. if torch.cuda.is_available():
  26. arr = arr.cuda()
  27. return arr
  28. def init_weights(m):
  29. classname = m.__class__.__name__
  30. if classname.startswith('Conv'):
  31. nn.init.orthogonal_(m.weight.data)
  32. m.bias.data.fill_(0)
  33. elif classname.find('Linear') != -1:
  34. nn.init.xavier_uniform_(m.weight)
  35. m.bias.data.fill_(0)
  36. elif classname.find('BatchNorm') != -1:
  37. m.weight.data.normal_(1.0, 0.02)
  38. m.bias.data.fill_(0)
  39. def num_params(model):
  40. pp=0
  41. for p in list(model.parameters()):
  42. nn=1
  43. for s in list(p.size()):
  44. nn = nn*s
  45. pp += nn
  46. return pp
  47. class Flatten(nn.Module):
  48. """
  49. Flatten layer, to flatten convolutional layer output
  50. """
  51. def forward(self, input):
  52. return input.view(input.size(0), -1)
  53. class Print(nn.Module):
  54. def forward(self, input):
  55. print(input.size())
  56. return input
  57. class Model(nn.Module):
  58. def __init__(self):
  59. super().__init__()
  60. self.layers = nn.Sequential(
  61. #Print(),
  62. nn.Conv2d(in_channels=3, out_channels=16, kernel_size=6, stride=2),
  63. nn.LeakyReLU(),
  64. nn.Conv2d(in_channels=16, out_channels=16, kernel_size=6, stride=2),
  65. nn.LeakyReLU(),
  66. nn.Conv2d(in_channels=16, out_channels=16, kernel_size=6, stride=2),
  67. nn.LeakyReLU(),
  68. nn.Conv2d(in_channels=16, out_channels=16, kernel_size=6, stride=2),
  69. nn.LeakyReLU(),
  70. #Print(),
  71. Flatten(),
  72. nn.Linear(144, 128),
  73. nn.LeakyReLU(),
  74. nn.Linear(128, 2),
  75. )
  76. self.apply(init_weights)
  77. def forward(self, obs):
  78. obs = obs / 16
  79. out = self.layers(obs)
  80. return out
  81. def present_prob(self, obs):
  82. obs = make_var(obs).unsqueeze(0)
  83. logits = self(obs)
  84. probs = F.softmax(logits, dim=-1)
  85. probs = probs.detach().cpu().squeeze().numpy()
  86. return probs[1]
  87. env = gym.make('BabyAI-GoToRedBall-v0')
  88. def sample_batch(batch_size=128):
  89. imgs = []
  90. labels = []
  91. for i in range(batch_size):
  92. obs = env.reset()['image']
  93. ball_visible = ('red', 'ball') in Grid.decode(obs)
  94. obs = env.get_obs_render(obs, tile_size=16, mode='rgb_array')
  95. #plt.imshow(obs)
  96. #plt.show()
  97. obs = obs.transpose([2, 0, 1])
  98. imgs.append(np.copy(obs))
  99. labels.append(ball_visible)
  100. imgs = np.stack(imgs).astype(np.float32)
  101. labels = np.array(labels, dtype=np.long)
  102. return imgs, labels
  103. print('Generating test set')
  104. test_imgs, test_labels = sample_batch(256)
  105. def eval_model(model):
  106. num_true = 0
  107. for idx in range(test_imgs.shape[0]):
  108. img = test_imgs[idx]
  109. label = test_labels[idx]
  110. p = model.present_prob(img)
  111. out_label = p > 0.5
  112. #print(out_label)
  113. if np.equal(out_label, label):
  114. num_true += 1
  115. #else:
  116. # if label:
  117. # print("incorrectly predicted as absent")
  118. # else:
  119. # print("incorrectly predicted as present")
  120. acc = 100 * (num_true / test_imgs.shape[0])
  121. return acc
  122. ##############################################################################
  123. batch_size = 64
  124. model = Model()
  125. model.cuda()
  126. print('Num params:', num_params(model))
  127. optimizer = optim.Adam(
  128. model.parameters(),
  129. lr=5e-4
  130. )
  131. criterion = nn.CrossEntropyLoss()
  132. running_loss = None
  133. for batch_no in range(1, 10000):
  134. batch_imgs, labels = sample_batch(batch_size)
  135. batch_imgs = make_var(batch_imgs)
  136. labels = make_var(labels)
  137. pred = model(batch_imgs)
  138. loss = criterion(pred, labels)
  139. optimizer.zero_grad()
  140. loss.backward()
  141. optimizer.step()
  142. loss = loss.data.detach().item()
  143. running_loss = loss if running_loss is None else 0.99 * running_loss + 0.01 * loss
  144. print('batch #{}, frames={}, loss={:.5f}'.format(
  145. batch_no,
  146. batch_no * batch_size,
  147. running_loss
  148. ))
  149. if batch_no % 25 == 0:
  150. acc = eval_model(model)
  151. print('accuracy: {:.2f}%'.format(acc))