今天要处理一些极坐标的数据,把他们转化成为直角坐标,本来这事不复杂,在python里面就一个命令的事情就可以把储存在极坐标下(矩阵大小M*N)的数据转化为任意大小的矩阵来储存直角坐标下的数据,其原理就是简单的线性插值。

http://liao.cpython.org/scipytutorial11/

但是坑爹的matlab的插值命令girddata(x,y,u,X,Y)要求x,y,u必须是一维数组,X,Y必须是二维数组,导致不像python一样直接能转换。写了一个简单的循环算法来把极坐标数据转换为直角坐标,用了循环嵌套和计数器,所以代码效率极低,但是也不是很想去优化,能用就行了。放在这里给大家借鉴,欢迎提出代码修改意见。

clear all
data = csvread('REP_2Dvelocity.csv');
r1 = reshape(data(:,1),[201,51])';
ux1 = reshape(data(:,3),[201,51])';
r1 = reshape(data(:,1),[201,51])';
z1 = reshape(data(:,2),[201,51])';
ux1 = reshape(data(:,3),[201,51])';
uz1 = reshape(data(:,4),[201,51])';
r2 = -flip(r1,2);
z2 = flip(z1,2);
ux2 = -flip(ux1,2);
uz2 = flip(uz1,2);
r = cat(2,r2,r1(:,2:201));
z = cat(2,z2,z1(:,2:201));
ux = cat(2,ux2,ux1(:,2:201));
uz = cat(2,uz2,uz1(:,2:201));r = r(41,:);
ux = ux(41,:);
n = 0;
for theta = 0:pi/360:pifor i = 1:401n = n+1;R = r(i);x(n) = r(i)*cos(theta);y(n) = r(i)*sin(theta);U(n) = ux(i)*cos(theta);V(n) = ux(i)*sin(theta);end
end
[X,Y] = meshgrid(-150:150,-150:150);
uq = griddata(x,y,U,X,Y);
vq = griddata(x,y,V,X,Y);
quiver(X,Y,uq,vq);

Matlab二维坐标下极坐标向直角坐标的转换相关推荐

  1. matlab二维离散/阶梯/极坐标图形绘制

    本博文源于matlab基础,主要针对离散图形.阶梯图形.极坐标图形进行绘制. 离散型图形 使用函数: stem(y) % 画火彩棍图,该图用线条显示数据点与x轴的距离,并在数据点处绘制小圆圈 stem ...

  2. matlab二维图形中确定一点的坐标系,在CAD中已知一个平面图中两个点的坐标如何建立坐标系以便求出图形中任一点的坐标?...

    在CAD中已知一个平面图中两个点的坐标如何建立坐标系以便求出图形中任一点的坐标?以下文字资料是由(历史新知网www.lishixinzhi.com)小编为大家搜集整理后发布的内容,让我们赶快一起来看一 ...

  3. Matlab快速将影像的二维坐标转换为经纬度坐标

    一般使用 [lat, lon] = pix2latlon(R,row,col) 逐个将影像的二维坐标转换为经纬度坐标,但这样用到双重循环,运行较慢,代码如下: [A, RA] = readgeoras ...

  4. matlab二维三维图形绘制和坐标轴范围设置

    matlab二维绘图 一. 二维图形(Two dimensional plotting) 1. 基本绘图函数(Basic plotting function):Plot, semilogx,   se ...

  5. matlab四维图形,MATLAB二维、三维、四维绘图

    二维图形绘制直角坐标极坐标统计图面域图填色图羽毛图.射线图Voronoi图彩带图二维半图三维图形绘制线状图数组图网格图面状图统计图瀑布图散点图函数图四维图形绘制切片图颜色权重图Reference 二维 ...

  6. matlab二维绘图部分

    matlab二维绘图部分 X,Y是向量,分别表示点集的横坐标和纵坐标 PLOT(X,Y,S) 符号函数(显函数.隐函数和参数方程)画图 (1) ezplot ezplot('f(x)',[a,b]) ...

  7. MATLAB二维图形绘制

    MATLAB二维图形绘制 数据点标记 数据点 颜色 线型 曲线图 一.plot函数 (一)最简单的plot函数调用格式: plot(x) (1)plot函数的参数 X 为普通向量, (2)plot函数 ...

  8. 基于栅格地图的路径规划(一)基于Matlab二维、三维栅格地图的构建

    基于栅格地图的路径规划(一)基于Matlab二维.三维栅格地图的构建 前言 1.二维栅格地图的创建 1.1.二维栅格地图构建原理 1.2.二维栅格地图构建例程 2.三维栅格地图的创建 2.1.三维栅格 ...

  9. 机器视觉学习之--贝叶斯学习 MATLAB二维正态分布二维图

    1.贝叶斯介绍 我个人一直很喜欢算法一类的东西,在我看来算法是人类智慧的精华,其中蕴含着无与伦比的美感.而每次将学过的算法应用到实际中,并解决了实际问题后,那种快感更是我在其它地方体会不到的. 一直想 ...

最新文章

  1. 比特币前核心开发者Mike Hearn三年前的预测一一应验
  2. 在开课吧学python的经历-做客李晨nic淘宝直播 胡海泉胡彦斌带货开课吧Python
  3. Ajax学习(6)---服务器端脚本和程序中用 JSON 进行响应和回复
  4. python不定长参数怎么相加_python函数不定长参数使用方法解析
  5. 网内病毒问题快速排除集锦
  6. 笔记-项目成本管理-制定预算的步骤与原则
  7. oracle开启/关闭归档模式
  8. cocos2d JS 源生js实现each方法
  9. React - 修改children(下)
  10. 新系统申请软件著作权详细操作说明
  11. OHSAS18000认证(转载)
  12. 50k热敏电阻温度对照表c语言,50K热敏电阻温度表..doc
  13. QAM调制解调的仿真实现
  14. 计算机关闭账号用户控制好吗,Win10系统彻底关闭用户帐户控制的方法
  15. 【案例分享】让新时代教育发展与“数”俱进
  16. rdo远程桌面管理快捷键在哪里?
  17. 单片机的基本构成及工作原理
  18. java初级工作总结_【Java初级程序员工作总结_Java初级程序员个人年终总结】-看准网...
  19. java awt 简单计算器,JAVA Swing 开发简易计算器(上)
  20. “听声避险”?这个AI可精确辩声 帮助防盗和规避交通事故

热门文章

  1. 进入一种全新的“区块链”思维方式,对我们的生活有什么影响?
  2. c语言 程序开机自启
  3. S.M.A.R.T 参数详解及推荐指标
  4. html中table设置滚动条
  5. WIF调试分享(ap6256)
  6. 传播模型——简单的元胞自动机(3)
  7. 基于python的管理系统_基于ssm的管理系统_基于python管理系统
  8. 浅谈const int *,int const *与int *const
  9. Main.obj : error LNK2019: 无法解析的外部符号 _Direct3DCreate9@4,该符号在函数 long __cdecl InitD3D(struct HWND__ *)
  10. 自己动手做一个绿色便携版的谷歌浏览器Chrome,可以放入U盘随意带走的