The original elevation has a spatial resolution of 0.0041666665°, with a total of 3124 columns and 2641 rows, as shown in the figure below:
Observing the above figure, we can see that the longitude range is from -102.0083363186 to -88.9916701726, and the latitude range is from 48.99791691515 to 60.00208314165. If the cmd file of the lapgrd process directly uses the mentioned range as the interpolation range, it will result in an error. The reason for this error is that:
1 2 | -102.0083363186+3124*0.0041666665=-88.9917 48.99791691515+2641*0.0041666665=60.0021 |
So, the range in the cmd file of lapgrd should be represented as follows:
1 2 | -102.0083363186 -88.9917 0.0041666665 48.99791691515 60.0021 0.0041666665 |
If the error message 'GETLIM - LIMITS EXCEED SURFACE COORDINATE LIMITS' is encountered, the root cause of the issue may be that the range set in the 'splina' command is smaller than the actual DEM (Digital Elevation Model) range. To resolve this, simply adjust the range in the 'splina' command to be slightly larger to cover the entire extent of the DEM. This should allow the interpolation process to proceed without encountering the error.
Share another successful case study:
Lapgrd:
1 2 | 0.29531831E+06 0.17793583E+07 30 0.24045270E+07 0.40906770E+07 30 |
When the spatial extent of the lapgrd's cmd file is difficult to determine, consider specifying the spatial range precisely in the splina's cmd file, and use the .sur file obtained from splina's computation to determine the spatial extent of lapgrd.
当插值温度(或者其他情况需要表面协变量的)时,需要同时输出多个表面结果,需要注意此时的DEM文件在第二步的参数中需要多次设置,设置数量要与输出数量一致。
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 | TEM11.sur 0 1 TEM11.cov 0 1 1 8871611.457651 10396111.457651 250 2 5131802.4564656 6137302.4564656 250 0 2 dem1.txt dem1.txt dem1.txt dem1.txt dem1.txt dem1.txt dem1.txt dem1.txt dem1.txt dem1.txt dem1.txt dem1.txt 2 -99.0 tem200001.grd tem200002.grd tem200003.grd tem200004.grd tem200005.grd tem200006.grd tem200007.grd tem200008.grd tem200009.grd tem200010.grd tem200011.grd tem200012.grd (100f10.3) |
不建议批量插值,可能会出现意想不到的错误,推荐是一次插值一张而已。
No comments:
Post a Comment