1

我正在尝试编写一个利用 sobel 滤波器检测图像边缘的程序。因此，首先，我在一些基本代码中写下了一些要求，例如将 x 和 y 方向过滤器作为数组以及尝试读取 pgm 图像：

program edges
implicit none
integer, dimension(:,:), allocatable :: inp, outim, GX, GY
integer, parameter :: dp = selected_real_kind(15,300)
integer :: ky, kx, x, y, out_unit = 10, M, N, sx, sy, i, j
real(kind=dp) :: G
M = 5
N = 5

allocate(inp(M,N))
allocate(outim(M-2,N-2))
allocate(GX(3,3))
allocate(GY(3,3))

open(file = 'clown.pgm',unit=out_unit,status= 'unknown') !opening file to         write to inp
read (out_unit,11) 'P2'      !pgm magic number
read (out_unit,12) 50,50     !width, height
read (out_unit,13) 1      !max gray value
do M=-25,25
    do N=-25,25
        read (out_unit,*) inp(M,N)
    end do
end do
11 format(a2)
12 format(i3,1x,i3)
13 format(i5)

这是我第一次在 FORTRAN 中使用图像处理，除了一次我将图像打印为 pbm 文件。读取图像的代码是我之前打印出来的代码的复制品，除了我将写入更改为读取。

所以我的问题是，如何将 pgm 格式的图像读入“inp”数组，以便应用 sobel 过滤器？当我运行我的尝试时，我收到以下错误：

read (out_unit,11) 'P2'      !pgm magic number
              1
Error: Expected variable in READ statement at (1)
sobel.f90:41:18:

 read (out_unit,12) 50,50     !width, height
              1
Error: Expected variable in READ statement at (1)
sobel.f90:42:18:

 read (out_unit,13) 1      !max gray value
              1
Error: Expected variable in READ statement at (1)

谢谢

4

1 回答 1

1

正如编译器非常清楚地指出的那样，第一个错误是在这一行：

read (out_unit,11) 'P2'

编译器期望，因为这是定义 Fortran 的方式，被告知从out_unit变量中读取一个值，但'P2'它是一个字符文字（或者无论标准如何称呼它们），它是一个字符串，它不是一个变量。

同样的错误也发生在下一行。编译器期望类似

integer :: p2  ! pgm magic number
...
read (out_unit,11) p2

等等。执行此版本的 read 语句后，变量p2保存从pgm文件中读取的幻数。

当我在写的时候，调用这个单元来读取out_unit是不正常的，最终会让你感到困惑。

于 2018-03-03T17:45:29.157 回答