从DEM中提取数据时一个点距无法正常输出

龙行天下

需要从DEM提取相应栅格的属性值，代码如下

[Fortran] 纯文本查看 复制代码

查看代码

打印?

01	PROGRAM MAIN

02	IMPLICIT NONE

03	INTEGER::I,J,K

04	INTEGER::DEMCOLS,DEMROWS    !DEM数据列数、行数

05	INTEGER::DEMNODATA          !DEM无效值

06	INTEGER,ALLOCATABLE,DIMENSION(:,:)::PROVALUE

07	REAL::CORNERLON,CORNERLAT,CELLSIZE

08	REAL,ALLOCATABLE,DIMENSION(:,:)::DEMLON,DEMLAT

09	CHARACTER(LEN=30)::DEMTEMP

10	OPEN(UNIT=1201503,FILE="DATA.txt")

11	READ(1201503,*)DEMTEMP,DEMCOLS

12	READ(1201503,*)DEMTEMP,DEMROWS

13	READ(1201503,*)DEMTEMP,CORNERLON

14	READ(1201503,*)DEMTEMP,CORNERLAT

15	READ(1201503,*)DEMTEMP,CELLSIZE

16	READ(1201503,*)DEMTEMP,DEMNODATA

17	ALLOCATE(PROVALUE(DEMROWS,DEMCOLS),DEMLON(DEMROWS,DEMCOLS),DEMLAT(DEMROWS,DEMCOLS))

18	DO I=1,DEMROWS

19	DEMLAT(I,:)=CORNERLAT+(DEMROWS-I)*CELLSIZE

20

END DO

21	DO I=1,DEMCOLS

22	DEMLON(:,I)=CORNERLON+(I-1)*CELLSIZE

23

END DO

24	DO I=1,DEMROWS

25	READ(1201503,*)(PROVALUE(I,J),J=1,DEMCOLS)

26

END DO

27	OPEN(UNIT=2201501,FILE="ResultsCAREPOINTS.txt")

28	DO J=1,DEMROWS

29	DO K=1,DEMCOLS

30	IF((DEMLON(J,K)==101.95.AND.DEMLAT(J,K)==27.05).OR.(DEMLON(J,K)==101.45.AND.DEMLAT(J,K)==26.95).OR.(DEMLON(J,K)==101.75.AND.DEMLAT(J,K)==26.65))THEN

31	WRITE(2201501,'(2F10.3,I10)')DEMLON(J,K),DEMLAT(J,K),PROVALUE(J,K)

32

END IF

33

END DO

34

END DO

35	END PROGRAM

数据如下： DATA.zip (31.26 KB, 下载次数: 1)

2015-4-13 15:06 上传

点击文件名下载附件

我需要提取经纬度坐标为（101.95,27.05），（101.45,26.95），（101.75，26.65）三个点的属性值
但是，结果只能提取出其中的两个点，如图

请各位大侠，帮我看一下
谢谢

vvt

编程时，时刻记住，浮点数是有误差的。

所以不做浮点数相等判断。

而改用差的绝对值小于某个很小的数。

[Fortran] 纯文本查看 复制代码

查看代码

打印?

1	Real , parameter :: eps = 1.0e-5

2	If ( (abs(demlon(j,k)-101.95)<eps.and.abs(demlat(j,k)-27.05)<eps) .or. &

3	(abs(demlon(j,k)-101.45)<eps.and.abs(demlat(j,k)-26.95)<eps) .or. &

4	(abs(demlon(j,k)-101.75)<eps.and.abs(demlat(j,k)-26.65)<eps) )then

5	WRITE(2201501,'(2F10.3,I10)')DEMLON(J,K),DEMLAT(J,K),PROVALUE(J,K)

6

END IF

龙行天下

脑子糊涂了，咋忘了这一条，还忙活了半天，判断其他地方呢

谢谢啦