diff --git a/doc/Projects/2023/Project1/DataFiles/SRTM_data_Norway_1.tif b/doc/Projects/2023/Project1/DataFiles/SRTM_data_Norway_1.tif
new file mode 100644
index 000000000..be81ff728
Binary files /dev/null and b/doc/Projects/2023/Project1/DataFiles/SRTM_data_Norway_1.tif differ
diff --git a/doc/Projects/2023/Project1/DataFiles/SRTM_data_Norway_2.tif b/doc/Projects/2023/Project1/DataFiles/SRTM_data_Norway_2.tif
new file mode 100644
index 000000000..90579e59a
Binary files /dev/null and b/doc/Projects/2023/Project1/DataFiles/SRTM_data_Norway_2.tif differ
diff --git a/doc/Projects/2023/Project1/DataFiles/test.py b/doc/Projects/2023/Project1/DataFiles/test.py
new file mode 100644
index 000000000..7c09c39e0
--- /dev/null
+++ b/doc/Projects/2023/Project1/DataFiles/test.py
@@ -0,0 +1,30 @@
+import numpy as np
+from imageio import imread
+import matplotlib.pyplot as plt
+from mpl_toolkits.mplot3d import Axes3D
+from matplotlib import cm
+import numpy as np
+from numpy.random import normal, uniform
+
+# Load the terrain
+terrain = imread('SRTM_data_Norway_1.tif')
+
+N = 1000
+m = 5 # polynomial order
+terrain = terrain[:N,:N]
+# Creates mesh of image pixels
+x = np.linspace(0,1, np.shape(terrain)[0])
+y = np.linspace(0,1, np.shape(terrain)[1])
+x_mesh, y_mesh = np.meshgrid(x,y)
+
+z = terrain
+X = create_X(x_mesh, y_mesh,m)
+
+
+# Show the terrain
+plt.figure()
+plt.title('Terrain over Norway 1')
+plt.imshow(terrain, cmap='gray')
+plt.xlabel('X')
+plt.ylabel('Y')
+plt.show()