日志
IDL之求二维数组极大值
|
x=dist(50)
;删选半径，极值判断的区域是r*2+1的正方形
r=3
for i=0+r,49-r do begin
for j=0+r,49-r do begin
;格子大小是r*2+1的正方形
zone=x[i-r:i+r,j-r:j+r]
;如果正中间这个数是格子中的最大值，那它就是极值：
if(x[i,j] ge max(zone))then begin
print,'极大值:',i,j,x[i,j]
endif
endfor
endfor
;通过调整半径r可以调整极值分布的分辨率
结果：
极大值: 25 25 35.3553