我正在尝试使用不同的5 * 5度瓷砖创建全局栅格。因此，我将标签+图像文件转换为带有gdal_translate的地理标记。之后我使用gdalwarp重新投影数据。这适用于0到180度经度之间的所有瓷砖，但会产生一个仅在-180到0度经度之间只有零的光栅。
数据来自Mars的CRISM摘要产品数据。
输入Proj4文件：
PROJCS["EQUIRECTANGULAR_MARS",GEOGCS["GCS_MARS",DATUM["D_MARS",SPHEROID["MARS_localRadius",3396000.0,0.0]],PRIMEM["Reference_Meridian",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Equidistant_Cylindrical"],PARAMETER["false_easting",0.0],PARAMETER["false_northing",0.0],PARAMETER["central_meridian",327.5],PARAMETER["standard_parallel_1",42.5],UNIT["Meter",1.0]]
输出Proj4文件：
GEOGCS["GCS_Mars_2000",DATUM["D_Mars_2000",SPHEROID["Mars_2000_IAU_IAG",3396190.0,169.8944472236118]],PRIMEM["Reference_Meridian",0.0],UNIT["Degree",0.0174532925199433],AUTHORITY["Esri",104905]]
有没有人知道为什么我的栅格中的数据仅在所描述的步骤之后由零组成？ （见下面的代码）
gdal_translate -a_nodata 65535
/home/kampsom/Documents/CRISM/Raw_data/$filename_lbl
/home/kampsom/Documents/CRISM/Converted/$filename_tif
gdalwarp -overwrite -t_srs
/home/kampsom/Documents/CRISM/Converted/Mars_2000.prj
/home/kampsom/Documents/CRISM/Converted/$filename_tif
/home/kampsom/Documents/CRISM/Converted/$filename_tif_repr