GP265 / EE355 Homework 7 (Final project 1) Solutions
|
|
- Monica Barrett
- 6 years ago
- Views:
Transcription
1 GP265 / EE355 Homework 7 (Final project 1) Solutions We form the interferogram by multiplying image S 1 and conjugated image S 2 together: S int = S 1 S 2. The multi-look version of the complex interferogram is computed by averaging 4 range pixels and 4 azimuth pixels, so it has dimensions 3072/4 = 768 range bins by 14336/4 = 3584 azimuth bins. We plot the phase of the multi-look interferogram, which ranges from π to π. We do not see any fringes; it looks like random noise, because we have not yet resampled the second image at the correct range and azimuth offsets. 2. Calculate the offsets over the original images by cross-correlating small regions of each image. 1
2 We can use the offset program executable downloaded from the EE 355 website under Software. I ran it on cardinal.stanford.edu with the following parameters: $./offset image.1 image n This will correlate the images over 3000 small regions: 30 across (in the range direction) and 100 down (in the azimuth direction). It is important that the regions are representative of the entire image, rather than just a subset of the image. We also need to specify the approximate offset of the images for the last two parameters, which we can estimate by manually picking the same point in image 1 and image 2 (for example, the bottom corner of the triangular dark spot) and finding the pixel shift between them. I estimated a shift of 17 pixels across (in range) and -20 pixels down (in azimuth). If we do not choose good starting values for these approximate offsets, the offset program will produce junk output. The offset program outputs a file called offset.out which we can read into MATLAB for the rest of the processing. Next, I fit a least squares line to the range offset data as a function of range bin. I first removed outliers from the data, where I defined outlier as any data point outside one standard deviation of the median of the range offset data - this definitely improved the quality of the fit to the data (and this gets rid of any correlation regions that were on the water region of the image). The least squares fit (after removing outliers) results in two parameters: the slope r slope = x 10 4 and the intercept r offset = pixels of the line relating range offset p r (in pixels) to range bin i r : p r = r offset + r slope i r (1) I follow similar steps to get the azimuth offset. The least squares fit (after removing outliers) results in two parameters: the slope az slope = x 10 5 and the intercept az offset = pixels of the line relating azimuth offset p az (in pixels) to azimuth bin i az : p az = az offset + az slope i az (2) The magnitude of the azimuth slope az slope is much smaller than that of the range slope r slope, so we can approximate the azimuth offset as constant for the entire image. The average azimuth offset is pixels. 2
3 3. Baseline calculation The equation for the range offset δ r (in meters) in terms of relative range from the first bin dr = i r r is: δ r = B par + B perp dr (3) r 0 tan θ 0 We can also express the range offset by multiplying both sides of Eq. (1) by r: where r is the slant range pixel spacing: r = δ r = p r r = r offset r + r slope i r r (4) c = m/s 2f s 2( Hz) = m. (5) Comparing Eq. (4) to Eq. (3), I can see that B par = r offset r = m. Also, I can see that r slope = B perp B perp = r slope r 0 tan θ 0 (6) r 0 tan θ 0 The range to the swath center is r 0 = m. The look angle to swath center, θ 0, is computed by this equation obtained using the law of cosines for a spherical earth 3
4 geometry: ( θ0 = cos 1 r02 + z 2 + 2re z 2r0 (z + re ) ) = rad = 32.53o (7) Plugging in the results from Eq. (7) into Eq. (6), along with rslope = x 10 4, I get Bperp = m. The baseline length B is: B= 2 + B2 Bpar perp = m The orientation angle α is: ( ) Bperp π 1 α = tan + θ0 = rad = 3.89o Bpar 2 (8) (9) 4. Resample image.2 to align with image.1 and form the new interferogram. Since the Doppler centroid is fdc = 0, there is no need to shift the images to baseband before doing the resampling. I used MATLAB s meshgrid() to compute the sampling indices with the offsets included, then used interp2() to resample image 2 to align with image 1 s coordinates. I form the interferogram by multiplying image S1 and resampled image R2 together: Tint = S1 R2. I display the phase of the multilook version of the interferogram (4 range looks, 4 azimuth looks), where I am now able to see the fringes. 4
5 We could also overlay the phase image on the amplitude image: 5
6 5. Calculate multi-look (4 range x 4 azimuth) correlation images I used this equation to compute the correlation between images P 1 and P 2 : C = i P 1,iP 2,i i P 1,iP1,i i P 2,iP 2,i (10) The correlation image from Problem 1 has very low correlation values, almost zero. This is not surprising, since image 2 has not been resampled to the same coordinates as image 1. In contrast, the correlation image from Problem 4 provides useful information about the data. The correlation is low over areas of the image with water (as well as some other parts of the image), and I do not see any fringes in the interferogram at these locations. However, the correlation is high in many other parts of the image, and I do see fringes here in the interferogram. % Final project: Part1 clear;clc;close all; set(0, defaultaxesfontsize,16); %% read data tic fid=fopen( image.1 ); data=fread(fid,[3072* ], float ); fclose(fid); 6
7 reald=data(1:2:end,:); imagd=data(2:2:end,:); img1=complex(reald,imagd); img1=img1. ; fid=fopen( image.2 ); data=fread(fid,[3072* ], float ); fclose(fid); reald=data(1:2:end,:); imagd=data(2:2:end,:); img2=complex(reald,imagd); img2=img2. ; display( images loaded ) toc nr=3072; naz=14336; clear data reald imagd %% 1 form an interferogram and take looks intgram=img1.*conj(img2); figure imagesc(angle(intgram)) colormap(gray) colorbar title( interferometric phase ) xlabel( range ) ylabel( azimuth ) saveas(gcf, 1_phase, png ) % average looks nlook=4; intgram_sm=cpxlooks(intgram,nr,naz,nlook,nlook); figure subplot(1,2,1) imagesc(abs(intgram_sm)); colormap gray cax=caxis; caxis(0.4*cax) xlabel( range ) ylabel( azimuth ) title( amplitude ) 7
8 subplot(1,2,2) imagesc(angle(intgram_sm)); colormap gray xlabel( range ) ylabel( azimuth ) title( phase ) colorbar saveas(gcf, 1_withlook, png ) clear intgram %% offset determination: done in computeoffset.m %% resample image.2 to image.1 using a linear interpolater load offsets_coef ind_azimuth=1:naz; ind_range=1:nr; azimi = coefaz(2)*ind_azimuth + coefaz(1); rangei = coefrg(2)*ind_range + coefrg(1); azimi = (azimi+(1:length(azimi))); rangei = (rangei+(1:length(rangei))); azimi(azimi > naz) = naz; azimi(azimi < 1) = 1; rangei(rangei > nr) = nr; rangei(rangei < 1) = 1; % Resample image 2 with new indices [meshrange, meshazim] = meshgrid(rangei, azimi); img2_new = interp2(img2, meshrange, meshazim); display( Resampling of image.2 completed ) %% form interferogram intgram2=img1.*conj(img2_new); % average looks nlook=4; intgram_sm2=cpxlooks(intgram2,nr,naz,nlook,nlook); amp=abs(intgram_sm2); phase=angle(intgram_sm2); scalepic=mycolormap(amp,phase,floor(nr/nlook),floor(naz/nlook)); 8
9 figure imagesc(angle(intgram_sm2)) colormap(jet) colorbar saveas(gcf, 4_phase, png ) figure imagesc(scalepic) axis off xlabel( range ) ylabel( azimuth ) saveas(gcf, 4_withlook, png ) save( intgrams, intgram_sm, intgram_sm2 ) %% correlation cc1=zeros(naz/4,nr/4); cc2=zeros(naz/4,nr/4); % average looks for image.1 and image.2 img1sm=powlook(img1,nr,naz,4,4); img2sm=powlook(img2,nr,naz,4,4); img2rssm=powlook(img2_new,nr,naz,4,4); h = waitbar(0, correlation calculating ); for i=1:naz/4 amp1=abs(intgram_sm(i,:)); amp2=abs(intgram_sm2(i,:)); slc1amp=(img1sm(i,:)).^0.5; slc2amp=(img2sm(i,:)).^0.5; slc3amp=(img2rssm(i,:)).^0.5; cc1(i,:)=amp1./slc1amp./slc2amp; cc2(i,:)=amp2./slc1amp./slc3amp; waitbar(i/naz*4) end close(h) pic1=mycolormap(amp,cc1,nr/4,naz/4); pic2=mycolormap(amp,cc2,nr/4,naz/4); figure subplot(1,2,1) imagesc(pic1) title( Before SLC coregisteration ) colormap jet colorbar subplot(1,2,2) imagesc(pic2) title( After SLC coregisteration ) 9
10 colormap jet colorbar saveas(gcf, correlation.png, png ) % save correlation data as byte files cc1=255*cc1; cc2=255*cc2; fid1=fopen( corr_1, w ); fid2=fopen( corr_2, w ); fwrite(fid1,cc1, uint8 ); fwrite(fid2,cc2, uint8 ); fclose(fid1); fclose(fid2); Function to compute multilooks: cpxlooks (average complex images) % multilook of complex images function imgsm=cpxlooks(imgbg,nr,naz,nlookrg,nlookaz) imgaz=zeros(floor(naz/nlookaz),nr); imgsm=zeros(floor(naz/nlookaz),floor(nr/nlookrg)); for i=1:floor(naz/nlookaz) imgaz(i,:)=sum(imgbg((i-1)*nlookaz+1:i*nlookaz,:)); end for i=1:floor(nr/nlookrg) imgsm(:,i)=sum(imgaz(:,(i-1)*nlookrg+1:i*nlookrg),2); end return Function to compute multilooks: powlooks (average power of complex images) % average power of complex images function imgsm=powlook(imgbg,nr,naz,nlookrg,nlookaz) imgaz=zeros(floor(naz/nlookaz),nr); imgsm=zeros(floor(naz/nlookaz),floor(nr/nlookrg)); powerimg=abs(imgbg).^2; 10
11 for i=1:floor(naz/nlookaz) imgaz(i,:)=sum(powerimg((i-1)*nlookaz+1:i*nlookaz,:)); end for i=1:floor(nr/nlookrg) imgsm(:,i)=sum(imgaz(:,(i-1)*nlookrg+1:i*nlookrg),2); end return Function to overlay phase image on the amplitude image function pic=mycolormap(amp,phase,nr,naz) amplow=prctile(amp(:),5); amphi=prctile(amp(:),95); amp(amp<amplow)=amplow; amp(amp>amphi)=amphi; scale=(amp-amplow)/(amphi-amplow); phlow=prctile(phase(:),1); phhi=prctile(phase(:),99); % create a color table colormap jet; map=colormap; stack=max(phase,phlow); stack=min(phase,phhi); colorstack=round((stack-phlow)/(phhi-phlow)*64); colorstack=max(colorstack,1); colorstack=min(colorstack,64); for k=1:naz for kk=1:nr red(k,kk)=map(colorstack(k,kk),1); green(k,kk)=map(colorstack(k,kk),2); blue(k,kk)=map(colorstack(k,kk),3); end end pic(:,:,1)=red.*scale; pic(:,:,2)=green.*scale; pic(:,:,3)=blue.*scale; Function to compute offset coefficients: % Estimate offsets between image.1 and image.2 11
12 clear; close all; clc; set(0, defaultaxesfontsize, 16); %% radar parameters c=3e8; fs=16e6; hgt=696000; re=6378e3; deltar=c/2/fs; rmid=844768; costheta0=((hgt+re)^2+rmid^2-re^2)/2/rmid/(hgt+re); theta0=acos(costheta0); nr=3072; naz=14336; %% calculate the offsets % use the offset.f program to get the offsets load offset.out rgpos=offset(:,1); rgoffset=offset(:,3); azpos=offset(:,2); azoffset=offset(:,4); meanazoffset=mean(azoffset); display(meanazoffset, Average azimuth offset is, ) meanrgoffset=mean(rgoffset); display(meanrgoffset, Average range offset is, ) % plot range offset vs. range figure subplot(1,2,1) dr=rgpos; scatter(dr,rgoffset) xlabel( range bins ) ylabel( range bin offsets ) subplot(1,2,2) daz=azpos; scatter(daz,azoffset) xlabel( azimuth bins ) ylabel( azimuth bin offsets ) saveas(gcf, offsetdata, tiff ) % get rid of offliers varrg=std(rgoffset); varaz=std(azoffset); num=length(rgpos); 12
13 index=find(rgoffset<meanrgoffset+0.5*varrg & rgoffset>meanrgoffset-0.5*varrg); newrgpos=rgpos(index); newrgoffset=rgoffset(index); index=find(azoffset<meanazoffset+0.5*varaz & azoffset>meanazoffset-0.5*varaz); newazpos=azpos(index); newazoffset=azoffset(index); figure subplot(1,2,1) dr=newrgpos; scatter(dr,newrgoffset) xlabel( range bins ) ylabel( range bin offsets ) subplot(1,2,2) daz=newazpos; scatter(daz,newazoffset) xlabel( azimuth bins ) ylabel( azimuth bin offsets ) % saveas(gcf, offsetdata.filtered, png ) %% estimate range and azimuth offsets as a liner function % get rid of more outliers % least square method YY=newrgoffset; dr=newrgpos; GG=[ones(size(newrgpos)), dr]; coefrg=gg\yy; est_rg=gg*coefrg; YY=newazoffset; daz=newazpos; GG=[ones(size(newazpos)), daz]; coefaz=gg\yy; est_az=gg*coefaz; % plot the estimated offsets figure subplot(1,2,1) scatter(rgpos,rgoffset) hold on plot(dr,est_rg, linewidth,2) hold off xlabel( range bins ) ylabel( range bin offsets ) subplot(1,2,2) scatter(azpos,azoffset) 13
14 hold on plot(daz,est_az, linewidth,2) hold off xlabel( azimuth bins ) ylabel( azimuth bin offsets ) saveas(gcf, offsetdata.filtered.fit.png, png ) %% calculate the baseline components B_par=coefrg(1)*deltar; B_perp=coefrg(2)*rmid*tan(theta0); B=sqrt(B_par^2+B_perp^2); sindiff=b_perp/b; alpha=atan(b_perp/b_par)+theta0-pi/2; alpha_dg=alpha/pi*180; save( offsets_coef, coefrg, coefaz ) 14
GP265 /EE 355 Homework 8 (Final project 2)
GP265 /EE 355 Homework 8 (Final project 2) March 14, 2018 1. Display interferogram phase 2. Calculate the curved earth fringe pattern and subtract it from the interferogram phase. First I compute the slant
More informationA Correlation Test: What were the interferometric observation conditions?
A Correlation Test: What were the interferometric observation conditions? Correlation in Practical Systems For Single-Pass Two-Aperture Interferometer Systems System noise and baseline/volumetric decorrelation
More informationScene Matching on Imagery
Scene Matching on Imagery There are a plethora of algorithms in existence for automatic scene matching, each with particular strengths and weaknesses SAR scenic matching for interferometry applications
More informationResults of UAVSAR Airborne SAR Repeat-Pass Multi- Aperture Interferometry
Results of UAVSAR Airborne SAR Repeat-Pass Multi- Aperture Interferometry Bryan Riel, Ron Muellerschoen Jet Propulsion Laboratory, California Institute of Technology 2011 California Institute of Technology.
More informationIndividual Interferograms to Stacks!
Individual Interferograms to Stacks! Piyush Agram! Jet Propulsion Laboratory!! Jun 29, 2015! @UNAVCO! Thanks to my colleagues from JPL, Caltech, Stanford University and from all over the world for providing
More informationBasic Simulation Lab with MATLAB
Chapter 3: Generation of Signals and Sequences 1. t = 0 : 0.001 : 1; Generate a vector of 1001 samples for t with a value between 0 & 1 with an increment of 0.001 2. y = 0.5 * t; Generate a linear ramp
More informationIndividual Interferograms to Stacks
Individual Interferograms to Stacks Piyush Agram Jet Propulsion Laboratory Aug 1, 2016 @UNAVCO Thanks to my colleagues from JPL, Caltech, Stanford University and from all over the world for providing images
More informationPSI Precision, accuracy and validation aspects
PSI Precision, accuracy and validation aspects Urs Wegmüller Charles Werner Gamma Remote Sensing AG, Gümligen, Switzerland, wegmuller@gamma-rs.ch Contents Aim is to obtain a deeper understanding of what
More informationSAR Interferometry. Dr. Rudi Gens. Alaska SAR Facility
SAR Interferometry Dr. Rudi Gens Alaska SAR Facility 2 Outline! Relevant terms! Geometry! What does InSAR do?! Why does InSAR work?! Processing chain " Data sets " Coregistration " Interferogram generation
More informationPrecise coregistration of Sentinel-1A TOPS data. Heresh Fattahi, Piyush Agram, Mark Simons
Precise coregistration of Sentinel-1A TOPS data Heresh Fattahi, Piyush Agram, Mark Simons Sentinel-1A TOPS Burst N Burst 3 Burst 2 Burst 1 [Prats-Iraola et al, 212] Swath1 Swath2 Swath3 [Sakar et al, 215]
More informationSentinel-1 Toolbox. Interferometry Tutorial Issued March 2015 Updated August Luis Veci
Sentinel-1 Toolbox Interferometry Tutorial Issued March 2015 Updated August 2016 Luis Veci Copyright 2015 Array Systems Computing Inc. http://www.array.ca/ http://step.esa.int Interferometry Tutorial The
More informationFirst TOPSAR image and interferometry results with TerraSAR-X
First TOPSAR image and interferometry results with TerraSAR-X A. Meta, P. Prats, U. Steinbrecher, R. Scheiber, J. Mittermayer DLR Folie 1 A. Meta - 29.11.2007 Introduction Outline TOPSAR acquisition mode
More informationOperational process interferometric for the generation of a digital model of ground Applied to the couple of images ERS-1 ERS-2 to the area of Algiers
Operational process interferometric for the generation of a digital model of ground Applied to the couple of images ERS-1 ERS-2 to the area of Algiers F. Hocine, M.Ouarzeddine, A. elhadj-aissa,, M. elhadj-aissa,,
More informationMotion compensation and the orbit restitution
InSA R Contents Introduction and objectives Pi-SAR Motion compensation and the orbit restitution InSAR algorithm DEM generation Evaluation Conclusion and future work Introduction and Objectives L-band
More informationInterferometry Tutorial with RADARSAT-2 Issued March 2014 Last Update November 2017
Sentinel-1 Toolbox with RADARSAT-2 Issued March 2014 Last Update November 2017 Luis Veci Copyright 2015 Array Systems Computing Inc. http://www.array.ca/ http://step.esa.int with RADARSAT-2 The goal of
More informationMore on Images and Matlab
More on Images and Matlab Prof. Eric Miller elmiller@ece.tufts.edu Fall 2007 EN 74-ECE Image Processing Lecture 3-1 Matlab Data Types Different means of representing numbers depending on what you want
More informationLateral Ground Movement Estimation from Space borne Radar by Differential Interferometry.
Lateral Ground Movement Estimation from Space borne Radar by Differential Interferometry. Abstract S.Sircar 1, 2, C.Randell 1, D.Power 1, J.Youden 1, E.Gill 2 and P.Han 1 Remote Sensing Group C-CORE 1
More informationRadar Coherent Backscatter!
Radar Coherent Backscatter! Pixels in a radar image are a complex phasor representation of the coherent backscatter from the resolution element on the ground and the propagation phase delay! Interferometric
More informationInterferometric processing. Rüdiger Gens
Rüdiger Gens Why InSAR processing? extracting three-dimensional information out of a radar image pair covering the same area digital elevation model change detection 2 Processing chain 3 Processing chain
More informationINSAR QUALITY CONTROL: ANALYSIS OF FIVE YEARS OF CORNER REFLECTOR TIME SERIES
INSAR QUALITY CONTROL: ANALYSIS OF FIVE YEARS OF CORNER REFLECTOR TIME SERIES Petar Marinkovic, Gini Ketelaar, Freek van Leijen, and Ramon Hanssen Delft University of Technology, Delft Institute of Earth
More informationfraction of Nyquist
differentiator 4 2.1.2.3.4.5.6.7.8.9 1 1 1/integrator 5.1.2.3.4.5.6.7.8.9 1 1 gain.5.1.2.3.4.5.6.7.8.9 1 fraction of Nyquist Figure 1. (top) Transfer functions of differential operators (dotted ideal derivative,
More informationBasic MATLAB Intro III
Basic MATLAB Intro III Plotting Here is a short example to carry out: >x=[0:.1:pi] >y1=sin(x); y2=sqrt(x); y3 = sin(x).*sqrt(x) >plot(x,y1); At this point, you should see a graph of sine. (If not, go to
More informationDo It Yourself 8. Polarization Coherence Tomography (P.C.T) Training Course
Do It Yourself 8 Polarization Coherence Tomography (P.C.T) Training Course 1 Objectives To provide a self taught introduction to Polarization Coherence Tomography (PCT) processing techniques to enable
More informationLAB 2: Resampling. Maria Magnusson, 2012 (last update August 2016) with contributions from Katarina Flood, Qingfen Lin and Henrik Turbell
LAB 2: Resampling Maria Magnusson, 2 (last update August 6) with contributions from Katarina Flood, Qingfen Lin and Henrik Turbell Computer Vision Laboratory, Dept. of Electrical Engineering, Linköping
More informationIn addition, the image registration and geocoding functionality is also available as a separate GEO package.
GAMMA Software information: GAMMA Software supports the entire processing from SAR raw data to products such as digital elevation models, displacement maps and landuse maps. The software is grouped into
More informationHow to learn MATLAB? Some predefined variables
ECE-S352 Lab 1 MATLAB Tutorial How to learn MATLAB? 1. MATLAB comes with good tutorial and detailed documents. a) Select MATLAB help from the MATLAB Help menu to open the help window. Follow MATLAB s Getting
More informationIntegrated Algebra 2 and Trigonometry. Quarter 1
Quarter 1 I: Functions: Composition I.1 (A.42) Composition of linear functions f(g(x)). f(x) + g(x). I.2 (A.42) Composition of linear and quadratic functions II: Functions: Quadratic II.1 Parabola The
More informationGAMMA Technical Report: Offset estimation programs update
GAMMA Technical Report: Offset estimation programs update Christophe Magnard, Charles Werner, and Urs Wegmüller Gamma Remote Sensing, Worbstrasse 225, CH-3073 Gümligen, Switzerland Nov 2017 TABLE OF CONTENTS
More informationConcept and methodology of SAR Interferometry technique
Concept and methodology of SAR Interferometry technique March 2016 Differen;al SAR Interferometry Young s double slit experiment - Construc;ve interference (bright) - Destruc;ve interference (dark) http://media-2.web.britannica.com/eb-media/96/96596-004-1d8e9f0f.jpg
More informationThe STUN algorithm for Persistent Scatterer Interferometry
[1/27] The STUN algorithm for Persistent Scatterer Interferometry Bert Kampes, Nico Adam 1. Theory 2. PSIC4 Processing 3. Conclusions [2/27] STUN Algorithm Spatio-Temporal Unwrapping Network (STUN) 4 1D
More informationDesign, Implementation and Performance Evaluation of Synthetic Aperture Radar Signal Processor on FPGAs
Design, Implementation and Performance Evaluation of Synthetic Aperture Radar Signal Processor on FPGAs Hemang Parekh Masters Thesis MS(Computer Engineering) University of Kansas 23rd June, 2000 Committee:
More informationPixel-Offset SBAS analysis: a tool for the investigation of deformation time-series with large dynamics
Pixel-Offset SBAS analysis: a tool for the investigation of deformation time-series with large dynamics F. Casu, A. Manconi, A. Pepe, M. Manzo, R. Lanari IREA-CNR Napoli, Italy casu.f@irea.cnr.it Summary
More informationSolutions For Homework #8
Solutions For Homework #8 Problem 1:[1 pts] (a) In this problem, the circular aperture functions representing Polyphemus and Odysseus eyes are expressed as functions of the unitless independent variable
More informationMATH 2221A Mathematics Laboratory II
MATH A Mathematics Laboratory II Lab Assignment 4 Name: Student ID.: In this assignment, you are asked to run MATLAB demos to see MATLAB at work. The color version of this assignment can be found in your
More informationSAOCOM 1A INTERFEROMETRIC ERROR MODEL AND ANALYSIS
SAOCOM A INTERFEROMETRIC ERROR MODEL AND ANALYSIS Pablo Andrés Euillades (), Leonardo Daniel Euillades (), Mario Azcueta (), Gustavo Sosa () () Instituto CEDIAC FI UNCuyo & CONICET, Centro Universitario,
More informationROI_PAC Documentation. Repeat Orbit Interferometry Package
ROI_PAC Documentation Repeat Orbit Interferometry Package Abridged version of Chapter 3 of a PhD thesis written by Sean Buckley CSR, UT Austin; edited by Paul Rosen and Patricia Persaud May 10, 2000 Version
More informationTopic 0b Graphics for Science & Engineering
Course Instructor Dr. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: rcrumpf@utep.edu Topic 0b Graphics for Science & Engineering EE 4386/5301 Computational Methods in EE Outline What are
More informationFiltering, unwrapping, and geocoding R. Mellors
Filtering, unwrapping, and geocoding R. Mellors or what to do when your interferogram looks like this correlation Haiti ALOS L-band (23 cm) ascend T447, F249 3/9/09-1/25/10 azimuth phase Bperp = 780 (gmtsar)
More informationECE251DN: Homework #3 Solutions
ECE251DN: Homework #3 Solutions Problem 3.7.2 (a) In this problem, we only have one null constraint. So N 1 j C = V(ψ ) = [e 2 ψ,..., 1,..., e j N 1 2 ψ ] T The weights of the least squares approximation
More informationInterferometric Synthetic-Aperture Radar (InSAR) Basics
Interferometric Synthetic-Aperture Radar (InSAR) Basics 1 Outline SAR limitations Interferometry SAR interferometry (InSAR) Single-pass InSAR Multipass InSAR InSAR geometry InSAR processing steps Phase
More informationmatlab_intro.html Page 1 of 5 Date: Tuesday, September 6, 2005
matlab_intro.html Page 1 of 5 % Introducing Matlab % adapted from Eero Simoncelli (http://www.cns.nyu.edu/~eero) % and Hany Farid (http://www.cs.dartmouth.edu/~farid) % and Serge Belongie (http://www-cse.ucsd.edu/~sjb)
More informationPC-MATLAB PRIMER. This is intended as a guided tour through PCMATLAB. Type as you go and watch what happens.
PC-MATLAB PRIMER This is intended as a guided tour through PCMATLAB. Type as you go and watch what happens. >> 2*3 ans = 6 PCMATLAB uses several lines for the answer, but I ve edited this to save space.
More informationMOTION COMPENSATION OF INTERFEROMETRIC SYNTHETIC APERTURE RADAR
MOTION COMPENSATION OF INTERFEROMETRIC SYNTHETIC APERTURE RADAR David P. Duncan Microwave Earth Remote Sensing Laboratory Brigham Young University Provo, UT 84602 PH: 801.422.4884, FAX: 801.422.6586 April
More informationUnit Maps: Grade 8 Math
Real Number Relationships 8.3 Number and operations. The student represents and use real numbers in a variety of forms. Representation of Real Numbers 8.3A extend previous knowledge of sets and subsets
More informationICE VELOCITY measurements are fundamentally important
102 IEEE GEOSCIENCE AND REMOTE SENSING LETTERS, VOL. 4, NO. 1, JANUARY 2007 Synergistic Fusion of Interferometric and Speckle-Tracking Methods for Deriving Surface Velocity From Interferometric SAR Data
More informationURBAN FOOTPRINT MAPPING WITH SENTINEL-1 DATA
URBAN FOOTPRINT MAPPING WITH SENTINEL-1 DATA Data: Sentinel-1A IW SLC 1SSV: S1A_IW_SLC 1SSV_20160102T005143_20160102T005208_009308_00D72A_849D S1A_IW_SLC 1SSV_20160126T005142_20160126T005207_009658_00E14A_49C0
More informationWIDE AREA DEFORMATION MAP GENERATION WITH TERRASAR-X DATA: THE TOHOKU-OKI EARTHQUAKE 2011 CASE
WIDE AREA DEFORMATION MAP GENERATION WITH TERRASAR-X DATA: THE TOHOKU-OKI EARTHQUAKE 2011 CASE Nestor Yague-Martinez (1), Michael Eineder (2), Christian Minet (2), Birgitt Schättler (2) (1) Starlab Barcelona
More informationNew Results on the Omega-K Algorithm for Processing Synthetic Aperture Radar Data
New Results on the Omega-K Algorithm for Processing Synthetic Aperture Radar Data Matthew A. Tolman and David G. Long Electrical and Computer Engineering Dept. Brigham Young University, 459 CB, Provo,
More informationINTERNATIONAL EDITION. MATLAB for Engineers. Third Edition. Holly Moore
INTERNATIONAL EDITION MATLAB for Engineers Third Edition Holly Moore 5.4 Three-Dimensional Plotting Figure 5.8 Simple mesh created with a single two-dimensional matrix. 5 5 Element,5 5 The code mesh(z)
More informationInterferometric Evaluation of Sentinel-1A TOPS data
Interferometric Evaluation of Sentinel-1A TOPS data N. Yague-Martinez, F. Rodriguez Gonzalez, R. Brcic, R. Shau Remote Sensing Technology Institute. DLR, Germany ESTEC/Contract No. 4000111074/14/NL/MP/lf
More informationInterferometry Module for Digital Elevation Model Generation
Interferometry Module for Digital Elevation Model Generation In order to fully exploit processes of the Interferometry Module for Digital Elevation Model generation, the European Space Agency (ESA) has
More informationSentinel-1 Toolbox. TOPS Interferometry Tutorial Issued May 2014
Sentinel-1 Toolbox TOPS Interferometry Tutorial Issued May 2014 Copyright 2015 Array Systems Computing Inc. http://www.array.ca/ https://sentinel.esa.int/web/sentinel/toolboxes Interferometry Tutorial
More informationGBT Commissioning Memo 11: Plate Scale and pointing effects of subreflector positioning at 2 GHz.
GBT Commissioning Memo 11: Plate Scale and pointing effects of subreflector positioning at 2 GHz. Keywords: low frequency Gregorian, plate scale, focus tracking, pointing. N. VanWey, F. Ghigo, R. Maddalena,
More informationAP Calculus Summer Review Packet
AP Calculus Summer Review Packet Name: Date began: Completed: **A Formula Sheet has been stapled to the back for your convenience!** Email anytime with questions: danna.seigle@henry.k1.ga.us Complex Fractions
More information2.2 Limit of a Function and Limit Laws
Limit of a Function and Limit Laws Section Notes Page Let s look at the graph y What is y()? That s right, its undefined, but what if we wanted to find the y value the graph is approaching as we get close
More informationMotion Compensation of Interferometric Synthetic Aperture Radar
Brigham Young University BYU ScholarsArchive All Theses and Dissertations 2004-07-07 Motion Compensation of Interferometric Synthetic Aperture Radar David P. Duncan Brigham Young University - Provo Follow
More informationIDENTIFICATION OF THE LOCATION PHASE SCREEN OF ERS-ENVISAT PERMANENT SCATTERERS
IDENTIFICATION OF THE LOCATION PHASE SCREEN OF ERS-ENVISAT PERMANENT SCATTERERS M. Arrigoni (1), C. Colesanti (1), A. Ferretti (2), D. Perissin (1), C. Prati (1), F. Rocca (1) (1) Dipartimento di Elettronica
More informationInSAR Processing. Sentinel 1 data Case study of subsidence in Mexico city. Marie-Pierre Doin, Cécile Lasserre, Raphaël Grandin, Erwan Pathier
1 InSAR Processing Sentinel 1 data Case study of subsidence in Mexico city Marie-Pierre Doin, Cécile Lasserre, Raphaël Grandin, Erwan Pathier NSBAS processing chain (based on ROI_PAC): ROI-PAC: Rosen et
More informationBasic plotting commands Types of plots Customizing plots graphically Specifying color Customizing plots programmatically Exporting figures
Basic plotting commands Types of plots Customizing plots graphically Specifying color Customizing plots programmatically Exporting figures Matlab is flexible enough to let you quickly visualize data, and
More informationRadiometric Calibration of S-1 Level-1 Products Generated by the S-1 IPF
Radiometric Calibration of S-1 Level-1 Products Generated by the S-1 IPF Prepared by Nuno Miranda, P.J. Meadows Reference ESA-EOPG-CSCOP-TN-0002 Issue 1 Revision 0 Date of Issue 21/05/2015 Status Final
More informationSignals and Systems Profs. Byron Yu and Pulkit Grover Fall Homework 1
18-290 Signals and Systems Profs. Byron Yu and Pulkit Grover Fall 2018 Homework 1 This homework is due in class on Thursday, September 6, 9:00am. Instructions Solve all non-matlab problems using only paper
More information34 Dead-Reckoning Error
34 DEAD-RECKONING ERROR 124 34 Dead-Reckoning Error An unmanned, untethered underwater vehicle operating near large metallic marine structures has to navigate without a magnetic compass, and in some cases
More informationENGR Fall Exam 1
ENGR 1300 Fall 01 Exam 1 INSTRUCTIONS: Duration: 60 minutes Keep your eyes on your own work! Keep your work covered at all times! 1. Each student is responsible for following directions. Read carefully..
More informationProblems of Plane analytic geometry
1) Consider the vectors u(16, 1) and v( 1, 1). Find out a vector w perpendicular (orthogonal) to v and verifies u w = 0. 2) Consider the vectors u( 6, p) and v(10, 2). Find out the value(s) of parameter
More informationAPPM 2460 PLOTTING IN MATLAB
APPM 2460 PLOTTING IN MATLAB. Introduction Matlab is great at crunching numbers, and one of the fundamental ways that we understand the output of this number-crunching is through visualization, or plots.
More informationECE 3793 Matlab Project 1
ECE 3793 Matlab Project 1 Spring 2017 Dr. Havlicek DUE: 02/04/2017, 11:59 PM Introduction: You will need to use Matlab to complete this assignment. So the first thing you need to do is figure out how you
More informationRepeat-pass SAR Interferometry Experiments with Gaofen-3: A Case Study of Ningbo Area
Repeat-pass SAR Interferometry Experiments with Gaofen-3: A Case Study of Ningbo Area Tao Zhang, Xiaolei Lv, Bing Han, Bin Lei and Jun Hong Key Laboratory of Technology in Geo-spatial Information Processing
More informationDepth estimation from stereo image pairs
Depth estimation from stereo image pairs Abhranil Das In this report I shall first present some analytical results concerning depth estimation from stereo image pairs, then describe a simple computational
More informationEE 168 Introduction to Digital Image Processing January 30, 2012 HOMEWORK 2 SOLUTIONS
EE 168 Introduction to Digital Image Processing January 3, 212 HOMEWORK 2 SOLUTIONS Problem 1: Image Scaling We determine the dimensions of the original image (read from the rodin.raw file) to be 348x261.
More informationUnit Maps: Grade 8 Math
Real Number Relationships 8.3 Number and operations. The student represents and use real numbers in a variety of forms. Representation of Real Numbers 8.3A extend previous knowledge of sets and subsets
More informationEE 301 Signals & Systems I MATLAB Tutorial with Questions
EE 301 Signals & Systems I MATLAB Tutorial with Questions Under the content of the course EE-301, this semester, some MATLAB questions will be assigned in addition to the usual theoretical questions. This
More informationInterferometric SAR Processing
Documentation - Theory Interferometric SAR Processing Version 1.0 November 2007 GAMMA Remote Sensing AG, Worbstrasse 225, CH-3073 Gümligen, Switzerland tel: +41-31-951 70 05, fax: +41-31-951 70 08, email:
More informationBackground and Accuracy Analysis of the Xfactor7 Table: Final Report on QuikScat X Factor Accuracy
Brigham Young University Department of Electrical and Computer Engineering 9 Clyde Building Provo, Utah 86 Background and Accuracy Analysis of the Xfactor7 Table: Final Report on QuikScat X Factor Accuracy
More information3 - SYNTHETIC APERTURE RADAR (SAR) SUMMARY David Sandwell, SIO 239, January, 2008
1 3 - SYNTHETIC APERTURE RADAR (SAR) SUMMARY David Sandwell, SIO 239, January, 2008 Fraunhoffer diffraction To understand why a synthetic aperture in needed for microwave remote sensing from orbital altitude
More informationLecture: Edge Detection
CMPUT 299 Winter 2007 Lecture: Edge Detection Irene Cheng Overview. What is a pixel in an image? 2. How does Photoshop, + human assistance, detect an edge in a picture/photograph? 3. Behind Photoshop -
More information1 >> Lecture 4 2 >> 3 >> -- Graphics 4 >> Zheng-Liang Lu 184 / 243
1 >> Lecture 4 >> 3 >> -- Graphics 4 >> Zheng-Liang Lu 184 / 43 Introduction ˆ Engineers use graphic techniques to make the information easier to understand. ˆ With graphs, it is easy to identify trends,
More informationLAB 1: Introduction to MATLAB Summer 2011
University of Illinois at Urbana-Champaign Department of Electrical and Computer Engineering ECE 311: Digital Signal Processing Lab Chandra Radhakrishnan Peter Kairouz LAB 1: Introduction to MATLAB Summer
More informationRADARGRAMMETRY AND INTERFEROMETRY SAR FOR DEM GENERATION
RADARGRAMMETRY AND INTERFEROMETRY SAR FOR DEM GENERATION Jung Hum Yu 1, Xiaojing Li, Linlin Ge, and Hsing-Chung Chang School of Surveying and Spatial Information Systems University of New South Wales,
More informationEE119 Homework 3. Due Monday, February 16, 2009
EE9 Homework 3 Professor: Jeff Bokor GSI: Julia Zaks Due Monday, February 6, 2009. In class we have discussed that the behavior of an optical system changes when immersed in a liquid. Show that the longitudinal
More informationPRECALCULUS MATH Trigonometry 9-12
1. Find angle measurements in degrees and radians based on the unit circle. 1. Students understand the notion of angle and how to measure it, both in degrees and radians. They can convert between degrees
More informationLogical Subscripting: This kind of subscripting can be done in one step by specifying the logical operation as the subscripting expression.
What is the answer? >> Logical Subscripting: This kind of subscripting can be done in one step by specifying the logical operation as the subscripting expression. The finite(x)is true for all finite numerical
More informationRESOLUTION enhancement is achieved by combining two
IEEE GEOSCIENCE AND REMOTE SENSING LETTERS, VOL. 3, NO. 1, JANUARY 2006 135 Range Resolution Improvement of Airborne SAR Images Stéphane Guillaso, Member, IEEE, Andreas Reigber, Member, IEEE, Laurent Ferro-Famil,
More informationAMS 27L LAB #2 Winter 2009
AMS 27L LAB #2 Winter 2009 Plots and Matrix Algebra in MATLAB Objectives: 1. To practice basic display methods 2. To learn how to program loops 3. To learn how to write m-files 1 Vectors Matlab handles
More informationA New Method for Correcting ScanSAR Scalloping Using Forests and inter SCAN Banding Employing Dynamic Filtering
A New Method for Correcting ScanSAR Scalloping Using Forests and inter SCAN Banding Employing Dynamic Filtering Masanobu Shimada Japan Aerospace Exploration Agency (JAXA), Earth Observation Research Center
More informationENGR Fall Exam 1 PRACTICE EXAM
ENGR 13100 Fall 2012 Exam 1 PRACTICE EXAM INSTRUCTIONS: Duration: 60 minutes Keep your eyes on your own work! Keep your work covered at all times! 1. Each student is responsible for following directions.
More informationPrecise orbits and accurate timing simplifies software and enables seamless mosaicing. Geometric validation of ERS, Envisat, and ALOS.
Geometric Calibration of GMTSAR Processing Software Using Corner Reflectors at Pinon Flat David Sandwell, UCSD/SIO CEOS Cal/Val Workshop, November 8, 2011 What is GMTSAR? Precise orbits and accurate timing
More informationIntegrated Math 1: Homework #2 Answers (Day 1)
Integrated Math 1: Homework #1 Answers 1-6. a: y = x 2 6 and then y = x 5 b: Yes, reverse the order of the machines ( y = x 5 and then y = x 2 6 ) and use an input of x = 6. 1-7. a: 54 b: 7 3 5 c: 2 d:
More informationPROGRAMMING WITH MATLAB WEEK 6
PROGRAMMING WITH MATLAB WEEK 6 Plot: Syntax: plot(x, y, r.- ) Color Marker Linestyle The line color, marker style and line style can be changed by adding a string argument. to select and delete lines
More informationThe 2017 InSAR package also provides support for the generation of interferograms for: PALSAR-2, TanDEM-X
Technical Specifications InSAR The Interferometric SAR (InSAR) package can be used to generate topographic products to characterize digital surface models (DSMs) or deformation products which identify
More informationCoherent Gradient Sensing Microscopy: Microinterferometric Technique. for Quantitative Cell Detection
Coherent Gradient Sensing Microscopy: Microinterferometric Technique for Quantitative Cell Detection Proceedings of the SEM Annual Conference June 7-10, 010 Indianapolis, Indiana USA 010 Society for Experimental
More informationChapter 11. Above: Principal contraction rates calculated from GPS velocities. Visualized using MATLAB.
Chapter 11 Above: Principal contraction rates calculated from GPS velocities. Visualized using MATLAB. We have used MATLAB to visualize data a lot in this course, but we have only scratched the surface
More informationData and Function Plotting with MATLAB (Linux-10)
Data and Function Plotting with MATLAB (Linux-10) This tutorial describes the use of MATLAB for general plotting of experimental data and equations and for special plots like histograms. (Astronomers -
More informationChapter 8 Complex Numbers & 3-D Plots
EGR115 Introduction to Computing for Engineers Complex Numbers & 3-D Plots from: S.J. Chapman, MATLAB Programming for Engineers, 5 th Ed. 2016 Cengage Learning Topics Introduction: Complex Numbers & 3-D
More informationCCSSM Curriculum Analysis Project Tool 1 Interpreting Functions in Grades 9-12
Tool 1: Standards for Mathematical ent: Interpreting Functions CCSSM Curriculum Analysis Project Tool 1 Interpreting Functions in Grades 9-12 Name of Reviewer School/District Date Name of Curriculum Materials:
More informationHP-35s Calculator Program Radiations 1
Calculate a Programmer: Dr. Bill Hazelton Date: March, 2008. Line Instruction Display User Instructions O001 LBL O LBL O O002 CLSTK CLEAR 5 O003 FS? 10 FLAGS 3.0 O004 GTO O008 O005 SF 1 FLAGS 1 1 O006
More informationGAMMA Interferometric Point Target Analysis Software (IPTA): Users Guide
GAMMA Interferometric Point Target Analysis Software (IPTA): Users Guide Contents User Handbook Introduction IPTA overview Input data Point list generation SLC point data Differential interferogram point
More informationDEM-BASED SAR PIXEL AREA ESTIMATION FOR ENHANCED GEOCODING REFINEMENT AND RADIOMETRIC NORMALIZATION.
DEM-BASED SAR PIXEL AREA ESTIMATION FOR ENHANCED GEOCODING REFINEMENT AND RADIOMETRIC NORMALIZATION Othmar Frey (1), Maurizio Santoro (2), Charles L. Werner (2), and Urs Wegmuller (2) (1) Gamma Remote
More informationBME I5000: Biomedical Imaging
BME I5000: Biomedical Imaging Lecture 1 Introduction Lucas C. Parra, parra@ccny.cuny.edu 1 Content Topics: Physics of medial imaging modalities (blue) Digital Image Processing (black) Schedule: 1. Introduction,
More informationChapter 37. Interference of Light Waves
Chapter 37 Interference of Light Waves Wave Optics Wave optics is a study concerned with phenomena that cannot be adequately explained by geometric (ray) optics These phenomena include: Interference Diffraction
More informationSentinel-1 processing with GAMMA software
Documentation User s Guide Sentinel-1 processing with GAMMA software Including an example of Sentinel-1 SLC co-registration and differential interferometry Version 1.1 May 2015 GAMMA Remote Sensing AG,
More information