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

相关文章

Python FastAPI+Celery+RabbitMQ实现分布式图片水印处理系统

《PythonFastAPI+Celery+RabbitMQ实现分布式图片水印处理系统》这篇文章主要为大家详细介绍了PythonFastAPI如何结合Celery以及RabbitMQ实现简单的分布式... 实现思路FastAPI 服务器Celery 任务队列RabbitMQ 作为消息代理定时任务处理完整

揭秘Python Socket网络编程的7种硬核用法

《揭秘PythonSocket网络编程的7种硬核用法》Socket不仅能做聊天室,还能干一大堆硬核操作,这篇文章就带大家看看Python网络编程的7种超实用玩法,感兴趣的小伙伴可以跟随小编一起... 目录1.端口扫描器:探测开放端口2.简易 HTTP 服务器:10 秒搭个网页3.局域网游戏:多人联机对战4.

C#使用SQLite进行大数据量高效处理的代码示例

《C#使用SQLite进行大数据量高效处理的代码示例》在软件开发中,高效处理大数据量是一个常见且具有挑战性的任务,SQLite因其零配置、嵌入式、跨平台的特性,成为许多开发者的首选数据库,本文将深入探... 目录前言准备工作数据实体核心技术批量插入:从乌龟到猎豹的蜕变分页查询:加载百万数据异步处理:拒绝界面

Java并发编程必备之Synchronized关键字深入解析

《Java并发编程必备之Synchronized关键字深入解析》本文我们深入探索了Java中的Synchronized关键字,包括其互斥性和可重入性的特性,文章详细介绍了Synchronized的三种... 目录一、前言二、Synchronized关键字2.1 Synchronized的特性1. 互斥2.

Springboot处理跨域的实现方式(附Demo)

《Springboot处理跨域的实现方式(附Demo)》:本文主要介绍Springboot处理跨域的实现方式(附Demo),具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不... 目录Springboot处理跨域的方式1. 基本知识2. @CrossOrigin3. 全局跨域设置4.

python+opencv处理颜色之将目标颜色转换实例代码

《python+opencv处理颜色之将目标颜色转换实例代码》OpenCV是一个的跨平台计算机视觉库,可以运行在Linux、Windows和MacOS操作系统上,:本文主要介绍python+ope... 目录下面是代码+ 效果 + 解释转HSV: 关于颜色总是要转HSV的掩膜再标注总结 目标:将红色的部分滤

Python实现自动化接收与处理手机验证码

《Python实现自动化接收与处理手机验证码》在移动互联网时代,短信验证码已成为身份验证、账号注册等环节的重要安全手段,本文将介绍如何利用Python实现验证码的自动接收,识别与转发,需要的可以参考下... 目录引言一、准备工作1.1 硬件与软件需求1.2 环境配置二、核心功能实现2.1 短信监听与获取2.

Python使用date模块进行日期处理的终极指南

《Python使用date模块进行日期处理的终极指南》在处理与时间相关的数据时,Python的date模块是开发者最趁手的工具之一,本文将用通俗的语言,结合真实案例,带您掌握date模块的六大核心功能... 目录引言一、date模块的核心功能1.1 日期表示1.2 日期计算1.3 日期比较二、六大常用方法详

利用Go语言开发文件操作工具轻松处理所有文件

《利用Go语言开发文件操作工具轻松处理所有文件》在后端开发中,文件操作是一个非常常见但又容易出错的场景,本文小编要向大家介绍一个强大的Go语言文件操作工具库,它能帮你轻松处理各种文件操作场景... 目录为什么需要这个工具?核心功能详解1. 文件/目录存javascript在性检查2. 批量创建目录3. 文件

Python异步编程中asyncio.gather的并发控制详解

《Python异步编程中asyncio.gather的并发控制详解》在Python异步编程生态中,asyncio.gather是并发任务调度的核心工具,本文将通过实际场景和代码示例,展示如何结合信号量... 目录一、asyncio.gather的原始行为解析二、信号量控制法:给并发装上"节流阀"三、进阶控制