²é¿´: 5165  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 0703 µ÷¼Á +3 ÎÒ¿ÉÒÔÉϰ¶µÄ¶Ô 2026-03-16 6/300 2026-03-23 08:52 by ×íÔÚ·çÀï
[¿¼ÑÐ] 0854µç×ÓÐÅÏ¢Çóµ÷¼Á +3 ¦Á____ 2026-03-22 3/150 2026-03-22 21:28 by zhq0425
[¿¼²©] ÕÐÊÕ²©Ê¿1-2ÈË +3 QGZDSYS 2026-03-18 4/200 2026-03-22 10:25 by QGZDSYS
[¿¼ÑÐ] ²ÄÁÏÇóµ÷¼Á +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[¿¼ÑÐ] Ò»Ö¾Ô¸ÄÏ´ó£¬0703»¯Ñ§£¬·ÖÊý336£¬Çóµ÷¼Á +3 ÊÕµ½VS 2026-03-21 3/150 2026-03-21 18:42 by ѧԱ8dgXkO
[¿¼ÑÐ] 0703»¯Ñ§297Çóµ÷¼Á +3 Daisy¡î 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[¿¼ÑÐ] 265Çóµ÷¼Á +12 ÁºÁºÐ£Ð£ 2026-03-19 14/700 2026-03-21 13:38 by lature00
[¿¼ÑÐ] 279·ÖÇóµ÷¼Á Ò»Ö¾Ô¸211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[¿¼ÑÐ] 316Çóµ÷¼Á +6 ÁºÜçö© 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666£¡
[¿¼ÑÐ] 332Çóµ÷¼Á +4 ydfyh 2026-03-17 4/200 2026-03-21 02:20 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸Äϲý´óѧ£¬327·Ö£¬²ÄÁÏÓ뻯¹¤085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[¿¼ÑÐ] 294Çóµ÷¼Á²ÄÁÏÓ뻯¹¤×¨Ë¶ +15 ݤÎÉ­ÁÖ 2026-03-18 15/750 2026-03-20 23:28 by JourneyLucky
[¿¼ÑÐ] 330Çóµ÷¼Á +4 С²Ä»¯±¾¿Æ 2026-03-18 4/200 2026-03-20 23:13 by JourneyLucky
[¿¼ÑÐ] Ò»Ö¾Ô¸ Î÷±±´óѧ £¬070300»¯Ñ§Ñ§Ë¶£¬×Ü·Ö287£¬Ë«·ÇÒ»±¾£¬Çóµ÷¼Á¡£ +4 ³¿»èÏßÓëÐǺ£ 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[¿¼ÑÐ] ÖÐÄÏ´óѧ»¯Ñ§Ñ§Ë¶337Çóµ÷¼Á +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[¿¼ÑÐ] Ò»Ö¾Ô¸Î÷°²½»Í¨´óѧ ѧ˶ 354Çóµ÷¼Á211»òÕß˫һÁ÷ +3 ÎÒÏëÒª¶ÁÑо¿Éú 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[¿¼ÑÐ] Çóµ÷¼Á +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
[¿¼ÑÐ] 085410È˹¤ÖÇÄÜר˶317Çóµ÷¼Á£¨0854¶¼¿ÉÒÔ£© +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by ²»168
[¿¼ÑÐ] 286Çóµ÷¼Á +6 lemonzzn 2026-03-16 10/500 2026-03-19 14:31 by lemonzzn
[¿¼ÑÐ] 070300»¯Ñ§Ñ§Ë¶Çóµ÷¼Á +6 Ì«Ïë½ø²½ÁË0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û