円を配列で表す

質問失礼します。

J = repmat(6.5,710,710);

この一律強度6.5の四角の上に強度650の円（中心,半径=(335,335),135）を描きたいです。

その後ガウス関数をかける予定です。

どの様にすれば描けますか。

どなたか教えてください。

[Hernia Baby]

円の中も塗りつぶしますか？

[英寿 野口]

全て一律強度650にしたいです

[英寿 野口]

円内強度一律650です

Answer 1

[Akira Agata]

いろいろなやり方があると思いますが、たとえばMATLABの基本関数のみを使う以下の方法はいかがでしょうか？

% 初期配列
J = repmat(6.5, 710, 710);

% 中心 (335, 335), 半径 135 の円内のグリッド点を示すインデックスを作成
[xg, yg] = meshgrid(1:710, 1:710);
idx = sqrt((xg-335).^2 + (yg-335).^2) <= 135;

% インデックス上の配列要素を650に置き換える
J(idx) = 650;

% 確認
imshow(J, [])

[英寿 野口]

ありがとうございます。出来ました。大変助かりました。