IDL之集合运算(自定义)

15195775117

学习IDL至今约6年了，我还不知道它有没有集合运算，反正我没找到。

由于编程中经常遇到集合运算的问题，或者有的问题可以通过集合运算简化，所以我索性自己写了一个，
见附件set_cal(即集合运算set_Calculation之意)
set_cal函数参数有2个:x,y，同时支持标量和2维向量，返回值是一个结构体，有4个结果，提供了x和y的交集、并集、2个差集(x-y,y-x)。


写PyQt、JS...这种软件开发的脚本，我觉得自己超聪明的，
但每次写这种数学性的算法，我就觉得自己是个弱智；
这个函数写了半天之久。


后来发现运行的时候很慢，我突然意识到自己的愚蠢：
明明只需要计算一个差集，我却连交集、并集和另一个差集一起算，不慢才怪！
于是把差集分出来，见附件set_dif函数。


以前我还试过从头开始写Krigging插值法和EMD分解，
前者我很快放弃了，因为IDL有现成的krig2D函数，后来发现python也有这个包，
而EMD分解我陆续写了很多天还没调整好，主要卡在端点问题上，后来发现python有现成的，网上说matlab也有这个算法，那我还费那个劲干嘛，吃力不讨好。

15195775117

python的集合运算

集合是python的基本数据结构，这比IDL简单多了
显然，python集合支持的数据类型是很丰富的，标量、多维矢量、字符等等都行！


x=[(1,2),(3,4),(2,3),(3,4)]
y=[(1,2),(2,3),(2,3),(2,4)]
print('x=',x)
print('y=',y)
print('x的集合=',set(x))
print('y的集合=',set(y))
print('交集=',set(x)&set(y))
print('并集=',set(x)|set(y))
print('并集=',set(x)-set(y))


x={'a','b','c'}
y={'b','c','d'}
print('交集=',set(x)&set(y))
print('并集=',set(x)|set(y))
print('并集=',set(x)-set(y))


x={1,2}
y={3,4}
if not x&y:
    print('空集')



输出：
x= [(1, 2), (3, 4), (2, 3), (3, 4)]
y= [(1, 2), (2, 3), (2, 3), (2, 4)]
x的集合= {(1, 2), (3, 4), (2, 3)}
y的集合= {(1, 2), (2, 3), (2, 4)}
交集= {(1, 2), (2, 3)}
并集= {(1, 2), (2, 3), (3, 4), (2, 4)}
并集= {(3, 4)}
交集= {'c', 'b'}
并集= {'c', 'a', 'd', 'b'}
并集= {'a'}
空集

15195775117

数组转集合

IDL剔除重复元素的方法：
function element,x
;x=[[1,2,1,3,3],[1,2,1,3,3]]
p=sort(x)
x2=x[p]
return,x2[uniq(x2)];一定要排序，uniq函数才能把不同的元素识别出来
end



python的方法：
x=[1,2,2,3,3,4]
print(set(x))

结果： {1, 2, 3, 4}