²é¿´: 5172  |  »Ø¸´: 3

norbie2008

гæ (³õÈëÎÄ̳)

[ÇóÖú] ÈçºÎ¸øÍ¼Ïñ¼ÓÒ»¸ö´óС¿É¿ØµÄ²´ËÉÔëÉù ÒÑÓÐ1È˲ÎÓë

ÈçÌ⣬Çë½Ì¸ßÊÖ
matlabÖÐÓÃimnoise²úÉúµÄpoissonÔëÉù£¬µ«ÔëÉù´óСÓÉÊäÈëͼÏñ×Ô¶¯Éú³É£¬²»¿É¿ØÖÆ¡£
ÎÒÏëÒª¿ÉÓɲÎÊý¿ØÖƵÄpoissonÔëÉù£¬´Óδ²úÉúº¬ÔëͼÏñ£¬ÓÃÓÚͼÏñ½µÔë¡£
лл´ó¼ÒÁË
»Ø¸´´ËÂ¥

» ²ÂÄãϲ»¶

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

feixiaolin

ÈÙÓþ°æÖ÷ (ÎÄ̳¾«Ó¢)

ÓÅÐã°æÖ÷

¡¾´ð°¸¡¿Ó¦Öú»ØÌû

¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï ¡ï
¸Ðл²ÎÓ룬ӦÖúÖ¸Êý +1
norbie2008(nono2009´ú·¢): ½ð±Ò+10 2017-07-18 21:53:44
P_noise0= numda *sqrt(-2* log( 1-x) );
numda  ×Ô¼º¶¨£¬xΪ [0, 1] ¾ù²¼Ëæ»úÔëÉù¡£
if(P_noise > 255)
      P_noise =255;      
if(P_noise < 0 )
      P_noise=0;
2Â¥2014-04-16 21:34:01
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

¹âÐÅСÉú

½ð³æ (СÓÐÃûÆø)

¿ÉÒÔÓÃpoissrnd£¬²»¹ýÕâÑùȫͼ¶¼ÊÇͬһ¸ö²ÎÊý

·¢×ÔСľ³æAndroid¿Í»§¶Ë
3Â¥2017-07-17 10:35:46
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû

ÀîÕð±±¾©Àí¹¤

Í­³æ (СÓÐÃûÆø)

ͼÏñ²´ËÉÔëÉùÊÇû·¨¿ØÖÆ´óСµÄ£¬²´ËÉÔëÉùÊǺÍÄãͼÏñÁÁ¶È¹Ò¹³µÄ£¬Í¼ÏñÔ½ÁÁ£¬²´ËÉÔëÉù¾Í½üËÆ¸ß˹·Ö²¼ÁË£¬Í¼ÏñÔ½°µ£¬ÆäÓ°ÏìÔ½ÑÏÖØ£¬Òò´Ëµ¥ÂÛÒ»ÕÅͼÀ´Ëµ£¬¿ØÖÆÆä²´ËÉÔëÉùµÄ´óСÊÇûÓÐÒâÒåµÄ¡£
Èç¹ûÄã·ÇÒª¸Ä±ä²´ËÉÔëÉù£¬¿ÉÒÔÓÐÈçÏ´úÂë
function [ b ] = poisson_fun( a,lamda )
    sizeA = size(a);
    a = double(a();

    %  (Monte-Carlo Rejection Method) Ref. Numerical
    %  Recipes in C, 2nd Edition, Press, Teukolsky,
    %  Vetterling, Flannery (Cambridge Press)
   
    b = zeros(size(a),'like', a);
    idx1 = find(a<50); % Cases where pixel intensities are less than 50 units
    if (~isempty(idx1))
      g = exp(-a(idx1)*lamda);
      em = -ones(size(g));
      t = ones(size(g));
      idx2 = (1:length(idx1))';
      while ~isempty(idx2)
        em(idx2) = em(idx2) + 1;
        t(idx2) = t(idx2) .* rand(size(idx2));
        idx2 = idx2(t(idx2) > g(idx2));
      end
      b(idx1) = em;
    end

    % For large pixel intensities the Poisson pdf becomes
    % very similar to a Gaussian pdf of mean and of variance
    % equal to the local pixel intensities. Ref. Mathematical Methods
    % of Physics, 2nd Edition, Mathews, Walker (Addison Wesley)
    idx1 = find(a >= 50); % Cases where pixel intensities are at least 50 units
    if (~isempty(idx1))
      b(idx1) = round(a(idx1) + sqrt(lamda*a(idx1)) .* randn(size(idx1)));
    end
   
    b =uint8(reshape(b,sizeA));


end
You_are_my_quest
4Â¥2018-05-02 11:19:26
ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢TA ¸øTA·¢ÏûÏ¢ ËÍTAºì»¨ TAµÄ»ØÌû
Ïà¹Ø°æ¿éÌø×ª ÎÒÒª¶©ÔÄÂ¥Ö÷ norbie2008 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 291 Çóµ÷¼Á +4 »¯¹¤2026½ì±ÏÒµÉ 2026-03-21 4/200 2026-03-23 09:59 by bingchuan
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤ +10 °²È«Éϰ¶£¡ 2026-03-16 10/500 2026-03-23 09:31 by Jy?
[¿¼ÑÐ] 291Çóµ÷¼Á +5 ‹üÈA 2026-03-22 5/250 2026-03-23 09:20 by haoshis
[¿¼ÑÐ] 317Çóµ÷¼Á +12 Éê×ÓÉêÉê 2026-03-19 18/900 2026-03-22 22:23 by luoyongfeng
[¿¼ÑÐ] 310Çóµ÷¼Á +4 baibai1314 2026-03-16 4/200 2026-03-22 20:19 by edmund7
[¿¼ÑÐ] 307Çóµ÷¼Á +11 ÀäóÏ123 2026-03-17 11/550 2026-03-22 20:16 by edmund7
[¿¼ÑÐ] 315·Ö£¬³ÏÇóµ÷¼Á£¬²ÄÁÏÓ뻯¹¤085600 +3 13756423260 2026-03-22 3/150 2026-03-22 20:11 by edmund7
[¿¼ÑÐ] ÉúÎïѧһ־Ը985£¬·ÖÊý349Çóµ÷¼Á +4 zxts12 2026-03-21 7/350 2026-03-22 09:57 by zxts12
[¿¼ÑÐ] ²ÄÁϹ¤³Ìר˶ 348·ÖÇóµ÷¼Á +3 ¶¬´Ç. 2026-03-17 5/250 2026-03-21 18:47 by ѧԱ8dgXkO
[¿¼ÑÐ] Ò»Ö¾Ô¸Î÷ÄϽ»´ó£¬Çóµ÷¼Á +5 ²Ä»¯ÖðÃÎÈË 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[¿¼ÑÐ] ÄϾ©´óѧ»¯Ñ§376Çóµ÷¼Á +3 hisfailed 2026-03-19 6/300 2026-03-20 23:43 by hisfailed
[¿¼ÑÐ] ¿¼Ñе÷¼ÁÇóÑ§Ð£ÍÆ¼ö +3 ²®ÀÖ29 2026-03-18 5/250 2026-03-20 22:59 by JourneyLucky
[¿¼ÑÐ] Çóµ÷¼Á +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[¿¼ÑÐ] 298-Ò»Ö¾Ô¸Öйúũҵ´óѧ-Çóµ÷¼Á +9 ÊÖ»úÓû§ 2026-03-17 9/450 2026-03-20 14:24 by ÎÞи¿É»÷111
[¿¼ÑÐ] ÉúÎïѧµ÷¼ÁÕÐÈË£¡£¡£¡ +3 ɽº£Ììá° 2026-03-17 4/200 2026-03-19 21:34 by ÔõôÊÍ»³
[¿¼ÑÐ] 288Çóµ÷¼Á£¬Ò»Ö¾Ô¸»ªÄÏÀí¹¤´óѧ071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[¿¼ÑÐ] 085600²ÄÁÏÓ뻯¹¤µ÷¼Á 324·Ö +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[¿¼ÑÐ] ¿¼ÑÐÇóµ÷¼Á +3 éÙËÌ. 2026-03-17 4/200 2026-03-17 21:43 by ÓÐÖ»ÀêÅ«
[ÂÛÎÄͶ¸å] ÓÐûÓдóÀз¢Ð¡ÂÛÎÄÄÜ´øÎÒ¸ö¶þ×÷ +3 ÔöÈñ©ÈË 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[¿¼ÑÐ] 11408 Ò»Ö¾Ô¸Î÷µç£¬277·ÖÇóµ÷¼Á +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û