妙用FFT之逆变换生成脸谱

2023-11-08 10:50

本文主要是介绍妙用FFT之逆变换生成脸谱,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

原创:项道德(daode3056,daode1212)

        快速傅立叶变换(FFT)有正向变换与反向变换,在极坐标系中,当对下图(左)正向变换之后,取不完全的比例进行逆变换,可生成以下右边的图案: 

以下各图略去了原始曲线,直接出脸谱了:

 

//不完全FFT逆变换_花脸曲线(C#编程语言源代码):
        private void button3_Click(object sender, EventArgs e)
        {
            /*
                public static void FFT(double[] real, double[] imag, Accord.Math.FourierTransform.Direction direction)
                Accord.Math.Transforms.FourierTransform2 的成员
             */

            //画布,作图工具:
            Bitmap bmp = new Bitmap(pictureBox1.Width, pictureBox1.Height);
            int xOrg = bmp.Width / 2; int yOrg = bmp.Height / 2;
            Graphics g = Graphics.FromImage(bmp);
            g.Clear(Color.Black);
            //画刷,画笔:
            Brush bh = new SolidBrush(Color.FromArgb(220, 220, 220));
            Brush bh2 = new SolidBrush(Color.FromArgb(220, 220, 0));
            Brush bh3 = new SolidBrush(Color.FromArgb(250, 0, 250));
            Pen pen = new Pen(Color.FromArgb(255, 255, 0), 4);
            Pen pen2 = new Pen(Color.FromArgb(255, 0, 255), 1);

            //样本数据:
            double[] real = new double[501];// { -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5 };
            double[] imag = new double[501];// { 1, 1.5, 3, 4, 5, 4.5, 5, 4, 3, 1.5, 1 };
            Random RD = new Random();
            double u1 = 0.5 * RD.NextDouble() + 0.5;
            double u2 = 0.4 * RD.NextDouble() + 0.6;
            double u3 = 0.3 * RD.NextDouble() + 0.7;
            double u4 = 0.2 * RD.NextDouble() + 0.8;

            for (int i = 0; i <= 500; i++)
            {
                float m = (float)((i - 250) / 25f);
                real[i] = m;
                //imag[i] = 40 * (8 * Math.Cos(m / u4 / 6) - 2 * Math.Cos(m / u3) + Math.Cos(2*m / u2 )) - 100;//OK
                imag[i] = 40 * (8 * Math.Cos(m / u2 / 8) - 2 * Math.Cos(m / u3) + Math.Cos(2 * m / u4)) - 100; //OK               
                float q = (float)(PI * m / 10);
                g.FillEllipse(bh3, (int)(imag[i]*Math.Sin(q) + xOrg*0.5), (int)(imag[i] * Math.Cos(q) + yOrg), 4, 4);//加圆点
            }

            //FFT正向变换,结果存入real1[],imag1[]:
            double[] real1 = real;
            double[] imag1 = imag;
            Accord.Math.Transforms.FourierTransform2.FFT(real1, imag1, Accord.Math.FourierTransform.Direction.Forward);

            //FFT逆向变换,结果存入real2[],imag2[]:            
            int NUM = (int)(real1.Length * (u1));//不完全截取,进行逆变换.  最大是: real1.Length;//
            double[] real2 = new double[NUM];
            double[] imag2 = new double[NUM];
            for (int i = 0; i < NUM; i++)
            {
                real2[i] = real1[i];
                imag2[i] = imag1[i];
            }
            Accord.Math.Transforms.FourierTransform2.FFT(real2, imag2, Accord.Math.FourierTransform.Direction.Backward);

            //数据:
            List<Point> LP = new List<Point>();
            //绘制线段集与点集:
            for (int i = 0; i < NUM; i++)
            {
                double x = real2[i];
                double y = imag2[i];
                g.FillEllipse(bh2, (int)(x - 2 + xOrg*1.5), (int)(y - 2 + yOrg-100), 4, 4);//加圆点
                LP.Add(new Point((int)(x + xOrg*1.5), (int)(y + yOrg-100)));//+ yOrg
            }

            //显示数据:
            g.FillPolygon(bh3, LP.ToArray());//填充多边形
            g.DrawPolygon(pen, LP.ToArray());//封闭多边形
            //string txt = string.Format("N3算法参数:{0:X}{1:X}{2:X}{3:X}", (int)(u1 * 0xFF), (int)(u2 * 0xFF), (int)(u3 * 0xFF), (int)(u4 * 0xFF));
            //myGraphics.DrawString(txt, new Font("", 12), darkBlueBrush, boxWidth / 2 - 90, boxHeight - 25);
            //------------------
            //生成文本与保存图片: ===============================
            string txt = sender.ToString().Split(':')[1] + string.Format("_{0:X}`{1:X}`{2:X}`{3:X}", (int)(100 * u1), (int)(100 * u2), (int)(100 * u3), (int)(100 * u4));
            g.DrawString(txt, new Font("", 12), bh, boxWidth / 2 - 150, boxHeight - 25);
            pictureBox1.Image = bmp;
            bmp.Save(txt + ".png"); this.Text = txt + ".png  ---文件已经保存";
            this.pictureBox1.Refresh();
        }

作者:daode3056[毕业于杭州师范大学数学系] 

这篇关于妙用FFT之逆变换生成脸谱的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!


原文地址:https://blog.csdn.net/baidu_39430410/article/details/126449115
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.chinasem.cn/article/369344

相关文章

Python实现自动化Word文档样式复制与内容生成

《Python实现自动化Word文档样式复制与内容生成》在办公自动化领域,高效处理Word文档的样式和内容复制是一个常见需求,本文将展示如何利用Python的python-docx库实现... 目录一、为什么需要自动化 Word 文档处理二、核心功能实现:样式与表格的深度复制1. 表格复制(含样式与内容)2

python如何生成指定文件大小

《python如何生成指定文件大小》:本文主要介绍python如何生成指定文件大小的实现方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录python生成指定文件大小方法一(速度最快)方法二(中等速度)方法三(生成可读文本文件–较慢)方法四(使用内存映射高效生成

Maven项目中集成数据库文档生成工具的操作步骤

《Maven项目中集成数据库文档生成工具的操作步骤》在Maven项目中,可以通过集成数据库文档生成工具来自动生成数据库文档,本文为大家整理了使用screw-maven-plugin(推荐)的完... 目录1. 添加插件配置到 pom.XML2. 配置数据库信息3. 执行生成命令4. 高级配置选项5. 注意事

MybatisX快速生成增删改查的方法示例

《MybatisX快速生成增删改查的方法示例》MybatisX是基于IDEA的MyBatis/MyBatis-Plus开发插件,本文主要介绍了MybatisX快速生成增删改查的方法示例,文中通过示例代... 目录1 安装2 基本功能2.1 XML跳转2.2 代码生成2.2.1 生成.xml中的sql语句头2

使用Python自动化生成PPT并结合LLM生成内容的代码解析

《使用Python自动化生成PPT并结合LLM生成内容的代码解析》PowerPoint是常用的文档工具,但手动设计和排版耗时耗力,本文将展示如何通过Python自动化提取PPT样式并生成新PPT,同时... 目录核心代码解析1. 提取 PPT 样式到 jsON关键步骤:代码片段:2. 应用 JSON 样式到

SpringBoot实现二维码生成的详细步骤与完整代码

《SpringBoot实现二维码生成的详细步骤与完整代码》如今,二维码的应用场景非常广泛,从支付到信息分享,二维码都扮演着重要角色,SpringBoot是一个非常流行的Java基于Spring框架的微... 目录一、环境搭建二、创建 Spring Boot 项目三、引入二维码生成依赖四、编写二维码生成代码五

Android与iOS设备MAC地址生成原理及Java实现详解

《Android与iOS设备MAC地址生成原理及Java实现详解》在无线网络通信中,MAC(MediaAccessControl)地址是设备的唯一网络标识符,本文主要介绍了Android与iOS设备M... 目录引言1. MAC地址基础1.1 MAC地址的组成1.2 MAC地址的分类2. android与I

PyQt5+Python-docx实现一键生成测试报告

《PyQt5+Python-docx实现一键生成测试报告》作为一名测试工程师,你是否经历过手动填写测试报告的痛苦,本文将用Python的PyQt5和python-docx库,打造一款测试报告一键生成工... 目录引言工具功能亮点工具设计思路1. 界面设计:PyQt5实现数据输入2. 文档生成:python-

IDEA自动生成注释模板的配置教程

《IDEA自动生成注释模板的配置教程》本文介绍了如何在IntelliJIDEA中配置类和方法的注释模板,包括自动生成项目名称、包名、日期和时间等内容,以及如何定制参数和返回值的注释格式,需要的朋友可以... 目录项目场景配置方法类注释模板定义类开头的注释步骤类注释效果方法注释模板定义方法开头的注释步骤方法注

Python如何自动生成环境依赖包requirements

《Python如何自动生成环境依赖包requirements》:本文主要介绍Python如何自动生成环境依赖包requirements问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑... 目录生成当前 python 环境 安装的所有依赖包1、命令2、常见问题只生成当前 项目 的所有依赖包1、