Fly~~~的专栏

import arcpyfrom arcpy import envfrom arcpy.sa import *import osimport os.pathimport sysarcpy.env.workspace = "E:\\xuancheng\\remotesensing\\tif"rootdir = 'E:\\xuancheng\\remotesensing\\tif'for

目的 现有某地区土地利用类型图，图中不同的数值代表不同的土地利用类型，如像元值为20的，我们重新赋值为0，像元值为123的，我们重新赋值为15等等。 图1 原始图像 图2 结果图像 代码如下from gdalconst import *from osgeo import gdalimport osrimport sysimport copy#实现栅格计算器中的con函数功

目的 现有两幅栅格图像，一个是某地区道路栅格图，一个是某地区土地利用类型图，需要将道路叠加到土地利用类型图中，即叠加后，重合的像元值以道路图为准，其余的像元值仍是土地利用类型图原有的像元值。 图1 道路信息图 图2 土地利用类型图 图3 结果图 具体实现 from gdalconst import *from osge

1.单个文件的extractimport arcpyfrom arcpy import envfrom arcpy.sa import *arcpy.env.workspace = "E:\\Data_Ex\\FJJXtif"arcpy.env.extent = "E:\\Data_Ex\\FJJXtif\\20006.tif"arcpy.env.snapRaster = "E:\\Data

转自：http://bookshadow.com/weblog/2014/10/12/how-to-create-windows-executable-exe-from-python-script/Python是一种简单而强大的编程语言，适用于编写脚本，甚至于应用程序的开发。Python可用的各种GUI包使得利用Python编写全功能的应用变为可能。这很好，但你有没有想过将你编

参考自：http://www.liaoxuefeng.com/wiki/001374738125095c955c1e6d8bb493182103fac9270762a000/001411031239400f7181f65f33a4623bc42276a605debf6000>>> from collections import OrderedDict>>> ident = OrderedDict(

程序仅供参考。import sysimport gdalfrom gdalconst import *from osr import SpatialReferencefn = r'E:\work\Ex1016\xian_mask1.tif'print fnds = gdal.Open(fn,GA_ReadOnly)if ds is None: print 'cannot open

1.加法操作import arcpyfrom arcpy import envenv.workspace = "E:\\work\\cppservice\\RasterCalc" #工作空间可以设置也可不设置inRaster1 = "E:\\work\\cppservice\\finalResult\\propertyGower.tif"inRaster2 = "E:\\work\\cppse

本程序所采用的方法并不是最优方法，ARCGIS已经提供了相关的函数供调用。本程序仅供参考。程序说明： 文件夹E://work//EVI_Data_tif中存放的是某地区2000-2010年的EVI图像，其中每个年份共13幅。目的是将每年的13幅图像的每个栅格相加求均值，生成相应年份的tif。例如，将2000年的13幅图像相加求均值生成2000.tif，里面的每个栅格的值就是13幅

在安装好grass之后，只需设置如下的环境变量，便可在python中调用grass中的模块功能：GISBASE= D:\Program Files (x86)\GRASS GIS 7.0.svnPATH=%GISBASE%\etc;%GISBASE%\etc\python;%GISBASE%\lib;%GISBASE%\bin;%GISBASE%\extrabin;C:\Python27\ArcG