Remote sensing for agricultural applications: principles and methods (2013-2014) Instructor: Prof. Tao Cheng (tcheng@njau.edu.cn) Nanjing Agricultural University Lab 1 Introduction to ENVI April 1 st, 2014
Why ENVI? One of the best remote sensing software packages in the world Very popular in academic institutions Famous for hyperspectral analysis Developed using the programming language IDL and therefore easy to be extended with IDL
ENVI versions Latest version: ENVI 5.1 + IDL 8.3 ENVI 5.1 A significant change in the interface Changes in ENVI and IDL programming ENVI 5.1 Classic has the ENVI heritage interface. If you are a long time ENVI user since ENVI 4.0, you may want to use the Classic one. Run your old IDL scripts in ENVI Classic.
ENVI Start( 开始 ) -> Programs( 程序 )-> ENVI 5.1 -> Tools -> ENVI Classic + IDL After ENVI Classic is started, you should see a menu as below: Interactive Displays File management Basic Tools Classification Tools Transform Tools Filter Tools Spectral Tools Map Tools Vector Tools Topographic Tools Radar Tools
Test data: a Landsat Scene Collected by the OLI on Landsat-8 Time of acquisition: 2014-03-16, 05:07:50 A Landsat-8 false color composite (4-5-3) Data_type = L1T (Landsat-8 Level 1 Terrain corrected product) Satellite orbit location at acquisition: path=119, row=238 File type: each band is in a separate GeoTIFF file Data format: 16-bit unsigned integer LANDSAT_SCENE_ID = "LC81190382014075LGN00" Data download: http://earthexplorer.usgs.gov/ Data files 8 TIFF files and one txt file We only use the 7 multispectral bands.
Landsat -8 band designations Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS) Launched February 11, 2013 Bands Wavelength (µm) Band center (µm) Resolution (m) Band 1 - Coastal 0.43-0.45 0.433 30 Band 2 - Blue 0.45-0.51 0.482 30 Band 3 - Green 0.53-0.59 0.562 30 Band 4 - Red 0.64-0.67 0.655 30 Band 5 - Near Infrared (NIR) 0.85-0.88 0.865 30 Band 6 - SWIR 1 1.57-1.65 1.610 30 Band 7 - SWIR 2 2.11-2.29 2.200 30 Band 8 - Panchromatic 0.50-0.68 0.590 15 Band 9 - Cirrus 1.36-1.38 1.375 30 Band 10 - Thermal Infrared (TIRS) 1 Band 11 - Thermal Infrared (TIRS) 2 10.60-11.19 10.800 100 11.50-12.51 12.000 100
What to do with this Landsat scene in ENVI? 1. Image loading 2. Image displaying 3. Image stacking 4. Image cropping 5. Image calibration
All the file and sensor types that supported by ENVI 5.1 classic. These files can be directly opened in ENVI. We are going to open Landsat-8 image files that are downloaded in GeoTIFF format. Tips: to select multiple files, hold the Ctrl key and click.
Load your data Once your files are loaded into ENVI, you should see this window (Available Bands List). Display mode Gray Scale: individual bands RGB Color: band combinations (e.g. R-NIR-G) Try making different color composites Data info: Dims 7791 x 7941: # of columns x # of rows Unsigned Int: data type = Unsigned integer (16 bit) (data range?) BSQ: data storage method What is the file size of a 7-band Landsat-8 image?
Interactive displays Image window 1. Be familiar with the display group. Navigate on your image and focus on your target area. Scroll window Scroll window Zoom control
2. Double click in the Image window to display cursor info. Image window How to read these values: Disp #: (col #, row #) Scrn: screen values (0~255) Projection: coordinate system Map: coordinates in m LL: coordinates in Lat/Lon Data: DN for the RGB bands Scroll window Zoom control Scroll window
3. Right click in any window to bring up a box with commonly used functions. Image window Right click with your mouse: How to locate the pixel for Ruigao (32.26598ºN, 120.76102ºE)? Scroll window Scroll window Zoom control
Save data as ENVI standard format Step 1 Step 2 Load Landsat-8 bands 1~7 and stack them up to generate a multispectral image file. Step 3 Step 4
Edit band names Step 2 Step 4 Step 1 Assign meaningful band names. Find this file in your working directory. Check if the file size matches the number you calculated in slide 8. Step 3
Image cropping: create a subset image Step 1: draw a rectangle ROI in the Scroll window to include your study area (hold left key and drag, right click to select).
The subset image is in display.
Conversion of DN to TOA (Top of Atmosphere) radiance L = M L DN + A L L λ = TOA spectral radiance (W/m2*srad*um) M L = the band-specific multiplicative factor (RADIANCE_MULT_BAND_x) A L = the band-specific multiplicative factor (RADIANCE_ADD_BAND_x) DN = Landsat-8 level-1 product pixel values Step 1 Step 2
Input band-specific gain and offset values Gain = RADIANCE_MULT_BAND_x (in the *_MTL.txt file) Offset = RADIANCE_ADD_BAND_x
Compare DN data to radiance data Radiance data Unit: W/m2*srad*um DN data Unit: dimensionless
Resources ENVI help documentation http://www.exelisvis.com/company/pressroom/blogs/imagerysp eaks.aspx - Imagery Speaks Blog (a good blog site) http://www.exelisvis.com/learn/eventstraining/webinars.aspx - ENVI/IDL webinars http://www.exelisvis.com/productsservices/idl.aspx - Exelis VIS, the developer of IDL. http://www.exelisvis.com/company/pressroom/blogs/idldatapo int.aspx - IDL data point (a good Blog site) LearnIDL.com Dr. Ken Bowman s web site http://www.idlcoyote.com/ - Coyote s guide to IDL programming (IDL tips, training, consulting) http://gumley.com/ - Practical IDL programming (IDL book) http://michaelgalloy.com/ - resources for IDL developers
Class project Think about what to do in your project. Make a tem of two students for each project. Start developing ideas and collecting materials (making notes is necessary) for your project. Allocate reasonable workload.