GACOS.PhCor: GACOS Atmospheric Phase delay correction

In SubhadipDatta/GInSARCorW: GACOS InSAR Correction Workflow

Description

GACOS Atmospheric Phase delay correction

Usage

GACOS.PhCor(
  unw_pha,
  re_dztd,
  wavelength = "in meter",
  inc_ang = 90,
  ref_lat = NA,
  ref_lon = NA
)

Arguments

unw_pha	Un-wrapped InSAR tile/raster.
re_dztd	Resampled Di-ZTD.
wavelength	SAR wavelength in meter.
inc_ang	SAR incident angle (to get output in LOS direction, don't use if not needed).
ref_lat	A reference point for correction, If NA, It use the tile center latitude.
ref_lon	A reference point for correction, If NA, It use the tile center longitude.

Author(s)

Subhadip Datta

Examples

library(raster)
library(GInSARCorW)
library(circular)
noDataAsNA<-FALSE
i1m<-system.file("td","20170317.ztd.rsc",package = "GInSARCorW")
i2m<-system.file("td","20170410.ztd.rsc",package = "GInSARCorW")
GACOS_ZTD_T1<-GACOS.Import(i1m,noDataAsNA)
GACOS_ZTD_T2<-GACOS.Import(i2m,noDataAsNA)
dztd<-d.ztd(GACOS_ZTD_T1,GACOS_ZTD_T2)
unw_pha<-raster(system.file("td","Unw_Phase_ifg_17Mar2017_10Apr2017_VV.img",package = "GInSARCorW"))
crs(unw_pha)<-CRS("+proj=longlat +datum=WGS84 +no_defs")
re_dztd<-d.ztd.resample(unw_pha,dztd)
GACOS.PhCor(unw_pha,re_dztd,0.055463,inc_ang=39.16362,ref_lat=NA,ref_lon=NA)

SubhadipDatta/GInSARCorW documentation built on July 3, 2020, 6:33 p.m.