For an assignment I want to resize a .jpg image with a python code, but without using the pil.image.resize() function or another similar function. I want to write the code myself but I can't figure out how. The image is RGB. I have found this can be solved by nearest neighbor interpolation (as well as other methods but this one is fine for my specific assignment). The height and the width should both be able to be made bigger or smaller. So far I only have this:

import numpy as np
import scipy as sc
import matplotlib as plt
import math
import PIL
from PIL import Imageimg = np.array(Image.open("foto1.jpg"))          height = img.shape[0]
width = img.shape[1]
dim = img.shape[2]  new_h = int(input("New height: "))  
new_w = int(input("New width: "))  imgR = img[:,:,0] #red pixels
imgG = img[:,:,1] #green pixels
imgB = img[:,:,2] #blue pixelsnewR = np.empty([new_h, new_w])             
newG = np.empty([new_h, new_w])            
newB = np.empty([new_h, new_w]) 

So now all three colours have a new array of the right dimensions. Unfortunately on the web I can only find people who use resize() functions... Does anyone know?

Thank in advance!

The key to doing any image transformation like resizing is to have a mapping from output coordinates to input coordinates. Then you can simply iterate over the entire output and grab a pixel from the input. Nearest neighbor makes this particularly easy, because there's never a need to interpolate a pixel that doesn't lie exactly on integer coordinates - you simply round the coordinates to the nearest integer.

for new_y in range(new_h):old_y = int(round(new_y * (new_h - 1) / (height - 1)))if old_y < 0: old_y = 0if old_y >= height: old_y = height - 1for new_x in range(new_w):old_x = int(round(new_x * (new_w - 1) / (width - 1)))if old_x < 0: old_x = 0if old_x >= width: old_x = width - 1newR[new_y,new_x] = imgR[old_y,old_x]newG[new_y,new_x] = imgG[old_y,old_x]newB[new_y,new_x] = imgB[old_y,old_x]