小编Lar*_*son的帖子

我正在尝试将具有不同单位的辅助y轴添加到3D绘图中.

[m2_array, ~ , ~] = F_readBin('amb.bin');
amb = m2_array(:,:,lat);

surfc(light,'LineWidth',0.001); 
ylim([1 24]); xlim([1 size(light,2)]); title(['@ ',num2str(lat),'°N']);
xticks([0:50:size(m2_array,2)]);
labels=cellstr(num2str((day_start:50:day_end)')); xticklabels(labels);
xlabel('Season days'); ylabel('Daytime{[hours]}');zlabel('surface light 
[\mumol m^2 s^-^1]')
colormap winter;

但是,我能找到的所有解决方案,例如yyaxis似乎只适用于2D图.是否有冲浪,网状,冲浪地块的工作？

我有一张海冰地图作为"geotiff".最终的目标是在特定的纬度/经度坐标处提取海冰浓度.

地理信息可在以下网址找到:http://www.iup.uni-bremen.de:8084 /amsr2data /asi_daygrid_swath/n6250/2015/jan/as-AMSR2-n6250-20150101-v5.tif

我想要做的是使用raster()加载geotiff,然后通过我的位置覆盖它,然后使用extract()函数从特定位置的栅格文件中获取值.

但是我的纬度/经度积累在地图的中心.我哪里出错了？非常感谢任何帮助或输入!

library(raster)
library(sp)

r1 = raster("test.tif")

##check plot
plot(r1)

## check projection
projection(r1)

mydf <- structure(list(longitude = rep(22,7), latitude = seq(60,90,5)),.Names = c("longitude","latitude"), class = "data.frame", row.names = c(NA, -7L)) 

### Get long and lat from data.frame.
xy <- mydf[,c(1,2)]
spdf <- SpatialPointsDataFrame(coords = xy, data = mydf,
                           proj4string = CRS("+proj=stere +lat_0=90 +lat_ts=70 +lon_0=-45 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"))

 points(spdf, col="red", lwd=2)

 ##extract RGB values as …