GAMMA初学笔记一

GAMMA初学笔记一

20200721学习记录

记录一些出错的问题，都是比较基础的问题
1、生成初始编码查询列表
哨兵数据D-InSAR
NOTE: enter - as output filename to avoid creating the corresponding output file

gc_map 20190424.mli.par - GEO/mosaic.dem.par GEO/mosaic.dem DEM_seg_par DEM_seg mosaic.eqa.rdc 1 1 sim_sar

gc_map 20190424.mli.par -: command not found
名字写错了，复制一下文件的名字（GAMMA里面好多问题，是因为名称不一致而出错的，尤其是新手）。我就吃了很多次这样的亏，还是仔细点吧。

ERROR: cannot open input DEM parameter file: GEO/mosaic.dem.par GEO/mosaic.dem

这个部分总是提醒dem参数文件无法输入。看一下gc_map的定义。
指导手册中这样描述：

一旦 DEM 和 SAR 影像编码完成，就可以确定雷达和地图几何结构之间的转换关系，即，生成编码查找表。必须采用程序 gc_map 生成初始查找表。

usage: gc_map <MLI_par> <OFF_par> <DEM_par> <DEM> <DEM_seg_par> <DEM_seg> <lookup_table> [lat_ovr] [lon_ovr] [sim_sar] [u] [v] [inc] [psi] [pix] [ls_map] [frame] [ls_mode] [r_ovr]

input parameters:
  MLI_par         (input) ISP MLI or SLC image parameter file (slant range geometry)
  OFF_par         (input) ISP offset/interferogram parameter file (enter - if geocoding SLC or MLI data)
  DEM_par         (input) DEM/MAP parameter file
  DEM             (input) DEM data file (or constant height value)
  DEM_seg_par     (input/output) DEM/MAP segment parameter file used for output products

  NOTE: If the DEM_seg_par already exists, then the output DEM parameters will be read from this file
        otherwise they are estimated from the image data.

  DEM_seg         (output) DEM segment used for output products, interpolated if lat_ovr > 1.0  or lon_ovr > 1.0
  lookup_table    (output) geocoding lookup table (fcomplex)
  lat_ovr         latitude or northing output DEM oversampling factor (enter - for default: 1.0)
  lon_ovr         longitude or easting output DEM oversampling factor (enter - for default: 1.0)
  sim_sar         (output) simulated SAR backscatter image in DEM geometry
  u               (output) zenith angle of surface normal vector n (angle between z and n)
  v               (output) orientation angle of n (between x and projection of n in xy plane)
  inc             (output) local incidence angle (between surface normal and look vector)
  psi             (output) projection angle (between surface normal and image plane normal)
  pix             (output) pixel area normalization factor
  ls_map          (output) layover and shadow map (in map projection)
  frame           number of DEM pixels to add around area covered by SAR image (enter - for default = 8)
  ls_mode         output lookup table values in regions of layover, shadow, or DEM gaps (enter - for default)
                    0: set to (0.,0.)
                    1: linear interpolation across these regions
                    2: actual value (default)
                    3: nn-thinned
  r_ovr           range over-sampling factor for nn-thinned layover/shadow mode(enter - for default: 2.0)

20200729学习GAMMA记录

geocode
*** Forward geocoding transformation using a lookup table ***
*** Copyright 2018, Gamma Remote Sensing, v4.2 6-Nov-2018 uw/clw/of/cm ***

usage: geocode <lookup_table> <data_in> <width_in> <data_out> <width_out> [nlines_out] [interp_mode] [dtype] [lr_in] [lr_out] [n_ovr] [rad_max] [nintr]

input parameters:
  lookup_table  (input) lookup table containing pairs of real-valued output data coordinates
  data_in       (input) data file (format as specified by format_flag parameter)
  width_in      width of input data file and gc_map lookup table
  data_out      (output) output data file
  width_out     width of output data file
  nlines_out    number of lines for the output data file (enter - or 0 for default: all lines)
  interp_mode   resampling interpolation mode (enter - for default)
                  0: 1/dist (default)
                  1: nearest neighbor
                  2: SQR(1/dist)
                  3: constant
                  4: Gauss weighting
  dtype        input/output data type (enter - for default)
                  0: FLOAT (default)
                  1: FCOMPLEX
                  2: Sun/BMP/TIFF 8 or 24-bit raster image
                  3: UNSIGNED CHAR
                  4: SHORT
                  5: SCOMPLEX
                  6: DOUBLE
  lr_in         input  Sun/BMP/TIFF raster image flipped left/right (enter - for default: 1: not flipped (default), -1: flipped)
  lr_out        output Sun/BMP/TIFF raster image flipped left/right (enter - for default: 1: not flipped (default), -1: flipped)
  n_ovr         interpolation oversampling factor (enter - for default: 2)
  rad_max       maximum interpolation search radius (enter - for default 4*n_ovr: 8)
  nintr         number of points required for interpolation when not nearest neighbor (enter - for default: 4)

学习一个命令geocode，把虚拟SAR影像从地图结构转换到雷达结构。
geocode <初始查询列表> <模拟地图结构的SAR影像> <sar影像宽度> <输出文件> <输出文件宽度> [输出数据文件的行数(输入-或默认为0:所有行 )] [重新采样插值模式(enter -默认)] [输入/输出数据类型(enter -为默认值)] [输入Sun/BMP/TIFF光栅图像左/右翻转(enter - for default: 1: not翻转(default)， -1:翻转)] [输出Sun/BMP/TIFF光栅图像左/右翻转(enter - for default: 1: not翻转(default)， -1:翻转)] [插值过采样因子(enter -默认值:2)] [最大插值搜索半径(enter -默认为4*n_ovr: 8)] [当不是最近邻时需要的点数(enter -默认为4)]

width_in是指模拟地图结构的sar影像的宽度
width_out是指sar影像的宽度（之前生成的强度图像，打开.par文件可以直接看）

基础知识点：看经纬度范围：

SLC_corners sar.mli.par

计算偏移值和生成配准多项式

create_diff_par
create_diff_par <DIFF_par> [PAR_type] [iflg]

显示图像

dis2pwr 20190424.mli nsim_sar 5696 5696
dispwr 20190424.mli.eqa 7992