GP265 /EE 355 Homework 8 (Final project 2)
|
|
- Franklin Welch
- 6 years ago
- Views:
Transcription
1 GP265 /EE 355 Homework 8 (Final project 2) March 14, Display interferogram phase 2. Calculate the curved earth fringe pattern and subtract it from the interferogram phase. First I compute the slant range pixel spacing r: c r = nlooks r = m/s 2f s 2( Hz) = m. (1) Then the range ρ to each range bin ir is: ρ(ir) = ρ 0 + (ir 1) r, ir = 0, 1,..., nrange 1 (2) where ρ 0 = m is the range to the first bin and nrange = 1024 range bins. 1
2 The look angle θ z0 to each range bin for z = 0 (no topography) is: ( (h + θ z0 = cos 1 re ) 2 + ρ 2 (z + r e ) 2 ) ( (h + z=0 = cos 1 re ) 2 + ρ 2 r 2 ) e 2ρ(h + r e ) 2ρ(h + r e ) where h = m is the radar altitude and r e = m is the earth radius. Then the curved earth phase ϕ z0 for z = 0 (no topography) is: ϕ z0 = 4π λ ( ) B 2 2ρ B sin (θ z0 α) (3) (4) where λ = m is the radar wavelength, B = 307 m is the baseline length, and α is the orientation angle. To remove this curved earth phase from the complex interferogram S (displayed in Problem 1), I apply a phase correction T = Se iϕ z0, and here is the resulting phase arg(t ), which should have only topography effects. The fringe lines look more like elevation contour lines here: 3. Flatten the interferogram by further correcting for any tilts over the image. We assume that the topography is flat over the caldera at the top of the volcano, at range bins and azimuth bins I used least squares to find the slope of the fringes within the caldera (in radians/pixel), in both the range and azimuth directions, and apply 2
3 these phase corrections to the interferogram from Problem 2 to flatten the tilt from the interferogram. The flattened interferogram no longer has fringes within the caldera, and the fringes look more like elevation contour lines around the volcano. Note that the area I chose is small enough that we can avoid the problem of phase unwrapping. This is not the only way to correct for tilts - many other solutions are acceptable. For example, you can also try guessing different values of orientation angle α that remove the tilt within the caldera. 4. Determine the ambiguity height, ha, for this geometry. The ambiguity height is defined as the height dz where the topography phase ϕtopo = 2π ϕtopo = 2π = 4π (Bsin(θ α) Bsin(θf lat α)) λ (5) We can compute θf lat using the law of cosines using Eq. (3) and ha can be computed once we obtain the look angle θ from Eq. (5), again using the law of cosines. The ambiguity height deps on the range ρ, so here is a plot of the range depence: 3
4 At the reference (first) range bin, ha = m. At the last range bin, ha = m. 5. Unwrap the flattened interferogram using the snaphu program. I write out the flattened interferogram to a file called flatinterferogram and then run it through snaphu to unwrap the phase, writing it out to a file called snaphu.out. I then read in the snaphu.out file into MATLAB, and plot its phase, which is now unwrapped - no more fringes here, and the phase is no longer limited to the interval [ π, π]. 4
5 6. Using the ambiguity height h a from (4), convert the unwrapped phase values to heights. You have now created a digital elevation model (DEM) in radar coordinates. I use the ambiguity height h a from Eq. (5), which varies as a function of range ρ, to scale the unwrapped phase values ϕ flat from Problem 5 into height values z DEM. I choose several zero elevation reference pixels near the coast, and subtract their average phase ϕ ref from every point in the unwrapped interferogram. Since we are dealing with topography here, I choose to plot the contour of the DEM. z DEM = (ϕ flat ϕ ref ) h a (6) 2π The elevation of Mauna Loa is 4169 m, so our DEM looks reasonable. There are some heights below 0 m, but they were in the region over water so I would not expect to have useful height values here. 7. Map correlation into height error First, I display the correlation to check that I read it in properly. 5
6 Problem 7: Correlation azimuth bin range bin I use this equation to compute the height error σz as a function of correlation ρc (range is ρ): λρ sin θz0 1 ρc σz = (7) 4πB cos(θz0 α) 2ρc Here is a plot of the height error σz, ranging from 0 to 50 m. It makes sense that the height error is greatest at the locations where the correlation is lowest, such as the water in the top-right corner, as well as the slopes in the top-left and bottom-right corners. 6
7 8. Map DEM to ground coordinates To compute the ground range, which deps on both range and azimuth since I take the DEM heights into account, I first compute the look angle (compare to Eq. (3)): ( ) (h + re ) + ρ (zdem + re ) θz = cos (8) 2ρ(h + re ) Then I can compute the angle βz at the center of the earth, using the law of sines: ) ( ρ 1 βz = sin sin θz re (9) The ground range rg,z is then rg,z = re βz (10) I get the platform velocity v from the effective platform velocity vef f = m/s, radar altitude h = m, and earth radius re = m: h + re v = vef f = ( m/s) = m/s. (11) re Then the azimuth bin ground spacing azg is azg = nlooksaz v re = 12 = m prf h + re (12)
8 I generated a ground coordinates grid with a spacing of az g between adjacent ground range bins, in order to get square pixels. I resampled only in range, not in azimuth. There are some holes in the image from resampling. 9. Map DEM to ground coordinates: Nice image with interpolation I resampled the new DEM using linear interpolation. There are no longer holes in the image. 8
9 % Hw8 Final Pt.2 clear; close all; clc; %% parameters nr=1024; naz=1190; nrlook=3; nazlook=12; B=307; % baseline,m B_par=-170; B_perp=-263; hgt=696000; % altitude in orbit,m wvl= ; % wavelength, m prf= ; fs=16e6; c=3e8; veff=7179.4; r0=844768; % near range deltar=c/2/fs*nrlook; rmid=r0+(nr/2)*deltar+deltar/2; % mid range re=6378e3; % Earth radius,m costheta0=((hgt+re)^2+rmid^2-re^2)/2/rmid/(hgt+re); theta0=acos(costheta0); alpha=atan(b_perp/b_par)+theta0-pi/2; % baseline angle, rad alphadg=alpha/pi*180; % baseline angle, degree %% read files fid=fopen( interferogram ); data=fread(fid,[nr*2 naz], float ); fclose(fid); reald=data(1:2:,:); imagd=data(2:2:,:); intgram=complex(reald,imagd); intgram=intgram. ; phase=angle(intgram); fid=fopen( amplitude ); data=fread(fid,[nr*2 inf], float ); fclose(fid); amp1=data(1:2:,:); amp2=data(2:2:,:); amp=amp1. ; fid=fopen( correlation ); data=fread(fid,[nr*2 inf], float ); fclose(fid); col=data(nr+1:,:); col=col. ; % display the phase of the interferogram 9
10 imagesc(phase) axis image saveas(gcf, 1.png, png ) %% calculate the curved Earth pattern and substract it. phase_bg=zeros(naz,nr); for i=1:nr rg=r0+(i-1)*deltar; costhetaflat=((hgt+re)^2+rg^2-re^2)/2/rg/(hgt+re); thetaflat=acos(costhetaflat); sindiff=sin(thetaflat-alpha); phase_bg(:,i)=-4*pi/wvl*(b^2/2/rg-b*sindiff); intgram_bg=exp(1j*phase_bg); imagesc(angle(intgram_bg)) colormap(jet) colorbar saveas(gcf, 2_2.png, png ) topo=intgram.*conj(intgram_bg); imagesc(angle(topo)) saveas(gcf, 2.png, png ) % pic=mycolormap(amp,angle(topo),nr,naz); % % imagesc(pic) %% flatten the interferogram by further correcting for any tilts over the image % select a small flat area where we expect no change in phase % note that the area should be small enough that there is no phase wrapping testsite=topo(925:960,110:150); testphase=angle(testsite); imagesc(testphase) colorbar % fit a 2-D plane [valid_x,valid_y]=meshgrid(1:size(testsite,2),1:size(testsite,1)); valid_x=valid_x(:); valid_y=valid_y(:); valid_data=testphase(:); design_a = [ones(size(valid_data)) valid_x valid_y]; coef = design_a\valid_data; ramppart=design_a*coef; imagesc(testphase-reshape(ramppart,size(testsite))) % get the ramp for the whole interferogram [all_x,all_y] = meshgrid (1:nr, 1:naz); 10
11 all_x1=all_x(:); all_y1=all_y(:); full_a=[ones(nr*naz,1) all_x1 all_y1 ]; ramp=full_a*coef; ramp=reshape(ramp,naz,nr); ramp=exp(1j*ramp); topo_detilt=topo.*conj(ramp); pic=mycolormap(amp,angle(topo_detilt),nr,naz); imagesc(pic) xlabel( range ) ylabel( azimuth ) title( corrected interferogram ) saveas(gcf, 3.png, png ) %% determine the ambiguity height as a function of range ha=zeros(nr,1); for i=1:nr rg=r0+(i-1)*deltar; costhetaflat=((hgt+re)^2+rg^2-re^2)/2/rg/(hgt+re); thetaflat=acos(costhetaflat); dffsin=wvl/2/b; sindifftopo=dffsin+sin(thetaflat-alpha); thetatopo=asin(sindifftopo)+alpha; ha(i)=sqrt((hgt+re)^2+rg^2-2*rg*(hgt+re)*cos(thetatopo))-re; ha(1) ha() plot(1:nr,ha, linewidth,2) grid on xlabel( range bin ) ylabel( ambiguity height (m) ) saveas(gcf, 4.png, png ) %% unwrap the flattened interferogram fid=fopen( intgram.uw, r ); dat=fread(fid,[2*nr,inf], float, ieee-le ); phase_uw=dat(nr+1:,:) ; imagesc(phase_uw) pic=mycolormap(amp,phase_uw,nr,naz); imagesc(pic) colorbar phaselo=prctile(phase_uw(:),0); phasehi=prctile(phase_uw(:),99); 11
12 caxis([phaselo,phasehi]) saveas(gcf, 5.png, png ) %% construct a DEM model % choose a reference point where height is zero and correlation is good. % % imagesc(col(1:200,500:1000)); rgref=[550,693,792,974]; azref=[31,94,167,159]; phref=0; for i=1:4 phref=phase_uw(azref(i),rgref(i))+phref; phref=mean(phref); DEM=nan(size(phase_uw)); for i=1:nr DEM(:,i)=(phase_uw(:,i)-phref)/2/pi*ha(i); DEM(DEM<0)=0; contourf(dem,40) colorbar ch=colorbar; ylabel(ch, Height,m ); set(gca, Ydir, reverse ) axis image saveas(gcf, 6.png, png ) %% map of standard deviation of map stdphase=((1-col)./2./col).^0.5; stddem=nan(size(stdphase)); for i=1:nr stddem(:,i)=ha(i)*stdphase(:,i)/2/pi; imagesc(stddem) caxis([0 50]) axis image colorbar ch=colorbar; ylabel(ch, Height Uncertainty,m ); title( uncertainty map ) saveas(gcf, 7.png, png ) %% rectify your map into ground range and azimuth coordinate vel=veff*sqrt((hgt+re)/re); deltaaz=vel/prf*re/(re+hgt)*nazlook; 12
13 % determine the nearest ground range demline=dem(:,1); cosgama=((hgt+re)^2-r0^2+(re+demline).^2)./2/(hgt+re)./(demline+re); gama=acos(cosgama); rnear=min(re*gama); % determine the farest ground range rmax=r0+(nr-1)*deltar; demline=dem(:,nr-1); cosgama=((hgt+re)^2-rmax^2+(re+demline).^2)./2/(hgt+re)./(demline+re); gama=acos(cosgama); rfar=max(re*gama); nrg=round((rfar-rnear)/deltaaz); newdem=nan(naz,nrg); % determine the offsets in range (pixels) and reesample rgoffset=zeros(naz,nr); for i=1:nr rsl=r0+(i-1)*deltar; demline=dem(:,i); cosgama=((hgt+re)^2-rsl^2+(re+demline).^2)./2/(hgt+re)./(demline+re); gama=acos(cosgama); rgr=re*gama; % ground range idxrgnew=((rgr-rnear)/deltaaz+1); % range index in new coordinates rgoffset(:,i)=idxrgnew; h = waitbar(0, Range resampling ); for j=1:naz idxrgnew=rgoffset(j,:); for i=1:nrg idxintper=find((abs(idxrgnew-i)<1)); if ~isempty(idxintper) newdem(j,i)= mean(dem(j,idxintper)); waitbar(j/naz) close(h) contourf(newdem,40) colorbar ch=colorbar; ylabel(ch, Height,m ); axis image set(gca, Ydir, reverse ) axis off saveas(gcf, 8.png, png ) 13
14 %% interpolate to obtain good looking pictures % linear interpolation [idxaz,idxrg]=find(isnan(newdem)); h = waitbar(0, Interpolating ); cull_window=1; DEM_interp=newDEM; for i=1:length(idxaz) rneari=min(rgoffset(idxaz(i),:)); rfari=max(rgoffset(idxaz(i),:)); if idxrg(i)>rneari && idxrg(i)<rfari index_az = max(1, idxaz(i)-cull_window):min(idxaz(i)+cull_window, naz); index_rg= max(1, idxrg(i)-cull_window):min(idxrg(i)+cull_window, nr); patch=newdem(index_az,index_rg); patch=patch(:); patch=patch(isfinite(patch)); DEM_interp(idxaz(i),idxrg(i))=mean(patch(:)); waitbar(i/length(idxaz)) close(h) contourf(dem_interp,40) colorbar ch=colorbar; ylabel(ch, Height,m ); set(gca, Ydir, reverse ) axis image axis off saveas(gcf, 9.png, png ) 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 ; map=colormap; stack=max(phase,phlow); stack=min(phase,phhi); colorstack=round((stack-phlow)/(phhi-phlow)*64); 14
15 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); pic(:,:,1)=red.*scale; pic(:,:,2)=green.*scale; pic(:,:,3)=blue.*scale; 15
GP265 / EE355 Homework 7 (Final project 1) Solutions
GP265 / EE355 Homework 7 (Final project 1) Solutions 2018-03-10 1. 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
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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationSynthetic Aperture Radar Interferometry (InSAR)
CEE 6100 / CSS 6600 Remote Sensing Fundamentals 1 Synthetic Aperture Radar Interferometry (InSAR) Adapted from and the ESA Interferometric SAR overview by Rocca et al. http://earth.esa.int/workshops/ers97/program-details/speeches/rocca-et-al/
More informationTerrain correction. Backward geocoding. Terrain correction and ortho-rectification. Why geometric terrain correction? Rüdiger Gens
Terrain correction and ortho-rectification Terrain correction Rüdiger Gens Why geometric terrain correction? Backward geocoding remove effects of side looking geometry of SAR images necessary step to allow
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 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 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 informationDIGITAL ELEVATION MODEL GENERATION FROM INTERFEROMETRIC SYNTHETIC APERTURE RADAR USING MULTI-SCALE METHOD
DIGITAL ELEVATION MODEL GENERATION FROM INTERFEROMETRIC SYNTHETIC APERTURE RADAR USING MULTI-SCALE METHOD Jung Hum Yu 1, Linlin Ge, Chris Rizos School of Surveying and Spatial Information Systems University
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 informationGEOG 4110/5100 Advanced Remote Sensing Lecture 4
GEOG 4110/5100 Advanced Remote Sensing Lecture 4 Geometric Distortion Relevant Reading: Richards, Sections 2.11-2.17 Geometric Distortion Geometric Distortion: Errors in image geometry, (location, dimensions,
More informationInSAR Operational and Processing Steps for DEM Generation
InSAR Operational and Processing Steps for DEM Generation By F. I. Okeke Department of Geoinformatics and Surveying, University of Nigeria, Enugu Campus Tel: 2-80-5627286 Email:francisokeke@yahoo.com Promoting
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 informationCOORDINATE TRANSFORMATION. Lecture 6
COORDINATE TRANSFORMATION Lecture 6 SGU 1053 SURVEY COMPUTATION 1 Introduction Geomatic professional are mostly confronted in their work with transformations from one two/three-dimensional coordinate system
More information: Find the values of the six trigonometric functions for θ. Special Right Triangles:
ALGEBRA 2 CHAPTER 13 NOTES Section 13-1 Right Triangle Trig Understand and use trigonometric relationships of acute angles in triangles. 12.F.TF.3 CC.9- Determine side lengths of right triangles by using
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 information7-5 Parametric Equations
3. Sketch the curve given by each pair of parametric equations over the given interval. Make a table of values for 6 t 6. t x y 6 19 28 5 16.5 17 4 14 8 3 11.5 1 2 9 4 1 6.5 7 0 4 8 1 1.5 7 2 1 4 3 3.5
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 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 informationUnit 7: Trigonometry Part 1
100 Unit 7: Trigonometry Part 1 Right Triangle Trigonometry Hypotenuse a) Sine sin( α ) = d) Cosecant csc( α ) = α Adjacent Opposite b) Cosine cos( α ) = e) Secant sec( α ) = c) Tangent f) Cotangent tan(
More informationDERIVATION of the BACKSCATTERING COEFFICIENT σ o in ESA ERS SAR PRI PRODUCTS
ERS SAR CALIBRATION DERIVATION of the BACKSCATTERING COEFFICIENT σ o in ESA ERS SAR PRI PRODUCTS H. Laur 1, P. Bally 2, P. Meadows 3, J. Sanchez 4, B. Schaettler 5, E. Lopinto 6, D. Esteban 4 Document
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 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 informationMath12 Pre-Calc Review - Trig
Math1 Pre-Calc Review - Trig Multiple Choice Identify the choice that best completes the statement or answers the question. 1. Which of the following angles, in degrees, is coterminal with, but not equal
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 informationSection 5.4: Modeling with Circular Functions
Section 5.4: Modeling with Circular Functions Circular Motion Example A ferris wheel with radius 25 feet is rotating at a rate of 3 revolutions per minute, When t = 0, a chair starts at its lowest point
More informationMid-Chapter Quiz: Lessons 9-1 through 9-3
Graph each point on a polar grid. 1. A( 2, 45 ) 3. Because = 45, locate the terminal side of a 45 angle with the polar axis as its initial side. Because r = 2, plot a point 2 units from the pole in the
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 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 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 informationTrigonometric Ratios and Functions
Algebra 2/Trig Unit 8 Notes Packet Name: Date: Period: # Trigonometric Ratios and Functions (1) Worksheet (Pythagorean Theorem and Special Right Triangles) (2) Worksheet (Special Right Triangles) (3) Page
More informationSection Graphs of the Sine and Cosine Functions
Section 5. - Graphs of the Sine and Cosine Functions In this section, we will graph the basic sine function and the basic cosine function and then graph other sine and cosine functions using transformations.
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 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 informationNelson Functions 11 Errata
Nelson Functions 11 Errata 1: Introduction to Functions Location Question Correct Answer Getting Started 6a Graph is correct but vertex and axis of symmetry are not labelled. Add blue point labelled (in
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 informationMath For Surveyors. James A. Coan Sr. PLS
Math For Surveyors James A. Coan Sr. PLS Topics Covered 1) The Right Triangle 2) Oblique Triangles 3) Azimuths, Angles, & Bearings 4) Coordinate geometry (COGO) 5) Law of Sines 6) Bearing, Bearing Intersections
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 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 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 informationGraphics in MATLAB. Responsible teacher: Anatoliy Malyarenko. November 10, Abstract. Basic Plotting Commands
Graphics in MATLAB Responsible teacher: Anatoliy Malyarenko November 10, 2003 Contents of the lecture: Two-dimensional graphics. Formatting graphs. Three-dimensional graphics. Specialised plots. Abstract
More informationFall 2016 Semester METR 3113 Atmospheric Dynamics I: Introduction to Atmospheric Kinematics and Dynamics
Fall 2016 Semester METR 3113 Atmospheric Dynamics I: Introduction to Atmospheric Kinematics and Dynamics Lecture 5 August 31 2016 Topics: Polar coordinate system Conversion of polar coordinates to 2-D
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 informationUpcoming altimeter measurements : nadir altimetry from Ku/C to Ka-band, SAR mode, interferometric SAR
Upcoming altimeter measurements : nadir altimetry from Ku/C to Ka-band, SAR mode, interferometric SAR Jean-Claude Souyris CNES, Service Altimétrie & Radar Acknowledgments : Nathalie Steunou, Roger Fjortoft,
More informationPractical work on SAR Interferometry Data Processing using DORIS Software Y.S. Rao CSRE, IIT Bombay
Practical work on SAR Interferometry Data Processing using DORIS Software Y.S. Rao CSRE, IIT Bombay ysrao@csre.iitb.ac.in Tools : DORIS (Delft Object-Oriented Radar Interferometric Software) is freely
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 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 informationPolSARpro v4.03 Forest Applications
PolSARpro v4.03 Forest Applications Laurent Ferro-Famil Lecture on polarimetric SAR Theory and applications to agriculture & vegetation Thursday 19 April, morning Pol-InSAR Tutorial Forest Application
More informationFAST SOLUTION OF PHASE UNWRAPPING PARTIAL DIFFERENTIAL EQUATION USING WAVELETS
Tenth MSU Conference on Differential Equations and Computational Simulations. Electronic Journal of Differential Equations, Conference 23 (2016), pp. 119 129. ISSN: 1072-6691. URL: http://ejde.math.txstate.edu
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 informationIMPROVING DEMS USING SAR INTERFEROMETRY. University of British Columbia. ABSTRACT
IMPROVING DEMS USING SAR INTERFEROMETRY Michael Seymour and Ian Cumming University of British Columbia 2356 Main Mall, Vancouver, B.C.,Canada V6T 1Z4 ph: +1-604-822-4988 fax: +1-604-822-5949 mseymour@mda.ca,
More informationANGLES 4/18/2017. Surveying Knowledge FE REVIEW COURSE SPRING /19/2017
FE REVIEW COURSE SPRING 2017 Surveying 4/19/2017 Surveying Knowledge 4 6 problems Angles, distances, & trigonometry Area computations Earthwork & volume computations Closure Coordinate systems State plane,
More informationStudy Guide and Review - Chapter 10
State whether each sentence is true or false. If false, replace the underlined word, phrase, expression, or number to make a true sentence. 1. A triangle with sides having measures of 3, 4, and 6 is a
More informationStudy Guide and Review - Chapter 10
State whether each sentence is true or false. If false, replace the underlined word, phrase, expression, or number to make a true sentence. 1. A triangle with sides having measures of 3, 4, and 6 is a
More informationSentinel-1 Toolbox. Offset Tracking Tutorial Issued August Jun Lu Luis Veci
Sentinel-1 Toolbox Offset Tracking Tutorial Issued August 2016 Jun Lu Luis Veci Copyright 2016 Array Systems Computing Inc. http://www.array.ca/ http://step.esa.int Offset Tracking Tutorial The goal of
More informationCh. 2 Trigonometry Notes
First Name: Last Name: Block: Ch. Trigonometry Notes.0 PRE-REQUISITES: SOLVING RIGHT TRIANGLES.1 ANGLES IN STANDARD POSITION 6 Ch..1 HW: p. 83 #1,, 4, 5, 7, 9, 10, 8. - TRIGONOMETRIC FUNCTIONS OF AN ANGLE
More informationPolar Coordinates. Chapter 10: Parametric Equations and Polar coordinates, Section 10.3: Polar coordinates 27 / 45
: Given any point P = (x, y) on the plane r stands for the distance from the origin (0, 0). θ stands for the angle from positive x-axis to OP. Polar coordinate: (r, θ) Chapter 10: Parametric Equations
More informationSeismic Reflection Method
Seismic Reflection Method 1/GPH221L9 I. Introduction and General considerations Seismic reflection is the most widely used geophysical technique. It can be used to derive important details about the geometry
More informationA METHOD OF GENERATING DEM FROM DSM BASED ON AIRBORNE INSAR DATA
A METHOD OF GENERATING DEM FROM DSM BASED ON AIRBORNE INSAR DATA Wenqing Lu 1,2,, Jixian Zhang 1,2, Guochao Xue 1,2, Chenggang Wang 1,2 1 Shan dong University of Science and Technology, Qingdao 266590,
More informationMemorandum. Clint Slatton Prof. Brian Evans Term project idea for Multidimensional Signal Processing (EE381k)
Memorandum From: To: Subject: Date : Clint Slatton Prof. Brian Evans Term project idea for Multidimensional Signal Processing (EE381k) 16-Sep-98 Project title: Minimizing segmentation discontinuities in
More informationSpatial Interpolation & Geostatistics
(Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Lag Mean Distance between pairs of points 1 Tobler s Law All places are related, but nearby places are related more than distant places Corollary:
More informationSEA SURFACE SPEED FROM TERRASAR-X ATI DATA
SEA SURFACE SPEED FROM TERRASAR-X ATI DATA Matteo Soccorsi (1) and Susanne Lehner (1) (1) German Aerospace Center, Remote Sensing Technology Institute, 82234 Weßling, Germany, Email: matteo.soccorsi@dlr.de
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 informationCorrection and Calibration 2. Preprocessing
Correction and Calibration Reading: Chapter 7, 8. 8.3 ECE/OPTI 53 Image Processing Lab for Remote Sensing Preprocessing Required for certain sensor characteristics and systematic defects Includes: noise
More informationGeometric Accuracy Evaluation, DEM Generation and Validation for SPOT-5 Level 1B Stereo Scene
Geometric Accuracy Evaluation, DEM Generation and Validation for SPOT-5 Level 1B Stereo Scene Buyuksalih, G.*, Oruc, M.*, Topan, H.*,.*, Jacobsen, K.** * Karaelmas University Zonguldak, Turkey **University
More informationDifferential Interferometry and Geocoding Software DIFF&GEO
Documentation User s Guide Differential Interferometry and Geocoding Software DIFF&GEO Geocoding and Image Registration Version 1.6 May 2011 GAMMA Remote Sensing AG, Worbstrasse 225, CH-3073 Gümligen,
More informationDefinition of Basic Polar Data Product
Definition of Basic Polar Data Product OPERA project 1d3: working document WD 02 03 Iwan Holleman and Gianmario Galli Date: September 22, 2003 1 1 Introduction During the OPERA project an additional deliverable
More informationINTERFEROMETRIC MULTI-CHROMATIC ANALYSIS OF HIGH RESOLUTION X-BAND DATA
INTERFEROMETRIC MULTI-CHROMATIC ANALYSIS OF HIGH RESOLUTION X-BAND DATA F. Bovenga (1), V. M. Giacovazzo (1), A. Refice (1), D.O. Nitti (2), N. Veneziani (1) (1) CNR-ISSIA, via Amendola 122 D, 70126 Bari,
More informationSpatial Interpolation - Geostatistics 4/3/2018
Spatial Interpolation - Geostatistics 4/3/201 (Z i Z j ) 2 / 2 Spatial Interpolation & Geostatistics Lag Distance between pairs of points Lag Mean Tobler s Law All places are related, but nearby places
More informationCEEN Engineering Measurements Final Exam Fall 2001 Closed Book, Calculator Required 3 Hour Time Limit
NAME Score CEEN 113-1 Engineering Measurements Final Exam Fall 001 Closed Book, Calculator Required 3 Hour Time Limit 1. (10 pts) You are interested in determining the height of a building. You are unable
More informationUnit 5 Day 5: Law of Sines and the Ambiguous Case
Unit 5 Day 5: Law of Sines and the Ambiguous Case Warm Up: Day 5 Draw a picture and solve. Label the picture with numbers and words including the angle of elevation/depression and height/length. 1. The
More information1. Be sure to complete the exploration before working on the rest of this worksheet.
PreCalculus Worksheet 4.1 1. Be sure to complete the exploration before working on the rest of this worksheet.. The following angles are given to you in radian measure. Without converting to degrees, draw
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 informationCoherence Based Polarimetric SAR Tomography
I J C T A, 9(3), 2016, pp. 133-141 International Science Press Coherence Based Polarimetric SAR Tomography P. Saranya*, and K. Vani** Abstract: Synthetic Aperture Radar (SAR) three dimensional image provides
More informationYoungstown State University Trigonometry Final Exam Review (Math 1511)
Youngstown State University Trigonometry Final Exam Review (Math 1511) 1. Convert each angle measure to decimal degree form. (Round your answers to thousandths place). a) 75 54 30" b) 145 18". Convert
More informationPlotting using Matlab. Vytautas Astromskas
Plotting using Matlab Vytautas Astromskas Contents 1 Manipulating a plot using GUI 2 2 Manipulating a plot using command line 4 3 Adding multiple lines to the same plot 5 4 Logarithmic axes 6 5 Plotting
More informationForest Retrievals. using SAR Polarimetry. (Practical Session D3P2a)
Forest Retrievals using SAR Polarimetry (Practical Session D3P2a) Laurent FERRO-FAMIL - Eric POTTIER University of Rennes 1 Pol-InSAR Practical Forest Application PolSARpro SIM PolSARproSim is a rapid,
More informationRange Imaging Through Triangulation. Range Imaging Through Triangulation. Range Imaging Through Triangulation. Range Imaging Through Triangulation
Obviously, this is a very slow process and not suitable for dynamic scenes. To speed things up, we can use a laser that projects a vertical line of light onto the scene. This laser rotates around its vertical
More informationLATE AND ABSENT HOMEWORK IS ACCEPTED UP TO THE TIME OF THE CHAPTER TEST ON
Trig/Math Anal Name No LATE AND ABSENT HOMEWORK IS ACCEPTED UP TO THE TIME OF THE CHAPTER TEST ON HW NO. SECTIONS ASSIGNMENT DUE TT 1 1 Practice Set D TT 1 6 TT 1 7 TT TT 1 8 & Application Problems 1 9
More informationTo graph the point (r, θ), simply go out r units along the initial ray, then rotate through the angle θ. The point (1, 5π 6. ) is graphed below:
Polar Coordinates Any point in the plane can be described by the Cartesian coordinates (x, y), where x and y are measured along the corresponding axes. However, this is not the only way to represent points
More informationTree Height Estimation Methodology With Xband and P-band InSAR Data. Lijun Lu Guoman Huang Qiwei Li CASM
Tree Height Estimation Methodology With Xband and P-band InSAR Data Lijun Lu Guoman Huang Qiwei Li CASM Outline CASMSAR dataset and test area Height Estimation with RVoG method Height Estimation with dual-band
More informationMATLAB Laboratory 09/23/10 Lecture. Chapters 5 and 9: Plotting
MATLAB Laboratory 09/23/10 Lecture Chapters 5 and 9: Plotting Lisa A. Oberbroeckling Loyola University Maryland loberbroeckling@loyola.edu L. Oberbroeckling (Loyola University) MATLAB 09/23/10 Lecture
More informationForward Time-of-Flight Detector Efficiency for CLAS12
Forward Time-of-Flight Detector Efficiency for CLAS12 D.S. Carman, Jefferson Laboratory ftof eff.tex May 29, 2014 Abstract This document details an absolute hit efficiency study of the FTOF panel-1a and
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 informationTriangle Trigonometry
Honors Finite/Brief: Trigonometry review notes packet Triangle Trigonometry Right Triangles All triangles (including non-right triangles) Law of Sines: a b c sin A sin B sin C Law of Cosines: a b c bccos
More informationPersistent Scatterer InSAR for Crustal Deformation Analysis, with Application to Volcán Alcedo, Galápagos
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. XXX, XXXX, DOI:1.19/6JB4763, 7 Persistent Scatterer InSAR for Crustal Deformation Analysis, with Application to Volcán Alcedo, Galápagos A. Hooper Department of Geophysics,
More informationEfficacy of Numerically Approximating Pi with an N-sided Polygon
Peter Vu Brewer MAT66 Honors Topic Efficacy of umerically Approximating Pi with an -sided Polygon The quest for precisely finding the irrational number pi has been an endeavor since early human history.
More informationLANDSLIDE PHENOMENA IN SEVAN NATIONAL PARK - ARMENIA
LANDSLIDE PHENOMENA IN SEVAN NATIONAL PARK - ARMENIA Andon Dimitrov Lazarov (1), Dimitar Minchev (2), Gurgen Aleksanyan (3), (1) Bourgas Free University, 62 San Stefano Str., 8000 Bourgas, Bulgaria,Email:
More information