本文目录
0 前言
1 函数对应记录 --- 2019年10月20日
1.1 读取与输出遥感影像
1.2 取出多维数组中的某一维
1.3 取出数组中满足某一条件的索引
1.4 打印、显示字符
1.5 for循环
1.6 字符长度
1.7 列表增加一个元素
1.8 卷积
1.9 卷积

0 前言
 如果你也同时用过Matlab与Python，你应该知道：其中一个应该算是可以编程的软件工具包，另一个算是一门高级编程语言。
 在编程方面，他们是十分相似的，语法结构十分简答、易读。
 当我们使用大于两门编程语言的时候，有些函数名儿、用法等会出现混淆问题。所以，本文的目的也就是为了记录同一功能，Matlab、Python所对应的的不同函数，方便日后的查阅，也希望能给其他朋友提供帮助。
 系统及配置：
  系统: win7 x64

  python版本: 3.6.4 x64

  python包: numpy\GDAL

  matlab版本: R2016b
1 函数对应记录 — 2019年10月20日
1.1 读取与输出遥感影像

Matlab:

% input
[image,g]=geotiffread([path, 'input.tif']);  % read image
p=geotiffinfo([path, 'input.tif']);  % read projection
% output
geotiffwrite('output.tif',data,g, 'GeoKeyDirectoryTag', p.GeoTIFFTags.GeoKeyDirectoryTag);  % 输出数据到tif


Python:

 使用python读取及输出遥感图像

1.2 取出多维数组中的某一维

Matlab:

red = image(:,:,4);


Python:

red = image[:,:,4];

1.3 取出数组中满足某一条件的索引

Matlab:

idx_nbg=find(red>0);


Python:

import numpy as np
red = np.where(red>0);

1.4 打印、显示字符

Matlab:

fprintf('背景值%d个\n', 100);


Python:

print('背景值{}个'.format(100))

1.5 for循环

Matlab:

for i = 1:1000
	% ...
end


Python:

for i in range(1, 1000, 1):

1.6 字符长度

Matlab:

a = '123'
b = length(a)


Python:

a = '123'
b = len(a)

1.7 列表增加一个元素

Matlab:

a = 1;
list(1, i) = a;


Python:

a = 1
list.append(a)

1.8 卷积

Matlab:

fil = ones((width, width));  % 预设卷积核
rvi_conv = conv2(rvi_over,fil,'same');  % full valid


Python:

fil = np.ones((width, width));  % 预设卷积核
rvi_conv = cv2.filter2D(rvi_over, -1, fil);

1.9 卷积

Matlab:

fil = ones((width, width));  % 预设卷积核
rvi_conv = conv2(rvi_over,fil,'same');  % full valid


Python:

fil = np.ones((width, width));  % 预设卷积核
rvi_conv = cv2.filter2D(rvi_over, -1, fil);

matlab输出语句有两种，分别是：1、自由格式，语法如“disp(23+454-29*4)”；2、格式化输出，语法如“fprintf('The area is %8.5f\n',area)”。
MATLAB的输入与输出语句
1.MATLAB的输入语句input函数用于接收用户的输入:
a.输入数据>> x=input('please input a number:')
please input a number:22
x =
22
b.输入字符串>> x=input('please input a string:','s')
please input a string:this is a string
x =
this is a string
2.MATLAB输出语句包括自由格式(disp)和格式化输出(fprintf)两种>> disp(23+454-29*4)
361
>> disp([11 22 33;44 55 66;77 88 99])
11   22    33
44   55    66
77   88   99
>> disp('this is a string')
this is a string
>> area=12.56637889;
>> fprintf('The area is %8.5f\n',area)
The area is 12.56638
>>

print()函数概述
print()函数用于打印输出，是python中最常见的一个函数
语法用法如下：
objects --表示输出的对象。输出多个对象时，需要用 , (逗号)分隔。
sep -- 用来间隔多个对象。
end -- 用来设定以什么结尾。默认值是换行符 ，我们可以换成其他字符。
file -- 要写入的文件对象。
* 数字类型的输出方式
如：print(1)
输出结果为 ：1
数字类型的输出方式
* 字符串的输出方式
如：print("hello world")
输出结果为：hello world
字符串的输出方式
* 一次性输出多个对象的输出方式
如：a = 1
 b=“abc”
 print(a,b) # 输出对象用逗号分隔
一次性输出多个对象的输出方式
变量的输出方式(无论什么类型的数据，包括但不局限于：数值型，布尔型，列表变量，字典变量...都可以直接输出。)
如：num = 19
 print(num) # 输出数值型变量
 输出：19
 输出数值型变量
如：str = "lue dong ke ji"
 print(str) # 输出字符串变量
输出：lue dong ke ji 
输出字符串变量
如：list = [1,2,"a"]
 print(list) # 输出列表变量
输出：[1,2，a]
输出列表变量
如：tuple = (1,2,"a")
 print(tuple) # 输出元组变量
输出：(1,2,a)
输出元组变量
如：dict = {'a':1,'b':2}
 print(dict) # 输出字典变量
输出：｛'a':1,'b':2}
扩展：
输出字典变量
* 数据的格式化输出
如：a = "lue dong ke ji"
 b =len(a) # len 表示长度
 print("The length of %s is %d" %(a,b))
输出：The length of lue dong ke ji is 14
数据的格式化输出
扩展：
The length of %s is %d 这部分叫做：格式控制符
(a,b) 这部分叫做：转换说明符
% 字符，表示标记转换说明的开始
格式字符 说明 格式字符 说明
%s 字符串采用str()的显示 %x 十六进制整数
%r 字符串(repr())的显示 %e 指数(基底写e)
%c 单个字符 %E 指数(基底写E)
%b 二进制整数 %f,%F 浮点数
%d 十进制整数 %g 指数(e)或浮点数(根据显示长度)
%i 十进制整数 %G 指数(E)或浮点数(根据显示长度)
%o 八进制整数 %% 字符%
设定最小字段宽度和精度
最小字段宽度：转换后的字符串至少应该具有该值指定的宽度。如果是*(星号)，则宽度会从值元组中读出。
点(.)后跟精度值：如果需要输出实数，精度值表示出现在小数点后的位数。如果需要输出字符串，那么该数字就表示最大字段宽度。如果是*，那么精度将从元组中读出。
 RN = 5.3679015
 print('%8.5f'%RN) # 字段宽8，精度6 %8 字段宽8，5f 精度5倍
输出 : 5.36790
设定最小字段宽度和精度
注意：因为设定字符宽度为8，所以输出时有缩进，如：*f 表示没有指定宽度，所以输出时没有缩进
*f 表示没有指定宽度，所以输出时没有缩进
* 所处的位置不同，读取的内容也不同
* 换行输出和不换行输出
在python中，输出函数总是默认换行
如： for a in range(0,3): # 循环范围0 - 3 (范围不包含3)
 print(a)
输出：0
 1
 2
换行输出
如： for a in range(0,3) : # 循环范围0 - 3 (范围不包含3)
 print(a, end=' ')
输出：0 1 2
不换行输出
扩展：
换行符： 
如： for a in range(0,3) : # 循环范围0 - 3 (范围不包含3)
 print(a, end=' ')
 for a in range(5,9): # 循环范围5 - 9 (范围不包含3)
 print(a, end=' ')
输出：0 1 2 5 6 7 8
不换行输出
如： for a in range(0,3) : # 循环范围0 - 3 (范围不包含3)
 print(a, end=' ')
 print('')
 for a in range(5,9): # 循环范围5 - 9 (范围不包含3)
 print(a, end=' ')
输出：0 1 2 
 5 6 7 8
 换行输出
print 本身自带换行输出

Matlab函数总结1
1. X=sqrtm(A)

功能：计算矩阵A的平方根1/2，即X*X=A。若矩阵A的每个特征值都有非负实部，则X是唯一的，若矩阵A的特征值有负的实部，X返回负矩阵，若A为奇异矩阵，则X不存在。检测到A的奇异性时，将打印警告信息。
2. A=inv(B)

功能：inv是矩阵求逆的意思。具体用法其中B是输入的可逆矩阵，输出A就是B的逆矩阵，逆矩阵满足性质 AB=BA=E (E是单位阵)。如果输入的是不可逆矩阵会弹出警告，并返回inf。

举例：DINV=inv（sqrtm(D)）;

3. eig

(1)e=eig(A)

功能：得到矩阵A的特征值组成的一个向量

举例：​​​​​​

​​​​​​

(2)[V,D]=eig(A)

功能：得到矩阵A的特征值组成的对角矩阵D，并得到每个特征值对应的特征向量组成的矩阵V

举例：

4. sort

(1)B=sort(A)

功能：实现对向量A的从小到大的排列

举例：

(2)B=sort(A,dim)

功能：对矩阵A实现按列（1）或者按行（2）从小到大的排列

举例：

(3)B=sort(_,direction)

功能：对矩阵实现按一定方向的排序

举例：

（4)[B,I]=sort(_)

功能：实现对数据信息排序的同时用矩阵I返回起索引值

举例：



5. dir

（1）dir(’.’)

功能：列出当前目录下所有子文件夹和文件

（2）dir(‘E:\Matlab’)

功能：列出指定目录下所有子文件夹和文件

（3）dir(E:\Matlab,’*.bmp’)

功能： 列出指定目录下后缀为.bmp的文件

6. minus函数

（1）dir(A,B)

功能：实现数组的减法，但是很少使用，它为类启用运算符重载。

举例：

7. times

times(A,B)=A.B

功能：矩阵点乘的表示，但是很少使用 ，它为类启用运算符重载。

举例：


8. bsxfun
(1)C = bsxfun(fun,A,B)

功能：对启用单例扩展的两个数组应用逐元素二进制操作。在启用单例扩展的情况下，将handle fun函数指定的逐元素二进制操作应用于数组A和B。fun可以是以下内置功能之一：

举例：

（1）使用bsxfun从矩阵的相应列元素中减去列平均值

（2）通过指定函数的句柄，使用bsxfun调用自定义的二进制函数



9. reshape

(1)B = reshape(A,sz)

功能：把一个向量重塑成一个数组。



(2)B = reshape(A,sz1,sz2,…,szN)
功能：将矩阵重塑为具有指定列数。如： 将6乘6的幻方矩阵重塑为只有3列的矩阵。指定[]作为第一个维度大小，以便“整形”自动计算适当的行数。​​​​​​
10. cat

(1)C = cat(dim, A, B)

功能：将矩阵A和B按照dim维度连接到一起

(2)C = cat(dim, A1, A2, A3, A4, …)

功能：将矩阵A1，A2，A3…按照dim维度连接到一起