def main(args): in_channels = 3 in_height = 224 in_width = 224 print("Importing ONNX model") net = eddl.import_net_from_onnx_file(args.model_fn, [in_channels, in_height, in_width]) # Add a softmax layer to get probabilities directly from the model input_ = net.lin[0] # getLayer(net,"input_layer_name") output = net.lout[0] # getLayer(net,"output_layer_name") new_output = eddl.Softmax(output) net = eddl.Model([input_], [new_output]) eddl.build( net, eddl.adam(0.001), # not used for prediction ["softmax_cross_entropy"], # not used for prediction ["categorical_accuracy"], # not used for prediction eddl.CS_GPU() if args.gpu else eddl.CS_CPU(), False # Disable model initialization, we want to use the ONNX weights ) eddl.summary(net) image = Tensor.load(args.img_fn) image_preprocessed = preprocess_input_resnet34(image, [in_height, in_width]) outputs = eddl.predict(net, [image_preprocessed]) print("Reading class names...") with open(args.classes_fn, "rt") as f: class_names = [_.strip() for _ in f] print("Top 5 predictions:") print(eddl.get_topk_predictions(outputs[0], class_names, 5))
def main(args): eddl.download_cifar10() num_classes = 10 in_ = eddl.Input([3, 32, 32]) layer = in_ layer = eddl.MaxPool(eddl.ReLu(Normalization( eddl.Conv(layer, 32, [3, 3], [1, 1]) )), [2, 2]) layer = eddl.MaxPool(eddl.ReLu(Normalization( eddl.Conv(layer, 64, [3, 3], [1, 1]) )), [2, 2]) layer = eddl.MaxPool(eddl.ReLu(Normalization( eddl.Conv(layer, 128, [3, 3], [1, 1]) )), [2, 2]) layer = eddl.MaxPool(eddl.ReLu(Normalization( eddl.Conv(layer, 256, [3, 3], [1, 1]) )), [2, 2]) layer = eddl.GlobalMaxPool(layer) layer = eddl.Flatten(layer) layer = eddl.Activation(eddl.Dense(layer, 128), "relu") out = eddl.Softmax(eddl.Dense(layer, num_classes)) net = eddl.Model([in_], [out]) eddl.build( net, eddl.adam(0.001), ["soft_cross_entropy"], ["categorical_accuracy"], eddl.CS_GPU(mem=args.mem) if args.gpu else eddl.CS_CPU(mem=args.mem) ) eddl.summary(net) eddl.plot(net, "model.pdf") x_train = Tensor.load("cifar_trX.bin") y_train = Tensor.load("cifar_trY.bin") x_train.div_(255.0) x_test = Tensor.load("cifar_tsX.bin") y_test = Tensor.load("cifar_tsY.bin") x_test.div_(255.0) if args.small: x_train = x_train.select([":5000"]) y_train = y_train.select([":5000"]) x_test = x_test.select([":1000"]) y_test = y_test.select([":1000"]) for i in range(args.epochs): eddl.fit(net, [x_train], [y_train], args.batch_size, 1) eddl.evaluate(net, [x_test], [y_test], bs=args.batch_size) print("All done")
def main(args): eddl.download_mnist() in_ = eddl.Input([784]) target = eddl.Reshape(in_, [1, 28, 28]) layer = in_ layer = eddl.Reshape(layer, [1, 28, 28]) layer = eddl.ReLu(eddl.Conv(layer, 8, [3, 3])) layer = eddl.ReLu(eddl.Conv(layer, 16, [3, 3])) layer = eddl.ReLu(eddl.Conv(layer, 8, [3, 3])) out = eddl.Sigmoid(eddl.Conv(layer, 1, [3, 3])) net = eddl.Model([in_], []) eddl.build( net, eddl.adam(0.001), [], [], eddl.CS_GPU(mem=args.mem) if args.gpu else eddl.CS_CPU(mem=args.mem) ) eddl.summary(net) x_train = Tensor.load("mnist_trX.bin") if args.small: x_train = x_train.select([":6000"]) x_train.div_(255.0) mse = eddl.newloss(mse_loss, [out, target], "mse_loss") dicei = eddl.newloss(dice_loss_img, [out, target], "dice_loss_img") dicep = eddl.newloss(dice_loss_pixel, [out, target], "dice_loss_pixel") batch = Tensor([args.batch_size, 784]) num_batches = x_train.shape[0] // args.batch_size for i in range(args.epochs): print("Epoch %d/%d (%d batches)" % (i + 1, args.epochs, num_batches)) diceploss = 0.0 diceiloss = 0.0 mseloss = 0 for j in range(num_batches): print("Batch %d " % j, end="", flush=True) eddl.next_batch([x_train], [batch]) eddl.zeroGrads(net) eddl.forward(net, [batch]) diceploss += eddl.compute_loss(dicep) / args.batch_size print("diceploss = %.6f " % (diceploss / (j + 1)), end="") diceiloss += eddl.compute_loss(dicei) / args.batch_size print("diceiloss = %.6f " % (diceiloss / (j + 1)), end="") mseloss += eddl.compute_loss(mse) / args.batch_size print("mseloss = %.6f\r" % (mseloss / (j + 1)), end="") eddl.optimize(dicep) eddl.update(net) print() print("All done")
def main(args): eddl.download_imdb_2000() epochs = 2 if args.small else 10 length = 250 embdim = 33 vocsize = 2000 in_ = eddl.Input([1]) # 1 word layer = in_ layer = eddl.RandomUniform(eddl.Embedding(layer, vocsize, 1, embdim), -0.05, 0.05) layer = eddl.GRU(layer, 37) layer = eddl.ReLu(eddl.Dense(layer, 256)) out = eddl.Sigmoid(eddl.Dense(layer, 1)) net = eddl.Model([in_], [out]) eddl.build( net, eddl.adam(0.001), ["cross_entropy"], ["binary_accuracy"], eddl.CS_GPU(mem=args.mem) if args.gpu else eddl.CS_CPU(mem=args.mem)) eddl.summary(net) x_train = Tensor.load("imdb_2000_trX.bin") y_train = Tensor.load("imdb_2000_trY.bin") x_test = Tensor.load("imdb_2000_tsX.bin") y_test = Tensor.load("imdb_2000_tsY.bin") # batch x timesteps x input_dim x_train.reshape_([x_train.shape[0], length, 1]) x_test.reshape_([x_test.shape[0], length, 1]) y_train.reshape_([y_train.shape[0], 1, 1]) y_test.reshape_([y_test.shape[0], 1, 1]) if args.small: x_train = x_train.select([":64", ":", ":"]) y_train = y_train.select([":64", ":", ":"]) x_test = x_test.select([":64", ":", ":"]) y_test = y_test.select([":64", ":", ":"]) for i in range(epochs): eddl.fit(net, [x_train], [y_train], args.batch_size, 1) eddl.evaluate(net, [x_test], [y_test]) print("All done")
def main(args): size = 256 // 2 # Conv3D expects (B, C, dim1, dim2, dim3) in_ = eddl.Input([3, 10, size, size]) layer = in_ layer = eddl.MaxPool3D(eddl.ReLu(eddl.Conv3D( layer, 4, [1, 3, 3], [1, 1, 1], "same" )), [1, 2, 2], [1, 2, 2], "same") layer = eddl.MaxPool3D(eddl.ReLu(eddl.Conv3D( layer, 8, [1, 3, 3], [1, 1, 1], "same" )), [1, 2, 2], [1, 2, 2], "same") layer = eddl.MaxPool3D(eddl.ReLu(eddl.Conv3D( layer, 16, [1, 3, 3], [1, 1, 1], "same" )), [1, 2, 2], [1, 2, 2], "same") layer = eddl.GlobalMaxPool3D(layer) layer = eddl.Reshape(layer, [-1]) layer = eddl.LSTM(layer, 128) layer = eddl.Dense(layer, 100) layer = eddl.ReLu(layer) layer = eddl.Dense(layer, 2) out = eddl.ReLu(layer) net = eddl.Model([in_], [out]) eddl.build( net, eddl.adam(), ["mse"], ["mse"], eddl.CS_GPU(mem=args.mem) if args.gpu else eddl.CS_CPU(mem=args.mem) ) eddl.summary(net) seqImages = Tensor.randu([32, 10, 3, 10, size, size]) seqLabels = Tensor.randu([32, 7, 2]) eddl.fit(net, [seqImages], [seqLabels], 4, 2 if args.small else 10)
def main(args): num_classes = 1 size = [512, 512] # size of images thresh = 0.5 best_dice = -1 if args.out_dir: os.makedirs(args.out_dir, exist_ok=True) in_ = eddl.Input([1, size[0], size[1]]) out = SegNetBN(in_, num_classes) out_sigm = eddl.Sigmoid(out) net = eddl.Model([in_], [out_sigm]) eddl.build(net, eddl.adam(0.0001), ["cross_entropy"], ["mean_squared_error"], eddl.CS_GPU([1], mem='low_mem') if args.gpu else eddl.CS_CPU()) eddl.summary(net) eddl.setlogfile(net, "pneumothorax_segmentation_training") if args.ckpts and os.path.exists(args.ckpts): print("Loading checkpoints '{}'".format(args.ckpts)) eddl.load(net, args.ckpts, 'bin') training_augs = ecvl.SequentialAugmentationContainer([ ecvl.AugResizeDim(size), ecvl.AugMirror(0.5), ecvl.AugRotate([-10, 10]), ecvl.AugBrightness([0, 30]), ecvl.AugGammaContrast([0, 3]), ]) validation_augs = ecvl.SequentialAugmentationContainer( [ecvl.AugResizeDim(size)]) dataset_augs = ecvl.DatasetAugmentations( [training_augs, validation_augs, None]) print("Reading dataset") d = ecvl.DLDataset(args.in_ds, args.batch_size, dataset_augs, ecvl.ColorType.GRAY) # Prepare tensors which store batch x = Tensor([args.batch_size, d.n_channels_, size[0], size[1]]) y = Tensor([args.batch_size, d.n_channels_gt_, size[0], size[1]]) # Retrieve indices of images with a black ground truth # which are not include in a split train_split = d.GetSplit(ecvl.SplitType.training) val_split = d.GetSplit(ecvl.SplitType.validation) test_split = d.GetSplit(ecvl.SplitType.test) all_split = set(train_split + val_split + test_split) images_list = set(range(len(d.samples_))) # Obtain images with black ground truth black_images = images_list - all_split # Add a 25% of training samples with black ground truth. num_samples_train = math.floor(len(train_split) * 1.25) num_batches_train = num_samples_train // args.batch_size # Add a 25% of validation samples with black ground truth. num_samples_validation = math.floor(len(val_split) * 1.25) num_batches_validation = num_samples_validation // args.batch_size black_images = list(black_images) black_training = black_images[0:-(num_samples_validation - len(val_split))] black_validation = black_images[-(num_samples_validation - len(val_split)):] indices = list(range(args.batch_size)) evaluator = utils.Evaluator() print("Starting training") for e in range(args.epochs): print("Epoch {:d}/{:d} - Training".format(e + 1, args.epochs), flush=True) d.SetSplit(ecvl.SplitType.training) eddl.reset_loss(net) s = d.GetSplit() random.shuffle(s) d.split_.training_ = s random.shuffle(black_training) d.ResetAllBatches() # Indices to track mask and black vector in PneumothoraxLoadBatch m_i = 0 b_i = 0 for i, b in enumerate(range(num_batches_train)): d, images, labels, _, m_i, b_i = PneumothoraxLoadBatch( d, black_training, m_i, b_i) x, y = fill_tensors(images, labels, x, y) x.div_(255.0) y.div_(255.0) eddl.train_batch(net, [x], [y], indices) if i % args.log_interval == 0: print("Epoch {:d}/{:d} (batch {:d}/{:d}) - ".format( e + 1, args.epochs, b + 1, num_batches_train), end="", flush=True) eddl.print_loss(net, b) print() d.SetSplit(ecvl.SplitType.validation) evaluator.ResetEval() print("Epoch %d/%d - Evaluation" % (e + 1, args.epochs), flush=True) m_i = 0 b_i = 0 for b in range(num_batches_validation): n = 0 print("Epoch {:d}/{:d} (batch {:d}/{:d}) ".format( e + 1, args.epochs, b + 1, num_batches_validation), end="", flush=True) d, images, labels, names, m_i, b_i = PneumothoraxLoadBatch( d, black_validation, m_i, b_i) x, y = fill_tensors(images, labels, x, y) x.div_(255.0) y.div_(255.0) eddl.forward(net, [x]) output = eddl.getOutput(out_sigm) # Compute Dice metric and optionally save the output images for k in range(args.batch_size): pred = output.select([str(k)]) gt = y.select([str(k)]) pred_np = np.array(pred, copy=False) gt_np = np.array(gt, copy=False) # DiceCoefficient modifies image as a side effect dice = evaluator.DiceCoefficient(pred_np, gt_np, thresh=thresh) print("- Dice: {:.6f} ".format(dice), end="", flush=True) if args.out_dir: # Save original image fused together with prediction and # ground truth pred_np *= 255 pred_ecvl = ecvl.TensorToImage(pred) pred_ecvl.colortype_ = ecvl.ColorType.GRAY pred_ecvl.channels_ = "xyc" ecvl.ResizeDim(pred_ecvl, pred_ecvl, (1024, 1024), ecvl.InterpolationType.nearest) filename_gt = names[n + 1] gt_ecvl = ecvl.ImRead(filename_gt, ecvl.ImReadMode.GRAYSCALE) filename = names[n] # Image as BGR img_ecvl = ecvl.ImRead(filename) ecvl.Stack([img_ecvl, img_ecvl, img_ecvl], img_ecvl) img_ecvl.channels_ = "xyc" img_ecvl.colortype_ = ecvl.ColorType.BGR image_np = np.array(img_ecvl, copy=False) pred_np = np.array(pred_ecvl, copy=False) gt_np = np.array(gt_ecvl, copy=False) pred_np = pred_np.squeeze() gt_np = gt_np.squeeze() # Prediction summed in R channel image_np[:, :, -1] = np.where(pred_np == 255, pred_np, image_np[:, :, -1]) # Ground truth summed in G channel image_np[:, :, 1] = np.where(gt_np == 255, gt_np, image_np[:, :, 1]) n += 2 head, tail = os.path.splitext(os.path.basename(filename)) bname = "{}.png".format(head) filepath = os.path.join(args.out_dir, bname) ecvl.ImWrite(filepath, img_ecvl) print() mean_dice = evaluator.MeanMetric() if mean_dice > best_dice: print("Saving weights") eddl.save( net, "pneumothorax_segnetBN_adam_lr_0.0001_" "loss_ce_size_512_{}.bin".format(e + 1), "bin") best_dice = mean_dice print("Mean Dice Coefficient: {:.6g}".format(mean_dice))
epochs = 50 batch_size = 10 num_classes = 2 shape = (256, 256) data_yml_path = "../data/padchest.yml" ################## # MODEL CREATION # ################## in_ = eddl.Input((1, ) + shape) out = VGG16(in_, num_classes) net = eddl.Model([in_], [out]) eddl.build( net, eddl.adam(0.0001), # Optimizer ["cross_entropy"], # Losses ["categorical_accuracy"], # Metrics eddl.CS_GPU([1]) # Computing service ) # Show model architecture eddl.summary(net) ################ # ECVL DATASET # ################ train_augs = ecvl.SequentialAugmentationContainer([ ecvl.AugResizeDim(shape, ecvl.InterpolationType.nearest), ecvl.AugRotate([-10, 10]) ])
l = defblock(l, bn, 256, 4, initializer) l = defblock(l, bn, 512, 4, initializer) l = defblock(l, bn, 512, 4, initializer) l = eddl.Flatten(l) for i in range(2): l = initializer(eddl.Dense(l, 4096)) if (bn): l = eddl.BatchNormalization(l, 0.99, 0.001, True, "") l = eddl.ReLu(l) out = eddl.Softmax(initializer(eddl.Dense(l, num_classes))) net = eddl.Model([inp], [out]) eddl.plot(net, "model.pdf") eddl.build(net, eddl.adam(0.00001), ["soft_cross_entropy"], ["categorical_accuracy"], eddl.CS_GPU() if gpu else eddl.CS_CPU()) eddl.summary(net) x_train = Tensor.load("cifar_trX.bin") y_train = Tensor.load("cifar_trY.bin") x_train.div_(255) x_test = Tensor.load("cifar_tsX.bin") y_test = Tensor.load("cifar_tsY.bin") x_test.div_(255) from time import time import numpy as np
def main(args): eddl.download_eutrans() epochs = 1 if args.small else 5 ilength = 30 olength = 30 invs = 687 outvs = 514 embedding = 64 # Encoder in_ = eddl.Input([1]) # 1 word layer = in_ lE = eddl.RandomUniform( eddl.Embedding(layer, invs, 1, embedding, True), -0.05, 0.05 ) enc = eddl.LSTM(lE, 128, True) cps = eddl.GetStates(enc) # Decoder ldin = eddl.Input([outvs]) ld = eddl.ReduceArgMax(ldin, [0]) ld = eddl.RandomUniform( eddl.Embedding(ld, outvs, 1, embedding), -0.05, 0.05 ) layer = eddl.LSTM([ld, cps], 128) out = eddl.Softmax(eddl.Dense(layer, outvs)) eddl.setDecoder(ldin) net = eddl.Model([in_], [out]) # Build model eddl.build( net, eddl.adam(0.01), ["softmax_cross_entropy"], ["accuracy"], eddl.CS_GPU(mem=args.mem) if args.gpu else eddl.CS_CPU(mem=args.mem) ) eddl.summary(net) # Load dataset x_train = Tensor.load("eutrans_trX.bin") y_train = Tensor.load("eutrans_trY.bin") y_train = Tensor.onehot(y_train, outvs) # batch x timesteps x input_dim x_train.reshape_([x_train.shape[0], ilength, 1]) # batch x timesteps x ouput_dim y_train.reshape_([y_train.shape[0], olength, outvs]) x_test = Tensor.load("eutrans_tsX.bin") y_test = Tensor.load("eutrans_tsY.bin") y_test = Tensor.onehot(y_test, outvs) # batch x timesteps x input_dim x_test.reshape_([x_test.shape[0], ilength, 1]) # batch x timesteps x ouput_dim y_test.reshape_([y_test.shape[0], olength, outvs]) if args.small: sel = [f":{3 * args.batch_size}", ":", ":"] x_train = x_train.select(sel) y_train = y_train.select(sel) x_test = x_test.select(sel) y_test = y_test.select(sel) # Train model ybatch = Tensor([args.batch_size, olength, outvs]) eddl.next_batch([y_train], [ybatch]) for i in range(epochs): eddl.fit(net, [x_train], [y_train], args.batch_size, 1) print("All done")
def main(args): num_classes = 1 size = [192, 192] # size of images thresh = 0.5 if args.out_dir: os.makedirs(args.out_dir, exist_ok=True) in_ = eddl.Input([3, size[0], size[1]]) out = SegNet(in_, num_classes) out_sigm = eddl.Sigmoid(out) net = eddl.Model([in_], [out_sigm]) eddl.build(net, eddl.adam(0.0001), ["cross_entropy"], ["mean_squared_error"], eddl.CS_GPU([1]) if args.gpu else eddl.CS_CPU()) eddl.summary(net) eddl.setlogfile(net, "skin_lesion_segmentation") training_augs = ecvl.SequentialAugmentationContainer([ ecvl.AugResizeDim(size), ecvl.AugMirror(0.5), ecvl.AugFlip(0.5), ecvl.AugRotate([-180, 180]), ecvl.AugAdditivePoissonNoise([0, 10]), ecvl.AugGammaContrast([0.5, 1.5]), ecvl.AugGaussianBlur([0, 0.8]), ecvl.AugCoarseDropout([0, 0.3], [0.02, 0.05], 0.5) ]) validation_augs = ecvl.SequentialAugmentationContainer( [ecvl.AugResizeDim(size)]) dataset_augs = ecvl.DatasetAugmentations( [training_augs, validation_augs, None]) print("Reading dataset") d = ecvl.DLDataset(args.in_ds, args.batch_size, dataset_augs) x = Tensor([args.batch_size, d.n_channels_, size[0], size[1]]) y = Tensor([args.batch_size, d.n_channels_gt_, size[0], size[1]]) num_samples_train = len(d.GetSplit()) num_batches_train = num_samples_train // args.batch_size d.SetSplit(ecvl.SplitType.validation) num_samples_validation = len(d.GetSplit()) num_batches_validation = num_samples_validation // args.batch_size indices = list(range(args.batch_size)) evaluator = utils.Evaluator() print("Starting training") for e in range(args.epochs): print("Epoch {:d}/{:d} - Training".format(e + 1, args.epochs), flush=True) d.SetSplit(ecvl.SplitType.training) eddl.reset_loss(net) s = d.GetSplit() random.shuffle(s) d.split_.training_ = s d.ResetAllBatches() for b in range(num_batches_train): print("Epoch {:d}/{:d} (batch {:d}/{:d}) - ".format( e + 1, args.epochs, b + 1, num_batches_train), end="", flush=True) d.LoadBatch(x, y) x.div_(255.0) y.div_(255.0) tx, ty = [x], [y] eddl.train_batch(net, tx, ty, indices) eddl.print_loss(net, b) print() print("Saving weights") eddl.save(net, "isic_segmentation_checkpoint_epoch_%s.bin" % e, "bin") d.SetSplit(ecvl.SplitType.validation) evaluator.ResetEval() print("Epoch %d/%d - Evaluation" % (e + 1, args.epochs), flush=True) for b in range(num_batches_validation): n = 0 print("Epoch {:d}/{:d} (batch {:d}/{:d}) ".format( e + 1, args.epochs, b + 1, num_batches_validation), end="", flush=True) d.LoadBatch(x, y) x.div_(255.0) y.div_(255.0) eddl.forward(net, [x]) output = eddl.getOutput(out_sigm) for k in range(args.batch_size): img = output.select([str(k)]) gt = y.select([str(k)]) img_np = np.array(img, copy=False) gt_np = np.array(gt, copy=False) iou = evaluator.BinaryIoU(img_np, gt_np, thresh=thresh) print("- IoU: %.6g " % iou, end="", flush=True) if args.out_dir: # C++ BinaryIoU modifies image as a side effect img_np[img_np >= thresh] = 1 img_np[img_np < thresh] = 0 img_t = ecvl.TensorToView(img) img_t.colortype_ = ecvl.ColorType.GRAY img_t.channels_ = "xyc" img.mult_(255.) # orig_img orig_img = x.select([str(k)]) orig_img.mult_(255.) orig_img_t = ecvl.TensorToImage(orig_img) orig_img_t.colortype_ = ecvl.ColorType.BGR orig_img_t.channels_ = "xyc" tmp, labels = ecvl.Image.empty(), ecvl.Image.empty() ecvl.CopyImage(img_t, tmp, ecvl.DataType.uint8) ecvl.ConnectedComponentsLabeling(tmp, labels) ecvl.CopyImage(labels, tmp, ecvl.DataType.uint8) contours = ecvl.FindContours(tmp) ecvl.CopyImage(orig_img_t, tmp, ecvl.DataType.uint8) tmp_np = np.array(tmp, copy=False) for cseq in contours: for c in cseq: tmp_np[c[0], c[1], 0] = 0 tmp_np[c[0], c[1], 1] = 0 tmp_np[c[0], c[1], 2] = 255 filename = d.samples_[d.GetSplit()[n]].location_[0] head, tail = os.path.splitext(os.path.basename(filename)) bname = "%s.png" % head output_fn = os.path.join(args.out_dir, bname) ecvl.ImWrite(output_fn, tmp) if e == 0: gt_t = ecvl.TensorToView(gt) gt_t.colortype_ = ecvl.ColorType.GRAY gt_t.channels_ = "xyc" gt.mult_(255.) gt_filename = d.samples_[d.GetSplit()[n]].label_path_ gt_fn = os.path.join(args.out_dir, os.path.basename(gt_filename)) ecvl.ImWrite(gt_fn, gt_t) n += 1 print() print("MIoU: %.6g" % evaluator.MeanMetric())
def main(args): num_classes = 1 size = [512, 512] # size of images thresh = 0.5 if args.out_dir: os.makedirs(args.out_dir, exist_ok=True) in_ = eddl.Input([1, size[0], size[1]]) out = SegNetBN(in_, num_classes) out_sigm = eddl.Sigmoid(out) net = eddl.Model([in_], [out_sigm]) eddl.build( net, eddl.adam(0.0001), ["cross_entropy"], ["mean_squared_error"], eddl.CS_GPU([1]) if args.gpu else eddl.CS_CPU() ) eddl.summary(net) eddl.setlogfile(net, "pneumothorax_segmentation_inference") if not os.path.exists(args.ckpts): raise RuntimeError('Checkpoint "{}" not found'.format(args.ckpts)) eddl.load(net, args.ckpts, "bin") training_augs = ecvl.SequentialAugmentationContainer([ ecvl.AugResizeDim(size), ]) test_augs = ecvl.SequentialAugmentationContainer([ ecvl.AugResizeDim(size), ]) dataset_augs = ecvl.DatasetAugmentations([training_augs, None, test_augs]) print("Reading dataset") d = ecvl.DLDataset(args.in_ds, args.batch_size, dataset_augs, ecvl.ColorType.GRAY) x = Tensor([args.batch_size, d.n_channels_, size[0], size[1]]) print("Testing") d.SetSplit(ecvl.SplitType.test) num_samples_test = len(d.GetSplit()) num_batches_test = num_samples_test // args.batch_size evaluator = utils.Evaluator() evaluator.ResetEval() for b in range(num_batches_test): n = 0 print("Batch {:d}/{:d} ".format( b + 1, num_batches_test), end="", flush=True) d.LoadBatch(x) x.div_(255.0) eddl.forward(net, [x]) if args.out_dir: output = eddl.getOutput(out_sigm) for k in range(args.batch_size): img = output.select([str(k)]) img_I = ecvl.TensorToImage(img) img_I.colortype_ = ecvl.ColorType.GRAY img_I.channels_ = "xyc" ecvl.Threshold(img_I, img_I, thresh, 255) filename = d.samples_[d.GetSplit()[n]].location_[0] head, tail = os.path.splitext(os.path.basename(filename)) bname = "{}.png".format(head) output_fn = os.path.join(args.out_dir, bname) ecvl.ImWrite(output_fn, img_I) n += 1 print()
def main(args): eddl.download_flickr() epochs = 2 if args.small else 50 olength = 20 outvs = 2000 embdim = 32 # True: remove last layers and set new top = flatten # new input_size: [3, 256, 256] (from [224, 224, 3]) net = eddl.download_resnet18(True, [3, 256, 256]) lreshape = eddl.getLayer(net, "top") # create a new model from input output image_in = eddl.getLayer(net, "input") # Decoder ldecin = eddl.Input([outvs]) ldec = eddl.ReduceArgMax(ldecin, [0]) ldec = eddl.RandomUniform(eddl.Embedding(ldec, outvs, 1, embdim, True), -0.05, 0.05) ldec = eddl.Concat([ldec, lreshape]) layer = eddl.LSTM(ldec, 512, True) out = eddl.Softmax(eddl.Dense(layer, outvs)) eddl.setDecoder(ldecin) net = eddl.Model([image_in], [out]) # Build model eddl.build( net, eddl.adam(0.01), ["softmax_cross_entropy"], ["accuracy"], eddl.CS_GPU(mem=args.mem) if args.gpu else eddl.CS_CPU(mem=args.mem)) eddl.summary(net) # Load dataset x_train = Tensor.load("flickr_trX.bin", "bin") y_train = Tensor.load("flickr_trY.bin", "bin") if args.small: x_train = x_train.select([f"0:{2 * args.batch_size}", ":", ":", ":"]) y_train = y_train.select([f"0:{2 * args.batch_size}", ":"]) xtrain = Tensor.permute(x_train, [0, 3, 1, 2]) y_train = Tensor.onehot(y_train, outvs) # batch x timesteps x input_dim y_train.reshape_([y_train.shape[0], olength, outvs]) eddl.fit(net, [xtrain], [y_train], args.batch_size, epochs) eddl.save(net, "img2text.bin", "bin") print("\n === INFERENCE ===\n") # Get all the reshapes of the images. Only use the CNN timage = Tensor([x_train.shape[0], 512]) # images reshape cnn = eddl.Model([image_in], [lreshape]) eddl.build( cnn, eddl.adam(0.001), # not relevant ["mse"], # not relevant ["mse"], # not relevant eddl.CS_GPU(mem=args.mem) if args.gpu else eddl.CS_CPU(mem=args.mem)) eddl.summary(cnn) # forward images xbatch = Tensor([args.batch_size, 3, 256, 256]) # numbatches = x_train.shape[0] / args.batch_size j = 0 eddl.next_batch([x_train], [xbatch]) eddl.forward(cnn, [xbatch]) ybatch = eddl.getOutput(lreshape) sample = str(j * args.batch_size) + ":" + str((j + 1) * args.batch_size) timage.set_select([sample, ":"], ybatch) # Create Decoder non recurrent for n-best ldecin = eddl.Input([outvs]) image = eddl.Input([512]) lstate = eddl.States([2, 512]) ldec = eddl.ReduceArgMax(ldecin, [0]) ldec = eddl.RandomUniform(eddl.Embedding(ldec, outvs, 1, embdim), -0.05, 0.05) ldec = eddl.Concat([ldec, image]) lstm = eddl.LSTM([ldec, lstate], 512, True) lstm.isrecurrent = False # Important out = eddl.Softmax(eddl.Dense(lstm, outvs)) decoder = eddl.Model([ldecin, image, lstate], [out]) eddl.build( decoder, eddl.adam(0.001), # not relevant ["softmax_cross_entropy"], # not relevant ["accuracy"], # not relevant eddl.CS_GPU(mem=args.mem) if args.gpu else eddl.CS_CPU(mem=args.mem)) eddl.summary(decoder) # Copy params from trained net eddl.copyParam(eddl.getLayer(net, "LSTM1"), eddl.getLayer(decoder, "LSTM2")) eddl.copyParam(eddl.getLayer(net, "dense1"), eddl.getLayer(decoder, "dense2")) eddl.copyParam(eddl.getLayer(net, "embedding1"), eddl.getLayer(decoder, "embedding2")) # N-best for sample s s = 1 if args.small else 100 # sample 100 # three input tensors with batch_size = 1 (one sentence) treshape = timage.select([str(s), ":"]) text = y_train.select([str(s), ":", ":"]) # 1 x olength x outvs for j in range(olength): print(f"Word: {j}") word = None if j == 0: word = Tensor.zeros([1, outvs]) else: word = text.select(["0", str(j - 1), ":"]) word.reshape_([1, outvs]) # batch = 1 treshape.reshape_([1, 512]) # batch = 1 state = Tensor.zeros([1, 2, 512]) # batch = 1 input_ = [word, treshape, state] eddl.forward(decoder, input_) # outword = eddl.getOutput(out) vstates = eddl.getStates(lstm) for i in range(len(vstates)): vstates[i].reshape_([1, 1, 512]) state.set_select([":", str(i), ":"], vstates[i]) print("All done")
def main(args): batch_size = args.batch_size image_size = args.size, args.size thresh = 0.5 if args.weights: os.makedirs(args.weights, exist_ok=True) training_augs = ecvl.SequentialAugmentationContainer([ ecvl.AugResizeDim(image_size, ecvl.InterpolationType.cubic, gt_interp=ecvl.InterpolationType.nearest), ecvl.AugMirror(.5), ecvl.AugFlip(.5), ecvl.AugRotate([-180, 180]), ecvl.AugAdditivePoissonNoise([0, 10]), ecvl.AugGammaContrast([0.5, 1.5]), ecvl.AugGaussianBlur([0, 0.8]), ecvl.AugCoarseDropout([0, 0.03], [0.02, 0.05], 0.25), ecvl.AugToFloat32(255, divisor_gt=255), ecvl.AugNormalize([0.6681, 0.5301, 0.5247], [0.1337, 0.1480, 0.1595]), # isic stats ]) validation_test_augs = ecvl.SequentialAugmentationContainer([ ecvl.AugResizeDim(image_size, ecvl.InterpolationType.cubic, gt_interp=ecvl.InterpolationType.nearest), ecvl.AugToFloat32(255, divisor_gt=255), ecvl.AugNormalize([0.6681, 0.5301, 0.5247], [0.1337, 0.1480, 0.1595]), # isic stats ]) dataset_augs = ecvl.DatasetAugmentations( [training_augs, validation_test_augs, validation_test_augs]) print('Reading dataset') d = ecvl.DLDataset(args.in_ds, args.batch_size, dataset_augs, ctype=ecvl.ColorType.RGB) num_classes = len(d.classes_) or d.n_channels_gt_ size = d.n_channels_, args.size, args.size if args.ckpts: net = eddl.import_net_from_onnx_file(args.ckpts, size) else: in_ = eddl.Input(size) out = Unet(in_, num_classes) out_sigm = eddl.Sigmoid(out) net = eddl.Model([in_], [out_sigm]) # model_path = utils.DownloadModel(segmentation_zoo[args.model]['url'], f'{args.model}.onnx', 'model_onnx') # net = eddl.import_net_from_onnx_file(model_path, size) # eddl.removeLayer(net, segmentation_zoo[args.model]['to_remove']) # top = eddl.getLayer(net, segmentation_zoo[args.model]['top']) # # out = eddl.Sigmoid(eddl.Conv(top, num_classes, [3, 3], name='last_layer')) # data_input = eddl.getLayer(net, segmentation_zoo[args.model]['input']) # input of the onnx # net = eddl.Model([data_input], [out]) loss_name = 'binary_cross_entropy' metric_name = 'mean_squared_error' eddl.build( net, eddl.adam(args.learning_rate), [loss_name], [metric_name], eddl.CS_GPU(args.gpu, mem="low_mem") if args.gpu else eddl.CS_CPU(), True) out = eddl.getOut(net)[0] # if not args.ckpts: # eddl.initializeLayer(net, "last_layer") eddl.summary(net) eddl.setlogfile(net, 'skin_lesion_segmentation') x = Tensor([args.batch_size, *size]) y = Tensor([args.batch_size, d.n_channels_gt_, size[1], size[2]]) miou = 0. if args.train: num_samples_train = len(d.GetSplit()) num_batches_train = num_samples_train // args.batch_size num_samples_val = len(d.GetSplit(ecvl.SplitType.validation)) num_batches_val = num_samples_val // args.batch_size evaluator = utils.Evaluator() print('Starting training') for e in range(args.epochs): d.SetSplit(ecvl.SplitType.training) eddl.reset_loss(net) s = d.GetSplit() random.shuffle(s) d.split_.training_ = s d.ResetAllBatches() for b in range(num_batches_train): d.LoadBatch(x, y) # x_ = x.select(["0"]) # x_.normalize_(0, 1) # x_.mult_(255.) # x_.save(f'images/train_{e}_{b}.png') # # y_ = y.select(["0"]) # # y_.mult_(255.) # y_.save(f'images/train_gt_{e}_{b}.png') eddl.train_batch(net, [x], [y]) losses = eddl.get_losses(net) metrics = eddl.get_metrics(net) print( f'Train - epoch [{e + 1}/{args.epochs}] - batch [{b + 1}/{num_batches_train}]' f' - {loss_name}={losses[0]:.3f} - {metric_name}={metrics[0]:.3f}', flush=True) d.SetSplit(ecvl.SplitType.validation) evaluator.ResetEval() eddl.reset_loss(net) for b in range(num_batches_val): n = 0 print( f'Validation - epoch [{e + 1}/{args.epochs}] - batch [{b + 1}/{num_batches_val}]' ) d.LoadBatch(x, y) eddl.forward(net, [x]) output = eddl.getOutput(out) for bs in range(args.batch_size): img = output.select([str(bs)]) gt = y.select([str(bs)]) img_np = np.array(img, copy=False) gt_np = np.array(gt, copy=False) iou = evaluator.BinaryIoU(img_np, gt_np, thresh=thresh) print(f' - IoU: {iou:.3f}', end="", flush=True) if args.out_dir: # C++ BinaryIoU modifies image as a side effect img_np[img_np >= thresh] = 1 img_np[img_np < thresh] = 0 img_t = ecvl.TensorToView(img) img_t.colortype_ = ecvl.ColorType.GRAY img_t.channels_ = "xyc" img.mult_(255.) # orig_img orig_img = x.select([str(bs)]) orig_img.mult_(255.) orig_img_t = ecvl.TensorToImage(orig_img) orig_img_t.colortype_ = ecvl.ColorType.BGR orig_img_t.channels_ = "xyc" tmp, labels = ecvl.Image.empty(), ecvl.Image.empty() ecvl.CopyImage(img_t, tmp, ecvl.DataType.uint8) ecvl.ConnectedComponentsLabeling(tmp, labels) ecvl.CopyImage(labels, tmp, ecvl.DataType.uint8) contours = ecvl.FindContours(tmp) ecvl.CopyImage(orig_img_t, tmp, ecvl.DataType.uint8) tmp_np = np.array(tmp, copy=False) for cseq in contours: for c in cseq: tmp_np[c[0], c[1], 0] = 0 tmp_np[c[0], c[1], 1] = 0 tmp_np[c[0], c[1], 2] = 255 filename = d.samples_[d.GetSplit()[n]].location_[0] head, tail = os.path.splitext( os.path.basename(filename)) bname = "%s.png" % head output_fn = os.path.join(args.out_dir, bname) ecvl.ImWrite(output_fn, tmp) if e == 0: gt_t = ecvl.TensorToView(gt) gt_t.colortype_ = ecvl.ColorType.GRAY gt_t.channels_ = "xyc" gt.mult_(255.) gt_filename = d.samples_[d.GetSplit() [n]].label_path_ gt_fn = os.path.join(args.out_dir, os.path.basename(gt_filename)) ecvl.ImWrite(gt_fn, gt_t) n += 1 print() last_miou = evaluator.MIoU() print( f'Validation - epoch [{e + 1}/{args.epochs}] - Total MIoU: {last_miou:.3f}' ) if last_miou > miou: miou = last_miou eddl.save_net_to_onnx_file( net, os.path.join(args.weights, f'isic_segm_{args.model}_epoch_{e + 1}.onnx')) print('Weights saved') elif args.test: evaluator = utils.Evaluator() evaluator.ResetEval() d.SetSplit(ecvl.SplitType.test) num_samples_test = len(d.GetSplit()) num_batches_test = num_samples_test // batch_size for b in range(num_batches_test): n = 0 print(f'Test - batch [{b + 1}/{num_batches_test}]') d.LoadBatch(x, y) eddl.forward(net, [x]) output = eddl.getOutput(out) for bs in range(args.batch_size): img = output.select([str(bs)]) gt = y.select([str(bs)]) img_np, gt_np = np.array(img, copy=False), np.array(gt, copy=False) iou = evaluator.BinaryIoU(img_np, gt_np, thresh=thresh) print(f' - IoU: {iou:.3f}', end="", flush=True) if args.out_dir: # C++ BinaryIoU modifies image as a side effect img_np[img_np >= thresh] = 1 img_np[img_np < thresh] = 0 img_t = ecvl.TensorToView(img) img_t.colortype_ = ecvl.ColorType.GRAY img_t.channels_ = "xyc" img.mult_(255.) # orig_img orig_img = x.select([str(bs)]) orig_img.mult_(255.) orig_img_t = ecvl.TensorToImage(orig_img) orig_img_t.colortype_ = ecvl.ColorType.BGR orig_img_t.channels_ = "xyc" tmp, labels = ecvl.Image.empty(), ecvl.Image.empty() ecvl.CopyImage(img_t, tmp, ecvl.DataType.uint8) ecvl.ConnectedComponentsLabeling(tmp, labels) ecvl.CopyImage(labels, tmp, ecvl.DataType.uint8) contours = ecvl.FindContours(tmp) ecvl.CopyImage(orig_img_t, tmp, ecvl.DataType.uint8) tmp_np = np.array(tmp, copy=False) for cseq in contours: for c in cseq: tmp_np[c[0], c[1], 0] = 0 tmp_np[c[0], c[1], 1] = 0 tmp_np[c[0], c[1], 2] = 255 filename = d.samples_[d.GetSplit()[n]].location_[0] head, tail = os.path.splitext(os.path.basename(filename)) bname = "%s.png" % head output_fn = os.path.join(args.out_dir, bname) ecvl.ImWrite(output_fn, tmp) gt_t = ecvl.TensorToView(gt) gt_t.colortype_ = ecvl.ColorType.GRAY gt_t.channels_ = "xyc" gt.mult_(255.) gt_filename = d.samples_[d.GetSplit()[n]].label_path_ gt_fn = os.path.join(args.out_dir, os.path.basename(gt_filename)) ecvl.ImWrite(gt_fn, gt_t) n += 1 miou = evaluator.MIoU() print(f'Test - Total MIoU: {miou:.3f}')
def main(args): batch_size = args.batch_size image_size = args.size, args.size if args.weights: os.makedirs(args.weights, exist_ok=True) training_augs = ecvl.SequentialAugmentationContainer([ ecvl.AugResizeDim(image_size, ecvl.InterpolationType.cubic), ecvl.AugMirror(.5), ecvl.AugFlip(.5), ecvl.AugRotate([-180, 180]), ecvl.AugAdditivePoissonNoise([0, 10]), ecvl.AugGammaContrast([0.5, 1.5]), ecvl.AugGaussianBlur([0, 0.8]), ecvl.AugCoarseDropout([0, 0.03], [0.02, 0.05], 0.25), ecvl.AugToFloat32(255), ]) validation_test_augs = ecvl.SequentialAugmentationContainer([ ecvl.AugResizeDim(image_size), ecvl.AugToFloat32(255), ]) dataset_augs = ecvl.DatasetAugmentations( [training_augs, validation_test_augs, validation_test_augs]) print('Reading dataset') d = ecvl.DLDataset(args.in_ds, args.batch_size, dataset_augs, ctype=ecvl.ColorType.RGB) num_classes = len(d.classes_) size = d.n_channels_, args.size, args.size if args.ckpts: net = eddl.import_net_from_onnx_file(args.ckpts, size) else: model_path = utils.DownloadModel(classification_zoo[args.model]['url'], f'{args.model}.onnx', 'model_onnx') net = eddl.import_net_from_onnx_file(model_path, size) eddl.removeLayer(net, classification_zoo[args.model] ['to_remove']) # remove last Linear of resnet top = eddl.getLayer( net, classification_zoo[args.model]['top']) # get flatten of resnet out = eddl.Softmax(eddl.Dense(top, num_classes, True, 'classifier')) # true is for the bias data_input = eddl.getLayer( net, classification_zoo[args.model]['input']) # input of the onnx net = eddl.Model([data_input], [out]) eddl.build( net, eddl.adam(args.learning_rate), ['softmax_cross_entropy'], ['accuracy'], eddl.CS_GPU(args.gpu, mem="low_mem") if args.gpu else eddl.CS_CPU(), False) out = eddl.getOut(net)[0] if not args.ckpts: eddl.initializeLayer(net, "classifier") eddl.summary(net) eddl.setlogfile(net, 'skin_lesion_classification') x = Tensor([batch_size, *size]) y = Tensor([batch_size, num_classes]) metric_fn = eddl.getMetric('accuracy') best_accuracy = 0. if args.train: num_samples_train = len(d.GetSplit()) num_batches_train = num_samples_train // args.batch_size num_samples_val = len(d.GetSplit(ecvl.SplitType.validation)) num_batches_val = num_samples_val // args.batch_size print('Starting training') for e in range(args.epochs): if args.out_dir: current_path = os.path.join(args.out_dir, f'Epoch_{e}') for c in d.classes_: c_dir = os.path.join(current_path, c) os.makedirs(c_dir, exist_ok=True) d.SetSplit(ecvl.SplitType.training) eddl.reset_loss(net) s = d.GetSplit() random.shuffle(s) d.split_.training_ = s d.ResetAllBatches() for b in range(num_batches_train): d.LoadBatch(x, y) eddl.train_batch(net, [x], [y]) losses = eddl.get_losses(net) metrics = eddl.get_metrics(net) print( f'Train - epoch [{e + 1}/{args.epochs}] - batch [{b + 1}/{num_batches_train}]' f' - loss={losses[0]:.3f} - accuracy={metrics[0]:.3f}', flush=True) d.SetSplit(ecvl.SplitType.validation) values = np.zeros(num_batches_val) eddl.reset_loss(net) for b in range(num_batches_val): n = 0 d.LoadBatch(x, y) eddl.forward(net, [x]) output = eddl.getOutput(out) value = metric_fn.value(y, output) values[b] = value if args.out_dir: for k in range(args.batch_size): result = output.select([str(k)]) target = y.select([str(k)]) result_a = np.array(result, copy=False) target_a = np.array(target, copy=False) classe = np.argmax(result_a).item() gt_class = np.argmax(target_a).item() single_image = x.select([str(k)]) img_t = ecvl.TensorToView(single_image) img_t.colortype_ = ecvl.ColorType.BGR single_image.mult_(255.) filename = d.samples_[d.GetSplit()[n]].location_[0] head, tail = os.path.splitext( os.path.basename(filename)) bname = '{}_gt_class_{}.png'.format(head, gt_class) cur_path = os.path.join(current_path, d.classes_[classe], bname) ecvl.ImWrite(cur_path, img_t) n += 1 print( f'Validation - epoch [{e + 1}/{args.epochs}] - batch [{b + 1}/{num_batches_val}] -' f' accuracy={np.mean(values[:b + 1] / batch_size):.3f}') last_accuracy = np.mean(values / batch_size) print( f'Validation - epoch [{e + 1}/{args.epochs}] - total accuracy={last_accuracy:.3f}' ) if last_accuracy > best_accuracy: best_accuracy = last_accuracy print('Saving weights') eddl.save_net_to_onnx_file( net, f'isic_classification_{args.model}_epoch_{e + 1}.onnx') elif args.test: d.SetSplit(ecvl.SplitType.test) num_samples_test = len(d.GetSplit()) num_batches_test = num_samples_test // batch_size values = np.zeros(num_batches_test) eddl.reset_loss(net) for b in range(num_batches_test): d.LoadBatch(x, y) eddl.forward(net, [x]) output = eddl.getOutput(out) value = metric_fn.value(y, output) values[b] = value if args.out_dir: n = 0 for k in range(args.batch_size): result = output.select([str(k)]) target = y.select([str(k)]) result_a = np.array(result, copy=False) target_a = np.array(target, copy=False) classe = np.argmax(result_a).item() gt_class = np.argmax(target_a).item() single_image = x.select([str(k)]) img_t = ecvl.TensorToView(single_image) img_t.colortype_ = ecvl.ColorType.BGR single_image.mult_(255.) filename = d.samples_[d.GetSplit()[n]].location_[0] head, tail = os.path.splitext(os.path.basename(filename)) bname = "%s_gt_class_%s.png" % (head, gt_class) cur_path = os.path.join(args.out_dir, d.classes_[classe], bname) ecvl.ImWrite(cur_path, img_t) n += 1 print( f'Test - batch [{b + 1}/{num_batches_test}] - accuracy={np.mean(values[:b + 1] / batch_size):.3f}' ) print(f'Test - total accuracy={np.mean(values / batch_size):.3f}')
def main(args): freeze_epochs = 2 unfreeze_epochs = 5 num_classes = 10 # 10 labels in cifar10 eddl.download_cifar10() eddl.download_model("resnet18.onnx", "re7jodd12srksd7") net = eddl.import_net_from_onnx_file("resnet18.onnx", [3, 32, 32], DEV_CPU) names = [_.name for _ in net.layers] # Remove dense output layer eddl.removeLayer(net, "resnetv15_dense0_fwd") # Get last layer to connect the new dense layer = eddl.getLayer(net, "flatten_170") out = eddl.Softmax(eddl.Dense(layer, num_classes, True, "new_dense")) # Get input layer in_ = eddl.getLayer(net, "data") # Create a new model net = eddl.Model([in_], [out]) eddl.build( net, eddl.adam(0.0001), ["softmax_cross_entropy"], ["categorical_accuracy"], eddl.CS_GPU(mem=args.mem) if args.gpu else eddl.CS_CPU(mem=args.mem), False # do not initialize weights to random values ) eddl.summary(net) # Force initialization of new layers eddl.initializeLayer(net, "new_dense") x_train = Tensor.load("cifar_trX.bin") y_train = Tensor.load("cifar_trY.bin") x_test = Tensor.load("cifar_tsX.bin") y_test = Tensor.load("cifar_tsY.bin") if args.small: sel = [f":{2 * args.batch_size}"] x_train = x_train.select(sel) y_train = y_train.select(sel) x_test = x_test.select(sel) y_test = y_test.select(sel) x_train.div_(255.0) x_test.div_(255.0) # Freeze pretrained weights for n in names: eddl.setTrainable(net, n, False) # Train new layers eddl.fit(net, [x_train], [y_train], args.batch_size, freeze_epochs) # Unfreeze weights for n in names: eddl.setTrainable(net, n, True) # Train all layers eddl.fit(net, [x_train], [y_train], args.batch_size, unfreeze_epochs) # Evaluate eddl.evaluate(net, [x_test], [y_test], args.batch_size) print("All done")
def main(args): num_classes = 1 size = [192, 192] # size of images thresh = 0.5 if args.out_dir: os.makedirs(args.out_dir, exist_ok=True) in_ = eddl.Input([3, size[0], size[1]]) out = SegNet(in_, num_classes) out_sigm = eddl.Sigmoid(out) net = eddl.Model([in_], [out_sigm]) eddl.build(net, eddl.adam(0.0001), ["cross_entropy"], ["mean_squared_error"], eddl.CS_GPU([1]) if args.gpu else eddl.CS_CPU()) eddl.summary(net) eddl.setlogfile(net, "skin_lesion_segmentation_inference") if not os.path.exists(args.ckpts): raise RuntimeError('Checkpoint "{}" not found'.format(args.ckpts)) eddl.load(net, args.ckpts, "bin") training_augs = ecvl.SequentialAugmentationContainer([ ecvl.AugResizeDim(size), ]) test_augs = ecvl.SequentialAugmentationContainer([ ecvl.AugResizeDim(size), ]) dataset_augs = ecvl.DatasetAugmentations([training_augs, None, test_augs]) print("Reading dataset") d = ecvl.DLDataset(args.in_ds, args.batch_size, dataset_augs) x = Tensor([args.batch_size, d.n_channels_, size[0], size[1]]) y = Tensor([args.batch_size, d.n_channels_gt_, size[0], size[1]]) print("Testing") d.SetSplit(ecvl.SplitType.test) num_samples_test = len(d.GetSplit()) num_batches_test = num_samples_test // args.batch_size evaluator = utils.Evaluator() evaluator.ResetEval() for b in range(num_batches_test): n = 0 print("Batch {:d}/{:d} ".format(b + 1, num_batches_test), end="", flush=True) d.LoadBatch(x, y) x.div_(255.0) y.div_(255.0) eddl.forward(net, [x]) output = eddl.getOutput(out_sigm) for k in range(args.batch_size): img = output.select([str(k)]) gt = y.select([str(k)]) img_np, gt_np = np.array(img, copy=False), np.array(gt, copy=False) iou = evaluator.BinaryIoU(img_np, gt_np, thresh=thresh) print("- IoU: %.6g " % iou, end="", flush=True) if args.out_dir: # C++ BinaryIoU modifies image as a side effect img_np[img_np >= thresh] = 1 img_np[img_np < thresh] = 0 img_t = ecvl.TensorToView(img) img_t.colortype_ = ecvl.ColorType.GRAY img_t.channels_ = "xyc" img.mult_(255.) # orig_img orig_img = x.select([str(k)]) orig_img.mult_(255.) orig_img_t = ecvl.TensorToImage(orig_img) orig_img_t.colortype_ = ecvl.ColorType.BGR orig_img_t.channels_ = "xyc" tmp, labels = ecvl.Image.empty(), ecvl.Image.empty() ecvl.CopyImage(img_t, tmp, ecvl.DataType.uint8) ecvl.ConnectedComponentsLabeling(tmp, labels) ecvl.CopyImage(labels, tmp, ecvl.DataType.uint8) contours = ecvl.FindContours(tmp) ecvl.CopyImage(orig_img_t, tmp, ecvl.DataType.uint8) tmp_np = np.array(tmp, copy=False) for cseq in contours: for c in cseq: tmp_np[c[0], c[1], 0] = 0 tmp_np[c[0], c[1], 1] = 0 tmp_np[c[0], c[1], 2] = 255 filename = d.samples_[d.GetSplit()[n]].location_[0] head, tail = os.path.splitext(os.path.basename(filename)) bname = "%s.png" % head output_fn = os.path.join(args.out_dir, bname) ecvl.ImWrite(output_fn, tmp) gt_t = ecvl.TensorToView(gt) gt_t.colortype_ = ecvl.ColorType.GRAY gt_t.channels_ = "xyc" gt.mult_(255.) gt_filename = d.samples_[d.GetSplit()[n]].label_path_ gt_fn = os.path.join(args.out_dir, os.path.basename(gt_filename)) ecvl.ImWrite(gt_fn, gt_t) n += 1 print() print("MIoU: %.6g" % evaluator.MeanMetric())
l = defblock(l, bn, 512, 3) l = defblock(l, bn, 512, 3) l = eddl.Flatten(l) for i in range(2): l = eddl.GlorotUniform(eddl.Dense(l, 4096)) if(bn): l = eddl.BatchNormalization(l, 0.99, 0.001, True, "") l = eddl.ReLu(l) out = eddl.Softmax(eddl.GlorotUniform(eddl.Dense(l, num_classes))) net = eddl.Model([inp], [out]) eddl.plot(net, "model.pdf") eddl.build(net, eddl.adam(0.0001), ["soft_cross_entropy"], ["categorical_accuracy"], eddl.CS_GPU() if gpu else eddl.CS_CPU() ) eddl.summary(net) x_train = Tensor.load("cifar_trX.bin") y_train = Tensor.load("cifar_trY.bin") x_train.div_(255) x_test = Tensor.load("cifar_tsX.bin") y_test = Tensor.load("cifar_tsY.bin") x_test.div_(255)
def main(args): eddl.download_cifar10() num_classes = 10 in_ = eddl.Input([3, 32, 32]) layer = in_ layer = eddl.RandomHorizontalFlip(layer) layer = eddl.RandomCropScale(layer, [0.8, 1.0]) layer = eddl.RandomCutout(layer, [0.1, 0.5], [0.1, 0.5]) layer = eddl.MaxPool(eddl.ReLu(eddl.BatchNormalization( eddl.HeUniform(eddl.Conv(layer, 32, [3, 3], [1, 1], "same", False)), True)), [2, 2]) layer = eddl.MaxPool(eddl.ReLu(eddl.BatchNormalization( eddl.HeUniform(eddl.Conv(layer, 64, [3, 3], [1, 1], "same", False)), True)), [2, 2]) layer = eddl.MaxPool(eddl.ReLu(eddl.BatchNormalization( eddl.HeUniform(eddl.Conv(layer, 128, [3, 3], [1, 1], "same", False)), True)), [2, 2]) layer = eddl.MaxPool(eddl.ReLu(eddl.BatchNormalization( eddl.HeUniform(eddl.Conv(layer, 256, [3, 3], [1, 1], "same", False)), True)), [2, 2]) layer = eddl.Reshape(layer, [-1]) layer = eddl.Activation(eddl.BatchNormalization( eddl.Dense(layer, 128), True ), "relu") out = eddl.Softmax(eddl.BatchNormalization( eddl.Dense(layer, num_classes), True )) net = eddl.Model([in_], [out]) eddl.build( net, eddl.adam(0.001), ["softmax_cross_entropy"], ["categorical_accuracy"], eddl.CS_GPU(mem=args.mem) if args.gpu else eddl.CS_CPU(mem=args.mem) ) eddl.summary(net) eddl.plot(net, "model.pdf") x_train = Tensor.load("cifar_trX.bin") y_train = Tensor.load("cifar_trY.bin") x_train.div_(255.0) x_test = Tensor.load("cifar_tsX.bin") y_test = Tensor.load("cifar_tsY.bin") x_test.div_(255.0) if args.small: x_train = x_train.select([":5000"]) y_train = y_train.select([":5000"]) x_test = x_test.select([":1000"]) y_test = y_test.select([":1000"]) for i in range(args.epochs): eddl.fit(net, [x_train], [y_train], args.batch_size, 1) eddl.evaluate(net, [x_test], [y_test], bs=args.batch_size) eddl.setlr(net, [0.0001]) for i in range(args.epochs): eddl.fit(net, [x_train], [y_train], args.batch_size, 1) eddl.evaluate(net, [x_test], [y_test], bs=args.batch_size) print("All done")
def main(args): eddl.download_drive() in_1 = eddl.Input([3, 584, 584]) in_2 = eddl.Input([1, 584, 584]) layer = eddl.Concat([in_1, in_2]) layer = eddl.RandomCropScale(layer, [0.9, 1.0]) layer = eddl.CenteredCrop(layer, [512, 512]) img = eddl.Select(layer, ["0:3"]) mask = eddl.Select(layer, ["3"]) # DA net danet = eddl.Model([in_1, in_2], []) eddl.build(danet) if args.gpu: eddl.toGPU(danet, mem="low_mem") eddl.summary(danet) # SegNet in_ = eddl.Input([3, 512, 512]) out = eddl.Sigmoid(UNetWithPadding(in_)) segnet = eddl.Model([in_], [out]) eddl.build( segnet, eddl.adam(0.00001), # Optimizer ["mse"], # Losses ["mse"], # Metrics eddl.CS_GPU(mem=args.mem) if args.gpu else eddl.CS_CPU(mem=args.mem) ) eddl.summary(segnet) print("Reading training data") # x_train_f = Tensor.fromarray(np.load("drive_trX.npy").astype(np.float32)) x_train_f = Tensor.load("drive_trX.bin") x_train = x_train_f.permute([0, 3, 1, 2]) x_train.info() x_train.div_(255.0) print("Reading test data") # y_train = Tensor.fromarray(np.load("drive_trY.npy").astype(np.float32)) y_train = Tensor.load("drive_trY.bin") y_train.info() y_train.reshape_([20, 1, 584, 584]) y_train.div_(255.0) xbatch = Tensor([args.batch_size, 3, 584, 584]) ybatch = Tensor([args.batch_size, 1, 584, 584]) print("Starting training") for i in range(args.epochs): print("\nEpoch %d/%d" % (i + 1, args.epochs)) eddl.reset_loss(segnet) for j in range(args.num_batches): eddl.next_batch([x_train, y_train], [xbatch, ybatch]) # DA net eddl.forward(danet, [xbatch, ybatch]) xbatch_da = eddl.getOutput(img) ybatch_da = eddl.getOutput(mask) # SegNet eddl.train_batch(segnet, [xbatch_da], [ybatch_da]) eddl.print_loss(segnet, j) if i == args.epochs - 1: yout = eddl.getOutput(out).select(["0"]) yout.save("./out_%d.jpg" % j) print() print("All done")