地圖代數(shù)中的條件運(yùn)算:

CON的語(yǔ)法:

Usage:? (out_grid) CON (<condition>, <true_expression>, {<condition>, <true_expression>}, ...
??????????????? {<condition>, <true_expression>}, {false_expression})

實(shí)例如下,將河流中NoData值變?yōu)?

con(isnull(D:\DB\河流.img),0,D:\DB\河流.img)

自動(dòng)生成一個(gè)以輸入柵格圖像的外接邊框?yàn)榉秶臇鸥駡D像。如果需要賦值為0的NoData區(qū)域大于這個(gè)外接長(zhǎng)方形，可在Single?Output Map Algebra對(duì)話框Environment->General setting->Extent。

注:setnull的用法,作用是把柵格值設(shè)為NoData

Usage:? (out_grid) SETNULL (<condition>, {grid | scalar | number})

??????????? setnull([raster] == 0 , [raster])

????????????其中[raster]為要處理的柵格,后面的[raster]表示若前面條件不滿足則用這個(gè)值來代替.

con函數(shù)可以嵌套，如下所示：

con ( [1991rockyIndex.img] < 1.1 , 1.0 , con ( [1991rockyIndex.img] > 1.1 & [1991rockyIndex.img] < 1.45 , 2.0 , con ( [1991rockyIndex.img] > 1.45 & [1991rockyIndex.img] < 1.75 , 3.0 , con ([1991rockyIndex.img] > 1.75 , 4.0 , [1991rockyIndex.img]))))

?

另一例：

outraster = con( [Extract_LUCC_GX_MCD12Q1.A2001001.Albers.500m.tif] == 0 , 1.07, con( [Extract_LUCC_GX_MCD12Q1.A2001001.Albers.500m.tif] == 1 , 5.13, con( [Extract_LUCC_GX_MCD12Q1.A2001001.Albers.500m.tif] == 2 , 5.13, con( [Extract_LUCC_GX_MCD12Q1.A2001001.Albers.500m.tif] == 3 , 5.13, con( [Extract_LUCC_GX_MCD12Q1.A2001001.Albers.500m.tif] == 4 , 5.13, con( [Extract_LUCC_GX_MCD12Q1.A2001001.Albers.500m.tif] == 5 , 4.14, con( [Extract_LUCC_GX_MCD12Q1.A2001001.Albers.500m.tif] == 6 , 6.17, con( [Extract_LUCC_GX_MCD12Q1.A2001001.Albers.500m.tif] == 7 , 5.43, con( [Extract_LUCC_GX_MCD12Q1.A2001001.Albers.500m.tif] == 8 , 5.43, con( [Extract_LUCC_GX_MCD12Q1.A2001001.Albers.500m.tif] == 9 , 1.08, con( [Extract_LUCC_GX_MCD12Q1.A2001001.Albers.500m.tif] == 10 , 1.08, con( [Extract_LUCC_GX_MCD12Q1.A2001001.Albers.500m.tif] == 11 , 1.08, 1.08))))))))))))

?

惱人的CON