blob: b116e38c9fedacd08a34a5c906ab120421aafe30 [file] [log] [blame]
#-------------------------------------------------------------
#
# Licensed to the Apache Software Foundation (ASF) under one
# or more contributor license agreements. See the NOTICE file
# distributed with this work for additional information
# regarding copyright ownership. The ASF licenses this file
# to you under the Apache License, Version 2.0 (the
# "License"); you may not use this file except in compliance
# with the License. You may obtain a copy of the License at
#
# http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing,
# software distributed under the License is distributed on an
# "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
# KIND, either express or implied. See the License for the
# specific language governing permissions and limitations
# under the License.
#
#-------------------------------------------------------------
args = commandArgs(TRUE)
options(digits=22)
library("Matrix")
image_crop = function(img_in, w, h, x_offset, y_offset) {
orig_w = ncol(img_in)
orig_h = nrow(img_in)
start_h = (ceiling((orig_h - h) / 2)) + y_offset
end_h = (start_h + h - 1)
start_w = (ceiling((orig_w - w) / 2)) + x_offset
end_w = (start_w + w - 1)
if((start_h < 0) | (end_h > orig_h) | (start_w < 0) | (end_w > orig_w)) {
print("Offset out of bounds! Returning input.")
img_out = img_in
}
else {
mask = matrix(0, orig_h, orig_w)
temp_mask = matrix(1, h , w )
mask[start_h:end_h, start_w:end_w] = temp_mask
mask = matrix(mask, 1, orig_w * orig_h)
img_out = input[start_h:end_h, start_w:end_w]
}
}
input = as.matrix(readMM(paste(args[1], "A.mtx", sep="")))
size = as.double(args[3])
input = matrix(input, as.integer(args[4]), as.integer(args[5]))
new_w = floor(ncol(input) * size)
new_h = floor(nrow(input) * size)
print(paste("New w/h=",new_w,"/",new_h))
crop1 = image_crop(input, new_w, new_h, as.integer(args[6]), as.integer(args[7]));
writeMM(as(crop1, "CsparseMatrix"), paste(args[2], "B", sep=""))