计算二维数组坐标“岛”的重心问题-Fortran算法研究-专业Fortran论坛

新手求指点 发表于 2014-3-6 13:34:37

计算二维数组坐标“岛”的重心问题

我刚接触fortran没多久，现在要写一个程序来算出下面这个数组的重心
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 1 1 1 1 1 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

我想把这个数组其中所有的1的值用坐标表示出来,比如（14，15）这样子，但是我不太清楚该怎么筛选和表示，求大神指点

楚香饭 发表于 2014-3-6 13:37:33

Do i = 1 , n
Do j = 1 , m
if ( a(i,j) == 1 ) then
!此时的 i j 即为坐标
end if
End Do
End Do

新手求指点 发表于 2014-3-6 13:43:38

chuxf 发表于 2014-3-6 13:37
Do i = 1 , n
Do j = 1 , m
if ( a(i,j) == 1 ) then

0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 1 1 1 0
如果是这样的，我想用程序判断上面的那些1把他们-20变到下面来该怎么做呢？

楚香饭 发表于 2014-3-6 13:45:19

什么叫把他们 -20 变到下面来？

新手求指点 发表于 2014-3-6 13:55:14

chuxf 发表于 2014-3-6 13:45
什么叫把他们 -20 变到下面来？

如果把所有的1看成一个岛，我想算这个岛重心，我想把所有的1的坐标全部相加最后除以25就是这个岛的重心，但是这是一个周期性的岛，就像第二组数组，有一部分1在上面，比如(5,20)，但实际这个1是该在(5,0)这个坐标上，但是因为下面无法表示了所以这个1就被加了20跑到上面去了，但在计算重心的时候需要把这个1减掉20，就是从（5，20）变到（5，0）来计算重心

楚香饭 发表于 2014-3-6 14:07:32

本帖最后由 chuxf 于 2014-3-6 14:15 编辑

这个问题可以很简单，也可以很复杂。

这取决于：
第一：整个区域内，你是否可以确定只存在一个“岛”？还是有可能有多个？
比如：
00000000000001111111100000000000
00000000000000011110000000000000
00000000000000000000000000000000
00000000110000000000000000000000
00000001111000000000000000000000
00000000000000000000000000000000
第二：岛会不会出现无限循环。
比如：
00000000000001111000000000000000
00000000000001111000000000000000
00000000001111111000000000000000
00000000000001111000000000000000
第三：岛中间是否可能出现空洞？
比如：
00000000000000000000000000000000
00000000000001111000000000000000
00000000000001111000000000000000
00000000001111001000000000000000
00000000000001111000000000000000
00000000000000000000000000000000

新手求指点 发表于 2014-3-6 14:26:03

chuxf 发表于 2014-3-6 14:07
这个问题可以很简单，也可以很复杂。

这取决于：

只有一个岛，就是加了周期或者减了周期后就只有一个岛了，不会有2个，岛中间没有空洞，无限循环那个不太懂。

楚香饭 发表于 2014-3-6 15:24:54

本帖最后由 chuxf 于 2014-3-6 19:22 编辑

我这里有个死办法，那就是如果有上下左右的联通，则移动整个区域，直到不再联通为止。（记录下 X Y 的偏移量）
最后算出重心，再把偏移量加回去。

Program www_fcode_cn
Implicit None
Integer , parameter :: N = 4
Integer , parameter :: M = 4
real :: gx , gy
Integer :: a( M , N )
Open( 12 , File = '文件名' )
Read( 12 , * ) a
Close( 12 )
call IslandGravity( a , N , M , gx , gy )
write(*,*) gx , gy

contains

Subroutine IslandGravity( d , n , m , gx , gy )
Integer :: n , m
Integer :: d( m , n ) , a( m , n )
integer :: i , j , c , iShiftX , iShiftY
Real :: gx , gy
Logical :: blink
a = d !// 复制一份 a

iShiftX = 0
iShiftY = 0
Do!// 向上滚动，直到没有上下链接为止
   blink = Check_Y_Link( a , n , m )
   if ( blink ) then
   call RollUp( a , n , m )
   iShiftY = iShiftY + 1
   else
   Exit
   end if
End Do
Do!// 向左滚动，直到没有左右链接为止
   blink = Check_X_Link( a , n , m )
   if ( blink ) then
   call RollLeft( a , n , m )
   iShiftX = iShiftX + 1
   else
   Exit
   end if
End Do
gx = 0.
gy = 0.
c = 0
Do i = 1 , N
   Do j = 1 , M
   if ( a(j,i) == 1 ) then
      gx = gx + j
      gy = gy + i
      c = c + 1
   end if
   End Do
End Do
gx = ( gx / c ) + iShiftX !// 最后再把偏移量加上去
if ( gx > m ) gx = gx - m
gy = ( gy / c ) + iShiftY
if ( gy > n ) gy = gy - n
End Subroutine IslandGravity

Subroutine RollUp( a , n , m ) !// 区域向上滚动
Integer :: n , m
Integer :: a( m , n ) , t( m )
t(:) = a( : , 1 )
a(:,1:n-1) = a(:,2:n)
a(:,n) = t(:)
End Subroutine RollUp

Subroutine RollLeft( a , n , m ) !// 区域向左滚动
Integer :: n , m
Integer :: a( m , n ) , t( n )
t(:) = a( 1 , : )
a(1:m-1,:) = a(2:m,:)
a(m,:) = t(:)
End Subroutine RollLeft

Logical Function Check_Y_Link( a , n , m ) !// 检查是否有上下链接
Integer :: n , m
Integer :: a( m , n )
integer :: i
Check_Y_Link = .false.
Do i = 1 , m
   if ( ( a(i,1) == 1 ) .and. ( a(i,n) == 1 ) ) then
   Check_Y_Link = .true.
   return
   end if
End Do
End Function Check_Y_Link

Logical Function Check_X_Link( a , n , m ) !// 检查是否有左右链接
Integer :: n , m
Integer :: a( m , n )
integer :: i
Check_X_Link = .false.
Do i = 1 , n
   if ( ( a(1,i) == 1 ) .and. ( a(m,i) == 1 ) ) then
   Check_X_Link = .true.
   return
   end if
End Do
End Function Check_X_Link

End Program www_fcode_cn

joejimming 发表于 2014-3-6 19:01:55

不清楚这里重心的意义为何，感觉有一点点问题：
例如：
1101
1101
0000
1101
重心按个人理解应该是1即坐标(1,1)，算出来的结果应该是(5,5)，超出了范围。

楚香饭 发表于 2014-3-6 19:23:39

joejimming 发表于 2014-3-6 19:01
不清楚这里重心的意义为何，感觉有一点点问题：
例如：
1101

感谢楼上指正。

已增加两行代码解决这一问题

gx = ( gx / c ) + iShiftX !// 最后再把偏移量加上去
if ( gx > m ) gx = gx - m
gy = ( gy / c ) + iShiftY
if ( gy > n ) gy = gy - n

页: [1] 2

Fortran Coder's Archiver

计算二维数组坐标“岛”的重心问题