我正在尝试使用R Raster软件包校正旋转的栅格。 以下是mnf_filename是我尝试读取的光栅文件的实现。
mnf <- brick(mnf_filename)
引发以下警告：
此文件轮播 对此类文件的支持是有限的，数据处理的结果可能是错误的。谨慎操作并考虑使用“纠正”功能
运行整流功能
mnf_r <- rectify(mnf)
引发以下错误：
.intersectExtent（x，y，validate = TRUE）中的错误： 对象不相交 另外：警告消息： 1：在.local（x，...）中：min y大于max y 2：在.local（x，...）中：min y大于max y
当像这样仅校正栅格堆栈中的一层时，
mnf_r <- rectify(subset(mnf, 1)
产生以下错误：
validityMethod（object）中的错误：无效范围：ymin> = ymax
我不知道怎么了。以下是“ extent（mnf）”的输出：
class：范围
xmin：720770.2
xmax：752353
ymin：4328909
ymax：4361487
以下是gdalinfo(mnf_filename)的输出：
[1]“驱动程序：GTiff / GeoTIFF”
[2]“文件：X：\ gcwu \ AVIRIS \ f150429t01p00r07_refl_hpc18_v1_topo_American_MNF2.tif”
[3]“尺寸为701，1821”
[4]“坐标系为：”
[5]“ PROJCS [\” WGS 84 / UTM区域10N \“，”
[6]“ GEOGCS [\” WGS 84 \“，”
[7]“ DATUM [\” WGS_1984 \“，”
[8]“ SPHEROID [\” WGS 84 \“，6378137,298.257223563，”
[9]“ AUTHORITY [\” EPSG \“，\” 7030 \“]]，”
[10]“ AUTHORITY [\” EPSG \“，\” 6326 \“]]，”
[11]“ PRIMEM [\”格林威治\“，0，”
[12]“ AUTHORITY [\” EPSG \“，\” 8901 \“]]，”
[13]“ UNIT [\” degree \“，0.0174532925199433，”
[14]“ AUTHORITY [\” EPSG \“，\” 9122 \“]]，”
[15]“ AUTHORITY [\” EPSG \“，\” 4326 \“]]，”
[16]“ PROJECTION [\” Transverse_Mercator \“]，”
[17]“ PARAMETER [\” latitude_of_origin \“，0]，”
[18]“ PARAMETER [\” central_meridian \“，-123]，”
[19]“ PARAMETER [\” scale_factor \“，0.9996]，”
[20]“ PARAMETER [\” false_easting \“，500000]，”
[21]“ PARAMETER [\” false_northing \“，0]，”
[22]“ UNIT [\” metre \“，1，”
[23]“ AUTHORITY [\” EPSG \“，\” 9001 \“]]，”
[24]“ AXIS [\” Easting \“，EAST]，”
[25]“ AXIS [\” Northing \“，NORTH]，”
[26]“ AUTHORITY [\” EPSG \“，\” 32610 \“]]”
[27]“ GeoTransform =”
[28]“ 720770.204，13.16436670051244，-12.27597040459283”
[29]“ 4352881.37，-12.27597040459283，-13.16436670051244”
[30]“元数据：”
[31]“ AREA_OR_POINT =面积”
[32]“ TIFFTAG_XRESOLUTION = 1”
[33]“ TIFFTAG_YRESOLUTION = 1”
[34]“图像结构元数据：”
[35]“ INTERLEAVE = BAND”
[36]“转角坐标：”
[37]“左上（（720770.204，4352881.370）（120d26'23.54 \” W，39d17'49.87 \“ N）”
[38]“左下（（698415.662，4328909.058）（120d42'21.35 \” W，39d 5'12.32 \“ N）”
[39]“右上（（729998.425，4344275.915）（120d20'9.24 \” W，39d13'2.37 \“ N）”
[40]“右下（（707643.883，4320303.603）（120d36'6.99 \” W，39d 0'25.66 \“ N）”
[41]“居中（714207.043，4336592.486）（120d31'15.89 \“ W，39d 9'7.89 \” N）“