Few times ago, I worked with daily land surface temperature satellite images of an area with the size of Alaska.
I had to process around one thousand images. I downloaded them, mosaiced, projected, clipped to shapefile boundaries, extracted means out of the clipped rasters, and plotted the means on a time-temperature graph. The time-temperature graph was the final output. I did all these in ArcGIS, and even though I created an ArcPy script which helped automating a part of the work, it was still painstaking.
Now I want to build a program that does the same thing, but 100% automatically. I merely want to input to the program the coordinates of the corners in order to define the area of interest and get the final output generated which is the time-temperature graph.
Do you think this is possible using Python with the GDAL library and PostGIS as a raster data store?