SAR成像处理编程——FFTW

2024-01-14 04:48
文章标签 编程 处理 成像 sar fftw

本文主要是介绍SAR成像处理编程——FFTW,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

FFTW概述

        FFTW ( the Faster Fourier Transform in the West) 是一个快速计算离散傅里叶变换的标准C语言程序集,其由MIT的M.Frigo 和S. Johnson 开发。可计算一维或多维实和复数据以及任意规模的DFT。

        FFTW 还包含对共享和分布式存储系统的并行变换,它可自动适应你的机器, 缓存,存储器大小,寄存器个数。

        FFTW 通常比其它开源Fourier变换程序都要快,最新版本为fftw-3.3.8。

        FFTW 的编码生成器采用面向对象设计技术和面向对象语言Caml 编写;它能自动适应系统硬件,因而可移植性很强。FFTW2.1.5 支持共享存储多线程并行和分布式存储 MPI 并行。FFTW 的运算性能远远领先于已有的其它 FFT 软件。

        FFTW 为任意大小的模式生成一个计划(plan),通过对该计划施行各种运算完成各种模式的转换;内部结构及其复杂性对用户透明;速度快 (适合各种机器的内部编译器、代码生成器利用 AST 在运行时生成代码并自我优化,而且不占用编译时间,采用分层存储技术)。

        FFTW 受到越来越多的科学研究和工程计算工作者的普遍青睐,并为量子物理、光谱分析、音视频流信号处理、石油勘探、地震预报、天气预报、概率论、编码理论、医学断层诊断等领域提供切实可行的大规模 FFT 计算。

       在配有Intel处理器的机器上,安装MKL数学库后,可以直接使用fftw,且速度比较快。但是没有intel处理器的机器上,不能安装MKL数学库,若想使用fftw,必须在自己机器上编译和安装。

Linux系统编译与安装

安装包:fftw-3.3.8.tar.gz

在安装包所在的路径下,一次执行以下命令,即可完成fftw的编译与安装

tar -xzvf fftw-3.3.8.tar.gz

cd fftw

./configure --prefix=/home/FFTW LDFLAGS=-L/home/FFTW/lib FC=ifort F77=ifort CC=gcc CXX=g++ --enable-mpi -enable-shared --enable-float

make

make install

根据实际需要,可以修改./configure的参数

使用示例

FFTW安装完成后,就可以使用了,下面给出了SAR成像的距离向FFT的例子。

相关代码和并行加速方案可以留言或私信

/************************************************************************************
*	距离向FFT
************************************************************************************/
int RangeFFT(MKL_Complex* sig, long Na, long Nr)
{long i, j, k, l, m, kk, ll, mm, nn;long lBlock_Size = 64;                        unsigned char   fftw_Flags = FFTW_MEASURE;#pragma omp parallel shared(sig) private(i,j){fftwf_plan  sS3_Plan;fftwf_complex *sS3_fftwf;sS3_fftwf = (fftwf_complex*)fftwf_malloc( Nr * sizeof( fftwf_complex ) );sS3_Plan = fftwf_plan_dft_1d( Nr , sS3_fftwf , sS3_fftwf , FFTW_FORWARD , fftw_Flags );				#pragma omp for schedule(static,1)  for ( i = 0; i < Na; i++ ){			for ( j = 0; j < Nr; j++ )                           {sS3_fftwf[j][0] = sig[i*Nr+j].re;sS3_fftwf[j][1] = sig[i*Nr+j].im;}	fftwf_execute( sS3_Plan );for ( j = 0; j < Nr; j++ )                           {sig[i*Nr+j].re = sS3_fftwf[j][0];sig[i*Nr+j].im = sS3_fftwf[j][1];}				} fftwf_free(sS3_fftwf);		fftwf_destroy_plan(sS3_Plan);		}return 0;
}

这篇关于SAR成像处理编程——FFTW的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/603992

相关文章

如何使用celery进行异步处理和定时任务(django)

《如何使用celery进行异步处理和定时任务(django)》文章介绍了Celery的基本概念、安装方法、如何使用Celery进行异步任务处理以及如何设置定时任务,通过Celery,可以在Web应用中... 目录一、celery的作用二、安装celery三、使用celery 异步执行任务四、使用celery

SpringBoot操作spark处理hdfs文件的操作方法

《SpringBoot操作spark处理hdfs文件的操作方法》本文介绍了如何使用SpringBoot操作Spark处理HDFS文件,包括导入依赖、配置Spark信息、编写Controller和Ser... 目录SpringBoot操作spark处理hdfs文件1、导入依赖2、配置spark信息3、cont

C#反射编程之GetConstructor()方法解读

《C#反射编程之GetConstructor()方法解读》C#中Type类的GetConstructor()方法用于获取指定类型的构造函数,该方法有多个重载版本,可以根据不同的参数获取不同特性的构造函... 目录C# GetConstructor()方法有4个重载以GetConstructor(Type[]

MyBatis延迟加载的处理方案

《MyBatis延迟加载的处理方案》MyBatis支持延迟加载(LazyLoading),允许在需要数据时才从数据库加载,而不是在查询结果第一次返回时就立即加载所有数据,延迟加载的核心思想是,将关联对... 目录MyBATis如何处理延迟加载?延迟加载的原理1. 开启延迟加载2. 延迟加载的配置2.1 使用

Android WebView的加载超时处理方案

《AndroidWebView的加载超时处理方案》在Android开发中,WebView是一个常用的组件,用于在应用中嵌入网页,然而,当网络状况不佳或页面加载过慢时,用户可能会遇到加载超时的问题,本... 目录引言一、WebView加载超时的原因二、加载超时处理方案1. 使用Handler和Timer进行超

Python中处理NaN值的技巧分享

《Python中处理NaN值的技巧分享》在数据科学和数据分析领域,NaN(NotaNumber)是一个常见的概念,它表示一个缺失或未定义的数值,在Python中,尤其是在使用pandas库处理数据时,... 目录NaN 值的来源和影响使用 pandas 的 isna()和 isnull()函数直接比较 Na

详解Python中通用工具类与异常处理

《详解Python中通用工具类与异常处理》在Python开发中,编写可重用的工具类和通用的异常处理机制是提高代码质量和开发效率的关键,本文将介绍如何将特定的异常类改写为更通用的ValidationEx... 目录1. 通用异常类:ValidationException2. 通用工具类:Utils3. 示例文

无人叉车3d激光slam多房间建图定位异常处理方案-墙体画线地图切分方案

墙体画线地图切分方案 针对问题:墙体两侧特征混淆误匹配,导致建图和定位偏差,表现为过门跳变、外月台走歪等 ·解决思路:预期的根治方案IGICP需要较长时间完成上线,先使用切分地图的工程化方案,即墙体两侧切分为不同地图,在某一侧只使用该侧地图进行定位 方案思路 切分原理:切分地图基于关键帧位置,而非点云。 理论基础:光照是直线的,一帧点云必定只能照射到墙的一侧,无法同时照到两侧实践考虑:关

Linux 网络编程 --- 应用层

一、自定义协议和序列化反序列化 代码: 序列化反序列化实现网络版本计算器 二、HTTP协议 1、谈两个简单的预备知识 https://www.baidu.com/ --- 域名 --- 域名解析 --- IP地址 http的端口号为80端口,https的端口号为443 url为统一资源定位符。CSDNhttps://mp.csdn.net/mp_blog/creation/editor

【Python编程】Linux创建虚拟环境并配置与notebook相连接

1.创建 使用 venv 创建虚拟环境。例如,在当前目录下创建一个名为 myenv 的虚拟环境: python3 -m venv myenv 2.激活 激活虚拟环境使其成为当前终端会话的活动环境。运行: source myenv/bin/activate 3.与notebook连接 在虚拟环境中,使用 pip 安装 Jupyter 和 ipykernel: pip instal