Calculate the Sum of the nightlights data from TIFF images

This program read the nightlights data from the satellite images (TIFF) to calculate the sum of the nightlights inside each geometry polygon and write the results in a CSV file.
I was working on this project during my training in KAPSARC (King Abdullah Petroleum Studies and Research Center).
With my supervisor Hector G. LOPEZ-RUIZ
from May- July 2018.

The two main algorithms that I have used in this program are:
1- [A,x,y,I]=geoimread(filename,varargin) by Aslak Grinsted
from : https://ww2.mathworks.cn/matlabcentral/fileexchange/46904-geoimread
2- insidepoly(gridsX,gridsY,shapeX,shapeY) by Bruno Luong
from: https://ww2.mathworks.cn/matlabcentral/fileexchange/27840-2d-polygon-interior-detection?s_tid=answers_rc2-1_p4_MLT