需求
一些MODIS EVI數(shù)據(jù)需要轉(zhuǎn)換為跟自己數(shù)據(jù)一致的投影巴帮。
定義坐標(biāo)系信息
可以先定義一個Shpfile的投影信息腹忽，然后把它的.prj文件用記事本打開，記得修改里面的符號“為'。
柵格圖像的投影定義
import arcpy
# 設(shè)置工作空間
arcpy.env.workspace="E:\\GPP\\test2"
# 得到所有tif格式影像
rasterlist = arcpy.ListRasters("*", "tif")
#定義坐標(biāo)系統(tǒng)
Coordinate_System = "PROJCS['WGS_1984_Albers',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Albers'],PARAMETER['false_easting',0.0],PARAMETER['false_northing',0.0],PARAMETER['central_meridian',105.0],PARAMETER['standard_parallel_1',25.0],PARAMETER['standard_parallel_2',47.0],PARAMETER['latitude_of_origin',0.0],UNIT['Meter',1.0]]"
#循環(huán)
for raster in rasterlist:
print str(raster)
arcpy.DefineProjection_management(raster, Coordinate_System)
print("OK!") #運(yùn)行完畢，顯示OK守伸。
要素類的投影定義
import arcpy
# 設(shè)置工作空間
arcpy.env.workspace="E:\\GPP\\test2"
# 得到所有要素
featurelist = arcpy.ListFeatureClasses()
# define coordinate system
Coordinate_System = "PROJCS['WGS_1984_Albers',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Albers'],PARAMETER['false_easting',0.0],PARAMETER['false_northing',0.0],PARAMETER['central_meridian',105.0],PARAMETER['standard_parallel_1',25.0],PARAMETER['standard_parallel_2',47.0],PARAMETER['latitude_of_origin',0.0],UNIT['Meter',1.0]]"
#循環(huán)
for ft in featurelist:
print str(ft)
arcpy.DefineProjection_management(ft, Coordinate_System)
print("OK!")
