Python+IDLE批量值提取至点/定义投影

# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
arcpy.env.workspace = r"H:\Graduatu student\TuhaBasin_tasks\ssrd\tiqu" #栅格路径
# Set local variables
inPointFeatures = "H:\数据集\辐射实测站点\西北辐射站点插值_矢量点\xxx.shp" #shp文件路径及名称
inRasterList = arcpy.ListRasters("*", "tif")
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute ExtractValuesToPoints
ExtractMultiValuesToPoints(inPointFeatures, inRasterList, "NONE")
print("finish!")
#路径中尽量别有中文

#别人分享的代码，不知道具体的引用来源，如有侵权请留言！

# -*- coding: UTF-8 -*-
#需要自行更改的如下：inws：输入路径；outws：输出路径；Coordinate_System：目标投影坐标系
import glob
import os
import arcpy

#输入路径  应该注意，中文路径，会导致读不出文件
inws = r"H:\Graduatu student\Resample\calulator"

# 参考文件路径    使用栅格数据集（从其导入方形像元大小）的路径
Coordinate_System="PROJCS['CGCS_2000_UTM_Zone_46N',GEOGCS['GCS_China_Geodetic_Coordinate_System_2000',DATUM['D_China_2000',SPHEROID['CGCS2000',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',93.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]"
# 输出路径
outws = r"H:\Graduatu student\Resample\1990-2020/"

# 利用glob包，将inws下的所有tif文件读存放到rasters中
rasters = glob.glob(os.path.join(inws, "*.tif"))

# 循环rasters中的所有影像，进行按掩模提取操作
for raster in rasters:
    print(str(raster))
    nameT = os.path.basename(raster).split(".")[0] + ".tif"    
    outname = os.path.join(outws, nameT)  # 合并输出文件名+输出路径

    arcpy.ProjectRaster_management(raster, outname, Coordinate_System, "NEAREST","#","#","#","#")

print("OK!")