- 博客(45)
- 资源 (3)
- 收藏
- 关注
RSS订阅原创 GIS 相关基础知识
本文介绍了GIS中的基础概念，包括地球模型、坐标系统和投影方法。地球被近似为椭球体，常用WGS-84坐标系描述，其参数包括半长轴和扁率。文中区分了天文纬度、地心纬度和大地纬度，以及海拔高与大地高的概念，并解释了它们之间的转换关系。此外，还介绍了投影坐标系，重点讲解了墨卡托投影、Web-Mercator投影和UTM投影的特点及应用场景。墨卡托投影保持角度不变但面积变形，Web-Mercator是网络地图标准，而UTM采用分度带投影适用于不同区域。这些概念为GIS空间数据处理提供了理论基础。
2025-09-30 20:25:14
777
原创 Tutorial 05 generate SAR RPC
本教程介绍了基于陆探一号（LT1）SAR数据生成RPC模型的方法，用于替代RD模型进行后向地理编码。RPC模型通过三次多项式表达像方与物方坐标关系，包含78个待求参数。求解过程包括：1）在三维空间均匀采样获取物方坐标点；2）通过RD模型计算对应像点坐标；3）确定归一化参数，将坐标转换至[-1,1]范围；4）将非线性模型转化为线性形式，采用SVD分解进行稳定求解。该方法可有效提高SAR影像后向地理编码的精度和效率。
2025-09-27 19:32:41
134
4
原创 Tutorial 04 Read SAR RPC
本文介绍了遥感影像中的有理函数模型（RPC），该模型通过比值多项式关联地面坐标与影像坐标，具有通用性强、计算效率高等特点。RPC模型采用三次有理多项式，适用于光学和SAR等多种遥感数据。文章详细说明了RPC的数学表达式、参数构成及文件格式（通常为.rpc或.rpb），并提供了代码实现功能，包括文件读写和坐标投影转换。RPC模型因其独立于传感器参数的特性，成为遥感数据处理中的重要工具。
2025-09-25 11:21:47
945
原创 Tutorial 03 SAR Geocode
本教程介绍了SAR影像的正向地理编码（geocode）过程，即从像方二维坐标和指定椭球高求解地物三维坐标。通过建立包含距离方程、多普勒方程和地球椭球方程的非线性方程组，采用牛顿迭代法进行求解。其中详细推导了各方程对目标位置分量的偏导数，构建雅可比矩阵，通过计算位置纠正量逐步逼近真实解。本文是Tutorial 02 SAR Backgeocode的逆过程，为LT1卫星SAR数据处理提供了完整的坐标转换方法。
2025-09-24 15:35:26
1061
原创 Tutorial 02 SAR Backgeocode
本文介绍了陆探一号(LT1)卫星SAR数据的后向地理编码(backgeocode)实现方法。主要内容包括：(1)距离-多普勒方程(RD方程)的数学表达及其几何意义，其中距离方程对应球面，多普勒方程对应圆锥面；(2)后向地理编码的处理流程，通过牛顿迭代法求解多普勒时间，并转换为SAR影像坐标；(3)提供MATLAB代码实现，包括参数读取、坐标转换及迭代求解过程。该技术可实现从地理坐标到SAR影像坐标的精确转换。
2025-09-23 15:42:46
965
原创 Tutorial 01 Obtain SAR Parameters
本教程将以陆探一号卫星（简称为"LT1"）的 SAR 数据产品为例，介绍如何使用 MATLAB 获取 SAR 定位模型所需的参数。
2025-09-23 00:15:35
1391
原创 MATLAB 函数签名器
MATLAB 函数签名器 ，在规范注释格式的基础上为函数文件或类文件自动生成函数签名，并保存至functionSignatures.json，函数签名能够允许使用者在编辑器调用这些（自定义）函数的时候具备代码提示和自动填充功能，提升编程体验。
2023-10-04 21:10:31
1723
原创 PDF目录书签生成器 PDG
基于 wxpython 和 PyPDF2 开发的一个简易的 GUI 程序 , 主要功能是给未添加书签的 PDF 添加目录书签。
2022-10-05 17:31:22
2173
1
原创 Eigen：eigen 与 opencv 矩阵的相互转换
eigen 与 opencv 矩阵的相互转换需要包含以下头文件：#include <Eigen/Dense>#include <Eigen/Core>// eigen 头文件必须放在 opencv 前面#include <opencv2/opencv.hpp>#include <opencv2/core/eigen.hpp>然后使用 cv::cv2eigen 和 cv::eigen2cv 进行转换。参考代码#include <iost
2022-02-21 12:26:50
4160
原创 Eigen：实现 sort
eigen 实现 sort使用标准库 std::sort()实现，可以用 lambda 表达式定义排序的规则。实现void sortVec(Eigen::ArrayXd & vec){ std::sort(vec.data(), vec.data()+vec.size());}void sortVec(const Eigen::ArrayXd & vec, Eigen::ArrayXd & dst, Eigen::ArrayXi & ind){
2022-02-21 12:26:10
1159
原创 Eigen： 实现 meshgrid
eigen 实现 meshgridvoid meshgrid(const Eigen::ArrayXd & x, const Eigen::ArrayXd & y, Eigen::ArrayXXd &xx, Eigen::ArrayXXd &yy){ Eigen::MatrixXd m_x = x.matrix(); Eigen::MatrixXd m_y = y.matrix(); int Nx = m_x.rows(); int Ny
2022-02-21 12:25:33
866
原创 Eigen：binaryExpr
binaryExpr字面意思就是“二元表达式”，主要的作用是实现”两个矩阵/数组的元素之间的运算“，如计算 atan2，实现按元素相加减的运算。MatrixBase::binaryExpr(const MatrixBase<OtherDerived> &,const CustomBinaryOp &) constCustomBinaryOp 可以是 lambda 表达式 或者是 ptr_fun 的仿函数。代码#include <iostream>#inc
2022-02-21 12:24:39
1105
原创 连续时间与离散时间傅里叶变换、傅里叶变换与傅里叶级数的关系
连续时间傅里叶变换(CTFT)与离散时间傅里叶变换(DTFT)，以及 傅里叶变换(FT)与傅里叶级数(FS)的关系
2022-01-13 12:52:13
6398
原创 使用 wx.tools.img2py（4）
使用 wx.tools.img2pyimg2py.py 将图像转换为PNG格式，并将其嵌入Python模块，以便在运行时将其加载到程序中。这样做的好处是，可以以*.pyc或使用freeze, py2exe等编译到程序中的形式来交付。具体用法：wx.tools.img2py — wxPython Phoenix 4.1.2a1 documentation下面的脚本实现将当前所在目录下的 png 文件夹的图片转换为 images.pyfrom wx.tools import img2pyimport
2022-01-10 19:52:41
630
原创 拉普拉斯变换的性质 - 对查表
【注意】初值定理要求：f(t)f(t)f(t) 连续可导；不包含任何阶次的冲激函数；F(s)F(s)F(s) 是真有理分式终值定理要求： x(t)x(t)x(t) 的终值存在，即 X(s)X(s)X(s) 的极点在左半 sss 平面点击查看 常见的拉普拉斯变换对 - 对查表...
2021-10-03 15:01:26
3696
2
原创 常见的拉普拉斯变换对 - 对查表
对于有理分式，求解拉氏逆变换最常用的方式是部分分式分解法。一个有理分式可以表示为H(s)=B(s)A(s)=∑n=0Nbnsn∑m=0MamsmH(s) = \frac{B(s)}{A(s)} = \frac{\displaystyle\sum_{n=0}^{N} b_n s^n}{\displaystyle\sum_{m=0}^{M} a_m s^m}H(s)=A(s)B(s)=m=0∑Mamsmn=0∑Nbnsn部分分式分解建立在极点分解的基础。极点即是分母 A(s)A(s)A(s).
2021-10-03 14:59:57
19796
原创 信号与系统05 拉普拉斯变换
1. 拉普拉斯变换文章目录1. 拉普拉斯变换1.1. 定义1.1.1. 计算公式1.1.2. 收敛域的计算1.1.3. 拉氏变换与傅氏变换的关系1.2. 性质1.3. 常见的拉氏变换对1.1. 定义1.1.1. 计算公式F(s)=∫−∞∞f(t)e−stdtf(t)=12πj∫σ−j∞σ+j∞F(s)estds\begin{aligned}F(s) &= \int_{-\infty}^{\infty} f(t) e^{-st} dt\\f(t) &= \frac{1}{2\p
2021-10-03 12:52:34
2789
原创 Matlab 手动实现plot——输入参数数可变的函数
Matlab 输入参数数可变的函数varargin 是函数定义语句中的一个输入变量，允许函数接受任意数量的输入参数。在执行函数时，varargin 是一个 1×N1×N1×N 元胞数组，其中 NNN 是函数在显式声明的输入后收到的输入项数。但是，如果该函数在显式声明的输入后未收到任何输入，则 varargin 是空元胞数组。例子 1重新包装 plot，使得绘制的线条的默认颜色为红色。function varargout = redplot(varargin) [varargout{1:na
2021-09-28 14:20:35
2545
原创 Matlab 结构体之间相同字段的赋值
Matlab 结构体之间相同字段的赋值假设有两个结构体 main_subject 和 all_subject 。其中 main_subject 中包含了某位同学的主修成绩和学号姓名信息，需要将这些信息地填到包含所有科目信息的另一个结构体 all_subject 对应的字段中（不存在的字段保持默认值），如何比较便利地实现呢?all_subject = 包含以下字段的 struct: Name: '' ID: 0 Chinese: 0 Math
2021-09-27 15:00:05
3398
原创 信号与系统04 离散时间傅里叶变换
1. 离散时间傅里叶变换1. 离散时间傅里叶变换1.1. 周期序列的离散傅里叶级数 DFS1.1.1. 计算公式1.1.2. 离散傅里叶级数的性质1.2. 离散时间傅里叶变换 DTFT1.2.1. 计算公式1.2.2. 性质1.2.2.1. 唯一性1.2.2.2. 奇偶不变性1.2.2.3. 周期性1.2.2.4. 线性1.2.2.5. 共轭对称性1.2.2.6. 时移、频移特性1.2.2.7. 尺度变换特性1.2.2.8. 差分、求和特性1.2.2.9. 频域
2021-09-26 17:13:18
598
原创 什么是吉布斯现象
1. 什么是吉布斯现象1. 什么是吉布斯现象1.1. 什么吉布斯现象？1.2. 吉布斯现象形成的原因？1.3. 如何减小吉布斯现象？1.1. 什么吉布斯现象？矛盾性：在时域描述一个不连续的信号要求信号的有无穷的频率成分，但实际情况中不可能采样到无穷的频率成分。实际中的信号采样系统只能采样一定的频率范围，对不连续信号（或有无穷频率成分的信号）采样将会存在频率截断。频率截断会引起时域信号在不连续处产生“振铃效应”，这个现象成为吉布斯现象。吉布斯现象：由于频率截断现象，具有无穷频率分
2021-09-17 22:08:58
21623
2
原创 部分分式展开
部分分式展开部分分式展开的步骤主要为：判断有理分式是否为假分式，若是则将其化为真分式。有理分式F(x)=N(x)D(x)=bmxm+bm−1xm−1+⋯b1x1+b0x0xn+an−1xn−1+⋯a1x1+a0x0\def\MY#1#2{ #1_{#2} x^{#2}}F(x) = \frac{N(x)}{D(x)}= \frac{ \MY{b}{m} + \MY{b}{m-1} + \cdots \MY{b}{1} + \MY{b}{0}}{x^n + \MY{a}{n-1} + \cd
2021-09-17 20:11:08
14173
2
原创 常见的傅里叶变换对
1. 常见的傅里叶变换对1. 常见的傅里叶变换对1.1. 矩形脉冲相关1.2. 阶跃信号相关1.3. 冲激信号相关1.4. 直流信号1.5. 指数信号1.6. 符号函数相关1.1. 矩形脉冲相关矩形脉冲信号Gτ(t)↔τSa(τ2w)G_\tau(t) \leftrightarrow \tau \mathrm{Sa} (\frac{\tau}{2} w)Gτ(t)↔τSa(2τw)采样信号Sa(wct)↔πwcG2wc(w)\mathrm{Sa}(w_c t)
2021-09-17 12:33:57
38659
3
原创 微分运算的时域扩展
微分运算的时域扩展【问题】设 g(t)g(t)g(t) 可导，令 f(t)=dg(t)dtf(t) = \frac{dg(t)}{dt}f(t)=dtdg(t)，求 f(2t)=dg(2t)dtf(2t)=\frac{dg(2t)}{dt}f(2t)=dtdg(2t) 是否成立?【答】否。将导数写成极限的形式，即f(t)=dg(t)dt=limΔt→0g(t+Δt)−g(t)Δt\begin{aligned}f(t) &= \frac{dg(t)}{dt}\\&= \li
2021-09-16 22:33:19
398
原创 信号与系统03 连续时间的傅里叶变换
1. 连续时间的傅里叶变换1. 连续时间的傅里叶变换1.1. 周期信号的傅里叶级数 CTFS1.1.1. 展开的条件1.1.2. 计算公式三角形式的傅里叶级数指数形式的傅里叶级数1.1.3. 周期信号的频谱分析波形对称性与谐波特性的关系频谱结构与波形参数的关系周期信号的平均功率1.2. 非周期信号的傅里叶变换 CTFT1.2.1. 计算公式1.2.2. 性质1.2.3. 常见的傅里叶变换对1.1. 周期信号的傅里叶级数 CTFS1.1.1.
2021-09-12 16:24:15
807
原创 信号与系统02 系统知识点
1. 系统知识点1. 系统知识点1.1. 系统的表示1.2. 系统的分类1.2.1. 连续/离散时间系统1.2.2. 线性/非线性系统1.2.3. 时变/时不变系统1.2.4. 因果/非因果系统1.2.5. 稳定/非稳定系统1.2.6. 记忆/无记忆系统1.3. 系统的互联1.1. 系统的表示箭头/方框表示 f(t)→y(t)f(t)\to y(t)f(t)→y(t)，f(t)→H→y(t)f(t)\to \boxed{H} \to y(t)f(t)→H→y(t
2021-09-10 23:56:22
721
原创 信号与系统01 信号知识点
1. 信号知识点1. 信号知识点1.1. 信号的分类1.1.1. 确定信号和随机信号1.1.2. 连续时间信号和离散时间信号1.1.3. 周期信号和非周期信号1.1.4. 对称信号和非对称信号1.1.5. 能量有限信号，功率有限信号，能量功率均无限信号1.1.6. (反)因果信号、非因果信号1.1.7. 左边信号、右边信号和双边信号1.2. 典型信号的特点1.2.1. 实指数信号1.2.2. 虚指数信号1.2.3. 一般的复指数信号1.2.4. 离散时间的单位脉冲信号
2021-09-10 23:54:26
1035
原创 连续时间单位冲激信号δ(t)的基本性质
连续时间单位冲激信号 δ(t)\delta(t)δ(t) 的基本性质筛选特性：x(t)δ(t−t0)=x(t0)δ(t−t0)x(t)\delta(t-t_0) = x(t_0)\delta(t-t_0)x(t)δ(t−t0)=x(t0)δ(t−t0)取样特性：∫−∞+∞x(t)δ(t−t0)dt=x(x0)\displaystyle\int_{-\infty}^{+\infty}x(t)\delta(t-t_0) dt = x(x_0)∫−∞+∞x(t)δ(t−t0)dt=x(x0
2021-09-10 23:51:14
7415
1
Tutorial 05 generate SAR RPC
2025-09-27
Tutorial 04 Read SAR RPC
2025-09-25
Tutorial 03 SAR Geocode
2025-09-24
Tutorial 02 SAR Backgeocode
2025-09-23
Tutorial 01 Obtain SAR Parameters
2025-09-23
MyGitRepository.zip
2021-08-17
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人