GNR401 Principles of Satellite Image Processing
|
|
- Jessie Phelps
- 5 years ago
- Views:
Transcription
1 Principles of Satellite Image Processing Instructor: Prof. CSRE, IIT Bombay Slot 5 Guest Lecture PCT and Band Arithmetic November 07, AM AM
2 IIT Bombay Slide 1 November 07, 2012 PCT and Band Arithmetic Contents of the ure Principal Component Transform Alternative Interpretations Inverse PCT Band Arithmetic Motivation for band arithmetic Band ratio Vegetation indices
3 IIT Bombay Slide 2 Multiband Image Operations Operations performed by combining gray levels recorded in different bands for the same pixel Applications Data reduction through decorrelation Highlighting specific features with significant difference in response in different bands The transformed data may be viewed like enhanced versions compared to originals
4 IIT Bombay Slide 3 Principal Component Transform Highlights the redundancy in the data sets due to similar response in some of the wavelengths Original bands variables represented along different coordinate axes, redundancy implies variables are correlated, not independent Gray level in a band at a pixel can be predicted from the knowledge of the pixel gray level in other bands
5 IIT Bombay Slide 4 Example of Redundancy in Data Example: Highly correlated data Values along band b1 leads to knowledge along band b2 of the data element Linear variation (nearly) between b1 and b2 Often true in case of visible bands b2 b1
6 IIT Bombay Slide 5 Example of Redundancy in Data Points projected onto the line a small error in the position of the point. Points represented by only one coordinate b1 half data reduced For highly correlated data, this error will be minimal b2 b2 b1 b1
7 IIT Bombay Slide 6 Decorrelating Multispectral Remotely Sensed Data How do we identify the optimum axes along which the remotely sensed data should be projected so that the transformed data would be uncorrelated? What should be the way to rank the new axes so that we can discard the least important dimensions of the transformed data? Invertibility of the transformation?
8 IIT Bombay Slide 7 Useful band statistics M N i 1 j 1 Mean i M M N 1 j 1 k l ( g )( g ) ij k ij l M. N. N g k i j Variance M Covariance i 1 j 1 k ( g ) M. N N ij k 2
9 IIT Bombay Slide 8 Covariance Matrix C = {C kl k = 1,, K, l = 1,, K} K is the number of bands in which the multispectral dataset was generated C is a symmetric matrix C kl = C lk Diagonal elements of C are the intra-band variances Off-diagonal elements are the inter-band covariances
10 IIT Bombay Slide 9 Relation between correlation and Correlation R kl = covariance It can be shown that R kl = C kl + m k m l M i N 1 j 1 M For data with zero-mean, correlation and co-variance will be equal. g N k i j g l i j
11 IIT Bombay Slide 10 Principal Component Transformation Problem to solve: Find a transformation to be applied to the input multispectral image such that the covariance matrix of the result is reduced to a diagonal matrix Further, we should find an axis v k such that the variance of the projected coordinates (z k = v kt x) is maximum.
12 IIT Bombay Slide 11 Solution Given the transformed vector z k = v k t x The variance s z2 = This simplifies to s z2 = v t Cv (Dropping subscript k for a moment!) C, the covariance matrix is a positive, semi-definite, real symmetric matrix. M N i 1 j 1 t t v ( x )( x ) v ij k ij l M. N
13 IIT Bombay Slide 12 Finding vector v To maximize the projected variance s z2, find a v such that v t Cv is maximum, subject to the constraint v t v = 1. Combining the maximization function with the constraint, we can write v t Cv l(v t v 1) = maximum Differentiating w.r.t. v, t t C ( 1) 0 v v v v v
14 IIT Bombay Slide 13 Finding v The derivative results in Cv = lv (Verify!) Therefore, v is an eigenvector of C Given that s z2 = v t Cv v t Cv = v t (lv) = lv t v = l = s z 2 This implies that v is the eigenvector of C with the largest eigenvalue Therefore all the eigenvectors with decreasing eigenvalues lead to axes with decreasing variance along them.
15 IIT Bombay Slide 14 Alternative Explanation to PCT Let the transformed pixel vector y = D t x Covariance matrix of y = S y = D t S x D (Note that S y = E{(y m y )(y-m y ) t } = E{(D t x D t m x )(D t x D t m x ) t } This simplifies to S y = D t E(x m x )(x m x ) t D D is a set of vectors independent of x)
16 IIT Bombay Slide 15 Alternative Explanation Covariance matrix of y = S y = D t S x D It is desired that S y be diagonal, i.e., the data in the transformed domain is uncorrelated Let S y = n
17 IIT Bombay Slide 16 Let S y = Alternative Explanation n Then S y = D t S x D is a similarity transformation with D containing eigenvectors of S x We can order l i in such a way that they are in descending order. Given that y = D t x, y 1 corresponds to direction given by e 1, that is the first row of D t, Each transformed pixel vector y is obtained from scalar products of eigenvectors of S x and x
18 IIT Bombay Slide 17 Sample Eigenvectors and Eigenvalues Covariance Matrix Eigenvalues Eigenvectors
19 IIT Bombay Slide 18 Sample Eigenvectors and Eigenvalues
20 IIT Bombay Slide 19 Transformation New component value = dot product of eigenvector and pixel vector (i,j) pixel position n eigenvectors for n principal components 1 st principal component dot product of pixel vector with eigenvector corresponding to largest eigenvalue
21 IIT Bombay Slide 20 Principal Components For n input bands, n principal components are computed The utility of the principal components gradually decreases from 1 st towards the last e.g., For Landsat TM, last three PCs are generally of very little value
22 IIT Bombay Slide 21 Visualization of PCT From J.R. Jensen s ure notes at Univ. South Carolina; used with permission
23 IIT Bombay Slide 22 Comments on PCT For IRS / IKONOS images, out of four bands, 2-3 principal components capture most of the useful information. The last 1-2 bands are redundant. Advantages Smaller data volume to handle Principal components appear to be enhanced versions of the originals, having contributions from all the four input bands Application scientists use composites of PC for interpretation of various features such as geology
24 IIT Bombay Slide 23 Band 1 (Blue)
25 IIT Bombay Slide 24 Band 2 (Green)
26 IIT Bombay Slide 25 Band 3 (Red)
27 IIT Bombay Slide 26 Band 4 (NIR)
28 IIT Bombay Slide 27 Band 5 (SWIR)
29 IIT Bombay Slide 28 Band 7 (SWIR)
30 IIT Bombay Slide 29 PC1
31 IIT Bombay Slide 30 PC2
32 IIT Bombay Slide 31 PC3
33 IIT Bombay Slide 32 PC6
34 IIT Bombay Slide 33 Input Image FCC
35 IIT Bombay Slide 34 Decorrelation Stretch
36 IIT Bombay Slide 35 Inverse PCT Inverse PCT is used to generate the bands in the original domain If ALL PCTs are retained, inverse will give back the original bands If any PCTs are dropped, inverse will give new bands in the original domain that may be close to the original bands depending on how many PCTs are discarded
37 IIT Bombay Slide 36 Inverse PCT From the principle of PCT, we have y = D t x D t contains eigenvectors of S x, covariance matrix from the original image Since D t is an orthonormal matrix, (D t ) t = (D t ) -1 From each pixel vector in PC domain, x = (D t ) t y
38 IIT Bombay Slide 37 Inverse PCT For k band image, matrix D is square, of size k x k If m principal components are dropped, we are left with a matrix (D 1 ) of size k x (k-m) The vector y is reduced to y 1 of size k-m x 1 Therefore the modified vector x 1 is given by x 1 = D 1 y 1 The difference between x and x 1 is a measure of the loss of information due to removal of some of the PCs
39 IIT Bombay Slide 38 Comments on PCT One of the other important applications of PCT is data fusion Replace first PC by the image from another sensor Apply inverse PC
40 Band Arithmetic
41 IIT Bombay Slide 39 Motivation
42 IIT Bombay Slide 40 Multiband Arithmetic In a given pair of bands the response of two objects is generally different. Pixel by pixel comparison between images can highlight pixels that have very high difference in ref ance in those bands Operations like band difference and band ratio or combinations of them are popularly used for this purpose
43 IIT Bombay Slide 41 Band Ratio Very common operation Ratio i,j (m,n) = Band i (m,n) / Band j (m,n) If Band j (m,n) = 0, suitable adjustment has to be made (e.g., add +1 to the denom.) Minimum ratio will be 0; Maximum ratio will be 255
44 IIT Bombay Slide 42 Input Image
45 IIT Bombay Slide 43 Input Image FCC
46 IIT Bombay Slide 44 IR/R
47 IIT Bombay Slide 45 Band Ratio For fast computing, approximations can be made such as: 0 Ratio i,j (m,n) 1, Ratio i,j (m,n) scaled = Round [Ratio i,j (m,n)x127] 1 < Ratio i,j (m,n) 255, Ratio i,j (m,n) scaled = Round [127 + Ratio i,j (m,n)/2] Advantage in one pass image is generated in range 0-255
48 IIT Bombay Slide 46 Band Difference Similar to band ratio, band difference can also be used to account for difference in ref ance by objects in two wavelengths Band ratio - more popular in practical applications such as geological mapping Topographic effects on the images are reduced by ratioing.
49 IIT Bombay Slide 47 Band Multiplication Pixel by pixel multiplication of two images Not used to multiply gray levels in one band with corresponding gray levels in another band Used in practice to mask some part of the image and retain the rest of it by preparing a mask image and performing image to image multiplication of pixels
50 IIT Bombay Slide 48 Band Addition Similar to Band Multiplication, band addition has no direct practical application in adding gray levels of two bands of an image This method too can be used to mask a portion of the image and retain the remaining part.
51 IIT Bombay Slide 49 Specialized Indices Combination of band differences, ratios and additions can result in useful outputs that can highlight features like green vegetation One such feature is Normalized Difference Vegetation Index (NDVI) NDVI(m,n) = BandIR ( m, n) BandR ( m, n) Band ( m, n) Band ( m, n) IR R
52 IIT Bombay Slide 50 NDVI NDVI results in high values where IR dominates red wavelength. This happens where vegetation is present Range of NDVI is [-1 +1] NDVI has been widely used in a wide ranging of agricultural, forestry and biomass estimation applications It is also used to measure the length of crop growth and dry-down periods by comparing NDVI computed from multidate images
53 IIT Bombay Slide 51 Input Image
54 IIT Bombay Slide 52 NIR
55 IIT Bombay Slide 53 RED
56 IIT Bombay Slide 54 NDVI
57 IIT Bombay Slide 55 Other Vegetation Indices Simple Ratio = Red/NIR NDVI6 = (Band 6 Band 5)/(Band 6 + Band 5) NDVI7 = (Band 7 Band 5)/(Band 7 + Band 5) Standard NDVI TM = (TM4 TM3)/(TM4 + TM3) These are applicable when seven band data like Landsat Thematic Mapper data are available For IRS LISS3 imagery, NDVI IRS = Band ( m, n) Band ( m, n) 4 3 Band ( m, n) Band ( m, n) 4 3
58 IIT Bombay Slide 56 IRS L4- NDVI
59 IIT Bombay Slide 57 Fast Computation of NDVI Range of NDVI [-1, +1] Scale suitably to generate an NDVI image For example, NDVI scaled =127(1+NDVI) This ensures that the resultant NDVI has a range of [0 254]
60 IIT Bombay Slide 58 Se ed Ref ance Curves From J.R. Jensen s ure notes at Univ. South Carolina Used with permission
61 IIT Bombay Slide 59 Time Series of 1984 and 1988 NDVI Measurements Derived from AVHRR Global Area Coverage (GAC) Data Region around El Obeid, Sudan, in Sub-Saharan Africa From J.R. Jensen s ure notes at Univ. South Carolina Used with permission
62 IIT Bombay Slide 60 Simple Ratio v/s NDVI From J.R. Jensen s ure notes at Univ. South Carolina Used with permission
63 IIT Bombay Slide 61 Infrared Index Traditional NDVI does not work very well when the soil is moist, as in case of wetlands. The Infrared Index (II) can tackle this situation better II NIR NIR MIR TM 4 TM 5 MIR TM 4 TM 5 Several bands needed in the infrared region, as in case of Landsat TM
64 IIT Bombay Slide 62 Soil Line From J.R. Jensen s ure notes at Univ. South Carolina Used with permission
65 IIT Bombay Slide 63 Perpendicular Vegetation Index PVI is defined as 2 2 PVI S, R V, R S, NIR V, NIR A vegetation index that assumes that the ref ance in the NIR and red varies with increasing vegetation density (such as leaf area index) and that these variations are parallel to the soil baseline. Therefore, the perpendicular distance from the baseline in a NIR-red plot determines the vegetation density. See for more definitions of various indices in remote sensing including PVI
66 IIT Bombay Slide 64 Soil Adjusted Vegetation Index The soil adjusted vegetation index (SAVI) introduces a soil calibration factor, L, to the NDVI equation to minimize soil background influences resulting from first order soilplant spectral interactions: SAVI (1 L) NIR red NIR red L Ref:A. R. Huete, A soil-adjusted vegetation index (SAVI), Rem. Sens. Env., vol. 25, pp , 1988.
67 IIT Bombay Slide 65 Atmospherically Adjusted Vegetation Index (ARVI) The atmospheric effects are accounted for in ARVI ARVI p* nir p* rb p* nir p* rb p* rb p* red p* blue p* red p* indicates the atmospherically corrected versions of NIR, Red and Blue bands for molecular scattering and ozone absorption (Ref. J.R. Jensen s notes)
68 IIT Bombay Slide 66 Enhanced Vegetation Index EVI is a mixture of SAVI and ARVI, in that both atmospheric effects and soil effects are accounted for. EVI p* nir p * red p * nir C p* red C p* blue L 1 2 C1 and C2 describe the use of the blue band in correction of the red band for atmospheric aerosol scattering. The coefficients, C1, C2,, and L are empirically determined as 6.0, 7.5, and 1.0, respectively for MODIS. This algorithm has improved sensitivity to high biomass regions and improved vegetation monitoring through a de-coupling of the canopy background signal and a reduction in atmospheric influences Source: Yoram J. Kaufman and Didier Tanre, Atmospherically Resistant Vegetation Index (ARVI) for EOS-MODIS, IEEE Trans. GERS, vol. 30, no. 2, pp , 1992
69 IIT Bombay Slide 67 From J.R. Jensen s ure notes at Univ. South Carolina; used with permission
70 Contd
GEOG 4110/5100 Advanced Remote Sensing Lecture 15
GEOG 4110/5100 Advanced Remote Sensing Lecture 15 Principal Component Analysis Relevant reading: Richards. Chapters 6.3* http://www.ce.yildiz.edu.tr/personal/songul/file/1097/principal_components.pdf *For
More informationEnvironmental Remote Sensing GEOG 2021
Environmental Remote Sensing GEOG 2021 Lecture 3 Spectral information in remote sensing Spectral Information 2 Outline Mechanisms of variations in reflectance Optical Microwave Visualisation/analysis Enhancements/transforms
More informationMany of remote sensing techniques are generic in nature and may be applied to a variety of vegetated landscapes, including
Remote Sensing of Vegetation Many of remote sensing techniques are generic in nature and may be applied to a variety of vegetated landscapes, including 1. Agriculture 2. Forest 3. Rangeland 4. Wetland,
More informationIntroduction to Satellite Derived Vegetation Indices
Introduction to the Use of Geospatial Information Technology for Drought Risk Management 13-17 November, 2017 Tonle Bassac II Restaurant, Phnom Penh, Cambodia Introduction to Satellite Derived Vegetation
More informationWhat is a vegetation index?
Vegetation indexes What is a vegetation index? A tool: a vegetation index is a mathematical formula used to estimate the likelihood of vegetation presence in remotely sensed data A product: the result
More informationVegetation Change Detection of Central part of Nepal using Landsat TM
Vegetation Change Detection of Central part of Nepal using Landsat TM Kalpana G. Bastakoti Department of Geography, University of Calgary, kalpanagb@gmail.com Abstract This paper presents a study of detecting
More informationAssessing Drought in Agricultural Area of central U.S. with the MODIS sensor
Assessing Drought in Agricultural Area of central U.S. with the MODIS sensor Di Wu George Mason University Oct 17 th, 2012 Introduction: Drought is one of the major natural hazards which has devastating
More informationCHAPTER 4 PRINCIPAL COMPONENT ANALYSIS-BASED FUSION
59 CHAPTER 4 PRINCIPAL COMPONENT ANALYSIS-BASED FUSION 4. INTRODUCTION Weighted average-based fusion algorithms are one of the widely used fusion methods for multi-sensor data integration. These methods
More informationMODULE 5 LECTURE NOTES 5 PRINCIPAL COMPONENT ANALYSIS
MODULE 5 LECTURE NOTES 5 PRINCIPAL COMPONENT ANALYSIS. (PCA) Principal component analysis (PCA), also known as Karhunen-Loeve analysis, transforms the information inherent in multispectral remotely sensed
More informationThe Principal Component Analysis
The Principal Component Analysis Philippe B. Laval KSU Fall 2017 Philippe B. Laval (KSU) PCA Fall 2017 1 / 27 Introduction Every 80 minutes, the two Landsat satellites go around the world, recording images
More informationSPECTRAL DISCRIMINATION OF ROCK TYPES IN THE ARAVALLI MOUNTAIN RANGES OF RAJASTHAN (INDIA) USING LANDSAT THEMATIC MAPPER DATA
SPECTRAL DISCRIMINATION OF ROCK TYPES IN THE ARAVALLI MOUNTAIN RANGES OF RAJASTHAN (INDIA) USING LANDSAT THEMATIC MAPPER DATA Dr. Nilanchal Patel Reader, Department of Remote Sensing Birla Institute of
More informationGMES: calibration of remote sensing datasets
GMES: calibration of remote sensing datasets Jeremy Morley Dept. Geomatic Engineering jmorley@ge.ucl.ac.uk December 2006 Outline Role of calibration & validation in remote sensing Types of calibration
More informationIntroduction to Machine Learning
10-701 Introduction to Machine Learning PCA Slides based on 18-661 Fall 2018 PCA Raw data can be Complex, High-dimensional To understand a phenomenon we measure various related quantities If we knew what
More informationEstimation of Wavelet Based Spatially Enhanced Evapotranspiration Using Energy Balance Approach
Estimation of Wavelet Based Spatially Enhanced Evapotranspiration Using Energy Balance Approach Dr.Gowri 1 Dr.Thirumalaivasan 2 1 Associate Professor, Jerusalem College of Engineering, Department of Civil
More informationMathematical foundations - linear algebra
Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar
More informationLearning Objectives. Thermal Remote Sensing. Thermal = Emitted Infrared
November 2014 lava flow on Kilauea (USGS Volcano Observatory) (http://hvo.wr.usgs.gov) Landsat-based thermal change of Nisyros Island (volcanic) Thermal Remote Sensing Distinguishing materials on the ground
More informationCoastal Characterization Using EO-1 Hyperion Data
Coastal Characterization Using EO-1 Hyperion Data Dr. Hsiao-hua K. Burke EO-1 SVT Meeting 18-21 November 2002 Sponsor: NOAA NESDIS GOES 2002-1 Channel Positions of Various Ocean- Color Sensors, 1978-2000*
More informationHYPERSPECTRAL IMAGING
1 HYPERSPECTRAL IMAGING Lecture 9 Multispectral Vs. Hyperspectral 2 The term hyperspectral usually refers to an instrument whose spectral bands are constrained to the region of solar illumination, i.e.,
More informationOverview on Land Cover and Land Use Monitoring in Russia
Russian Academy of Sciences Space Research Institute Overview on Land Cover and Land Use Monitoring in Russia Sergey Bartalev Joint NASA LCLUC Science Team Meeting and GOFC-GOLD/NERIN, NEESPI Workshop
More informationFundamentals of Remote Sensing
Division of Spatial Information Science Graduate School Life and Environment Sciences University of Tsukuba Fundamentals of Remote Sensing Prof. Dr. Yuji Murayama Surantha Dassanayake 10/6/2010 1 Fundamentals
More information12.2 Dimensionality Reduction
510 Chapter 12 of this dimensionality problem, regularization techniques such as SVD are almost always needed to perform the covariance matrix inversion. Because it appears to be a fundamental property
More informationLecture Topics. 1. Vegetation Indices 2. Global NDVI data sets 3. Analysis of temporal NDVI trends
Lecture Topics 1. Vegetation Indices 2. Global NDVI data sets 3. Analysis of temporal NDVI trends Why use NDVI? Normalize external effects of sun angle, viewing angle, and atmospheric effects Normalize
More informationIR-MAD Iteratively Re-weighted Multivariate Alteration Detection
IR-MAD Iteratively Re-weighted Multivariate Alteration Detection Nielsen, A. A., Conradsen, K., & Simpson, J. J. (1998). Multivariate Alteration Detection (MAD) and MAF Postprocessing in Multispectral,
More informationPercentage of Vegetation Cover Change Monitoring in Wuhan Region Based on Remote Sensing
Available online at www.sciencedirect.com Procedia Environmental Sciences 10 (2011 ) 1466 1472 2011 3rd International Conference on Environmental Science and Information Conference Application Title Technology
More informationEECS490: Digital Image Processing. Lecture #26
Lecture #26 Moments; invariant moments Eigenvector, principal component analysis Boundary coding Image primitives Image representation: trees, graphs Object recognition and classes Minimum distance classifiers
More informationVegetation Remote Sensing
Vegetation Remote Sensing Huade Guan Prepared for Remote Sensing class Earth & Environmental Science University of Texas at San Antonio November 2, 2005 Outline Why do we study vegetation remote sensing?
More informationLECTURE 16: PCA AND SVD
Instructor: Sael Lee CS549 Computational Biology LECTURE 16: PCA AND SVD Resource: PCA Slide by Iyad Batal Chapter 12 of PRML Shlens, J. (2003). A tutorial on principal component analysis. CONTENT Principal
More information1 Principal Components Analysis
Lecture 3 and 4 Sept. 18 and Sept.20-2006 Data Visualization STAT 442 / 890, CM 462 Lecture: Ali Ghodsi 1 Principal Components Analysis Principal components analysis (PCA) is a very popular technique for
More informationISO MODIS NDVI Weekly Composites for Canada South of 60 N Data Product Specification
ISO 19131 MODIS NDVI Weekly Composites for South of 60 N Data Product Specification Revision: A Data specification: MODIS NDVI Composites for South of 60 N - Table of Contents - 1. OVERVIEW... 3 1.1. Informal
More informationDrought Estimation Maps by Means of Multidate Landsat Fused Images
Remote Sensing for Science, Education, Rainer Reuter (Editor) and Natural and Cultural Heritage EARSeL, 2010 Drought Estimation Maps by Means of Multidate Landsat Fused Images Diego RENZA, Estíbaliz MARTINEZ,
More informationGEO 874 Remote Sensing. Zihan Lin, Nafiseh Hagtalab, Ranjeet John
GEO 874 Remote Sensing Zihan Lin, Nafiseh Hagtalab, Ranjeet John http://onlinelibrary.wiley.com/doi/10.1002/wat2.1147/full Landscape Heterogeneity S-I: homogeneous landscape using only cropland for the
More informationDigital Change Detection Using Remotely Sensed Data for Monitoring Green Space Destruction in Tabriz
Int. J. Environ. Res. 1 (1): 35-41, Winter 2007 ISSN:1735-6865 Graduate Faculty of Environment University of Tehran Digital Change Detection Using Remotely Sensed Data for Monitoring Green Space Destruction
More informationA high spectral resolution global land surface infrared emissivity database
A high spectral resolution global land surface infrared emissivity database Eva E. Borbas, Robert O. Knuteson, Suzanne W. Seemann, Elisabeth Weisz, Leslie Moy, and Hung-Lung Huang Space Science and Engineering
More informationPrincipal Components Analysis (PCA)
Principal Components Analysis (PCA) Principal Components Analysis (PCA) a technique for finding patterns in data of high dimension Outline:. Eigenvectors and eigenvalues. PCA: a) Getting the data b) Centering
More informationPrincipal Component Analysis-I Geog 210C Introduction to Spatial Data Analysis. Chris Funk. Lecture 17
Principal Component Analysis-I Geog 210C Introduction to Spatial Data Analysis Chris Funk Lecture 17 Outline Filters and Rotations Generating co-varying random fields Translating co-varying fields into
More informationDETECTION OF HOTSPOTS IN NOAA/AVHRR IMAGES USING PRINCIPAL COMPONENT ANALYSIS AND INFORMATION FUSION TECHNIQUE
DETECTION OF HOTSPOTS IN NOAA/AVHRR IMAGES USING PRINCIPAL COMPONENT ANALYSIS AND INFORMATION FUSION TECHNIQUE R. S. Gautam, D. Singh, A. Mittal {rsgcsdec@iitr.ernet.in, dharmfec@iitr.ernet.in, ankumfec@iitr.ernet.in}
More informationCHANGE DETECTION USING REMOTE SENSING- LAND COVER CHANGE ANALYSIS OF THE TEBA CATCHMENT IN SPAIN (A CASE STUDY)
CHANGE DETECTION USING REMOTE SENSING- LAND COVER CHANGE ANALYSIS OF THE TEBA CATCHMENT IN SPAIN (A CASE STUDY) Sharda Singh, Professor & Programme Director CENTRE FOR GEO-INFORMATICS RESEARCH AND TRAINING
More informationLanduse and Landcover change analysis in Selaiyur village, Tambaram taluk, Chennai
Landuse and Landcover change analysis in Selaiyur village, Tambaram taluk, Chennai K. Ilayaraja Department of Civil Engineering BIST, Bharath University Selaiyur, Chennai 73 ABSTRACT The synoptic picture
More informationPrincipal Components Theory Notes
Principal Components Theory Notes Charles J. Geyer August 29, 2007 1 Introduction These are class notes for Stat 5601 (nonparametrics) taught at the University of Minnesota, Spring 2006. This not a theory
More informationElementary Linear Algebra Review for Exam 2 Exam is Monday, November 16th.
Elementary Linear Algebra Review for Exam Exam is Monday, November 6th. The exam will cover sections:.4,..4, 5. 5., 7., the class notes on Markov Models. You must be able to do each of the following. Section.4
More informationA Tutorial on Data Reduction. Principal Component Analysis Theoretical Discussion. By Shireen Elhabian and Aly Farag
A Tutorial on Data Reduction Principal Component Analysis Theoretical Discussion By Shireen Elhabian and Aly Farag University of Louisville, CVIP Lab November 2008 PCA PCA is A backbone of modern data
More informationCross-Sensor Continuity Science Algorithm
Cross-Sensor Continuity Science Algorithm - Long Term Vegetation Index and Phenology Workshop - Javzan Tsend-Ayush and Tomoaki Miura Department of Natural Resources and Environmental Management University
More informationLand cover/land use mapping and cha Mongolian plateau using remote sens. Title. Author(s) Bagan, Hasi; Yamagata, Yoshiki. Citation Japan.
Title Land cover/land use mapping and cha Mongolian plateau using remote sens Author(s) Bagan, Hasi; Yamagata, Yoshiki International Symposium on "The Imp Citation Region Specific Systems". 6 Nove Japan.
More informationA GLOBAL ANALYSIS OF URBAN REFLECTANCE. Christopher SMALL
A GLOBAL ANALYSIS OF URBAN REFLECTANCE Christopher SMALL Lamont Doherty Earth Observatory Columbia University Palisades, NY 10964 USA small@ldeo.columbia.edu ABSTRACT Spectral characterization of urban
More informationidentify tile lines. The imagery used in tile lines identification should be processed in digital format.
Question and Answers: Automated identification of tile drainage from remotely sensed data Bibi Naz, Srinivasulu Ale, Laura Bowling and Chris Johannsen Introduction: Subsurface drainage (popularly known
More information7. Variable extraction and dimensionality reduction
7. Variable extraction and dimensionality reduction The goal of the variable selection in the preceding chapter was to find least useful variables so that it would be possible to reduce the dimensionality
More informationChapter 6 Spectral Domain Image Transforms
Chapter 6 Spectral Domain Image Transforms 6.1 Introduction Because of the multispectral or vector character of remote sensing image data we can devise transformations that generate new pixel vector descriptions
More informationPrinciple Components Analysis (PCA) Relationship Between a Linear Combination of Variables and Axes Rotation for PCA
Principle Components Analysis (PCA) Relationship Between a Linear Combination of Variables and Axes Rotation for PCA Principle Components Analysis: Uses one group of variables (we will call this X) In
More informationVectors To begin, let us describe an element of the state space as a point with numerical coordinates, that is x 1. x 2. x =
Linear Algebra Review Vectors To begin, let us describe an element of the state space as a point with numerical coordinates, that is x 1 x x = 2. x n Vectors of up to three dimensions are easy to diagram.
More informationLecture 16: Small Sample Size Problems (Covariance Estimation) Many thanks to Carlos Thomaz who authored the original version of these slides
Lecture 16: Small Sample Size Problems (Covariance Estimation) Many thanks to Carlos Thomaz who authored the original version of these slides Intelligent Data Analysis and Probabilistic Inference Lecture
More informationPrincipal Component Analysis
CSci 5525: Machine Learning Dec 3, 2008 The Main Idea Given a dataset X = {x 1,..., x N } The Main Idea Given a dataset X = {x 1,..., x N } Find a low-dimensional linear projection The Main Idea Given
More informationCHAPTER 6 HYDROTHERMAL ALTERATION ZONES MAPPING USING PRINCIPAL COMPONENT ANALYSIS AND BAND RATIO TECHNIQUES
CHAPTER 6 HYDROTHERMAL ALTERATION ZONES MAPPING USING PRINCIPAL COMPONENT ANALYSIS AND BAND RATIO TECHNIQUES 6.1 INTRODUCTION Hydrothermal alteration zones mapping is one of the most common applications
More informationMULTI-CHANNEL REMOTE SENSING DATA AND ORTHOGONAL TRANSFORMATIONS FOR CHANGE DETECTION 1
MULTI-CHANNEL REMOTE SENSING DATA AND ORTHOGONAL TRANSFORMATIONS FOR CHANGE DETECTION 1 Allan Aasbjerg Nielsen IMM, Department of Mathematical Modelling Technical University of Denmark, Building 321 DK
More informationPrincipal Component Analysis -- PCA (also called Karhunen-Loeve transformation)
Principal Component Analysis -- PCA (also called Karhunen-Loeve transformation) PCA transforms the original input space into a lower dimensional space, by constructing dimensions that are linear combinations
More informationMachine Learning 2nd Edition
INTRODUCTION TO Lecture Slides for Machine Learning 2nd Edition ETHEM ALPAYDIN, modified by Leonardo Bobadilla and some parts from http://www.cs.tau.ac.il/~apartzin/machinelearning/ The MIT Press, 2010
More informationStructure in Data. A major objective in data analysis is to identify interesting features or structure in the data.
Structure in Data A major objective in data analysis is to identify interesting features or structure in the data. The graphical methods are very useful in discovering structure. There are basically two
More informationPrincipal Component Analysis!! Lecture 11!
Principal Component Analysis Lecture 11 1 Eigenvectors and Eigenvalues g Consider this problem of spreading butter on a bread slice 2 Eigenvectors and Eigenvalues g Consider this problem of stretching
More informatione 2 e 1 (a) (b) (d) (c)
2.13 Rotated principal component analysis [Book, Sect. 2.2] Fig.: PCA applied to a dataset composed of (a) 1 cluster, (b) 2 clusters, (c) and (d) 4 clusters. In (c), an orthonormal rotation and (d) an
More informationThermal Infrared (TIR) Remote Sensing: Challenges in Hot Spot Detection
Thermal Infrared (TIR) Remote Sensing: Challenges in Hot Spot Detection ASF Seminar Series Anupma Prakash Day : Tuesday Date : March 9, 2004 Time : 2.00 pm to 2.30 pm Place : GI Auditorium Geophysical
More informationL3: Review of linear algebra and MATLAB
L3: Review of linear algebra and MATLAB Vector and matrix notation Vectors Matrices Vector spaces Linear transformations Eigenvalues and eigenvectors MATLAB primer CSCE 666 Pattern Analysis Ricardo Gutierrez-Osuna
More informationVectors and Matrices Statistics with Vectors and Matrices
Vectors and Matrices Statistics with Vectors and Matrices Lecture 3 September 7, 005 Analysis Lecture #3-9/7/005 Slide 1 of 55 Today s Lecture Vectors and Matrices (Supplement A - augmented with SAS proc
More informationMETRIC tm. Mapping Evapotranspiration at high Resolution with Internalized Calibration. Shifa Dinesh
METRIC tm Mapping Evapotranspiration at high Resolution with Internalized Calibration Shifa Dinesh Outline Introduction Background of METRIC tm Surface Energy Balance Image Processing Estimation of Energy
More informationMonitoring Vegetation Growth of Spectrally Landsat Satellite Imagery ETM+ 7 & TM 5 for Western Region of Iraq by Using Remote Sensing Techniques.
Monitoring Vegetation Growth of Spectrally Landsat Satellite Imagery ETM+ 7 & TM 5 for Western Region of Iraq by Using Remote Sensing Techniques. Fouad K. Mashee, Ahmed A. Zaeen & Gheidaa S. Hadi Remote
More informationMULTI-SOURCE IMAGE CLASSIFICATION
MULTI-SOURCE IMAGE CLASSIFICATION Hillary Tribby, James Kroll, Daniel Unger, I-Kuai Hung, Hans Williams Corresponding Author: Daniel Unger (unger@sfasu.edu Arthur Temple College of Forestry and Agriculture
More informationEigenvalues, Eigenvectors, and an Intro to PCA
Eigenvalues, Eigenvectors, and an Intro to PCA Eigenvalues, Eigenvectors, and an Intro to PCA Changing Basis We ve talked so far about re-writing our data using a new set of variables, or a new basis.
More informationEigenvalues, Eigenvectors, and an Intro to PCA
Eigenvalues, Eigenvectors, and an Intro to PCA Eigenvalues, Eigenvectors, and an Intro to PCA Changing Basis We ve talked so far about re-writing our data using a new set of variables, or a new basis.
More informationASSESSMENT OF NDVI FOR DIFFERENT LAND COVERS BEFORE AND AFTER ATMOSPHERIC CORRECTIONS
Bulletin of the Transilvania University of Braşov Series II: Forestry Wood Industry Agricultural Food Engineering Vol. 7 (56) No. 1-2014 ASSESSMENT OF NDVI FOR DIFFERENT LAND COVERS BEFORE AND AFTER ATMOSPHERIC
More informationImage Registration Lecture 2: Vectors and Matrices
Image Registration Lecture 2: Vectors and Matrices Prof. Charlene Tsai Lecture Overview Vectors Matrices Basics Orthogonal matrices Singular Value Decomposition (SVD) 2 1 Preliminary Comments Some of this
More informationFundamental Interactions with Earth Surface
Fundamental Interactions with Earth Surface 1 / 96 http://speclab.cr.usgs.gov/papers/tetracorder/ GEO 827 - Digital Image Processing and Analysis 2 / 96 http://www.markelowitz.com/hyperspectral.html 3
More informationMaths for Signals and Systems Linear Algebra in Engineering
Maths for Signals and Systems Linear Algebra in Engineering Lecture 18, Friday 18 th November 2016 DR TANIA STATHAKI READER (ASSOCIATE PROFFESOR) IN SIGNAL PROCESSING IMPERIAL COLLEGE LONDON Mathematics
More informationIndian National (Weather) SATellites for Agrometeorological Applications
Indian National (Weather) SATellites for Agrometeorological Applications Bimal K. Bhattacharya Agriculture-Terrestrial Biosphere- Hydrology Group Space Applications Centre (ISRO) Ahmedabad 380015, India
More informationDETECTION OF CHANGES IN FOREST LANDCOVER TYPE AFTER FIRES IN PORTUGAL
DETECTION OF CHANGES IN FOREST LANDCOVER TYPE AFTER FIRES IN PORTUGAL Paulo M. Barbosa, Mário R. Caetano, and Teresa G. Santos Centro Nacional de Informação Geográfica (CNIG), Portugal barp@cnig.pt, mario@cnig.pt,
More informationGreening of Arctic: Knowledge and Uncertainties
Greening of Arctic: Knowledge and Uncertainties Jiong Jia, Hesong Wang Chinese Academy of Science jiong@tea.ac.cn Howie Epstein Skip Walker Moscow, January 28, 2008 Global Warming and Its Impact IMPACTS
More informationPrincipal component analysis
Principal component analysis Angela Montanari 1 Introduction Principal component analysis (PCA) is one of the most popular multivariate statistical methods. It was first introduced by Pearson (1901) and
More informationPRINCIPAL COMPONENT ANALYSIS
PRINCIPAL COMPONENT ANALYSIS Dimensionality Reduction Tzompanaki Katerina Dimensionality Reduction Unsupervised learning Goal: Find hidden patterns in the data. Used for Visualization Data compression
More informationMachine Learning (Spring 2012) Principal Component Analysis
1-71 Machine Learning (Spring 1) Principal Component Analysis Yang Xu This note is partly based on Chapter 1.1 in Chris Bishop s book on PRML and the lecture slides on PCA written by Carlos Guestrin in
More informationSingular Value Decomposition. 1 Singular Value Decomposition and the Four Fundamental Subspaces
Singular Value Decomposition This handout is a review of some basic concepts in linear algebra For a detailed introduction, consult a linear algebra text Linear lgebra and its pplications by Gilbert Strang
More informationTesting the Sensitivity of Vegetation Indices for Crop Type Classification Using Rapideye Imagery
Testing the Sensitivity of Vegetation Indices for Crop Type Classification Using Rapideye Imagery Mustafa USTUNER and Fusun BALIK SANLI, Turkey Keywords: Remote Sensing, Agriculture, Geoinformation SUMMARY
More informationPrincipal Component Analysis (PCA) CSC411/2515 Tutorial
Principal Component Analysis (PCA) CSC411/2515 Tutorial Harris Chan Based on previous tutorial slides by Wenjie Luo, Ladislav Rampasek University of Toronto hchan@cs.toronto.edu October 19th, 2017 (UofT)
More informationData Mining and Analysis: Fundamental Concepts and Algorithms
Data Mining and Analysis: Fundamental Concepts and Algorithms dataminingbook.info Mohammed J. Zaki 1 Wagner Meira Jr. 2 1 Department of Computer Science Rensselaer Polytechnic Institute, Troy, NY, USA
More informationThe Effects of Haze on the Accuracy of. Satellite Land Cover Classification
Applied Mathematical Sciences, Vol. 9, 215, no. 49, 2433-2443 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ams.215.52157 The Effects of Haze on the Accuracy of Satellite Land Cover ification
More informationMAPPING LAND USE/ LAND COVER OF WEST GODAVARI DISTRICT USING NDVI TECHNIQUES AND GIS Anusha. B 1, Sridhar. P 2
MAPPING LAND USE/ LAND COVER OF WEST GODAVARI DISTRICT USING NDVI TECHNIQUES AND GIS Anusha. B 1, Sridhar. P 2 1 M. Tech. Student, Department of Geoinformatics, SVECW, Bhimavaram, A.P, India 2 Assistant
More informationSATELLITE REMOTE SENSING
SATELLITE REMOTE SENSING of NATURAL RESOURCES David L. Verbyla LEWIS PUBLISHERS Boca Raton New York London Tokyo Contents CHAPTER 1. SATELLITE IMAGES 1 Raster Image Data 2 Remote Sensing Detectors 2 Analog
More informationPCA, Kernel PCA, ICA
PCA, Kernel PCA, ICA Learning Representations. Dimensionality Reduction. Maria-Florina Balcan 04/08/2015 Big & High-Dimensional Data High-Dimensions = Lot of Features Document classification Features per
More informationPCA & ICA. CE-717: Machine Learning Sharif University of Technology Spring Soleymani
PCA & ICA CE-717: Machine Learning Sharif University of Technology Spring 2015 Soleymani Dimensionality Reduction: Feature Selection vs. Feature Extraction Feature selection Select a subset of a given
More informationPRINCIPAL COMPONENT ANALYSIS
PRINCIPAL COMPONENT ANALYSIS 1 INTRODUCTION One of the main problems inherent in statistics with more than two variables is the issue of visualising or interpreting data. Fortunately, quite often the problem
More informationSeasonal and interannual relations between precipitation, soil moisture and vegetation in the North American monsoon region
Seasonal and interannual relations between precipitation, soil moisture and vegetation in the North American monsoon region Luis A. Mendez-Barroso 1, Enrique R. Vivoni 1, Christopher J. Watts 2 and Julio
More informationMATH 20F: LINEAR ALGEBRA LECTURE B00 (T. KEMP)
MATH 20F: LINEAR ALGEBRA LECTURE B00 (T KEMP) Definition 01 If T (x) = Ax is a linear transformation from R n to R m then Nul (T ) = {x R n : T (x) = 0} = Nul (A) Ran (T ) = {Ax R m : x R n } = {b R m
More informationPrincipal Components Analysis. Sargur Srihari University at Buffalo
Principal Components Analysis Sargur Srihari University at Buffalo 1 Topics Projection Pursuit Methods Principal Components Examples of using PCA Graphical use of PCA Multidimensional Scaling Srihari 2
More informationData Preprocessing Tasks
Data Tasks 1 2 3 Data Reduction 4 We re here. 1 Dimensionality Reduction Dimensionality reduction is a commonly used approach for generating fewer features. Typically used because too many features can
More informationEMPIRICAL ESTIMATION OF VEGETATION PARAMETERS USING MULTISENSOR DATA FUSION
EMPIRICAL ESTIMATION OF VEGETATION PARAMETERS USING MULTISENSOR DATA FUSION Franz KURZ and Olaf HELLWICH Chair for Photogrammetry and Remote Sensing Technische Universität München, D-80290 Munich, Germany
More informationLecture: Face Recognition and Feature Reduction
Lecture: Face Recognition and Feature Reduction Juan Carlos Niebles and Ranjay Krishna Stanford Vision and Learning Lab Lecture 11-1 Recap - Curse of dimensionality Assume 5000 points uniformly distributed
More informationMathematical foundations - linear algebra
Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar
More informationQuantitative Understanding in Biology Principal Components Analysis
Quantitative Understanding in Biology Principal Components Analysis Introduction Throughout this course we have seen examples of complex mathematical phenomena being represented as linear combinations
More informationComputer Assisted Image Analysis
Computer Assisted Image Analysis Lecture 0 - Object Descriptors II Amin Allalou amin@cb.uu.se Centre for Image Analysis Uppsala University 2009-04-27 A. Allalou (Uppsala University) Object Descriptors
More informationTHE LAND-SAF SURFACE ALBEDO AND DOWNWELLING SHORTWAVE RADIATION FLUX PRODUCTS
THE LAND-SAF SURFACE ALBEDO AND DOWNWELLING SHORTWAVE RADIATION FLUX PRODUCTS Bernhard Geiger, Dulce Lajas, Laurent Franchistéguy, Dominique Carrer, Jean-Louis Roujean, Siham Lanjeri, and Catherine Meurey
More informationLinear Dimensionality Reduction
Outline Hong Chang Institute of Computing Technology, Chinese Academy of Sciences Machine Learning Methods (Fall 2012) Outline Outline I 1 Introduction 2 Principal Component Analysis 3 Factor Analysis
More informationAFRI aerosol free vegetation index
Remote Sensing of Environment 77 (2001) 10 21 www.elsevier.com/locate/rse AFRI aerosol free vegetation index Arnon Karnieli a, *Yoram J. Kaufman b Lorraine Remer b Andrew Wald c a The Remote Sensing Laboratory,
More informationReflectivity in Remote Sensing
Reflectivity in Remote Sensing The amount of absorbance and reflection of white light by a substance is dependent upon the molecular makeup of the substance. Humans have used dyes for years to obtain colors-
More informationPrincipal Component Analysis (PCA)
Principal Component Analysis (PCA) Salvador Dalí, Galatea of the Spheres CSC411/2515: Machine Learning and Data Mining, Winter 2018 Michael Guerzhoy and Lisa Zhang Some slides from Derek Hoiem and Alysha
More information