Imread で　imcrop を使った画像　本文のプログラムで　エラーは出ないが画像が出ない

Question

Yoshifumi Harada 2019-7-22

　次のプログラムでグレーイすケールの画像変換をしたプログラムの画像が出ない。

切り抜きの、幅と高さは、わかりますが、Xmin, Ymin の決定方法を教えて下さい。

　　TP_Blackhol_2019.jpg の画像が表示されません　切り抜きの設定方法を教えて下さい

% -- 読み込み、切り出し --------------------

I = imread('VECTOR_FIELD_KP002.JPG'); % 画像の読み込み

I = rgb2gray(I); % グレースケールの画像に変換

I_cr = imcrop(I, [403, 151, 470, 279]); % 切り抜き

imshow(I_cr) % 切り抜いた箇所を表示

% -- 処理 ---------------------------------

I_db = double(I_cr); % FFTを行うため、uint8型からdouble型に変換

Y = fft2(I_db); % 2D-FFT

Ys = fftshift(Y); % 直流成分を中心にシフト

Yx = 20 * log10(abs(Ys)); % 対数に変換

% -- 表示 ----------------------------------

figure

imagesc(Yx) % 強度データを表示

colormap(gray)

axis square

xlim([515,815]) % X軸の表示範囲の設定

ylim([515,815]) % Y軸の表示範囲の設定

% ---------------- 追加 ---------------------

figure

plot(trapz(Yx,2))

title('X方向の積分')

figure

plot(trapz(Yx,1))

title('Y方向の積分')

Imread で imcrop を使った画像 本文の​プログラムで エラー​は出ないが画像が出な​い