gdal读取栅格数据,其中行列号和返回的数组索引是相反的,所以赋值的时候一定是行列互换,否则读取的时候就会越界 array_band raster_band.ReadAsArray(0, 0, raster_ds.RasterXSize, raster_ds.RasterYSize).astype(np.float64)for i in rang…
grep得到的内容用sed处理
-E or Extended Regular Expressions
I mentioned extended regular expressions earlier. FreeBSD (and Mac OS X) uses “-E” to enable this. However, FreeBSD later added the -r command to be compatible with GNU sed. -e <script>或…