24СʱÈÈÃŰæ¿éÅÅÐаñ    

Znn3bq.jpeg
²é¿´: 5185  |  »Ø¸´: 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 µÄÖ÷Ìâ¸üÐÂ
×î¾ßÈËÆøÈÈÌûÍÆ¼ö [²é¿´È«²¿] ×÷Õß »Ø/¿´ ×îºó·¢±í
[¿¼ÑÐ] È˹¤ÖÇÄÜ320µ÷¼Á08¹¤À໹Óлú»áÂð +13 Õñ¡ªTZ 2026-04-10 14/700 2026-04-13 11:25 by wp06
[¿¼ÑÐ] Çóµ÷¼Á +7 ºÎÆøÕý 2026-04-13 8/400 2026-04-13 09:33 by ÁãÁã¶þ
[¿¼ÑÐ] ²ÄÁÏ¿¼Ñе÷¼Á +28 ÔÆÄ¾´ï´ï 2026-04-11 30/1500 2026-04-13 08:51 by Sammy2
[½Ìʦ֮¼Ò] ɽ¶«Ë«·ÇԺУ¿¼ºË³¬¼¶ÎÞµ×Ïߣ¬Áìµ¼ÐÒÔÖÀÖ»ö£¬½ÌʦÔâÑê¿Ö +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[ÕÒ¹¤×÷] ɽ¶«¸ßУ½Ìʦ¿¼ºË³¬¼¶ÎÞµ×Ïߣ¬Ô±¹¤¹ý²»ÏÂÈ¥À² +4 qut2026 2026-04-09 9/450 2026-04-12 00:54 by qut2026
[¿¼ÑÐ] 086003µ÷¼ÁÇóÖú +21 ËÕß®Íò 2026-04-09 22/1100 2026-04-11 20:25 by dongdian1
[¿¼ÑÐ] Çóµ÷¼Á +11 ôæôæÒ»ÊéÉú 2026-04-09 11/550 2026-04-11 19:57 by ÄæË®³Ë·ç
[¿¼ÑÐ] 283Çóµ÷¼Á +22 ÄǸöàà×Ó 2026-04-09 22/1100 2026-04-11 10:41 by ÄæË®³Ë·ç
[¿¼ÑÐ] 22408µ÷¼ÁÇóÖú +7 ì±12 2026-04-09 9/450 2026-04-11 09:23 by ŶŶ123
[¿¼ÑÐ] ²ÄÁÏÀà284µ÷¼Á +40 Ïë»»ÊÖ»ú²»Ïë½âÊ 2026-04-08 48/2400 2026-04-10 23:28 by 314126402
[¿¼ÑÐ] µ÷¼Á +19 СÕÅZA 2026-04-10 20/1000 2026-04-10 22:08 by Öí»á·É
[¿¼ÑÐ] ¼ÆËã»úÀàÇóµ÷¼Á£¬22408-274·Ö +7 Éϰ¶deС³æ 2026-04-09 8/400 2026-04-10 19:56 by fxue1114
[¿¼ÑÐ] 298Çóµ÷¼Á +13 ¶¤¶£ß˶¬¹Ï 2026-04-09 13/650 2026-04-10 15:49 by jiajinhpu
[¿¼ÑÐ] »úе»¹Óл¹ÓÐÃû¶îÂð£¿Ì«ÄÑÁË +6 ЦЦԬ 2026-04-10 6/300 2026-04-10 11:54 by ¸ßά´º
[¿¼ÑÐ] µ÷¼ÁÉêÇë086000Ò»Ö¾Ô¸Î÷±±Å©ÁֿƼ¼´óѧÉúÎïÓëÒ½Ò©320·Ö-±¾¿ÆÆë³¹¤Òµ´óѧ +3 ÃÀÃÀŮʿ 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[¿¼ÑÐ] 086000ÉúÎïÓëÒ½Ò©µ÷¼Á +7 awwwwwooooo 2026-04-09 7/350 2026-04-09 13:31 by ±±¼«159263
[¿¼ÑÐ] 286Çóµ÷¼Á +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by ŶŶ123
[¿¼ÑÐ] 304Çóµ÷¼Á +10 ËØÄê¼ÀÓï 2026-04-06 17/850 2026-04-08 09:05 by À¶ÔÆË¼Óê
[¿¼ÑÐ] Çóµ÷¼Á +15 ÐܶþÏëÉϰ¶ 2026-04-06 15/750 2026-04-08 04:53 by Î޼ʵIJÝÔ­
[¿¼ÑÐ] ÇóÖú +3 ¿¨¿¨¶«88 2026-04-06 4/200 2026-04-06 15:28 by going home
ÐÅÏ¢Ìáʾ
ÇëÌî´¦ÀíÒâ¼û