数学在绘图中的应用可视化（判断点是否在多边形内，矩阵的二维旋转，矩阵的三维旋转）

最新推荐文章于 2021-03-17 04:46:54 发布
最新推荐文章于 2021-03-17 04:46:54 发布
阅读量291 收藏 1
点赞数 1
分类专栏： 数据结构
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xiaobaobao123321/article/details/108611047
版权

前言：代码全部实现可视化，能够带来更直观的体验，java语言。

1：判断点是否在多边形内

这个问题最初遇到是我在毕设的游戏设计地图时候，如何防止玩家走出边界。

思想可以参考   https://blog.csdn.net/u011722133/article/details/52813374

就是通过计算一条射线和多边形的边相交点数判断。

效果图如下，当然也可以实现任意的多边形。

package Tools;

import java.awt.*;
import java.util.concurrent.TimeUnit;

import javax.swing.*;

import org.apache.commons.lang.math.RandomUtils;

public class RandomPoint extends JPanel {

	public static JFrame jFrame;
	public static int JFRAME_WIDTH = 1300;
	public static int JFRAME_HEIGHT = 950;
	static boolean first = true;
	int[][] map = {{-1561, 162}, {-891, 619}, {-351, 277}, {-1013, -195}};

	public RandomPoint() {

		jFrame = new JFrame();
		Dimension ScreenSize = Toolkit.getDefaultToolkit().getScreenSize();
		int JFrame_X = (int) (ScreenSize.getWidth() - JFRAME_WIDTH) / 2;
		int JFrame_Y = (int) (ScreenSize.getHeight() - JFRAME_HEIGHT) / 2;
		jFrame.setLocation(JFrame_X, JFrame_Y);
		jFrame.setSize(JFRAME_WIDTH, JFRAME_HEIGHT);
		jFrame.setLayout(null);
		this.setBounds(0, 0, JFRAME_WIDTH, JFRAME_HEIGHT);
		this.setLayout(null);

		jFrame.add(this);
		jFrame.setDefaultCloseOperation(JFrame.EXIT_ON_CLOSE);
		jFrame.setVisible(true);

		new Thread(() -> {
			for (; ; ) {
				try {
					TimeUnit.MILLISECONDS.sleep(1);
				} catch (InterruptedException e) {
					e.printStackTrace();
				}

				this.repaint();
			}
		}).start();
	}

	public static void main(String[] args) {
		new RandomPoint();
	}

	@Override
	public void paint(Graphics g) {
		int M = 1600;
		int N = 200;

		if (first) {
			g.setColor(Color.GREEN);
			((Graphics2D) g).setStroke(new BasicStroke(5.0f));
			for (int i = 0, j = 3; i < 4; j = i++) {
				g.drawLine(map[i][0] + M, map[i][1] + N, map[j][0] + M, map[j][1] + N);
			}
			first = false;
		}
		g.setColor(Color.RED);

		int w = RandomUtils.nextInt(JFRAME_WIDTH);
		int h = RandomUtils.nextInt(JFRAME_HEIGHT);

		while (!map_collision(w - M, h - N)) {
			w = RandomUtils.nextInt(JFRAME_WIDTH);
			h = RandomUtils.nextInt(JFRAME_HEIGHT);
		}
		g.fillOval(w, h, 20, 20);

	}

	public boolean map_collision(float testx, float testy) {
		int[][] ver = map;
		int nvert = map.length;
		int i, j;
		boolean c = false;
		for (i = 0, j = nvert - 1; i < nvert; j = i++) {
			if (((ver[i][1] > testy) != (ver[j][1] > testy)) && (testx < (ver[j][0] - ver[i][0]) * (testy - ver[i][1]) / (ver[j][1] - ver[i][1]) + ver[i][0]))
				c = !c;
		}
		return c;
	}

}

2：二维矩阵的旋转

思想可以参考：https://baike.baidu.com/item/%E6%97%8B%E8%BD%AC%E7%9F%A9%E9%98%B5 里的二维空间矩阵计算。

高中知识体系

效果如下

package Tools;

import java.awt.*;
import java.util.concurrent.TimeUnit;

import javax.swing.*;

/**
 * 旋转的矩阵
 */
public class RotationMatrix extends JPanel {

	public static JFrame jFrame;

	public static int JFRAME_WIDTH = 818;
	public static int JFRAME_HEIGHT = 848;

	double degrees = 0.0;

	public RotationMatrix() {

		jFrame = new JFrame();
		Dimension ScreenSize = Toolkit.getDefaultToolkit().getScreenSize();
		int JFrame_X = (int) (ScreenSize.getWidth() - JFRAME_WIDTH) / 2;
		int JFrame_Y = (int) (ScreenSize.getHeight() - JFRAME_HEIGHT) / 2;
		jFrame.setLocation(JFrame_X, JFrame_Y);
		jFrame.setSize(JFRAME_WIDTH, JFRAME_HEIGHT);
		jFrame.setLayout(null);
		this.setBounds(0, 0, JFRAME_WIDTH, JFRAME_HEIGHT);
		this.setLayout(null);

		jFrame.add(this);
		jFrame.setDefaultCloseOperation(JFrame.EXIT_ON_CLOSE);
		jFrame.setVisible(true);

		new Thread(() -> {
			for (; ; ) {
				try {
					TimeUnit.MILLISECONDS.sleep(1);
				} catch (InterruptedException e) {
					e.printStackTrace();
				}

				degrees += 0.05;

				this.repaint();
			}
		}).start();
	}

	public static void main(String[] args) {
		new RotationMatrix();
	}

	@Override
	public void paint(Graphics g) {
		super.paint(g);
		g.setColor(Color.RED);
		int m = 10;
		for (int i = 0; i < m; i++) {
			for (int j = 0; j < m; j++) {
				double x = i * 30.0;
				double y = j * 30.0;

				double _x = x * Math.cos(Math.toRadians(degrees)) - y * Math.sin(Math.toRadians(degrees));
				double _y = x * Math.sin(Math.toRadians(degrees)) + y * Math.cos(Math.toRadians(degrees));

				g.fillOval(400 + (int) _x, 400 + (int) _y, 10, 10);
			}
		}
	}

}

3：比较难的三维矩阵在二维上的投影，可旋转

投影的计算步骤是：

1：三维坐标的旋转 参考：https://baike.baidu.com/item/%E6%97%8B%E8%BD%AC%E7%9F%A9%E9%98%B5 中的三维矩阵

2：计算三维坐标与平面的交点，求得投影平面的三维坐标 参考：https://blog.csdn.net/abcjennifer/article/details/6688080

3：投影平面的三维坐标转换成二维坐标，已知投影平面的三维中心点，两个坐标可以算得，在xy方向的距离

4：单得到距离还差方向，就是进行投影后的xy方向的比较。

说明

代码实现上尽可能减少重复的计算，设计了很多了缓存计算量的变量，包括枚举类的设计，在视图上没有处理好不可见点和线的隐藏。

效果如下：

看代码吧

package Tools;

import java.awt.*;
import java.util.concurrent.TimeUnit;

import javax.swing.*;

public class Matrix_3D extends JPanel {

	public static JFrame jFrame;

	static double length = 300.0;
	static double[][] points = {{-length, -length, length}, {length, -length, length}, {length, length, length}, {-length, length, length},
			{-length, -length, -length}, {length, -length, -length}, {length, length, -length}, {-length, length, -length}};
	//投影平面的中心点
	static double[] canvas = {0, 500, 500};
	//投影平面的法向量
	static double[] canvas_c = {0, 100, 100};
	//投影平面的上的一点，此点到canvas视为投影平面的最上方
	static double[] canvas_p = {0, 400, 600};
	//摄像机位置
	static double[] eye = {0, 600, 600};

	//计算结果缓存
	double[][] d = new double[8][3];
	int[][] b = new int[8][2];

	int center_x = 500;
	int center_y = 500;
	int JFRAME_WIDTH = 1000;
	int JFRAME_HEIGHT = 1000;

	enum XYZ {

		X(0, 0),
		Y(0, 0),
		Z(0, 0.1);

		private final double change;

		private double degree;
		double sin;
		double cos;

		XYZ(double init, double change) {
			this.degree = init;
			this.change = change;
		}

		public void init() {
			double radians = Math.toRadians(degree);
			sin = Math.sin(radians);
			cos = Math.cos(radians);
		}

		public void change() {
			degree = (degree + change) % 360;
			init();
		}
	}

	public Matrix_3D() {

		jFrame = new JFrame();
		Dimension ScreenSize = Toolkit.getDefaultToolkit().getScreenSize();
		int JFrame_X = (int) (ScreenSize.getWidth() - JFRAME_WIDTH) / 2;
		int JFrame_Y = (int) (ScreenSize.getHeight() - JFRAME_HEIGHT) / 2;
		jFrame.setLocation(JFrame_X, JFrame_Y);
		jFrame.setSize(JFRAME_WIDTH, JFRAME_HEIGHT);
		jFrame.setLayout(null);
		this.setBounds(0, 0, JFRAME_WIDTH, JFRAME_HEIGHT);
		this.setLayout(null);

		jFrame.add(this);
		jFrame.setDefaultCloseOperation(JFrame.EXIT_ON_CLOSE);
		jFrame.setVisible(true);

		XYZ.X.init();
		XYZ.Y.init();
		XYZ.Z.init();

		new Thread(() -> {
			for (; ; ) {
				try {
					TimeUnit.MILLISECONDS.sleep(1);
				} catch (InterruptedException e) {
					e.printStackTrace();
				}
				XYZ.X.change();
				XYZ.Y.change();
				XYZ.Z.change();
				this.repaint();
			}
		}).start();
	}

	public static void main(String[] args) {
		new Matrix_3D();
	}

	@Override
	public void paint(Graphics g) {
		super.paint(g);

		for (int i = 0; i < 8; i++) {
			get(i);
		}

		paint___(g, 0, 1);
		paint___(g, 1, 2);
		paint___(g, 2, 3);
		paint___(g, 3, 0);

		paint___(g, 0, 4);
		paint___(g, 1, 5);
		paint___(g, 2, 6);
		paint___(g, 3, 7);

		paint___(g, 4, 5);
		paint___(g, 5, 6);
		paint___(g, 6, 7);
		paint___(g, 7, 4);

	}

	public void paint___(Graphics g, int i, int j) {
		g.drawLine(center_x + b[i][0], center_y + b[i][1], center_x + b[j][0], center_y + b[j][1]);
	}

	public void get(int index) {
		d[index][0] = points[index][0];
		d[index][1] = points[index][1] * XYZ.X.cos - points[index][2] * XYZ.X.sin;
		d[index][2] = points[index][1] * XYZ.X.sin + points[index][2] * XYZ.X.cos;

		double qwe = d[index][0] * XYZ.Y.cos + d[index][2] * XYZ.Y.sin;
		d[index][2] = -d[index][0] * XYZ.Y.sin + d[index][2] * XYZ.Y.cos;

		double rty = qwe * XYZ.Z.cos - d[index][1] * XYZ.Z.sin;
		d[index][1] = qwe * XYZ.Z.sin + d[index][1] * XYZ.Z.cos;
		d[index][0] = rty;

		double[] t = aaaaCalPlaneLineIntersectPoint(d[index]);
		if (t == null) {
			return;
		}

		double o = gerDistance(t, canvas_p);
		double p = gerDistance(canvas_p, canvas);
		double i = gerDistance(t, canvas);
		double ii = Math.sqrt(i);

		double cosA = Math.abs((i + p - o) / (2 * ii * Math.sqrt(p)));
		double sina = Math.sqrt((1 - Math.pow(cosA, 2)));

		double xx = ii * sina * (t[0] >= canvas[0] ? 1 : -1);
		double yy = ii * cosA * (t[1] >= canvas[1] ? 1 : -1);

		b[index][0] = (int) xx;
		b[index][1] = (int) yy;
	}

	private double[] aaaaCalPlaneLineIntersectPoint(double[] linePoint) {
		double[] returnResult = new double[3];
		double vp1, vp2, vp3, n1, n2, n3, v1, v2, v3, m1, m2, m3, t, vpt;
		// 平面的法线向量
		vp1 = canvas_c[0];
		vp2 = canvas_c[1];
		vp3 = canvas_c[2];
		// 平面经过的一点坐标
		n1 = canvas[0];
		n2 = canvas[1];
		n3 = canvas[2];
		// 直线的方向向量
		v1 = eye[0] - linePoint[0];
		v2 = eye[1] - linePoint[1];
		v3 = eye[2] - linePoint[2];
		// 直线经过的一点坐标
		m1 = eye[0];
		m2 = eye[1];
		m3 = eye[2];
		vpt = v1 * vp1 + v2 * vp2 + v3 * vp3;
		if (vpt == 0) {
			returnResult = null;
		} else {
			t = ((n1 - m1) * vp1 + (n2 - m2) * vp2 + (n3 - m3) * vp3) / vpt;
			returnResult[0] = m1 + v1 * t;
			returnResult[1] = m2 + v2 * t;
			returnResult[2] = m3 + v3 * t;
		}
		return returnResult;
	}

	public double gerDistance(double[] a, double[] b) {
		double i = 0;
		for (int j = 0; j < 3; j++) {
			i += Math.pow(Math.abs(b[j] - a[j]), 2);
		}
		return i;
	}

}

 

小包包007
关注
专栏目录
canvas判断坐标点是否在给定的多边形内
w643462659的博客
06-02 357
canvas判断坐标点是否在给定的多边形内。
实对称矩阵的特征值求法_旋转之三 - 旋转矩阵
weixin_39557373的博客
11-03 1008
绕坐标轴的旋转绕着z轴旋转的3d矩阵写出来(推导过程可以见之前的文章)： 类似可以写出绕x，y轴旋转的矩阵： 绕向量u旋转假设物体绕单位向量 旋转 θ, 可以写出旋转矩阵： 具体推导过程可以参见这篇文章9.2部分.同时 R 也可以展开写成： 其中： 可以把这个矩阵看成外积或者转置形式。而 则是一个反对称矩阵（反对称矩阵为 ，即 ），可以找到它对应的向量 , 通常我们会这样表示，如...
旋转矩阵简单可视化与平移向量方向问题探讨
一点儿也不萌的萌萌
12-06 1582
接触了SLAM好久了，但是对于变换矩阵和平移矩阵始终觉得没有完全理解透彻，刚好碰到了因为理解导致的错误，就借此好好总结一下！ 直接开始! 1.只有平移时： O′−x′y′z′O'-x'y'z'O′−x′y′z′坐标系就是将O−xyzO-xyzO−xyz坐标系沿着yyy轴向右平移10个单位。 PwP_wPw​在O−xyzO-xyzO−xyz坐标系下的坐标值为[0,0,1]T[0,0,1]^T[0,...
判断点是否在多边形内(MFC)
06-03
用射线法判断点是否在多边形中，编译环境VC6.0，鼠标左键实现绘制多边形，右键进行判断。
图像旋转、伸缩的自写matlab实现
weixin_33840661的博客
03-12 1918
一、图像的旋转 　　今天的代码不是自己写的，缺少一些时间。但是认认真真推导了一下旋转的公式，代码的思想与原博博主一致，致敬！ 　　愚以为，自己来实现图像旋转算法的关键点有二：其一，确定旋转后的图像边界。其二，确定图像新旧坐标系下的图像变换矩阵。 　　下面重点介绍图像变换矩阵（算法）。由于我们常见的旋转算法是根据图像中心来进行的，并非是ROC坐标轴左上角的位置。如果我们直接进行变换，那么将难以...
C#二维三维图形绘制工程实例
05-28
在C#编程环境中，二维和三维图形的绘制是一项重要的技能，尤其对于游戏开发、科学可视化、数据呈现等领域至关重要。这个名为“C#二维三维图形绘制工程实例”的资源为我们提供了丰富的学习和实践素材，帮助开发者深入...
java二维数组货柜,装箱布局三维可视化控件开发<sup>*</sup>
weixin_29664819的博客
03-10 1053
杨淑德21　翟晓庆22　华顺刚21(1大连理工大学机械工程学院　2山东理工大学机械工程学院)摘　要:结合OpenGL和ActiveX技术，研究了三维装箱的可视化问题，在VisualStudio环境下，进行面向装箱布局的三维显示控件开发。通过接口方法获得货物配载信息，实现货物配载结果的可视化；为了观察集装箱内货物布局，根据用户观察方向与箱壁平面法线方向，判别集装箱箱壁平面是否绘制。以线框形...
matlab二维三维裂隙.zip
10-23
在MATLAB中，二维和三维随机裂隙程序主要用于地质建模、材料科学或者工程计算等领域。这些程序通常涉及几何形状的生成，特别是裂隙网络，它们可以帮助研究人员模拟裂隙对流体流动、物质传输以及结构稳定性的影响。...
winfromGraphics2维绘制三维图效果
04-24
在Windows Forms（WinForms）开发中，我们常常需要在应用程序中展示各种图形，包括二维和三维图。本文将深入探讨如何使用WinForms内置的Graphics类来实现看似三维的图形效果，尽管Graphics本质上只支持二维绘图。 ...
判断点是否在给定三角形内的matlab程序
08-09
- **地理信息系统（GIS）**：在GIS中，常常需要对地形进行三维建模，这时会用到Delaunay三角网来构建地形表面。 - **计算机图形学**：在游戏开发和电影特效制作中，需要创建逼真的地形模型，Delaunay插值技术可以...
旋转矩阵求旋转角度_（加餐）欧拉角及矩阵旋转
weixin_39633455的博客
12-18 3200
(1 本文约2000字，大约需要阅读20分钟)(2 本文适合有一定雷达及信号与系统知识基础的读者)乡亲们大家好，今天的这期内容就像有钱人的生活，比较枯燥，不过能完整读下来也确能有所收获，本文是分析微多普勒的数学基础，或者说研究微多普勒的必备数学工具，为后续聊微多普勒做铺垫，也是就是(如何做好一款4D 高分辨毫米波雷达)中的Micro Doppler部分，咱们开始。Preliminarie...
矩阵内寻找多边形面
Sur的工作日子
06-09 801
<br /><br />先查找点e和f是否在4边形a,b,c,d内，是！线段过这个4边形。<br />如果线段的两端分别在19,28,37,46内线段过这个4边形。<br />如果线段在26,68,48,24内，以26为例线段ef与dc的交点比过<br />g,而g如果小于d那么线段在4边形内。<br /> <br /><br />从横向和纵向扫描矩阵，扫描过程中碰到的第一点作为起始或结束点，在起始和结束点中间的<br />矩阵sell都设置成为多边形面所在的区域。进过n*n次扫描后多变形所在在矩阵内所覆盖
已知旋转矩阵，求关于各个轴的角度
gbz3300255的博客
10-09 1968
eigen库里就有现成的函数。 Eigen::Matrix3d R;//定义旋转矩阵的eigen形式 cv2eigen(R_opt,R);//R_opt是opencv的函数计算得到的旋转矩阵 是mat类型 做个转换 Eigen::Matrix3d Rotation = R.block<3, 3>(0, 0); Eigen::Vector3d eulerAngle = Rotat...
已知旋转矩阵求角度_旋转的知识点都在这儿了，搞懂了，轻松应对旋转的题目...
weixin_39534002的博客
11-05 1601
考点分析：(1)旋转与三角形、四边形相结合，求某个角的度数或者某条边的长度(多在选择题、填空题、解答题出现)(2)识别中心对称图形与中心对称，根据中心对称来画图(多在选择题、填空题、画图题中出现)(3)旋转作图，根据旋转求坐标(多在选择题、填空题、画图题出现)图形的旋转(1)把一个图形绕着某一个点O转动一个角度的图形变换叫做旋转，点O叫做旋转中心，转动的角度叫做旋转角。(2)特征： 图形旋转过程中...
根据旋转矩阵（roational matrix）计算三个坐标轴的旋转角度
热门推荐
mxs30443的专栏
05-29 1万+
通过当前的旋转矩阵计算出当前这个旋转矩阵对应的三个轴的旋转角度： 参考连接：http://stackoverflow.com/questions/15022630/how-to-calculate-the-angle-from-roational-matrix 相关资料下载：http://download.csdn.net/detail/mxs30443/9855408 或者http://
判断点在2D多边形内的一个算法
shichen35的专栏
05-30 576
test 123 321 end test
python人工智能之：多边形矩阵热图程序实战篇（二）
顾鸢的博客
11-01 1199
前言 这是一个六边形组件热图程序，大概花了我将近一个多月的时间才完成的，主要用到的知识有matplotlib、Numpy、Pandas,也就是python中的科学计算、数据分析以及可视化绘图库。 六边形组件热图 下面是这个程序运行效果图： 六边形热图实现 废话不说了，直接上代码。。。 # -*- coding: utf-8 -*- &quot;&quot;&quot; Created on Thu Jun 21 14:32:...
python 空间法向量可视化_利用空间法向量求二面角具体方法
weixin_33558072的博客
12-29 1299
如果已经求得各点坐标,或者说我们说的,能够建系,就用“法向量法”,所谓法向量,是指垂直于一个平面的直线,根据向量可在平面内任意平移,我们可以知道,一个平面的法向量有无数多条.以上是理论知识简介,因不知道你懂不,所以只得在此阐述下,不然可能会对下面的问题的理解不透产生障碍.具体做法：1. 设分别设出两个平面的法向量,n1=(x1, y1, z1); n2=(x2, y2, z2)2. 求出平面内线段...
matlab 多边形转矩阵,matlab – 如何在矩阵或二进制图像中绘制多边形以便使用图像处理功能？...
weixin_39615596的博客
03-17 217
我正在开发一个matlab程序,其中我使用多边形(凹面或凸面).我需要在多边形上使用imdilate或imerode等图像处理功能.为此,我应该将我的多边形转换为图像.我想知道是否有一种方法可以直接在二进制矩阵中绘制多边形(1为前景,0为背景)？目前,我使用’getframe’,然后’frame2im’,然后使用’im2bw’函数来执行此操作.但它的缺点是我无法控制最终图像(=矩阵)的大小(即,将...
MATLAB 三维多边形绘图详解
"这篇资源主要讲解了如何在MATLAB中进行三维多边形的绘制，以及二维绘图的基础知识，特别是plot命令的使用。" 在MATLAB中，三维多边形的绘制可以通过`fill3`函数实现，其使用方法与二维多边形的填充基本一致。例如...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值