栅格面积统计方法_计算机软件及应用_IT计算机_专业资料

1、栅格面积统计统计栅格影像的面积，要求栅格数据的坐标必须为平面直角坐标；如不是则 需要进行投影转换。经纬度到平面直角坐标转换仓u建一个本地数据库 file geodatabase,右键菜单 new一file geodatabasecontents preview i metadataq arccatalog arclnfo d:(5±1e-f申审etalo6d:vudddnud3be:f:mdaqintsecaa® 田es s e e s s ® gcagl fsem-zio folderfile geodatabasepersonal geodatabasenam

2、e_| srecycle.bin_i downloads_|envi二envi.f oracle二system volume information vod.cache.datahtpgfolderfolderfolderfolderfolderfolderfolderpersonal geodatabaseraster dataset layer.q group layer口 shapefile.貉 coverage relati on ship class.3 lurn feature class.0 toolboxo arcinfo workspace13 dbase table回 in

3、fo table.为创建的数据库输入名称，如图屮的rasterdata，也可使用默认名称© catalog由g c:曰 d:w _j 1recycle.binffl _| downloadsm enviw i envi.f1*1 _i oracle田 _| system volume informationiq rasterdata iki ij ext0906.1005pl亠e：f：database connectionsdatabase serversgis serverscontents | pr«vi«v metedtta |name_isrecycle.

4、bin_j downloadsi envi_jenv1.fij oraclesystem volume information vod cache da怙j rasterdata |cjtestihl ext0906.1005pli typefolderfolderfolderfolderfolderfolderfolderfile geodatabase personal geodatabase raster dataset:gqj冒紋 ma-& 田 0 ainteroperability conn ections search results导入要处理的栅格影像。在上一步创建的ra

5、sterdata数据库节点上鼠标右键import一rasterdatasetsj catalog 由；bu 曰6s q srecycle.bins o downloads田 q envis cj envi.fs cd orades lj system volume information ffl cj vod.cache.dataa ext(je： xf：data bi 0copyctrlfcnamei typedeleterenamerefreshf2生怙匕 new gis seinteropsearchimport ifeature £las$ (single)十exportfe

6、ature class (multiple)”compress file geodatabase.tafele (single)yncompress rle geodatabase.iable (multiple)q compact databaseraster datasets.raster datasets (mosaic).e search.q publish to arcgis server为 xml workspace document.distributed geodatabase卤 properties.t选择要导入的影像数据，然后点击“0k”按钮导入数据。弟 raster to

7、 geodatabase (multiple)output geoaatabase|c： rasterdata. gdbconfiguration keyword (optional)cancelenvironments.show help »设置栅格影像的空间投影。在导入的影像数据节点上鼠标右键，选择propertiesstylesheet： |fgdc esri二1 首 筍国空ul- xcontents preview | metadatanametfl catalog 宙 0 c： 目 gl d:srecycle.bindow nloads0-:eb-.0a-:ea-envie

8、nvij oraclevod.cache.data:b(3 rasterdata3 ® ext0906 1005d1 |testi 田翻 ext0906_10 宙 g* e: 审3 f: 0 g database con3 database serx电 gis servers 直 interoperabilits 紋 search resultsrenamef2create layer.load屉 copyctrl+cexportbuild pyramids.calculate statistics.总 publish to arcgis server.properties.在 ra

9、seter dataset properties 窗口中点击 spatial reference 右侧的“edit” 按钮raster dataset propertiesgeneralpropertyvueno data value-3.4028235e*038colormapabsentpyramidslevel: 2. resamping nearest neighborbuid.compcetsionlz77s extenttop43 488304074left118 839202799right125 791202799bottom38 66s304074spatial refere

10、ncev undefinedsedit. |linear ixrmgularunte) statisticsoptions bandbudd parametersmin182 2750009765625max830.31298828125mannz辺tu唧q翊|i取消 |应用) |在弹出的 spatial reference properties 窗u 中，点击 “select”按钮;spatial reference properties在 browse for coordinate system 对话彳匡屮双击“geographic coordinate systems进入;browse

11、for coordinate systemlook in:coordinate systems|二医nt|ss|) geographic coordinate systems_l projected coordinate systemsname:addshow of type： |coordinate systems二|h 二-cancela双击“asia”进入;选择栅格影像数据的坐标系,如"xianl980.prj,点击“add”按钮完成添加在arcmap下加载导入到rasterdata数拯库中的栅格影像;cg* 由囱由由由曲曲曲曲田由甲申 甲申！eb-;由曲；(=!- o. 1k

12、b/s -r o.o3kb/s设置 投影。双击 arctoolbox一 datamanagement tools 一 projections andtransformations一raster一project raster30 analyst toolsanalysis toolscartography toolscorwersion toolsdata interoperability toolsdata monagemect toolsdata comparisondatabasedi connected editing distributed geoddtabasedorr>ain

13、sfeature classfeaturesfieldsfile geocjitibesegerterdgeneralizationindexesjoinslayers arid tbleprojections or»d trnsformationssfeaturesrastera hip；kirrorizminxuuu iinput raster：点击右侧的下拉箭头，选择加载到arcmap中的影像r«t«re - 旦團回團-+1xj21丄 3 -okcancelen vronments.show heb »project rasteroutput co

14、ordinate system：点击右侧的按钮;p project rasterinput raster |ext0906j005plinput coordinate system (optional)|gcsjxanj980output raster dataset |d:krtsterdta. <dbext0906_1005p1_projectr<sterloutput coordinate systemgo<rphic tr<n讦onhuion (oplioxul)canedenvironmentsshow help »在弹出的 spatial refe

15、rence properties 窗口中，点击"select”按钮;spatial reference propertiesxy coordinate syxtemname：details：|unknown确走 | 取消 | 匚用*|在 browse for coordinate system 对话框屮双击“projected coordinate systemsv进入;browse for coordinate systemlook in:j i projected coordnate systems3 z拥|匡制開i选择合适的椭球体，如“gauss kruger”高斯克吕格，双击

16、进入;brov/se for coordinate system_j arc (equal arc-second) _| world _j continentalcountysystemsj gauss kruger._i national gridslj polar 一j state planelj state systems utmname:show of type选择栅格影像的坐标系,如“xian 1980",双击进入;d刮ol鋼=1底iml弗|browse for coordinatw systemlook in:gauss kruger_ beijing 1954other

17、 gcs_j pulkovo 1942_l pulkovo 1995 xian 1980 jname:show of type: | coordinate systems根据山歌影像数据经纬度范围等信息，选择合适的坐标带。如选择“西安80,123度中央经线所在的3度分带”。browse for coordinate systemlook in: |!_j xan 1980oxian 1980 3 degree gk cm 102e.prj 'xian 1980 3 degree gk cm 105e.pr xian 1980 3 degree gk cm 108e.pr ®x

18、ian 1980 3 degree gk cm llle.pg ©xian 1980 3 degree gk cm 114e.pr 'xian 1980 3 degree gk cm 117e.pr 噬，幻刖 1980 3 dear” gk cm 120e.orixian 1980 3 degree gk cm 123e.prji&'xian 1980 3 degree gk cm 126e.pg ®xian 1980 3 degree gk cm 129e.prj ® xian 1980 3 degree gk cm 132e.prj 觀

19、xian 1980 3 degree gk cm 135e.prj ®xian 1980 3 degree gk cm 75e.prj ®xian 1980 3 degree gk cm 78e.prj i® xian 1980 3 degree gk cm 81e.prj ©xian 1980 3 degree gk cm 84e.prj ®xian 1980 3 degree gk cm 87e.pfj 癒xian 1980 3 degree gk cm 90e.prjname:rrri xian 1980 3 degree gkcm 12

20、3e.ptshow of typei coordinate systemscancel |点击“ok”按钮，完成栅格影像的投影转换。input rasterzjcancelenvir onmen ts.创xltj丄show help >>|ext0906j005plinput coordinate system (optional) |gcs-xiu980output raster dataset (darasterdat& <dbex2906005pl_projectrsterloutput coordinate system|xi 0.1980.3 jucr“.g

21、kj：il123rgo(rphic transformation (optional)投影转换后的数据被重新加载到arcmap中qrawing *ntrtkd ft -1 "j emt0906_1005p l_firotectr*ster 111valuehigh : 8233删话的 low : 1376hr格事0三 0 ext0906.1005plvaluehigh : 830.313lew; 182.276xg arctoc<>ox- 30 aral/tt tooh a gp tools 匍 c«rt09f«phy too = s co耐took

22、s dau lnterop«fab b oata mr-»9<mef oata compridisconnected distnboted ge domm feature chssfield/file geod«ub generalaaffifflffl回sewmdh83320424 4902662416 meters点击view菜单卜的data frame properties菜单项,q untitled - arcmap - arcinfofile edit view bookmarks insert selection ools window plts

23、 help话& data view 口 layout view1:4, 141z543editor eaturezoom datael自0zoom layoutioolbars7 status baroverflow annotationeh scrollbars图 data frame properties.在 data frame properties 窗口屮选择 “coordinate system” 选项栏,data frame propertiesfrweproduct librtry on gr 3 dz map cccho |annotation groxq>x|e

24、xtent b«ctonglessize poseoe 丨ginkgener1 da.ta ?rone i coordint systen i iuunm&ticurrent coordinate system:gcsjqanl980 datum: d xian 1960select a coordinate system: o favorttes lj predefined 电匚layers s- *£j <ojstom>gcs xian 39001选择z前投影转换时选屮的坐标带，如卜图所示。(predefinedprojectedcoordinate

25、 systems一gauss kruger一xian 1980一xian 1980 3degree gk cm 123e)data frame propertiesi ©data frame properties选择"general"选项卡，更改下面的display项为“meters”，点击“确定”按钮。此吋，地图的坐标系以更改为平面直角坐标系，见arcmap窗口右下角状态栏的坐标信息便可知。统计栅格不同颜色像素的面积更改栅格影像专题图样式。在投影转换后的栅格影像上点击鼠标右键，选择properties 菜单;exto9o6.1oo5p l_pie) layers

26、e) s2oom to lyerzfiom to raster resolution yisible scale rangelow : 182.276display nodata asr show dass breaks usng cd values n use hdlshade effect确定valuehigh : 823.3hh low : 182.276 e) (2 exto9o6.1oo5pl valuehigh : 830.313在弹出的layer properties对话框中,左侧列表选择"classified"项，口j在右侧的“classes”项选择级数，点

27、击“确定”按钮。选择spatial analyst il具栏下拉菜单的“reclassify”项;interpolate to rastersurface analysiscell statistics.neighborhood statistics.zonal statisticszonal histogram.geclassifyraster calculatorconvertoptions在弹出的窗口中点击“ok”按钮，按当前的分类方法重分类。重分类后的栅格影像，如下图所示s 自 layers曰 3 reclass of ext0906.1005pl.projeceditor 1认：|

28、create new featurespatial analystel 823ext0906_1005pl_prpjectra$terl< value > 182.276001 318.022257 318.0222571 383.3815655 383.3815656 443.7132348 443.7132349 531.6969193 531.6969194 823.2999878 e b ext0906j005plvaluehigh : 830.313| low : 182.276为重分类后的影像添加面积字段。在重分类后的影像上点击鼠标右键，选择"open att

29、ribute table”untitled arcmap ardnfofile edit view bookmarks insert selection tools windoweditor 4 task:create new feature二j q s x ® x me 炉 11:4,125,593layers0 1口 2 3口 4口 5口 823电 copy x remove團joins and relatesopen attribute table zoom to layer0 ext0906.1c* zoom '©makevisible <valui

30、液 zfiom to raster resolution 182.276visible scale range 318.022 383.381 443.713 531.696datasave as layer file.create lyer package.13 ext0906_lc 囹 properties.value high : 830.313 | low : 182.276属性表窗口中点击“options”按钮,在弹出的菜单中选择“add fieldx it?) ® x5 c i|1:4,12s.593创 x k?ask: create new featuress of e

31、xt0906_1005pl_bc_prqjectrth) attributes of reclass of ext0906_1005plk.projectr q4(3 ho 必氏色次spa n find & replace.|绘 select by attributes.o clear selectionif?i switch selection 匡1 select all 回l3(xrowidvalue coumt0x460412909：23：03883448624s19193823xadd field.-ium all fields onsho field aliasesresto

32、re default column widths2.8.3<2oin$ and relates related jables create graph. add table to layout0 reload cacfce © print.reportsexport.record: m| !1 » | >i| show: | al selected records (0 out of 6 sejected)appearance.输入面积字段的名称,如“area”,字段类型选择“double”，点击“ok”按钮添加字段add fieldname: lype：areadoublefield propertiesprecision0scale0ok | cancel i获得像素单元格大小。在重分类后的影像上点击鼠标右键，选择“properties”,editor | 4 | 夕曰自layers回reclass o|毎 copy x remove口 5 823叩 open attribute iable2oin$ and relates zoom to layere) ext0906.1l<vall 佥 zqom to raster resolution 182