Znn3bq.jpeg
ÉÇÍ·´óѧº£Ñó¿ÆÑ§½ÓÊܵ÷¼Á
²é¿´: 5187  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] 085404 298·ÖÇóµ÷¼Á +11 ºôÀ²ºôÀ²ºôºôºô 2026-04-10 12/600 2026-04-14 08:38 by wfj257
[¿¼ÑÐ] 0854µ÷¼Á +10 ³¤¹­°Á 2026-04-11 11/550 2026-04-13 10:38 by wp06
[¿¼ÑÐ] ¿¼ÑжþÂÖµ÷¼Á +10 ¹ÊÈË?? 2026-04-09 10/500 2026-04-13 09:55 by szhize
[¿¼ÑÐ] 22408 352·ÖÇóµ÷¼Á +5 ŬÁ¦µÄÏÄÄ© 2026-04-09 5/250 2026-04-12 19:17 by wj165256
[¿¼ÑÐ] ÉúÎïѧ308Çóµ÷¼Á +5 ÏàÐűػá¹ââÍòÕ 2026-04-11 5/250 2026-04-12 18:14 by zhouxiaoyu
[¿¼ÑÐ] 305Çóµ÷¼Á +7 Â꿨°Í¿¨boom 2026-04-11 7/350 2026-04-12 07:35 by zhouxiaoyu
[¿¼ÑÐ] µ÷¼Á +6 ÇàµÆ²»¸º 2026-04-09 6/300 2026-04-11 20:35 by dongdian1
[¿¼ÑÐ] 086003µ÷¼ÁÇóÖú +21 ËÕß®Íò 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[¿¼ÑÐ] 270Çóµ÷¼Á +14 ÑîÀÖ369 2026-04-11 14/700 2026-04-11 20:16 by À¶ÔÆË¼Óê
[¿¼ÑÐ] 11408¡£358Çóµ÷¼Á +3 TMYzds 2026-04-07 3/150 2026-04-11 17:10 by µªÆøÆøÆø
[¿¼ÑÐ] ÖпÆÔº×Ü·Ö315Çóµ÷¼Á +8 lallalh 2026-04-09 8/400 2026-04-10 19:30 by dick_runner
[¿¼ÑÐ] 071000ÉúÎïѧµ÷¼ÁÇóÖú +17 zzzzwww 2026-04-09 20/1000 2026-04-10 15:55 by Çóµ÷¼Ázz
[¿¼ÑÐ] 085800 ÄÜÔ´¶¯Á¦Çóµ÷¼Á +6 °¢biu°¡°¡°¡°¡°¡ 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[¿¼ÑÐ] Ò»Ö¾Ô¸ÖÐÄÏ´óѧÎïÀíѧ£¬Ó¢Ò»66£¬Çóµ÷¼Á +4 ³¤ÑÌì½ì» 2026-04-08 5/250 2026-04-10 10:31 by Ó±¹û¶ù
[¿¼ÑÐ] ½­ËÕ´óѧ ¹¤¿Æµ÷¼Á ¼ñ© +3 Evan_Liu 2026-04-09 5/250 2026-04-10 10:22 by Evan_Liu
[¿¼²©] ²©Ê¿×Ô¼ö +7 ¿É¿ÉСÅÖ 2026-04-08 7/350 2026-04-10 08:28 by kimhero
[¿¼ÑÐ] 332£¬085601Çóµ÷¼Á +12 ydfyh 2026-04-09 14/700 2026-04-09 17:28 by wp06
[¿¼ÑÐ] Çóµ÷¼Á +13 Æâluck 2026-04-07 13/650 2026-04-08 22:46 by Öí»á·É
[¿¼ÑÐ] 331Çóµ÷¼Á +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
[¿¼ÑÐ] 307Çóµ÷¼Á +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by СºÚ²»ÅÂÄÑ
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û