我在netcdf中有ERA5气候数据,在将ERA5数据从每小时转换为每天,并将它们计算在一起以得到一个每日产品后,我需要在该气候数据和一组GeoTIFF中相同地理区域的相应像素之间进行一些计算,具有更高的空间分辨率和不同的投影。netCDF在EPSG 4326中进行了投影,并配有lat/lon瓷砖。geotiff的分辨率更高(3.125km2),并在南极投影中投影。我已经在QGIS中一起打开了netcdf和GeoTIFF,并验证了一切正常,并且它们匹配
我一直不确定如何实现这一点,但我的感觉是,我需要重新投影原始netCDF以匹配Geotiff,同时重新填充/上采样netCDF数据,以便像素匹配
有没有办法简单地“匹配”这两个文件
(我曾尝试在gdal中重新投影,并在python中手动设置新投影,这不会导致任何错误,但也不会创建任何文件!-只是什么都没有发生:
import gdal
path = '//Users/gj/Documents/Sea Ice Project/AutoSeaIceProduc/adaptor.mars.internal-1618515503.3594973-11959-19-93fbbedb-7e77-4565-ba2b-05ef0e8b1c71.nc'#Path to file
out = '//Users/gj/Documents/Sea Ice Project/AutoSeaIceProduc/ReprojectedTest.nc'#Path to file
gdal.Warp(out,path, srcSRS='EPSG:4326',dstSRS='EPSG:3976')
NetCDF的Gdal信息:
Driver: netCDF/Network Common Data Format
Files: adaptor.mars.internal-1618515503.3594973-11959-19-93fbbedb-7e77-4565-ba2b-05ef0e8b1c71.nc
adaptor.mars.internal-1618515503.3594973-11959-19-93fbbedb-7e77-4565-ba2b-05ef0e8b1c71.nc.aux.xml
Size is 512, 512
Metadata:
NC_GLOBAL#Conventions=CF-1.6
NC_GLOBAL#history=2021-04-15 19:39:59 GMT by grib_to_netcdf-2.16.0: /opt/ecmwf/eccodes/bin/grib_to_netcdf -S param -o /cache/data8/adaptor.mars.internal-1618515503.3594973-11959-19-93fbbedb-7e77-4565-ba2b-05ef0e8b1c71.nc /cache/tmp/93fbbedb-7e77-4565-ba2b-05ef0e8b1c71-adaptor.mars.internal-1618515503.360069-11959-2-tmp.grib
Subdatasets:
SUBDATASET_1_NAME=NETCDF:"adaptor.mars.internal-1618515503.3594973-11959-19-93fbbedb-7e77-4565-ba2b-05ef0e8b1c71.nc":u10
SUBDATASET_1_DESC=[720x13x85] u10 (16-bit integer)
SUBDATASET_2_NAME=NETCDF:"adaptor.mars.internal-1618515503.3594973-11959-19-93fbbedb-7e77-4565-ba2b-05ef0e8b1c71.nc":v10
SUBDATASET_2_DESC=[720x13x85] v10 (16-bit integer)
SUBDATASET_3_NAME=NETCDF:"adaptor.mars.internal-1618515503.3594973-11959-19-93fbbedb-7e77-4565-ba2b-05ef0e8b1c71.nc":d2m
SUBDATASET_3_DESC=[720x13x85] d2m (16-bit integer)
SUBDATASET_4_NAME=NETCDF:"adaptor.mars.internal-1618515503.3594973-11959-19-93fbbedb-7e77-4565-ba2b-05ef0e8b1c71.nc":t2m
SUBDATASET_4_DESC=[720x13x85] t2m (16-bit integer)
SUBDATASET_5_NAME=NETCDF:"adaptor.mars.internal-1618515503.3594973-11959-19-93fbbedb-7e77-4565-ba2b-05ef0e8b1c71.nc":sp
SUBDATASET_5_DESC=[720x13x85] surface_air_pressure (16-bit integer)
SUBDATASET_6_NAME=NETCDF:"adaptor.mars.internal-1618515503.3594973-11959-19-93fbbedb-7e77-4565-ba2b-05ef0e8b1c71.nc":ssrd
SUBDATASET_6_DESC=[720x13x85] surface_downwelling_shortwave_flux_in_air (16-bit integer)
SUBDATASET_7_NAME=NETCDF:"adaptor.mars.internal-1618515503.3594973-11959-19-93fbbedb-7e77-4565-ba2b-05ef0e8b1c71.nc":strd
SUBDATASET_7_DESC=[720x13x85] strd (16-bit integer)
Corner Coordinates:
Upper Left ( 0.0, 0.0)
Lower Left ( 0.0, 512.0)
Upper Right ( 512.0, 0.0)
Lower Right ( 512.0, 512.0)
Center ( 256.0, 256.0)
样本geotiff的相关gdal信息
Driver: GTiff/GeoTIFF
Files: clip_asi-AMSR2-s3125-20200301-v5.4.tif
Size is 117, 163
Coordinate System is:
PROJCRS["IDL GeoTIFF Suport
Projection = Polar Stereographic
True scale: -70.0000deg
Gunnar Spreen, Apr 2004",
BASEGEOGCRS["WGS 84",
DATUM["World Geodetic System 1984",
ELLIPSOID["WGS 84",6378137,298.257223563,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
ANGLEUNIT["degree",0.0174532925199433]],
ID["EPSG",4326]],
CONVERSION["Polar Stereographic (variant B)",
METHOD["Polar Stereographic (variant B)",
ID["EPSG",9829]],
PARAMETER["Latitude of standard parallel",-70,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8832]],
PARAMETER["Longitude of origin",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8833]],
PARAMETER["False easting",0,
LENGTHUNIT["metre",1],
ID["EPSG",8806]],
PARAMETER["False northing",0,
LENGTHUNIT["metre",1],
ID["EPSG",8807]]],
CS[Cartesian,2],
AXIS["(E)",north,
MERIDIAN[90,
ANGLEUNIT["degree",0.0174532925199433,
ID["EPSG",9122]]],
ORDER[1],
LENGTHUNIT["metre",1]],
AXIS["(N)",north,
MERIDIAN[0,
ANGLEUNIT["degree",0.0174532925199433,
ID["EPSG",9122]]],
ORDER[2],
LENGTHUNIT["metre",1]]]
Data axis to CRS axis mapping: 1,2
Origin = (-1810937.500000000000000,-554687.500000000000000)
Pixel Size = (3125.000000000000000,-3125.000000000000000)
Metadata:
AREA_OR_POINT=Area
TIFFTAG_DATETIME=2020:03:02 05:15:07
TIFFTAG_DOCUMENTNAME=/ssmi/www/htdocs/data/amsr2/asi_daygrid_swath/s3125/2020/mar/Amundsen/asi-AMSR2-s3125-20200301-v5.4.tif
TIFFTAG_IMAGEDESCRIPTION=IDL TIFF file
TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
TIFFTAG_SOFTWARE=IDL 8.1, ITT Visual Information Solutions
TIFFTAG_XRESOLUTION=100
TIFFTAG_YRESOLUTION=100
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left (-1810937.500, -554687.500) (107d 1'47.10"W, 72d38'40.82"S)
Lower Left (-1810937.500,-1064062.500) (120d26'14.82"W, 70d47' 7.34"S)
Upper Right (-1445312.500, -554687.500) (110d59'45.57"W, 75d46'46.47"S)
Lower Right (-1445312.500,-1064062.500) (126d21'39.97"W, 73d32'29.83"S)
Center (-1628125.000, -809375.000) (116d25'58.56"W, 73d19'46.09"S)
目前没有回答
相关问题 更多 >
编程相关推荐