²é¿´: 5169  |  »Ø¸´: 3
µ±Ç°Ö»ÏÔʾÂú×ãÖ¸¶¨Ìõ¼þµÄ»ØÌû£¬µã»÷ÕâÀï²é¿´±¾»°ÌâµÄËùÓлØÌû

norbie2008

гæ (³õÈëÎÄ̳)

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

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

» ²ÂÄãϲ»¶

ÒÑÔÄ   »Ø¸´´ËÂ¥   ¹Ø×¢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µÄ»ØÌû
²é¿´È«²¿ 4 ¸ö»Ø´ð

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µÄ»ØÌû
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 306Çó0703µ÷¼ÁÒ»Ö¾Ô¸»ªÖÐʦ·¶ +6 Ö½Óãly 2026-03-21 6/300 2026-03-23 12:08 by Iveryant
[¿¼ÑÐ] ²ÄÁÏÓ뻯¹¤¿¼Ñе÷¼Á +3 ‹üÈA 2026-03-22 3/150 2026-03-23 11:20 by barlinike
[¿¼ÑÐ] 08¹¤Ñ§µ÷¼Á +6 Óû§573181 2026-03-20 10/500 2026-03-23 11:04 by Óû§573181
[¿¼ÑÐ] 291 Çóµ÷¼Á +4 »¯¹¤2026½ì±ÏÒµÉ 2026-03-21 4/200 2026-03-23 09:59 by bingchuan
[¿¼ÑÐ] 323Çóµ÷¼Á +6 ÍÝСͰ 2026-03-18 6/300 2026-03-23 00:29 by king123£¡
[¿¼ÑÐ] 311Çóµ÷¼Á +3 26ÑÐ0 2026-03-20 3/150 2026-03-22 14:46 by ColorlessPI
[¿¼ÑÐ] 318Çóµ÷¼Á +4 plumÀî×Ó 2026-03-21 7/350 2026-03-22 14:17 by ColorlessPI
[¿¼ÑÐ] ²ÄÁÏÇóµ÷¼Á +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[¿¼ÑÐ] Ò»Ö¾Ô¸Éî´ó£¬0703»¯Ñ§£¬×Ü·Ö302£¬Çóµ÷¼Á +4 ÆßÔÂ-ÆßÆß 2026-03-21 4/200 2026-03-21 18:20 by ѧԱ8dgXkO
[¿¼ÑÐ] 316Çóµ÷¼Á +6 ÁºÜçö© 2026-03-19 6/300 2026-03-21 06:32 by Ecowxq666£¡
[¿¼ÑÐ] ³õʼ318·ÖÇóµ÷¼Á£¨Óй¤×÷¾­Ñ飩 +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[¿¼ÑÐ] Çóµ÷¼Á +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[¿¼ÑÐ] Çóµ÷¼Á£¬Ò»Ö¾Ô¸:ÄϾ©º½¿Õº½Ìì´óѧ´óѧ £¬080500²ÄÁÏ¿ÆÑ§Ó빤³Ìѧ˶£¬×Ü·Ö289·Ö +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[¿¼ÑÐ] ÖÐÄÏ´óѧ»¯Ñ§Ñ§Ë¶337Çóµ÷¼Á +3 niko- 2026-03-19 6/300 2026-03-20 21:58 by luoyongfeng
[¿¼ÑÐ] 290Çóµ÷¼Á +7 ^O^Ø¿ 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[¿¼ÑÐ] 085410È˹¤ÖÇÄÜר˶317Çóµ÷¼Á£¨0854¶¼¿ÉÒÔ£© +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by ²»168
[¿¼ÑÐ] ÊÕ¸´ÊÔµ÷¼ÁÉú +4 ÓêºóÇïºÉ 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[¿¼ÑÐ] 0703»¯Ñ§µ÷¼Á +3 ÄÝÄÝninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[¿¼ÑÐ] 11408 Ò»Ö¾Ô¸Î÷µç£¬277·ÖÇóµ÷¼Á +3 zhouzhen654 2026-03-16 3/150 2026-03-17 07:03 by laoshidan
[¿¼ÑÐ] [µ¼Ê¦ÍƼö]Î÷ÄϿƼ¼´óѧ¹ú·À/²ÄÁϵ¼Ê¦ÍƼö +3 ¼â½ÇСºÉ 2026-03-16 6/300 2026-03-16 23:21 by ¼â½ÇСºÉ
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û