24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2756  |  回复: 22

小孩一样在

新虫 (初入文坛)

引用回帖:
9楼: Originally posted by nickandmiles at 2015-02-01 12:18:41
找到这个论文的程序了。这个血管的程序挺多人用的。期望以后和你多交流。我也做过血管增强方面的研究。

这个程序下不了,能不能再发一份? 谢谢了  很需要
21楼2016-11-25 15:46:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小孩一样在

新虫 (初入文坛)

引用回帖:
18楼: Originally posted by 小明7034 at 2015-12-05 09:47:38
我也在做这个  想请问一下您找到相关的程序了吗(就是使用高斯模板对不同方向的微血管图像增强)  能不能给点建议 ...

一年过去了  你做的怎么样了?  能不能分享点经验给 我?
22楼2016-11-25 15:48:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

蓝青Ching

新虫 (初入文坛)

【答案】应助回帖

老师,请问能不能帮我解决一下python的问题?
考虑一个三层对称平板波导,包层材料为二氧化硅,芯层材料为硅。计算波导TE基膜的模场宽度(E场的半高全宽)随波导厚度变化的趋势曲线。
这是我写的,希望指点。。。。。


from camfr import *
from matplotlib.pyplot import *
from numpy import *

set_lambda(1.55)
set_N(1)
set_polarisation(TE)

sio2=Material(1.45)
silicon=Material(3.46)
air=Material(1.0)

w2=arange(0.01,6,0.8)
transmission=[]

for tw in w2:
    slab=Slab(sio2(4.995)+silicon(0.01)+sio2(4.995))
    slab=Slab(sio2(5-w2/2)+silicon(w2)+sio2(5-w2/2)

x=arrange(0,tw,0.8)
              
field1=[slab.mode(0).field((coord(tx,0,0).E1() for tx in x]
field2=[slab.mode(0).field((coord(tx,0,0)).E1()*(slab.n(Coord(tx,0,0)))**2 for tx in x]
index=[slab.n((coord(tx,0,0)) for tx in x]

M=max(field)
y=M/2
[m,n]=find(abs(M-y)<0.00001)
Q=m-n
      
plot(x,real(Q))
plot(x,absolute(index))
show()
23楼2016-11-27 07:39:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 emily1227 的主题更新
信息提示
请填处理意见