Inpolygon returning a NAN
조회 수: 1 (최근 30일)
이전 댓글 표시
data=netcdf.open('gisst_full.nc','NC_NOWRITE')
lat=netcdf.getVar(data,0);
lon=netcdf.getVar(data,1);
time=netcdf.getVar(data,2);
sst=netcdf.getVar(data,3);
sst(sst==-32768)=NaN;
ti=datetime(1800,1,1)+days(time);
n=length(ti);
sst1=[sst(181:360,:,:);sst(1:180,:,:)];
lon1=[lon(181:360);lon(1:180)]
sst2=permute(sst1,[2,1,3]);
lonn=wrapTo360(lon1);
[X,Y]=meshgrid(lonn,lat);
idx=inpolygon(X,Y,[160 170],[-35 -25]);
data=sst2(:,:,1333:12:1584);
data1=mean(data,3),'omitnan';
Jan=mean(nonzeros(data1(idx))),'omitnan';
When I run the above script I get as far as data1 but the last line to calculate Jan value I am getting 'NaN' for Jan instead of the mean value of my inpolygon. The gisst data is monthly SST values and I just want the mean for an area 160-170E and 25-35S as specified in my idx. What am I missing?
댓글 수: 0
채택된 답변
Steven Lord
2022년 12월 15일
The inpolygon function can't return a NaN value. The problem is later on in your code. If you ask for the mean of an empty array the answer is NaN.
mean([])
That implies that neither of the points in your inpolygon call are inside your polygon. You can check this by calling any. If it returns false, neither of the points are inside.
any([false false])
any([false true])
I can't offer any guidance for how to correct this problem since it's not clear to me from your code what your goal is. If you describe in words not code what the variables in your code represent (some of the names seem suggestive of their purposes, but suggestions can be deceiving) and describe what you intend this code to do (again in words not code) we may be able to offer some suggestions.
추가 답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Creating and Concatenating Matrices에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!