遥感图像ENVI水体提取步骤

admin管理员组

文章数量:1531475

2024年7月25日发(作者：)

数据要求：

1 .下载的影像数据，尽量为同日期或者尽量靠近，不能相差时间太长，提供的影像为

2004

年第259天，1994年第295天，2004年第268天。其中1994年的影像肯定不行

LT5124tM02.004^5 9 BJCW

LT51

240391994

ag 5

EJCOO

LT5

123 04020042 &8 EJCOO

2 .下载的影像数据，尽量没有云层覆盖

类似这种研究区域中水体部分存在云层时，该影像不能用，需用接近该日期的影像替代。

水体提取步骤如下

(一)7个单波段合并成一个文件

1 .ENVI软件中File-Open Image File,弹出以下对话框，选择文件夹下b1-b7影像并打开，如

下:

Enter Data I ilenames

祖祝， 新叔坤共

T

谒 做

Elem

!!S

+ LT5123BJCM

蔻至

LT3

爵时

.，

声

S

画 /

L 31

片

;刘

B2 ■

国和

/I la^nfstacking

•

括阿

」

g

3

咨

F

素

st

白

tking

J

营不

碣

生施封

L

LT51230402004

。犯

JC0CB1

LT5 ⑵ 040200426

犯

JC00 B2

LT5 ⑵ 04020042M

犯北

00

胎

LT5 ⑵ 040200426

犯

LT5 ⑵ 04020042M

犯

JC00 B5

LT5 ⑵ 040200426

犯

JC00 B6

LT5 ⑵ 04020042M

犯

JC00 B7

LT5123(H02OM268BJCOO_GCP

LT5123(M02OM268BJCOO_Mn.

J|LT5123iM02(XH26gBJCOO_VER

LT5123M20M26

犯

JCCKLVER

计宾机

& Windmn3_D5 |

■-就暗

［6

)

r拼加卷

(E:)

r-a Elements (EJ

交

i+SONJ：

B

LT51l23B4D2D042&0O_Br -LT51Z3D /

2 .将7个波段合成一个影像文件，操作如下图:

Basic Took Clas sifi cation Transfo-rm Fi It

-

Rs-dze Data (^patisl/Spectral)

Subset Daiia via ROIs

Rotate/Flip Data

L

1

3

Layer £ta clcin g

Li”.

E

「

LI

r?in nrm

3.点击Import File,选择所有波段

5.点击 Reorder Files

Layer Stacking Parameters

Selected. Pnles for Layer Stickmc：

LT51 ^CK02flM2BaiCCm J7 7IF [Emd 1]

LT5123D4020M2BaHjrm^Ba 7IF L：B

皿[

1]

LTS1 E^0qi：iti：i0q2efiBJCi：i：iJ5 TIF [H

皿

d 1 ：l LTS]230q0200q2£SBJC[(IJ L TIP

[Biml ]]

LT5]£3MOZO[MZESlLfCLDJ3. ILF [j

皿

4

1151£330501268

隧口。』

£

」于 口 皿

A

Output. H-ap Froj act ion U

AV

.

LT51HtM

口£

13

空

EHBJCmjEL TIE L：B

皿

R 1

:Reorder Files .

Output File Buiie,

Delete

Inclusive

：

rance ezicorpasse-s all the files

i

' Ewc'.uiivn

：

A^corip bSE4s f

L

1

A

awLip

Dutput Rgsult to '"* ?il« ■"■ Iknory

居“

isr Output FaLciisiiis Choose

吨

Xearest jfle: fihtor

二三

I

二

］

匚工”二

确保波段

鼠标拖动，

I

1-7序号，从b1-b7,排序如下:

6.右边窗口设置坐标系如下：UTM,WGS-84,49N

7.定义文件名后，生成一个整的影像文件

同理，依次将其他文件夹下的7个波段合并成各自文件。

（二）多个文件镶嵌拼接成一个整的文件

注意：该步操作比较复杂，拼接文件可能存在色差不均衡问题，具体请多网上查些资料;

1.基于地理坐标进行拼接，操作如下：

File

Essie

T

DO

I

C

Cla^^ificstiDn TranEform Filter Spectral Map Vectoir T opo gra]

粘贴

更

D

Resize Data [Spatial/Speclral] Subset

Data via ROIs Rotaae/Fli p Data

Layer Stacking

启口/①缨叵]三三.

嚏。

A / A ③== = == = IB

二

1

：■— 口

Layer

田

and

：■— O L4year

小皿

3

：■— □

Layer (Sand

Comeirt Data (BSQj BIL BIP)

Stretch Data

Stali«ic&

Spatisl Statistics

匚

hwrigu Decertion

■■■ □ Layer

(BanH

：■— n L*yer

$皿

8

□ Layer (Band 0 ™i M«p

Life

卜

卜

►

Me-asurErnc-na

T

OD

I

Band Meth

<

,* Cray Scale

「

RJ

J

]

Spectral Math

Segmentation Image

Region Of Interest

卜

卜

卜

Selected I

L

叼叱

(E

皿

d 1:LT51ZS>

Preprocessing

Files

将上步生成的三个文件导入进来

Mwai eking

Masking

PizeJ Eas&d

G cDrerferenced

Tiled QuickBird Product

卜 卜

Mosaic Input Files

#1 1 ayer

E

t aclang3. img [R

电

[Gr e

#2 layerstacki

d]

enj

#3 layerEtagki KLue

# w XU * Y]

3 .分别右键文件名，选择Edit Entry （三个文件操作一致）

#2 1 ay er s t ^cki . img [Green]

#3 layer sticking, img [Elue ]

1a1

r

且

rztacking3. img： [Rgd]

Rai&e Image-10 Top

Rai&e Image- one Position

Remove Entry

X

0

4 .设置Data Value to Ignore背景值为0,羽化距离根据需要设置（不固定）；Color

Balancing（颜

色平衡参数，其中Fixed为以该文件为标准，其他影像进行调整，可对其中一个文件设置

为

Fixed，其他两个文件设置为Adjust）

5 .File-Apply，影像拼接

> Mosaic 12SS1 x 123^1 (UT

File Import O p-ti o-n s Help

九

ppb

Ca n

GE

I

拼接结果如下:

（三）水体区域提取

中波段运算，如下:

Easic Tools ClaEsificaaion rranEform Fili

Resize Data (Spatial/Specira Ij

Subset 0

日也

via ROIs

Rctats/Flip Data

Lajffiir Stacking

Convert Data (BSQ. BIL BFPj

Stretch Dets

Sestictirs

Spatial Statistics

Change Deteoicn

M&asurement Tool

Band Math

Specfcral Math

Segmentstion Image

Region Cif Interest

Mosaicking

Masking

PrEprocessing

2 .输入以下表达式

(

b2*1.0-b4)/(b2+b4) gt 0 (可用其他方法，依实际情况而定)

Ba nd Math

3 .分别设置算法中各个变量对应的波段，b2表示第3个波段，b4为第5个波段

4 .根据研究区域进行裁剪，并统计其中为1的像元个数，影像加载显示后，加载矢量文

件:

0 #1

Band Math (3

万

1.0-b4)/(b2+

一

File Oue ri ay Enhance rook Window

Ar r i...

CI

JZZ

■ cuticn..

'ontoir I in

，

三,

HH|

Zeniit/

二『

id Lire£...

>cyiLii L- Jilcrct：t,

■IH

vcctzrc

：

..

MUmH—i

5 .加载区域shp文件，第一次加载时后缀选择.shp

-si

Impart Vector Files Parameters Bifl

会自动生成一个evf文件，下次打开直接加载evf即

可。

6.

将文件转换成ROI兴趣区

7 .基于Roi文件对影像进行裁剪;

Basic Tools Classification Transfomi Fil

Resize Data (Spatial/SpectraI) Subset Data via R.015 R.&tate/Flip Data

Layer Stacking

Convert Data (BSQ BIL, BIP)

8 .裁剪后范围如下:

9.对裁剪后影像中所有水体进行统计（值为1的像元个数）

Basic Tools ClassificatiD-n Transfo-rm Filler SpectriaI Map Ve-ctor Topogtra

Re-size Data (Spatia l/Sp cctra I)

£ub

E

et Data via ROIs Rotate/Flip

Data Layer Stacking

Conveiirt Data (ESQ BIL, BIP)

Str&tch Data

Statistics

Spatial Statistics

Change Detection Measurement

Tool

Band Math

Spectral Math

Segmentation Image

Region Of Interest Mcsaidcing

Masking

Preprocessing

Spati al Si

Jutpnt Resu

Eni. er Dutp

-

Compute Statistics

View Statistics File

Sum Data Bands

Generate Random Sample

nr

勾选"Histogram”,如下:

。

Compute Statistics Parameters

Samples Resize Factor 1- UO'J--

IIIII IV

Lines Eesi ze Factor L 0

口

000

# 0input to the Screen

Oaiput to

i

Stitiaiica File

-

Oatput to & Teit Report File

OK QiifiUb C uical Raj

av

t Pr ac i

E L

an.

Stat'sties Results: [Memory2'

I File Options

- 匚

Select Plot

T

|

Clear Plot 1

Select Stat,

，

]

istcgrcift

and 1

H

LD

更二.器

1

Mcc.n

匚

t dev

0.2125

：

3

Npts 1 ota_ tercel

0 0. J47i^7

151445351 151445351

7

匚

43E7n 1 nFQRaq?" 4

n

4.

ilenane: [MenorvZ] (12t

31

K

12

IV - Fi: 1 I

G

IEHW

(1 nF 9FR 921 1nirtv

：

F5.25!

)

其中

1

的像元个数为水体，

Npts*30*30

(

30

为分辨率，表示

个像元的大小)为水体面积值

本文标签： 文件影像水体进行波段