利用arcpy和arctoolbox中的定义投影工具编写代码为pixelsValue.shp定义投影： WGS84 UTM 50度带；然后利用ExtractMultiValuesToPoints工具编写代码将GF1_PMS2_fuse2m_clip.tif所有波段的值提取到pixelsValue.shp点要素文件中；利用ExportXYv （Export Feature Attribute to ASCII）工具编写代码将pixelsValue.shp中各个波段的值导出到文本文件gf1_sample_points.txt中

以下是利用arcpy和arctoolbox中的定义投影工具编写代码为pixelsValue.shp定义投影为WGS84 UTM 50度带的示例：

import arcpy

# 设置工作空间和输出路径
arcpy.env.workspace = "C:/data"
output_shp = "C:/data/pixelsValue.shp"

# 定义投影参数
projection = arcpy.SpatialReference(32650)  # WGS84 UTM 50度带的投影参数

# 定义投影工具
arcpy.DefineProjection_management(output_shp, projection)

# 利用ExtractMultiValuesToPoints工具将GF1_PMS2_fuse2m_clip.tif所有波段的值提取到pixelsValue.shp点要素文件中
input_raster = "C:/data/GF1_PMS2_fuse2m_clip.tif"
arcpy.ExtractMultiValuesToPoints(output_shp, [[input_raster, "BAND1"], [input_raster, "BAND2"], [input_raster, "BAND3"]],
                                "NONE")

# 利用ExportXYv工具将pixelsValue.shp中各个波段的值导出到文本文件gf1_sample_points.txt中
output_txt = "C:/data/gf1_sample_points.txt"
arcpy.ExportXYv_stats(output_shp, [["BAND1", "MEAN"], ["BAND2", "MEAN"], ["BAND3", "MEAN"]],
                      "SPACE", output_txt, "ADD_FIELD_NAMES")

请注意，上述示例中的路径和文件名是示意性的，请根据实际情况修改为你的数据路径和文件名。