最近的一个项目需要在电子海图中下载已知水深点,导出点的地理坐标(经纬度).然后在arcgis中打开这些地理坐标输出为shp,利用GDAL读取不同波段的点对应的像元值,从而构建水深和像元值的对应关系. 其实想要根据经纬度得到像元值有两种方法. 1.通过arcpy选择点输出像元值,这个方法简单快捷,依赖于arcgis. 2.通过GDAL读取经纬度(度分秒),转换成小数形式.使用GDAL的CoordinateTransformation方法将经纬度转换成投影坐标,再转换成行列坐标,最后读取像元值.这里