1. 程式人生 > >matlab實現生物細胞數目統計方案原理

matlab實現生物細胞數目統計方案原理

學生物和醫學最煩人的可能統計細胞數目是排名很前的了,因此有款軟體叫imageJ的在這方面就給予了很大的幫助,本著學習的態度,本人花了兩天的時間將細胞統計的功能做了復現,接下來就是方案解釋與實現的過程,不想看我廢話的可以直接跳到文末有介面圖和整體實現流程。

本人使用的平臺是MATLAB2017b,下載地址我將放在文末,而matlab這個軟體不需要我多說,它是一款很強大的軟體,因此我也將用它來實現模仿imageJ的功能軟體設計,這一章節主要是對遇到的問題做下分析與探討解決方案,具體實現過程在下一篇。

我們看圖說話,第一張據說是血細胞的圖片,我也不是學生物的,不太懂,第二很張據說是什麼細胞小核啥的。大家不難發現它們的共同點都是細胞巨多,一個個數要數到天黑,此時如果使用軟體將會大大減少工作量,動動滑鼠1

分鐘就搞定。


血細胞



細胞小核

    我來介紹下如何實現這個功能,本文的方法比較粗暴,沒有用到什麼分水嶺的方法,而是簡單的使用了二值化然後通過計算連通域的數量最終得出細胞數目。但是有個很棘手的問題,細胞連在一起怎麼辦?那還數什麼,不用急,細胞連在一起但是細胞核不連在一起呀,所以本文是通過二值化,只提取細胞核。

    拿下圖細胞分裂的圖片來舉例說明,可以看出染色較深的地方是細胞核和染色體,但是有的細胞染色較深有的較淺,很難分辨出來。

原圖

     接下來將三通道的彩色圖轉換為單通道的灰度圖,如下圖,此時問題就更加明顯了,可以看出有些地方的細胞核其實和其他深色地帶的細胞質的顏色是一樣的了,這樣會導致使用全域性二值化時沒辦法分辨出來,如下圖。

    而此時本人使用一種方法,就是通過USM銳化,還真別說,用了之後感覺這是神技啊,對於這些的銳化提高對比度的效果也太好了吧,如下圖就是通過銳化後的灰度圖,他下面那張是全域性二值化的圖,我敢用人格擔保,真的不是為了誇大這個方法而故意把上面那個二值化圖弄得好像一坨米田共那樣的。

    USM銳化其實理解起來是真的非常簡單,無非就是通過高斯模糊,然後通過一條小學生都可以理解的公式將原灰度圖和模糊後的圖片融合在一起,就可以實現銳化的效果,雖然公式看著簡單,但是我真的對提出這個公式的那位大神五體投地,它真的好簡潔好美。

(灰度影象-w*高斯模糊影象)/(1-w) 其中w是銳化的程度。通過這條公式,圖片看起來會更清晰,而且細胞核和細胞質的分界也越來越明顯。其中高斯模糊的方法可以參看本人的另一篇文章

最終通過MATLAB自帶的連通區域檢測函式檢測二值化後的影象中的連通區域就可以得到細胞數,如下圖

具體實現流程:

1、載入圖片

2、將rgb彩色圖轉換成灰度圖

3、通過高斯模糊將原影象模糊

4、通過USM銳化將原影象對比度提高

5、通過全域性二值化的方法獲取細胞核

6、通過連通區域檢測獲取細胞核的數目

matlab2017b下載地址:https://pan.baidu.com/s/1m9qo4Sbk9q_lyB-9vfx61w

有任何問題或者有技術探討可以發郵件給我,大家一起學習,郵箱:[email protected]

相關推薦

matlab實現生物細胞數目統計方案原理

學生物和醫學最煩人的可能統計細胞數目是排名很前的了,因此有款軟體叫imageJ的在這方面就給予了很大的幫助,本著學習的態度,本人花了兩天的時間將細胞統計的功能做了復現,接下來就是方案解釋與實現的過程,不想看我廢話的可以直接跳到文末有介面圖和整體實現流程。本人使用的平臺是MAT

生物細胞計數軟體MATLAB程式碼實現

    之前因為工作原因一直沒將程式碼貼上來,這次我比較詳細介紹下MATLAB的程式碼實現,有時間的話可能會出個視訊教程。     首先說下介面設計,我這裡只是說下都用了什麼元件以及一些引數,具體的GUI怎麼設計的我這裡就不詳細介紹,又需要介紹的可以給我留言我可以抽出點時間

線性規劃中的單純形法與內點法(原理、步驟以及matlab實現)(三)

應用 最大化 round 並不是 兩個 生產 陰影 3.3 ima 在本系列的第三篇博客中,筆者討論對偶單純形法的相關理論和應用 2.3 Dual Simplex Method(對偶單純形法) Contents   2.3.1 對偶問題產生的原因   2.3.2 對偶問題的

描述性統計matlab實現

pre tool nes http 調用 一點 bsp log 文件 理論講的再多不會做也白弄 直接上手 一.針對接近正態分布的(均值,方差,標準差,極差,變異系數,偏度,峰度) 這裏我必須提前說明一點就是,你在寫好函數後,函數的名是dts,你保存的文件名也必須是dts.m

Excel值班表方案:使用Excel進行不同班次排班並實現多功能自動化統計及檢視轉換

標題 本文介紹了用Excel進行多班次排班的一種方案,由於格式轉換問題,本文采用了圖片展示的形式,有關表格所用到的公式,煩請大家自己手動敲打了,這樣一來,還更方便大家加深記憶哦n_n. (本文省略掉了如何根據“班次轉置”工作表的內容上傳到釘釘考勤排班系統,因為操

卷積神經網路CNN原理——結合例項matlab實現

卷積神經網路CNN是深度學習的一個重要組成部分,由於其優異的學習效能(尤其是對圖片的識別)。近年來研究異常火爆,出現了很多模型LeNet、Alex net、ZF net等等

KNN分類演算法原理及其Matlab實現

KNN演算法原理 K近鄰演算法是一種簡單的監督學習演算法。對於給定測試樣本,直接計算該樣本和訓練集的距離,將距離最近的k個“鄰居”點的類別作為參考,作為預測結果返回。 測試資料 程式碼 matlab版

高斯消元法原理Matlab實現

直接法解線性方程組-高斯消元法 1.高斯消元法思想 設有線性方程組如下所示: ⎧⎩⎨⎪⎪⎪⎪⎪⎪a11x1+a12x2+⋯+a1nxn=b1,a21x1+a22x2+⋯+a2nxn=b2,⋮an1x1+an2x2+⋯+annxn=bn,

Android中實現「類方法指令抽取方式」加固方案原理解析

一、前言Android中加固方案一直在進步,因為新的加固方案出來就會被人無情的破解脫殼了,從第一代加固方案落地加密dex檔案,第二代加固方案不落地加密dex檔案,在到第三代加固方案類方法抽取,以後後面的

區域生長演算法原理MATLAB實現

% Segment based on area, Region Growing; clear all; close all; clc [fileName,pathName] = uigetfile('*.*','Please select an image');%檔案筐,選擇檔案 if(fileName)

Android的延遲實現的幾種解決方案以及原理分析

寫這篇文章的目的,是看到群裡有人在實現延遲的時候,用如下的第四種方法,個人感覺有點不妥,為了防止更多的人有這種想法,所以自己抽空深入分析,就分析的結果,寫下此文,希望對部分人有啟示作用。 1.實現延遲的幾種方法? 答: 1.java.util.Time

直方圖均衡原理Matlab實現

直方圖均衡(Histogram Equalization) 預備知識(PMF、CDF) 概述 實現 Matlab程式碼 預備知識(PMF、CDF) 1.PMF 概率質量函式(Probability Mass Function)是離散隨機變數在各特

K-means原理Matlab實現

前言 作為勵志在機器學習界闖出一片天地的小女子,在此整理學習到的機器學習方法,並使用Matlab及Python實現。希望可以和大家互相交流和探討。 聯絡方式:[email protected] K-means原理 k-means應該是入

傅立葉變換的原理、意義以及如何用Matlab實現快速傅立葉變換

本帖最後由 xiaoliu 於 2011-7-28 21:00 編輯一、傅立葉變換的由來關於傅立葉變換,無論是書本還是在網上可以很容易找到關於傅立葉變換的描述,但是大都是些故弄玄虛的文章,太過抽象,盡是一些讓人看了就望而生畏的公式的羅列,讓人很難能夠從感性上得到理解,最近,

方差分析的數學原理matlab實現

方差分析的概念 方差分析(Analysis of Variance,ANOVA),又稱“變異數分析”,1920年由波蘭數學家R.A.Fisher發明,通常用於兩個以上的樣本引數差別的顯著性檢驗。 為什麼要使用方差分析 統計學中的假設檢驗理論已經為樣本引數為兩個

邏輯迴歸原理matlab實現

個人部落格文章連結:  http://www.huqj.top/article?id=163 對於某些分類問題,自變數可能是連續的,但是因變數卻可能是離散的,例如:根據腫瘤大小判斷該腫瘤是否是良性。這種問題不適合用線性迴歸來解決,雖然可以將連續的因變數值對映到離散的分類上,但

BP神經網路原理matlab實現

個人部落格文章連結: http://www.huqj.top/article?id=166 BP神經網路,即“反向傳播”神經網路,是一種被廣泛運用的人工神經網路,它的思想就是通過反向傳播不斷調整模型權重,最終使得模型輸出與預期輸出的誤差控制在一個小範圍內。其中反向傳播的演算法(BP演算

吉布斯取樣——原理matlab實現

MCMC(馬爾可夫鏈蒙特卡洛方法):the Gibbs Sampler(吉布斯取樣)         在之前的部落格中,我們對比了在一個多元概率分佈p(x)中,分別採用分組(block-wise)和分成分(component-wise)實現梅特羅波利斯哈斯廷斯演算法。對

網站統計中的資料收集原理實現_埋點統計

網站資料統計分析工具是網站站長和運營人員經常使用的一種工具,比較常用的有谷歌分析、百度統計和騰訊分析等等。所有這些統計分析工具的第一步都是網站訪問資料的收集。目前主流的資料收集方式基本都是基於javascript的。本文將簡要分析這種資料收集的原理,並一步一步實際搭建一個

用Burg法估計AR模型的引數原理詳解及matlab實現

用Burg法估計AR模型的引數。藉助如圖所示的格型預測誤差濾波器,伯格法通過求出前向預測誤差和後向預測誤差的平均功率來選取最佳的反射係數k,使誤差的平均功率取得最小值,進而通過反饋求出模型係數和噪聲方差,該演算法可直接通過分析觀測資料得到需要的模型引數,而不用取求解計算量較大