I am doing remote sensing work and have a dataset containing a few hundred images. Each image has several bands, which contain information such as reflectance, incidence angle, longitude and latitude information. My goal is to create a database of reflectance and other information contained in the other bands on a pixel-by-pixel basis.

1. To reduce the large amount of data this would generate, I would like to create a region of interest (ROI) around a specified location of the surface and apply it to all bands of all images. How can I do this in IDL?

2. Then I would like to collect reflectance and other information from the bands to create the database. How can I accomplish this?

Thank you