我下载的数据包含以下几个集:
[1] "navigation_data/latitude" "navigation_data/longitude"
[3] "geophysical_data/Chl" "geophysical_data/flag"
Chl是我需要的数据,本来应该Chl是有坐标的,但是这个卫星的数据经度、纬度是单独的两个数据集,和Chl是分开的。我想实现的结果是数据坐标值提取对应的Chl值,但是现在我提出来的Chl没有坐标,是一个2780*2780的栅格。
我用R的raster包提取GOCI2的chl数据,代码如下
nc2chl <- raster("Wavelet/GOCI/GK2B_GOCI2_L2_20220707_061530_LA_S009_Chl.nc",
varname="geophysical_data/Chl")
output:
[1] "vobjtovarid4: **** WARNING **** I was asked to get a varid for dimension named pixels_per_line BUT this dimension HAS NO DIMVAR! Code will probably fail at this point"
[1] "vobjtovarid4: **** WARNING **** I was asked to get a varid for dimension named number_of_lines BUT this dimension HAS NO DIMVAR! Code will probably fail at this point"
使用crs查看nc2chl的投影
crs(nc2chl)
Coordinate Reference System:
Deprecated Proj.4 representation: NA
Warning message:
In wkt(x) : CRS object has no comment
想请大佬帮忙解决一下这个问题!不胜感激!
http://www.khoa.go.kr/nosc/satellite/searchL2.do
数据16M大于10M上传不了附件,这个网站已经到下载的界面了,右上角切换英文,Output选择Chlorophyll Concentration,点Search就可以下载nc文件了,随便下一个就好。
相似问题