Files

205 lines
9.6 KiB
Plaintext

# *********************************** PyINT Template *******************************************#
# Template for PyINT v2.1
# Please fill in the parameters according to your project requirements
# ========================= Basic parameters ========================
# DEM file path (if not provided, SRTM will be downloaded automatically from OpenTopography)
#DEM = /path/to/your/dem.dem
# Master date for interferometric processing
masterDate = 20250711 # Master date in YYYYMMDD format
# Swath and burst parameters (for Sentinel-1 TOPS mode only)
start_swath = 1 # Starting swath number
end_swath = 3 # Ending swath number
start_burst = 1 # Starting burst position (relative to master image)
end_burst = 9 # Ending burst position (relative to master image)
# Multilooking parameters
range_looks = 10 # Multi-look numbers in range direction
azimuth_looks = 2 # Multi-look numbers in azimuth direction
# ========================= Download Data ========================
# Data source and search parameters
Data_Hub = ASF # Data download source: SSARA, Scihub, ASF
sensor = Sentinel-1A # Envisat, ERS-1/2, Sentinel-1A/B, etc.
track = 171 # Track number of the SAR image
frame = 96 # Frame number [defined according to ASF]
start_time = 2024-11-01 # Start date for data search
end_time = 2026-02-04 # End date for data search
# ASF download specific parameters
asf_platform = Sentinel-1A # Platform: Sentinel-1A, Sentinel-1B
asf_beam_mode = IW # Beam mode: IW (Interferometric Wide)
asf_polarization = VV # Polarization: VV, VH, HH, HV
asf_flight_direction = ascending # Flight direction: ascending, descending
asf_bbox = 120.1,30.7,120.2,30.9 # Bounding box: West,South,East,North
# ========================= DEM Generation ========================
# DEM download parameters (for OpenTopography)
opentopo_api_key = 09ad77d34545607fdf5cb182b64ac64e # OpenTopography API key
opentopo_dem_type = SRTMGL1 # DEM type: SRTMGL1(30m), SRTMGL3(90m), NASADEM(30m), COP30, COP90
opentopo_auto_tile = 1 # Auto-tile large regions (recommended: 1)
# DEM processing parameters
dem_lat_ovr = 1 # Latitude oversampling (as to SRTM-1: 0.5 for 60m, 2 for 15m)
dem_lon_ovr = 1 # Longitude oversampling
# ========================= Interferometry Parameters ========================
# Spectral filtering
Igram_Spsflg = 1 # Implement range spectral filtering
rbw_min = 1 # Minimum range bandwidth fraction (0.1 --> 1.0)
# Coregistration parameters
rwin4cor = 256 # Range window length for coregistration
azwin4cor = 256 # Azimuth window length for coregistration
rsample4cor = 32 # Range samples for fitting coregistration parameters
azsample4cor = 32 # Azimuth samples for fitting coregistration parameters
thresh4cor = 0.3 # Coherence threshold (GAMMA 2016+, SNR for older versions)
# Coherence estimation
Igram_Cor_rwin = 5 # Range window length for cc_wave coherence estimation
Igram_Cor_awin = 5 # Azimuth window length for cc_wave coherence estimation
Igram_Cor_Win = 5 # Window for cc estimation in adf
adf_alpha = 0.4 # Alpha for Gold-Stein filtering
# ========================= Simulation Phase ========================
Igram_Flag_TDM = N # Y for Tandem-X
Simphase_rpos = -
Simphase_azpos = -
Simphase_rwin = 256
Simphase_azwin = 256
Simphase_thresh = -
# ========================= Unwrapping ========================
mcf_triangular = 1 # Triangular type of mcf [0: regular; 1: delaunay]
unwrap_patr = 1 # Unwrap patches in range direction
unwrap_pataz = 1 # Unwrap patches in azimuth direction
unwrapThreshold = 0.05 # Coherence threshold for unwrapping
auto_unw = 1 # Auto or manual reference point [1: auto; 0: manual]
make_mask = 0 # Make mask based on coherence [0: don't mask; 1: make mask]
init_flag = 1 # Flag to set phase at reference point [0: use initial; 1: set to 0.0]
r_refer = - # Phase reference point range offset
a_refer = - # Phase reference point azimuth offset
# ========================= Geocoding ========================
geo_interp = 0 # Interpolation method [0: nearest; 1: bicubic spline]
# ========================= Satellite Configuration ========================
satelite = S1A # Satellite: S1A, S1B, CSK, TSX, ALOS2, ALOS, ENVISAT
# ========================= Pair Selection ========================
network_method = sbas # Network method: sbas, sequential, delaunay, stars
endDate = YYYYMMDD # Exclude dates after this date
startDate = YYYYMMDD # Exclude dates before this date
conNumb = 15 # Connect number for sequential network
max_tb = 365 # Maximum temporal baseline (days)
max_sb = 100 # Maximum spatial baseline (meters)
min_tb = 1 # Minimum temporal baseline (days)
min_sb = 0 # Minimum spatial baseline (meters)
# ========================= Phase Gradient (Optional) ========================
dataformat = cpxfloat32 # Data format of interferogram
scale = 10000 # Original phase/scale
step_windows = 1 # Step length for spatial phase gradient
filter_wins = [7, 7] # Windows of median filter
directions = ['grad_east', 'grad_north', 'grad_northeast', 'grad_southeast']
# Mapping parameters for phase gradient
in_low_high = [0.001, 1] # Input range for mapping
out_low_high = [0.001, 1] # Output range for mapping
gamma = 0.8 # Gamma for curve shape
choose_linear = 'n' # Linear or nonlinear mapping ['y': linear; 'n': nonlinear]
# ========================= Hyp3 Format Conversion (NEW) ========================
# Convert GAMMA outputs to Hyp3-compatible GeoTIFF format
hyp3format = 0 # Enable Hyp3 format conversion [0: skip; 1: process]
hyp3format_parallel = 4 # Number of parallel processors for Hyp3 conversion
hyp3_output_dir = - # Output directory for Hyp3 products [default: projectName/Hyp3Products/]
# ========================= Phase Bias Correction (NEW) ========================
# Apply phase bias correction based on loop closures algorithm
# Reference: 10.1016/j.remote.2022.100013
phasebias = 0 # Enable phase bias correction [0: skip; 1: process]
phasebias_interval = 12 # Temporal baseline in days [6, 12, or 24]
phasebias_nlook = 10 # Number of looks for multilooking
phasebias_num_a = 2 # Number of calibration parameters to estimate
phasebias_estimate_an = 0 # Estimate an from data [0: use defaults; 1: estimate from loop closures]
phasebias_max_con = 5 # Maximum number of connections
# Default an values (used if phasebias_estimate_an = 0)
a1_6_day = 0.50
a2_6_day = 0.36
a3_6_day = 0.299
a4_6_day = 0.2476
a1_12_day = 0.494
a2_12_day = 0.297
a3_12_day = 0.24
a4_12_day = 0.22
a1_24_day = 0.48
a2_24_day = 0.28
a3_24_day = 0.22
a4_24_day = 0.20
phasebias_start = - # Start date YYYYMMDD [default: auto-detect]
phasebias_end = - # End date YYYYMMDD [default: auto-detect]
# ========================= Atmospheric Correction ========================
# GACOS tropospheric correction
gacos_correction = 0 # Enable GACOS correction [0: skip; 1: process]
gacos_parallel = 4 # Number of parallel processors
# Ionospheric correction
ionosphere_correction = 0 # Enable ionospheric correction [0: skip; 1: process]
ionosphere_parallel = 4 # Number of parallel processors
# ========================= Processing Workflow Control ========================
download_data = 1 # Enable data download [0: skip; 1: process]
down_parallel = 4 # Number of parallel processors for download
raw2slc_all = 0 # Convert raw to SLC [0: skip; 1: process]
raw2slc_all_parallel = 4 # Number of parallel processors
extract_burst_all = 0 # Extract bursts [0: skip; 1: process]
extract_all_parallel = 4 # Number of parallel processors
coreg_all = 0 # Coregister SLCs [0: skip; 1: process]
coreg_all_parallel = 4 # Number of parallel processors
select_pairs = 0 # Select interferometric pairs [0: skip; 1: process]
diff_all = 0 # Generate differential interferograms [0: skip; 1: process]
diff_all_parallel = 4 # Number of parallel processors
unwrap_all = 0 # Unwrap interferograms [0: skip; 1: process]
unwrap_all_parallel = 4 # Number of parallel processors
geocode_all = 0 # Geocode products [0: skip; 1: process]
geocode_all_parallel = 4 # Number of parallel processors
hyp3format_all = 0 # Convert to Hyp3 format [0: skip; 1: process]
hyp3format_all_parallel = 4 # Number of parallel processors
phasebias_all = 0 # Apply phase bias correction [0: skip; 1: process]
load_data = 0 # Load data for MintPy processing [0: skip; 1: process]
# ========================= Good Luck =======================================